#!/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.'

var=hi
avgfs=''
trendfs=''
for i in hi_??.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 -ensmean -selname,$var $trendfs hi_trend_ensmean.nc
cdo -ensmean -selname,$var $avgfs hi_avg_ensmean.nc