filename1=RENAMED_NORESM_NOASIA_3D.nc
filename2=RENAMED_NORESM_SP2005_3D.nc
#filename1=RENAMED_ECHAM_NOASIA_3D.nc
#filename2=RENAMED_ECHAM_SP2005_3D.nc
#filename_kernel=ECHAM6_CTRL_kernel_remap.nc
filename_kernel=ECHAM6_CTRL_kernel_remap.nc
tag=noresm_asia_warming
#tag=echam_asia_warming
#
rm brac_*
rm *_RENA*
rm K_*
rm Delta*
rm temp*
rm A_*
rm diff_K_*
rm gamma*
rm mu*
rm 1.nc
#
# Total values
#
cdo selname,tas $filename1 tas_$filename1
cdo selname,tas $filename2 tas_$filename2
cdo selname,rsdt $filename1 rsdt_$filename1
cdo selname,rsdt $filename2 rsdt_$filename2
cdo selname,rsut $filename1 rsut_$filename1
cdo selname,rsut $filename2 rsut_$filename2
cdo selname,rlut $filename1 rlut_$filename1
cdo selname,rlut $filename2 rlut_$filename2
cdo selname,rsds $filename1 rsds_$filename1
cdo selname,rsds $filename2 rsds_$filename2
cdo selname,rsus $filename1 rsus_$filename1
cdo selname,rsus $filename2 rsus_$filename2
cdo selname,rlds $filename1 rlds_$filename1
cdo selname,rlds $filename2 rlds_$filename2
cdo selname,rlus $filename1 rlus_$filename1
cdo selname,rlus $filename2 rlus_$filename2
cdo selname,hfls $filename1 hfls_$filename1
cdo selname,hfls $filename2 hfls_$filename2
cdo selname,hfss $filename1 hfss_$filename1
cdo selname,hfss $filename2 hfss_$filename2
#cdo selname,pr   $filename1 pr_$filename1
#cdo selname,pr   $filename2 pr_$filename2
#cdo selname,prsn $filename1 prsn_$filename1
#cdo selname,prsn $filename2 prsn_$filename2
cdo selname,ts   $filename1 ts_$filename1   #surface temp
cdo selname,ts   $filename2 ts_$filename2   
cdo selname,st   $filename1 st_$filename1   #3D atm. temp
cdo selname,st   $filename2 st_$filename2 
cdo selname,lnq  $filename1 lnq_$filename1   #3D log humidity
cdo selname,lnq  $filename2 lnq_$filename2 
#
# Clearsky values
#
cdo selname,rsutcs $filename1 rsutcs_$filename1
cdo selname,rsutcs $filename2 rsutcs_$filename2
cdo selname,rlutcs $filename1 rlutcs_$filename1
cdo selname,rlutcs $filename2 rlutcs_$filename2
cdo selname,rsdscs $filename1 rsdscs_$filename1
cdo selname,rsdscs $filename2 rsdscs_$filename2
cdo selname,rsuscs $filename1 rsuscs_$filename1
cdo selname,rsuscs $filename2 rsuscs_$filename2
cdo selname,clt $filename1 clt_$filename1
cdo selname,clt $filename2 clt_$filename2
#cdo selname,clt $filename1 clt100_$filename1
#cdo selname,clt $filename2 clt100_$filename2
#cdo divc,100  clt100_$filename1  clt_$filename1 #clt should be given as a fraction (0-1), but here it was in percentage
#cdo divc,100  clt100_$filename2  clt_$filename2 
#
# ECHAM6 Kernel from https://swiftbrowser.dkrz.de/public/dkrz_0c07783a-0bdc-4d5e-9f3b-c1b86fac060d/Radiative_kernels/
# (see Block & Mauritsen, JAMES 2013)
# Kernel is currently used only for LW components to estimate cloud LW
# ECHAM6 trad0 nad traf0 directions DO NOT need to be reversed
#
cdo selname,Ts_trad0 $filename_kernel K_Tslw.nc
cdo selname,Ts_traf0 $filename_kernel K_Tslw_cl.nc
cdo selname,Ta_trad0 $filename_kernel K_Talw.nc
cdo selname,Ta_traf0 $filename_kernel K_Talw_cl.nc
cdo selname,Q_trad0 $filename_kernel K_Qlw.nc
cdo selname,Q_traf0 $filename_kernel K_Qlw_cl.nc
#
# Overcast values for Taylor
#
cdo div tas_$filename1 tas_$filename1 1.nc                                       # unitary matrix
cdo sub 1.nc clt_$filename1 temp1.nc
cdo mul temp1.nc rsutcs_$filename1 temp2.nc
cdo sub rsut_$filename1 temp2.nc temp3.nc
cdo div temp3.nc clt_$filename1 rsutoc_$filename1                                #rsutoc1
cdo mul temp1.nc rlutcs_$filename1 temp2.nc
cdo sub rlut_$filename1 temp2.nc temp3.nc
cdo div temp3.nc clt_$filename1 rlutoc_$filename1                                #rlutoc1
cdo mul temp1.nc rsdscs_$filename1 temp2.nc
cdo sub rsds_$filename1 temp2.nc temp3.nc
cdo div temp3.nc clt_$filename1 rsdsoc_$filename1                                #rsdsoc1
cdo mul temp1.nc rsuscs_$filename1 temp2.nc
cdo sub rsus_$filename1 temp2.nc temp3.nc
cdo div temp3.nc clt_$filename1 rsusoc_$filename1                                #rsusoc1
cdo sub 1.nc clt_$filename2 temp1.nc
cdo mul temp1.nc rsutcs_$filename2 temp2.nc
cdo sub rsut_$filename2 temp2.nc temp3.nc
cdo div temp3.nc clt_$filename2 rsutoc_$filename2                                #rsutoc1
cdo mul temp1.nc rlutcs_$filename2 temp2.nc
cdo sub rlut_$filename2 temp2.nc temp3.nc
cdo div temp3.nc clt_$filename2 rlutoc_$filename2                                #rlutoc2
cdo mul temp1.nc rsdscs_$filename2 temp2.nc
cdo sub rsds_$filename2 temp2.nc temp3.nc
cdo div temp3.nc clt_$filename2 rsdsoc_$filename2                                #rsdsoc2
cdo mul temp1.nc rsuscs_$filename2 temp2.nc
cdo sub rsus_$filename2 temp2.nc temp3.nc
cdo div temp3.nc clt_$filename2 rsusoc_$filename2                                #rsusoc2
#
cdo sub tas_$filename1 tas_$filename2 temp1.nc
cdo chname,tas,DeltaT temp1.nc DeltaT_$tag.nc                                    #DeltaT
#
# Calculate Raisanen factors
#
cdo mul tas_$filename1 tas_$filename1 temp2.nc
cdo mul temp2.nc temp2.nc temp3a.nc
cdo mulc,5.670373E-8 temp3a.nc temp4.nc
cdo div rlut_$filename1 temp4.nc eps_$filename1                                  #eps1
cdo div rlutcs_$filename1 temp4.nc epscs_$filename1                              #epscs1
cdo mul tas_$filename2 tas_$filename2 temp2.nc
cdo mul temp2.nc temp2.nc temp3a.nc
cdo mulc,5.670373E-8 temp3a.nc temp4.nc
cdo div rlut_$filename2 temp4.nc eps_$filename2                                  #eps2
cdo div rlutcs_$filename2 temp4.nc epscs_$filename2                              #epscs2
#
cdo add tas_$filename1 tas_$filename2 temp1.nc
cdo mulc,0.5 temp1.nc brac_tas_$tag.nc                                           #brac_tas
cdo add eps_$filename1 eps_$filename2 temp1.nc
cdo mulc,0.5 temp1.nc brac_eps_$tag.nc                                           #brac_eps
cdo mul brac_tas_$tag.nc brac_tas_$tag.nc temp1.nc
cdo mul brac_tas_$tag.nc temp1.nc temp2.nc
cdo mulc,5.670373E-8 temp2.nc temp3.nc
cdo mulc,4.0 temp3.nc temp4.nc
cdo mul temp4.nc brac_eps_$tag.nc D_$tag.nc                                      #D
#
# Raisanen LW components
#
cdo mul tas_$filename1 tas_$filename1 temp1.nc
cdo mul temp1.nc temp1.nc  temp2.nc
cdo mul tas_$filename2 tas_$filename2 temp1.nc
cdo mul temp1.nc temp1.nc  temp3.nc
cdo add temp2.nc temp3.nc temp4.nc
cdo mulc,0.5 temp4.nc brac_taspower4_$tag.nc                                     #brac_taspower4
cdo sub eps_$filename1 eps_$filename2 Delta_eps_$tag.nc                          #Delta_eps
cdo mulc,-5.670373E-8 Delta_eps_$tag.nc temp1.nc
cdo mul  temp1.nc brac_taspower4_$tag.nc temp2.nc
cdo div temp2.nc D_$tag.nc temp3.nc
cdo chname,rlut,DeltaT_LW temp3.nc DeltaT_LW_$tag.nc                             #DeltaT_LW
rm temp?.nc
cdo sub epscs_$filename1 epscs_$filename2 Delta_epscs_$tag.nc                    #Delta_epscs 
cdo mulc,-5.670373E-8 Delta_epscs_$tag.nc temp1.nc
cdo mul  temp1.nc brac_taspower4_$tag.nc temp2.nc
cdo div temp2.nc D_$tag.nc temp3.nc
cdo chname,rlutcs,DeltaT_LW_cs temp3.nc DeltaT_LWcs_$tag.nc                      #DeltaT_LWcs
cdo sub DeltaT_LW_$tag.nc DeltaT_LWcs_$tag.nc temp4.nc
cdo chname,DeltaT_LW,DeltaT_LW_cre temp4.nc DeltaT_LWcre_$tag.nc                      #DeltaT_LWcre
#
# Kernel LW
#
cdo sub K_Tslw.nc K_Tslw_cl.nc diff_K_Tslw1.nc
cdo sub K_Talw.nc K_Talw_cl.nc diff_K_Talw1.nc
cdo sub K_Qlw.nc K_Qlw_cl.nc diff_K_Qlw1.nc
for year in $(seq -w 40 1 100); do cdo setyear,$year diff_K_Tslw1.nc exp_diff_K_Tslw_${year}.nc ; done 
cdo mergetime exp_diff_K_Tslw_* diff_K_Tslw.nc
rm exp_diff_K_Tslw_*
for year in $(seq -w 40 1 100); do cdo setyear,$year diff_K_Talw1.nc exp_diff_K_Talw_${year}.nc ; done 
cdo mergetime exp_diff_K_Talw_* diff_K_Talw.nc
rm exp_diff_K_Talw_*
for year in $(seq -w 40 1 100); do cdo setyear,$year diff_K_Qlw1.nc exp_diff_K_Qlw_${year}.nc ; done 
cdo mergetime exp_diff_K_Qlw_* diff_K_Qlw.nc
rm exp_diff_K_Qlw_*
cdo sub ts_$filename1 ts_$filename2 Delta_ts_$tag.nc                             #Delta_ts 
#cp DeltaT_$tag.nc Delta_ts_$tag.nc
cdo sub st_$filename1 st_$filename2 Delta_st_$tag.nc                             #Delta_st
cdo sub lnq_$filename1 lnq_$filename2 Delta_lnq_$tag.nc                           #Delta_lnq
cdo mul diff_K_Tslw.nc Delta_ts_$tag.nc Delta_LWts_$tag.nc                         #Delta_LWts
cdo mul diff_K_Talw.nc Delta_st_$tag.nc temp_LWta_$tag.nc
cdo vertsum temp_LWta_$tag.nc Delta_LWta_$tag.nc                                 #Delta_LWta
cdo mul diff_K_Qlw.nc Delta_lnq_$tag.nc temp_LWq_$tag.nc
cdo vertsum temp_LWq_$tag.nc Delta_LWq_$tag.nc                                   #Delta_LWq
cdo add Delta_LWts_$tag.nc Delta_LWta_$tag.nc temp1.nc
cdo add temp1.nc Delta_LWq_$tag.nc Delta_LW_cor_$tag.nc                          #Delta_LW_cor
cdo div Delta_LW_cor_$tag.nc D_$tag.nc temp3.nc 
cdo chname,Ts_trad0,DeltaT_LW_cor temp3.nc DeltaT_LW_cor_$tag.nc                 #DeltaT_LW_cor
cdo sub DeltaT_LWcre_$tag.nc DeltaT_LW_cor_$tag.nc temp4.nc
cdo chname,DeltaT_LW_cre,DeltaT_LW_cld temp4.nc DeltaT_LW_cld_$tag.nc           #DeltaT_LW_cld 
cdo sub DeltaT_LW_$tag.nc DeltaT_LW_cld_$tag.nc temp5.nc
cdo chname,DeltaT_LW,DeltaT_LW_clr temp5.nc  DeltaT_LW_clr_$tag.nc            #DeltaT_LW_clr
cdo div Delta_LWts_$tag.nc D_$tag.nc temp1.nc
cdo chname,Ts_trad0,DeltaT_LWts temp1.nc DeltaT_LWts_$tag.nc
cdo div Delta_LWta_$tag.nc D_$tag.nc temp2.nc
cdo chname,Ta_trad0,DeltaT_LWta temp2.nc DeltaT_LWta_$tag.nc
cdo div Delta_LWq_$tag.nc D_$tag.nc temp3.nc
cdo chname,Q_trad0,DeltaT_LWq temp3.nc DeltaT_LWq_$tag.nc
#
cdo timmean DeltaT_LWts_$tag.nc DeltaT_LWts_yearmean_$tag.nc
cdo timmean DeltaT_LWta_$tag.nc DeltaT_LWta_yearmean_$tag.nc
cdo timmean DeltaT_LWq_$tag.nc DeltaT_LWq_yearmean_$tag.nc
cdo timmean DeltaT_LW_cor_$tag.nc DeltaT_LW_cor_yearmean_$tag.nc
cdo timmean DeltaT_LW_clr_$tag.nc DeltaT_LW_clr_yearmean_$tag.nc
cdo timmean DeltaT_LW_cld_$tag.nc DeltaT_LW_cld_yearmean_$tag.nc
#
#
#
# Raisanen SW without decomposition
#
cdo sub rsdt_$filename1 rsut_$filename1 temp1.nc
cdo sub rsdt_$filename2 rsut_$filename2 temp2.nc
cp temp1.nc NETTOASW_$filename1
cp temp2.nc NETTOASW_$filename2
cdo sub temp1.nc temp2.nc temp3.nc
cdo div temp3.nc D_$tag.nc temp1.nc
cdo chname,rsdt,DeltaT_SW temp1.nc DeltaT_SW_$tag.nc                             #DeltaT_SW
cdo sub rsdt_$filename1 rsutcs_$filename1 temp1.nc
cdo sub rsdt_$filename2 rsutcs_$filename2 temp2.nc
cp temp1.nc NETTOASWcs_$filename1
cp temp2.nc NETTOASWcs_$filename2
cdo sub temp1.nc temp2.nc temp3.nc
cdo div temp3.nc D_$tag.nc temp1.nc
cdo chname,rsdt,DeltaT_SW_cs temp1.nc DeltaT_SWcs_$tag.nc                        #DeltaT_SWcs
cdo sub DeltaT_SW_$tag.nc DeltaT_SWcs_$tag.nc temp2.nc
cdo chname,DeltaT_SW,DeltaT_SW_cre temp2.nc DeltaT_SWcre_$tag.nc                #DeltaT_SWcre
#
# Raisanen SURF
#
cdo add hfls_$filename1 hfss_$filename1 temp1.nc
cdo sub rsus_$filename1 rsds_$filename1 temp2.nc
cdo sub rlus_$filename1 rlds_$filename1 temp3.nc
cdo add temp2.nc temp3.nc temp4.nc
cdo add temp1.nc temp4.nc temp5.nc
cp  temp5.nc NETSURF_$filename1
cdo add hfls_$filename2 hfss_$filename2 temp1.nc
cdo sub rsus_$filename2 rsds_$filename2 temp2.nc
cdo sub rlus_$filename2 rlds_$filename2 temp3.nc
cdo add temp2.nc temp3.nc temp4.nc
cdo add temp1.nc temp4.nc temp6.nc
cp  temp6.nc NETSURF_$filename2
cdo sub temp5.nc temp6.nc temp7.nc
cdo div temp7.nc D_$tag.nc temp8.nc
cdo chname,hfls,DeltaT_SURF temp8.nc DeltaT_SURF_$tag.nc                         #DeltaT_SURF
#
# NETTOA (outgoing);  here rlut (LW up) equals net LW at TOA
cdo sub rlut_$filename1 NETTOASW_$filename1 NETTOA_$filename1 # This is now net outgoing
cdo sub rlut_$filename2 NETTOASW_$filename2 NETTOA_$filename2
#
# Convergent (heat transport) components. First calculate as Russotto and Ackermann (2018).
# Later on, calculate check values from residual DeltaT
#
# NETTOA up, NETSURF up
# MSE:Moist static energy;LE latent energy;DSE dry static energy transport
#
cdo sub  NETTOA_$filename1 NETSURF_$filename1 MSE_$filename1 #MSE now positive if the column gains heat
#cdo mulc,2257000 pr_$filename1   hfls_vap_$filename1
#cdo mulc,333550  prsn_$filename1 hfls_snow_$filename1
#cdo add hfls_vap_$filename1 hfls_snow_$filename1 hfls_MSE_$filename1 
#cdo sub hfls_$filename1 hfls_MSE_$filename1 LE_$filename1
#cdo sub MSE_$filename1 LE_$filename1 DSE_$filename1

