#!/bin/bash

## Just run code in recipes, use it carefully
## example: use cdo to make cache file. Since cdo is quicker than ncl.
#;DIAG_NORCPM; RUNTHESECODES: echo '    no codes run here.'

test -f PSL_trend_ensmean.nc && test -f PSL_avg_ensmean.nc && exit 0
var=PSL
avgfs=''
trendfs=''
for i in PSL_??.nc; do
    avgf=$(echo $i | sed 's/\.nc/_avg.nc/')
    trendf=$(echo $i | sed 's/\.nc/_trend.nc/')
    test -f $avgf && test -f $trendf && continue
    cdo -s -O -trend $i $avgf $trendf 
    trendfs="$trendfs $trendf"
    avgfs="$avgfs $avgf"
done
cdo -s -O -ensmean -selname,$var $trendfs PSL_trend_ensmean.nc
cdo -s -O -ensmean -selname,$var $avgfs PSL_avg_ensmean.nc