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


dir=/projects/NS9039K/shared/norcpm/cases/NorCPM/norcpm-cmip6_pacemaker_18500115
var=saln
lonlatbox=120,200,10,30
level=1

gridfile=$(find $dir -name \*pop.h.\*|head -n1)
griddes=griddes.txt
## make griddes first
#cdo griddes2 $gridfile > $griddes
#sed -i -e'1,5d' $griddes

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

outfile="npac_saln_fca_ts.nc"
if [ ! -f "$outfile" ]; then
    ifiles=$(ls $dir/RESULT/{1982..2018}_{01..12}//forecast_avg_1.nc)
    cdo -s  -O       \
        -fldmean  \
        -setctomiss,nan \
        -sellonlatbox,${lonlatbox}  \
        -setgrid,$griddes \
        -select,name=${var},level=${level} \
        ${ifiles} \
        ${outfile}
fi