#!/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/norcpm/cases/NorCPM/ana_19500115_bccrfast
outfile=NPac_ts.nc
var=sst
lonlatbox=120,200,10,30
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_NPac_$(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_NPac_*.nc avg_ts.nc
cdo -O -chname,${var},std -ensstd ts_NPac_*.nc std_ts.nc
cdo -O -chname,${var},min -ensmin ts_NPac_*.nc min_ts.nc
cdo -O -chname,${var},max -ensmax ts_NPac_*.nc max_ts.nc
rm -f ts_NPac_*.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