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

## ccipc2 env
module purge

dir=/trd-project4/NS9039K/shared/pgchiu/archive/norcpm_ana_f09_tn14/
var=temp
lonlatbox=-150,-90,-5,5
level=1

gridfile=/nird/home/pgchiu/scratch/norcpm_verify/grid_tnx1v4_20170622.nc
griddes=griddes.txt
## make griddes first
#sed -i -e'1,5d' $griddes

outfile="nino34_temp_ana_ts.nc"
if [ ! -f "$outfile" ]; then
    ifiles=$(ls $dir/RESULT/{1980..2020}_{01..12}//analysis_avg_1.nc)
    cdo              \
        -fldmean  \
        -sellonlatbox,${lonlatbox}  \
        -setgrid,$griddes \
        -select,name=${var},level=${level} \
        ${ifiles} \
        ${outfile} 
fi

outfile="nino34_temp_fca_ts.nc"
if [ ! -f "$outfile" ]; then
    ifiles=$(ls $dir/RESULT/{1980..2020}_{01..12}//forecast_avg_1.nc)
    cdo              \
        -fldmean  \
        -sellonlatbox,${lonlatbox}  \
        -setgrid,$griddes \
        -select,name=${var},level=${level} \
        ${ifiles} \
        ${outfile} 2>/dev/null
fi