#!/bin/bash -e

#run="f09_g16.B1850C5.TaiESM.exp2.7rhc"
#run="f09_g16.B1850C5.TaiESM.exp2.1nrhc.snap.testv2.tn15"
#un="f09_g16.CAM5.B1850-2000"
#run="f09.B_1850-2005.TAI.tn9"
#un="f09.B-hist.Tai.tn15"

run="historical"
isy=1900
iey=2014

#run="amip"
#isy=1979
#iey=2014

#run="piControl"
#isy=451
#iey=700

#run="ssp370"
#isy=2015
#iey=2100

#isy=20
#iey=50
export ISYY=$isy
export IEYY=$iey

#for run in  "historical" ##"amip" ##"hist-noSNAP" "hist-noTFZM" "hist-no3DRT" "hist-noGTSC"
#for run in "hist-noSNAP" "hist-noTFZM" "hist-no3DRT" "hist-noGTSC"
#for run in "hist-SNAP" "hist-TFZM" "hist-3DRT" "hist-GTSC"
#for run in "historical" "hist-aer" "hist-GHG" "hist-nat"
#for run in "hist-GHG" "hist-nat"
#for run in "hist-CAM5"
#for run in "control-2000"
#for run in "historical"
#do
export MODE=$run
#ndir="/lfs/archive/taiesm_ppost/$run/v20220321/"
#ndir="/lfs/archive/taiesm_ppost/$run/v20220329/"
#ndir="/lfs/archive/taiesm_ppost/$run/v20220824/"
#ndir="/lfs/archive/taiesm_ppost/control-2000/v20240911/r1i1p1f1/"  ## HALO
#ndir=" /lfs/archive/taiesm_ppost/control-2000/v20250608/r1i1p1f1/" ## HAHO
ndir=" /nird/datapeak/NS9034K/CMIP6/CMIP/NCC/NorESM2-MM/historical/r1i1p1f1/" ## NorESM
ndir=" /nird/datapeak/NS9034K/CMIP6/CMIP/NCC/NorESM2-MM/historical/r1i1p1f1/Amon/*/gn/latest/"
model="NorESM2-MM"
export model

#ndir="/archive/CMIP6/TaiESM1/$run/atmos/mon/r1i1p1f1/"
export path=$ndir

if [ ! -d $run ]; then
 mkdir $run
fi
#ncl output_cloud_W_lev_clm.ncl
#ncl output_ssh_ocn-remp.ncl
#exit

# ENSO index and case

ncl stp1_nino34_for_ERSST.ncl
ncl stp2_Dag_nino34_find-case.ncl

#ncl plot_nino34_STD_xy.ncl
#ncl Disp_TS_X-time_comp.ncl
#ncl plot_lag_corr_with_nino34.ncl

# Regression 
#ncl plot_nino34_reg_rsds.ncl
ncl plot_nino34_reg_ts_uv10_prec.ncl
  
#ncl plot_reg_SSH_taux_nino34_area.ncl
#ncl plot_reg_SSH_taux_nino4_area.ncl

#ncl plot_nino34_reg_taux.ncl

#ncl plot_nino34_duration.ncl
#ncl plot_Scatter_SSH_nino34.ncl
#ncl plot_SODA_SSH_EOF.ncl
#exit

# Case composties X-time

#ncl Disp_SSH_taux_X-time_comp.ncl
#ncl Disp_PRECT_U1000_X-time_comp.ncl
ncl Disp_TS_rsds_X-time_comp.ncl
#ncl Disp_TS_hfls_X-time_comp.ncl
#ncl Disp_TS_NET-flux_X-time_comp.ncl
#exit
#ncl Disp_TS_U1000_X-time_comp.ncl
#ncl Disp_Pr_rsds_X-time_comp.ncl

# Case composties cloud 

#ncl plot_cloud_lev_comp.ncl
#ncl plot_cloud_sig-lev_comp.ncl
#ncl plot_cloud_lev_comp_test-lama-file.ncl
#ncl plot_cloud_lev_comp_fun-contour-vc.ncl
    
# Case composties map 
#ncl plot_Ts_PSL_nino_case_cp.ncl
#ncl plot_PRECT_psi_nino_case_cp.ncl
#ncl output_Ts_Pr_nino_case_cp_zonal.ncl

# Case composties ocean x-Z

#cl plot_ocean_temp_lon_lev.ncl
#ncl plot_ocean_U_lon_lev.ncl
#ncl plot_ocean_temp_lon_lev_clm_diff.ncl

# mean state 
#ncl plot_ts_uv10_prec_diff.ncl
#ncl plot_Scatter_TS_W500.ncl
#ncl plot_Scatter_U850_W850.ncl
#ncl plot_Scatter_TS_tauu.ncl
#ncl plot_Scatter_tauu_TS.ncl
#ncl plot_Scatter_HMXL_tauu.ncl
ncl plot_sst_STD_map.ncl
#mv TaiESM_${run}*.* ${run}/
exit
#mv ${run} ../
#done