#!/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 a super ensemble mean among TBI pacemaker runs
## old plotCases:
##  cases1='norcpm1_tro10Atl5m_ng2l1dA_19800101 norcpm1_tro10Atl5m_ng2l1d_19800101'
##  cases2='norcpm2a_tro10Atl5m_ng2l1dA         norcpm2a_tro10Atl5m_ng2l1d'
cases1='TBI_NUG_Ano TBI_NUG_FF'
cases2='TBI_N2_Ano_w2 TBI_N2_FF_w2'
cases="$cases1 $cases2"

check_file_suf='sst_ensmean_ano_TBI_NUG_FF_19800101_fixed.nc sst_ensmean_ano_TBI_NUG_Ano_19800101_fixed.nc sst_ensmean_ano_TBI_N2_FF_w2_fixed.nc sst_ensmean_ano_TBI_N2_Ano_w2_fixed.nc'
skip=0
for i in $check_file_suf ;  do
    test -f $i && skip='1'
done
test $skip == '1' && exit

ocngridf1='/nird/projects/NS9039K/shared/pgchiu/diag_norcpm/grid_norcpm1_ocn.nc'
ocngridf2='/nird/projects/NS9039K/shared/pgchiu/diag_norcpm/grid_tnx1v4_20170622.nc'

## Make ensmean file for each run
dfpath='01_sst_regional-cor/sst_??.nc'
for cn in $cases; do 
    if [ ! -f sst_ensmean_${cn}.nc ] ; then
        ## ncea WARNING due to libnetcdf 1.8.1
        ncea -v sst ../../${cn}/$dfpath -o sst_ensmean_${cn}.nc
    fi
done
## Remove seasonal cycle with ensemble mean ( save to a file )
for cn in $cases; do 
    if [ ! -f sst_ensmean_ano_${cn}.nc ] ; then
        cdo ymonavg sst_ensmean_${cn}.nc sst_ensmean_cli_${cn}.nc
        cdo ymonsub sst_ensmean_${cn}.nc sst_ensmean_cli_${cn}.nc sst_ensmean_ano_${cn}.nc
    fi
done

## Interpolate to AMIP tos grid(fix grid 180x360)
## use 00_sst_cor/amip_tos.nc
## skip, use python to regrid

## Do ensemble mean ( with trend and without )
##    result two files, (t,y,x) SST anomaly
## skip, use python