cdo sub NETTOA_$filename2 NETSURF_$filename2 MSE_$filename2
#cdo mulc,2257000 pr_$filename2   hfls_vap_$filename2
#cdo mulc,333550  prsn_$filename2 hfls_snow_$filename2
#cdo add hfls_vap_$filename2 hfls_snow_$filename2 hfls_MSE_$filename2 
#cdo sub hfls_$filename2 hfls_MSE_$filename2 LE_$filename2
#cdo sub MSE_$filename2 LE_$filename2 DSE_$filename2
cdo sub MSE_$filename1 MSE_$filename2 temp1.nc
cdo div temp1.nc D_$tag.nc temp2.nc
cdo chname,rlut,DeltaT_MSE temp2.nc DeltaT_MSE_$tag.nc                           #DeltaT_MSE
#cdo sub LE_$filename1 LE_$filename2 temp1.nc
#cdo div temp1.nc D_$tag.nc temp2.nc 
#cdo chname,hfls,DeltaT_LE temp2.nc DeltaT_LE_$tag.nc                             #DeltaT_LE
#cdo sub DSE_$filename1 DSE_$filename2 temp1.nc
#cdo div temp1.nc D_$tag.nc temp2.nc
#cdo chname,rlut,DeltaT_DSE temp2.nc DeltaT_DSE_$tag.nc                           #DeltaT_DSE
#
# Check MSE (=CONV)
cdo add DeltaT_LW_$tag.nc DeltaT_SW_$tag.nc temp1.nc 
cdo add DeltaT_SURF_$tag.nc temp1.nc temp2.nc
cdo sub DeltaT_$tag.nc temp2.nc temp3.nc
cdo chname,DeltaT,DeltaT_CONV temp3.nc DeltaT_CONV_$tag.nc 
#
# Check DeltaT
cdo add DeltaT_MSE_$tag.nc DeltaT_SURF_$tag.nc temp1.nc
cdo add DeltaT_SW_$tag.nc DeltaT_LW_$tag.nc temp2.nc
cdo add temp1.nc temp2.nc temp3.nc
cdo chname,DeltaT_MSE,DeltaT_check temp3.nc DeltaT_check_$tag.nc                 #DeltaT_check
cdo timmean DeltaT_check_$tag.nc DeltaT_check_yearmean_$tag.nc
cdo timmean DeltaT_$tag.nc DeltaT_yearmean_$tag.nc
#