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

cases='TBI_NUG_Ano_detrend TBI_NUG_FF_detrend TBI_N2_Ano_detrend TBI_N2_FF_detrend'
for s in '_' '_MAM_' '_JJA_'; do
for p in 1980-2000 2000-2020 1980-2020; do
for i in $cases ;do
    ii=$(echo $i | sed 's/TBI_//')
    fig="21_prec/cor_ens_${ii}_PRECT-GPCP${s}${p}.png"
    src="../../${i}/${fig}"
    dst="cor_ens_${i}_PRECT-GPCP${s}${p}.png"
    tsrc=$(echo $src | sed -e's/.png/_thumb.png/')  
    tdst=$(echo $dst | sed -e's/.png/_thumb.png/')  
    cp -v "$src" "./$dst"
    cp "$tsrc" "./$tdst"
done
done
done
ln -sf ../../TBI_NUG_FF_detrend/21_prec/gpcp_prec_atm_grid_1980-2020.nc ./  ## for lower resolution
cdo -O -selyear,1980/2000  "gpcp_prec_atm_grid_1980-2020.nc" "gpcp_prec_atm_grid_1980-2000.nc"
cdo -O -selyear,2000/2020  "gpcp_prec_atm_grid_1980-2020.nc" "gpcp_prec_atm_grid_2000-2020.nc"