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

## env for cdo
module -q purge
#module load CDO/1.9.9-iompi-2020a
source /cluster/projects/nn9039k/people/pgchiu/conda/py3env.sh

## ensemble
dir=/cluster/shared/NS9039K/norcpm_ana_hindcast/archive/norcpm1_assim-i1_TroPac_19800115
avgfile=ens_TREFHT_avg.nc
anofile=ens_TREFHT_ano.nc
var=TREFHT
test -f ${avgfile} && test -f ${anofile} && exit 0
## loop for every member
memdirs=$(ls -d ${dir}/*mem??)
i=0
for d in $memdirs; do
    ifiles=$(ls ${d}/atm/hist/*.cam*.h0.{1982..2010}-??.nc)
    i=$(($i+1))
    o1ts=TREFHT_$(printf "%2.2d" $i).nc
    test -f "$o1ts" && continue
    ncrcat -O -h -v ${var} ${ifiles} ${o1ts}
    #cdo -s           \
        #-select,name=${var} \
        #${ifiles} \
        #${o1ts}
done
## ensemble variables
cdo -ensavg TREFHT_??.nc ${avgfile}
cdo -ymonsub ${avgfile} -ymonmean ${avgfile} ${anofile}
## cat to 1 file
#cdo --reduce_dim -O -merge avg_ts.nc std_ts.nc min_ts.nc max_ts.nc ${outfile}