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

vn=TGCLDLWP
ybe=1980/2014
var=CLWP

casefn=noresm2_lm-snow-rad02_${var}.nc
ctrlfn=CTRL_${var}.nc
if [ ! -f "$casefn" ] ;then
fns='/nird/datalake/NS9039K/users/pgchiu/n209_rm-snowrad02/atm/hist/*.h0.*'
cdo -select,name=${vn},year=${ybe} -cat "$fns" "$casefn"
fi

if [ ! -f "$ctrlfn" ] ;then
fns='/nird/datalake/NS9039K/projects/TRIATLAS/NorESM2_LM/atm/hist/NHIST_f19_tn14_20190710.cam.h0.*.nc'
cdo -select,name=${vn},year=${ybe} -cat "$fns" "$ctrlfn"
fi

## make avg
cdo -O --reduce_dim -timmean "$casefn" "mean_${casefn}"
cdo -O --reduce_dim -timmean "$ctrlfn" "mean_${ctrlfn}"

## make variance with seasonal cycle removed
cdo -O --reduce_dim -timstd -ymonsub "${casefn}" -ymonmean "${casefn}" "anovar_${casefn}"
cdo -O --reduce_dim -timstd -ymonsub "${ctrlfn}" -ymonmean "${ctrlfn}" "anovar_${ctrlfn}"