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

## ensemble
dir=/trd-project4/NS9039K/shared/pgchiu/archive/norcpm_ana_f09_tn14/
outfile=NINO34_ts.nc
var=sst
lonlatbox=-150,-90,-5,5
test -f ${outfile} && exit 0
## loop for every member
i=0
for d in $(ls model_sst_ano_mem*.nc); do
i=$(($i+1))
o1ts=ts_NINO34_$(printf "%2.2d" $i).nc
test -f "$o1ts" && continue
## NorESM (BLOM/MICOM)
cdo -s  -b f32         \
    -fldmean  \
    -sellonlatbox,${lonlatbox}  \
    -setgrid,griddes.txt \
    -select,name=${var} \
    ${d} \
    ${o1ts}
done
## ensemble variables
cdo -O -chname,${var},avg -ensavg ts_NINO34_*.nc avg_ts.nc
cdo -O -chname,${var},std -ensstd ts_NINO34_*.nc std_ts.nc
cdo -O -chname,${var},min -ensmin ts_NINO34_*.nc min_ts.nc
cdo -O -chname,${var},max -ensmax ts_NINO34_*.nc max_ts.nc
rm -f ts_NINO34_*.nc
## cat to 1 file
cdo --reduce_dim -O -merge avg_ts.nc std_ts.nc min_ts.nc max_ts.nc ${outfile}
rm -f  avg_ts.nc std_ts.nc min_ts.nc max_ts.nc