#!/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
#
#
outfile=era5_t2m_cor_ANN.nc
modfile=ens_TREFHT_ano.nc
var=t2m
outDJF=${outfile%_ANN.*}_DJF.nc
outMAM=${outfile%_ANN.*}_MAM.nc
outJJA=${outfile%_ANN.*}_JJA.nc
outSON=${outfile%_ANN.*}_SON.nc
test -f ${outfile} && exit 0
## loop for every member
ifiles=$(ls /cluster/shared/noresm/norcpm/Obs/AIR_TEM/ERA5/T2M/t2m.{1982..2018}.nc)
cdo -O    -b F32         \
    -select,name=${var} \
    ${ifiles} \
    ERA5_t2m_1982-2018.nc
    
# correlation, need regrid
cdo  -O  \
    -ymonsub ERA5_t2m_1982-2018.nc -ymonmean ERA5_t2m_1982-2018.nc  \
    ERA5_t2m_ano_1982-2018.nc
## Somehow remap is easily crash, should be run along
cdo  -O  \
    -remapbic,${modfile} \
    ERA5_t2m_ano_1982-2018.nc \
    tmp2.nc
cdo  -O -timcor    \
    -selvar,t2m tmp2.nc    \
    -select,name=TREFHT ${modfile} \
    ${outfile}
cdo  -O -timcor    \
    -selseason,DJF -selvar,t2m tmp2.nc    \
    -selseason,DJF -select,name=TREFHT ${modfile} \
    ${outDJF}
cdo  -O -timcor    \
    -selseason,MAM -selvar,t2m tmp2.nc    \
    -selseason,MAM -select,name=TREFHT ${modfile} \
    ${outMAM}
cdo  -O -timcor    \
    -selseason,JJA -selvar,t2m tmp2.nc    \
    -selseason,JJA -select,name=TREFHT ${modfile} \
    ${outJJA}
cdo  -O -timcor    \
    -selseason,SON -selvar,t2m tmp2.nc    \
    -selseason,SON -select,name=TREFHT ${modfile} \
    ${outSON}
rm -f ERA5_t2m_ano_1982-2018.nc tmp2.nc