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

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


## make model climatology
ofn=model_TS_cli.nc
if [ ! -f $ofn ];then
    #filelist=$(ls /cluster/shared/NS9039K/norcpm_ana_hindcast/archive/norcpm1_assim-i1_TroPac_19800115/*/atm/hist/*.h0.{1982..2018}-??.nc | sort )
    #echo $filelist | ncrcat -v TS -o ${ofn}_tmp.nc
    #ls /cluster/shared/NS9039K/norcpm_ana_hindcast/archive/norcpm1_assim-i1_TroPac_19800115/*/atm/hist/*.h0.{1982..2018}-??.nc | sort | ncrcat -O -v TS -o ${ofn}_tmp.nc
    find  /cluster/shared/NS9039K/norcpm_ana_hindcast/archive/norcpm1_assim-i1_TroPac_19800115/*/atm/hist/ -name \*.h0.????-??.nc | grep -e{1982..2018}-...nc | sort | ncrcat -H -O -v TS -o ${ofn}_tmp.nc
    cdo -s -ymonmean -selname,TS ${ofn}_tmp.nc $ofn 
fi
## make HadISST mean climatology
ofn=obs_sst_cli.nc
if [ ! -f $ofn ];then
    ##filelist=$(ls ~/work/HadISST_2.1.0.0/*.nc)
    ##cdo -s -ymonmean -vertavg -select,name=sst $filelist $ofn 
    filelist=$(ls /cluster/shared/noresm/norcpm/Obs/SST/HADISST/HadISST_sst.nc)
    cdo -s -ymonmean -setattribute,sst@units="degC" -addc,273.15 -setmissval,-1000 -selyear,1982/2018 -select,name=sst $filelist $ofn 
fi
## diff
ofn=dif_sst.nc
if [ ! -f $ofn ];then
    cdo -sub -selvar,TS model_TS_cli.nc -remapbic,model_TS_cli.nc -selvar,sst obs_sst_cli.nc $ofn
fi