#!/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=aice
avgfs=''
trendfs=''
for i in aice_??.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 aice_trend_ensmean.nc
cdo -ensmean -selname,$var $avgfs aice_avg_ensmean.nc