#!/bin/bash

# set variables

expin=itczmip008
expout=Land4xCO2

model=echam6-1

year1=0126
year2=0135

level='100000,92500,85000,70000,60000,50000,40000,30000,25000,20000,15000,10000,7000,5000,3000,2000,1000'


# merge BOT and ATM daily files
cdo -O mergetime BOT_dm_${expin}_012[6789] BOT_dm_${expin}_013[012345] BOT_dm_${expin}_${year1}_${year2}
cdo -O mergetime ATM_dm_${expin}_012[6789] ATM_dm_${expin}_013[012345] ATM_dm_${expin}_${year1}_${year2}

############## 2d data #####################################

# tas
cdo -r -f nc -t echam6 -setname,tas -selvar,temp2 BOT_dm_${expin}_${year1}_${year2} \
 	tas_Aday_${model}_${expout}_r1i1p1_${year1}01-${year2}12.nc
ncatted -O -a long_name,tas,o,c,air_temperature tas_Aday_${model}_${expout}_r1i1p1_${year1}01-${year2}12.nc

# ts
cdo -r -f nc -t echam6 -setname,ts -selvar,tsurf BOT_dm_${expin}_${year1}_${year2} \
 	ts_Aday_${model}_${expout}_r1i1p1_${year1}01-${year2}12.nc
ncatted -O -a long_name,ts,o,c,surface_temperature ts_Aday_${model}_${expout}_r1i1p1_${year1}01-${year2}12.nc

# tasmin
cdo -r -f nc -t echam6 -setname,tasmin -selvar,t2min BOT_dm_${expin}_${year1}_${year2} \
 	tasmin_Aday_${model}_${expout}_r1i1p1_${year1}01-${year2}12.nc
ncatted -O -a long_name,tasmin,o,c,air_temperature tasmin_Aday_${model}_${expout}_r1i1p1_${year1}01-${year2}12.nc

# tasmax
cdo -r -f nc -t echam6 -setname,tasmax -selvar,t2max BOT_dm_${expin}_${year1}_${year2} \
 	tasmax_Aday_${model}_${expout}_r1i1p1_${year1}01-${year2}12.nc
ncatted -O -a long_name,tasmax,o,c,air_temperature tasmax_Aday_${model}_${expout}_r1i1p1_${year1}01-${year2}12.nc

# psl
cdo -r -f nc -t echam6 -setname,psl -selvar,slp BOT_dm_${expin}_${year1}_${year2} \
 	psl_Aday_${model}_${expout}_r1i1p1_${year1}01-${year2}12.nc
ncatted -O -a long_name,psl,o,c,air_pressure_at_sea_level psl_Aday_${model}_${expout}_r1i1p1_${year1}01-${year2}12.nc

# ps
cdo -r -f nc -t echam6 -setname,ps -selvar,aps BOT_dm_${expin}_${year1}_${year2} \
 	ps_Aday_${model}_${expout}_r1i1p1_${year1}01-${year2}12.nc
ncatted -O -a long_name,ps,o,c,surface_air_pressure ps_Aday_${model}_${expout}_r1i1p1_${year1}01-${year2}12.nc

# uas
cdo -r -f nc -t echam6 -setname,uas -selvar,u10 BOT_dm_${expin}_${year1}_${year2} \
 	uas_Aday_${model}_${expout}_r1i1p1_${year1}01-${year2}12.nc
ncatted -O -a long_name,uas,o,c,eastward_wind uas_Aday_${model}_${expout}_r1i1p1_${year1}01-${year2}12.nc

# vas
cdo -r -f nc -t echam6 -setname,vas -selvar,v10 BOT_dm_${expin}_${year1}_${year2} \
 	vas_Aday_${model}_${expout}_r1i1p1_${year1}01-${year2}12.nc
ncatted -O -a long_name,vas,o,c,northward_wind vas_Aday_${model}_${expout}_r1i1p1_${year1}01-${year2}12.nc

# sfcWind
cdo -r -f nc -t echam6 -setname,sfcWind -selvar,wind10 BOT_dm_${expin}_${year1}_${year2} \
 	sfcWind_Aday_${model}_${expout}_r1i1p1_${year1}01-${year2}12.nc
ncatted -O -a long_name,sfcWind,o,c,wind_speed sfcWind_Aday_${model}_${expout}_r1i1p1_${year1}01-${year2}12.nc

# pr
cdo -r -f nc -t echam6 -setname,pr -selvar,precip BOT_dm_${expin}_${year1}_${year2} \
 	pr_Aday_${model}_${expout}_r1i1p1_${year1}01-${year2}12.nc
ncatted -O -a long_name,pr,o,c,precipitation_flux pr_Aday_${model}_${expout}_r1i1p1_${year1}01-${year2}12.nc

# prsn
cdo -r -f nc -t echam6 -setname,prsn -selvar,aprs BOT_dm_${expin}_${year1}_${year2} \
 	prsn_Aday_${model}_${expout}_r1i1p1_${year1}01-${year2}12.nc
ncatted -O -a long_name,prsn,o,c,snowfall_flux prsn_Aday_${model}_${expout}_r1i1p1_${year1}01-${year2}12.nc

# prc
cdo -r -f nc -t echam6 -setname,prc -selvar,aprc BOT_dm_${expin}_${year1}_${year2} \
 	prc_Aday_${model}_${expout}_r1i1p1_${year1}01-${year2}12.nc
ncatted -O -a long_name,prc,o,c,convective_precipitation_flux prc_Aday_${model}_${expout}_r1i1p1_${year1}01-${year2}12.nc

# evspsbl
cdo -r -f nc -t echam6 -mulc,-1 -setname,evspsbl -selvar,evap BOT_dm_${expin}_${year1}_${year2} \
 	evspsbl_Aday_${model}_${expout}_r1i1p1_${year1}01-${year2}12.nc
ncatted -O -a long_name,evspsbl,o,c,water_evaporation_flux evspsbl_Aday_${model}_${expout}_r1i1p1_${year1}01-${year2}12.nc

# tauu
cdo -r -f nc -t echam6 -setname,tauu -selvar,ustr BOT_dm_${expin}_${year1}_${year2} \
 	tauu_Aday_${model}_${expout}_r1i1p1_${year1}01-${year2}12.nc
ncatted -O -a long_name,tauu,o,c,surface_downward_eastward_stress tauu_Aday_${model}_${expout}_r1i1p1_${year1}01-${year2}12.nc

# tauv
cdo -r -f nc -t echam6 -setname,tauv -selvar,vstr BOT_dm_${expin}_${year1}_${year2} \
 	tauv_Aday_${model}_${expout}_r1i1p1_${year1}01-${year2}12.nc
ncatted -O -a long_name,tauv,o,c,surface_downward_northward_stress tauv_Aday_${model}_${expout}_r1i1p1_${year1}01-${year2}12.nc

# hfls
cdo -r -f nc -t echam6 -mulc,-1 -setname,hfls -selvar,ahfl BOT_dm_${expin}_${year1}_${year2} \
 	hfls_Aday_${model}_${expout}_r1i1p1_${year1}01-${year2}12.nc
ncatted -O -a long_name,hfls,o,c,surface_upward_latent_heat_flux hfls_Aday_${model}_${expout}_r1i1p1_${year1}01-${year2}12.nc

# hfss
cdo -r -f nc -t echam6 -mulc,-1 -setname,hfss -selvar,ahfs BOT_dm_${expin}_${year1}_${year2} \
 	hfss_Aday_${model}_${expout}_r1i1p1_${year1}01-${year2}12.nc
ncatted -O -a long_name,hfss,o,c,surface_upward_sensible_heat_flux hfss_Aday_${model}_${expout}_r1i1p1_${year1}01-${year2}12.nc

#rlds
cdo -r -f nc -t echam6 -setname,rlds -sub -selvar,trads BOT_dm_${expin}_${year1}_${year2} \
        -selvar,tradsu BOT_dm_${expin}_${year1}_${year2} rlds_Aday_${model}_${expout}_r1i1p1_${year1}01-${year2}12.nc
ncatted -O -a long_name,rlds,o,c,surface_downwelling_longwave_flux_in_air rlds_Aday_${model}_${expout}_r1i1p1_${year1}01-${year2}12.nc

#rlus
cdo -r -f nc -t echam6 -mulc,-1 -setname,rlus -selvar,tradsu BOT_dm_${expin}_${year1}_${year2} \
        rlus_Aday_${model}_${expout}_r1i1p1_${year1}01-${year2}12.nc
ncatted -O -a long_name,rlus,o,c,surface_upwelling_longwave_flux_in_air rlus_Aday_${model}_${expout}_r1i1p1_${year1}01-${year2}12.nc

#rsds
cdo -r -f nc -t echam6 -setname,rsds -sub -selvar,srads BOT_dm_${expin}_${year1}_${year2} \
        -selvar,sradsu BOT_dm_${expin}_${year1}_${year2} rsds_Aday_${model}_${expout}_r1i1p1_${year1}01-${year2}12.nc
ncatted -O -a long_name,rsds,o,c,surface_downwelling_shortwave_flux_in_air rsds_Aday_${model}_${expout}_r1i1p1_${year1}01-${year2}12.nc

#rsus
cdo -r -f nc -t echam6 -mulc,-1 -setname,rsus -selvar,sradsu BOT_dm_${expin}_${year1}_${year2} \
        rsus_Aday_${model}_${expout}_r1i1p1_${year1}01-${year2}12.nc
ncatted -O -a long_name,rsus,o,c,surface_upwelling_shortwave_flux_in_air rsus_Aday_${model}_${expout}_r1i1p1_${year1}01-${year2}12.nc

########### clear-sky sw surface fluxes (start) 
#need to calcule effective sfc albedo first
rm -f sfcalbedo
cdo div rsus_Aday_${model}_${expout}_r1i1p1_${year1}01-${year2}12.nc \
        rsds_Aday_${model}_${expout}_r1i1p1_${year1}01-${year2}12.nc \
        sfcalbedo
#rsdscs
cdo -r -f nc -t echam6 -setmisstoc,0 -setname,rsdscs -div -selvar,srafs BOT_dm_${expin}_${year1}_${year2} \
        -addc,1 -mulc,-1 sfcalbedo rsdscs_Aday_${model}_${expout}_r1i1p1_${year1}01-${year2}12.nc 
ncatted -O -a long_name,rsdscs,o,c,surface_downwelling_shortwave_flux_in_air_assuming_clear_sky rsdscs_Aday_${model}_${expout}_r1i1p1_${year1}01-${year2}12.nc

#rsuscs
cdo -r -f nc -t echam6 -mulc,-1 -setname,rsuscs -sub -selvar,srafs BOT_dm_${expin}_${year1}_${year2} \
        rsdscs_Aday_${model}_${expout}_r1i1p1_${year1}01-${year2}12.nc \
        rsuscs_Aday_${model}_${expout}_r1i1p1_${year1}01-${year2}12.nc
ncatted -O -a long_name,rsuscs,o,c,surface_upwelling_shortwave_flux_in_air_assuming_clear_sky rsuscs_Aday_${model}_${expout}_r1i1p1_${year1}01-${year2}12.nc

rm -f sfcalbedo
########### clear-sky sw surface fluxes (end)

#rldscs
cdo -r -f nc -t echam6 -setname,rldscs -sub -selvar,trafs BOT_dm_${expin}_${year1}_${year2} \
        -selvar,tradsu BOT_dm_${expin}_${year1}_${year2} rldscs_Aday_${model}_${expout}_r1i1p1_${year1}01-${year2}12.nc
ncatted -O -a long_name,rldscs,o,c,surface_downwelling_longwave_flux_in_air_assuming_clear_sky rldscs_Aday_${model}_${expout}_r1i1p1_${year1}01-${year2}12.nc

# rsdt 
cdo -r -f nc -t echam6 -setname,rsdt -selvar,srad0d BOT_dm_${expin}_${year1}_${year2} \
        rsdt_Aday_${model}_${expout}_r1i1p1_${year1}01-${year2}12.nc
ncatted -O -a long_name,rsdt,o,c,toa_incoming_shortwave_flux rsdt_Aday_${model}_${expout}_r1i1p1_${year1}01-${year2}12.nc

# rsut 
cdo -r -f nc -t echam6 -mulc,-1 -setname,rsut -selvar,srad0u BOT_dm_${expin}_${year1}_${year2} \
        rsut_Aday_${model}_${expout}_r1i1p1_${year1}01-${year2}12.nc
ncatted -O -a long_name,rsut,o,c,toa_outgoing_shortwave_flux rsut_Aday_${model}_${expout}_r1i1p1_${year1}01-${year2}12.nc

# rlut
cdo -r -f nc -t echam6 -mulc,-1 -setname,rlut -selvar,trad0 BOT_dm_${expin}_${year1}_${year2} \
        rlut_Aday_${model}_${expout}_r1i1p1_${year1}01-${year2}12.nc
ncatted -O -a long_name,rlut,o,c,toa_outgoing_longwave_flux rlut_Aday_${model}_${expout}_r1i1p1_${year1}01-${year2}12.nc

# rlutcs
cdo -r -f nc -t echam6 -mulc,-1 -setname,rlutcs -selvar,traf0 BOT_dm_${expin}_${year1}_${year2} \
        rlutcs_Aday_${model}_${expout}_r1i1p1_${year1}01-${year2}12.nc
ncatted -O -a long_name,rlutcs,o,c,toa_outgoing_longwave_flux_assuming_clear_sky rlutcs_Aday_${model}_${expout}_r1i1p1_${year1}01-${year2}12.nc

# rsutcs
cdo -r -f nc -t echam6 -mulc,-1 -setname,rsutcs -sub -selvar,sraf0 BOT_dm_${expin}_${year1}_${year2} \
        -selvar,srad0d BOT_dm_${expin}_${year1}_${year2} rsutcs_Aday_${model}_${expout}_r1i1p1_${year1}01-${year2}12.nc
ncatted -O -a long_name,rsutcs,o,c,toa_outgoing_shortwave_flux_assuming_clear_sky rsutcs_Aday_${model}_${expout}_r1i1p1_${year1}01-${year2}12.nc

# prw
cdo -r -f nc -t echam6 -setname,prw -selvar,qvi BOT_dm_${expin}_${year1}_${year2} \
 	prw_Aday_${model}_${expout}_r1i1p1_${year1}01-${year2}12.nc
ncatted -O -a long_name,prw,o,c,atmosphere_water_vapor_content prw_Aday_${model}_${expout}_r1i1p1_${year1}01-${year2}12.nc

# clt
cdo -r -f nc -t echam6 -mulc,100 -setname,clt -selvar,aclcov BOT_dm_${expin}_${year1}_${year2} \
 	clt_Aday_${model}_${expout}_r1i1p1_${year1}01-${year2}12.nc
ncatted -O -a long_name,clt,o,c,cloud_area_fraction clt_Aday_${model}_${expout}_r1i1p1_${year1}01-${year2}12.nc

# clwvi
cdo -r -f nc -t echam6 -setname,clwvi -selvar,xlvi BOT_dm_${expin}_${year1}_${year2} \
 	clwvi_Aday_${model}_${expout}_r1i1p1_${year1}01-${year2}12.nc
ncatted -O -a long_name,clwvi,o,c,atmosphere_cloud_condensed_water_content clwvi_Aday_${model}_${expout}_r1i1p1_${year1}01-${year2}12.nc

# clivi
cdo -r -f nc -t echam6 -setname,clivi -selvar,xivi BOT_dm_${expin}_${year1}_${year2} \
 	clivi_Aday_${model}_${expout}_r1i1p1_${year1}01-${year2}12.nc      
ncatted -O -a long_name,clivi,o,c,atmosphere_cloud_ice_content clivi_Aday_${model}_${expout}_r1i1p1_${year1}01-${year2}12.nc

############## 3d data #####################################

# cl
cdo -r -f nc -t echam6 -mulc,100 -setname,cl -selvar,aclcac -sellevel,${level} ATM_dm_${expin}_${year1}_${year2} \
        cl_Aday_${model}_${expout}_r1i1p1_${year1}01-${year2}12.nc
ncatted -O -a long_name,cl,o,c,cloud_area_fraction_in_atmosphere_layer cl_Aday_${model}_${expout}_r1i1p1_${year1}01-${year2}12.nc
ncrename -d lev,plev -v lev,plev cl_Aday_${model}_${expout}_r1i1p1_${year1}01-${year2}12.nc

# clw
cdo -r -f nc -t echam6 -mulc,100 -setname,clw -selvar,xl -sellevel,${level} ATM_dm_${expin}_${year1}_${year2} \
        clw_Aday_${model}_${expout}_r1i1p1_${year1}01-${year2}12.nc
ncatted -O -a long_name,clw,o,c,mass_fraction_of_cloud_liquid_water_in_air clw_Aday_${model}_${expout}_r1i1p1_${year1}01-${year2}12.nc
ncrename -d lev,plev -v lev,plev clw_Aday_${model}_${expout}_r1i1p1_${year1}01-${year2}12.nc

# cli
cdo -r -f nc -t echam6 -mulc,100 -setname,cli -selvar,xi -sellevel,${level} ATM_dm_${expin}_${year1}_${year2} \
        cli_Aday_${model}_${expout}_r1i1p1_${year1}01-${year2}12.nc
ncatted -O -a long_name,cli,o,c,mass_fraction_of_cloud_ice_in_air cli_Aday_${model}_${expout}_r1i1p1_${year1}01-${year2}12.nc
ncrename -d lev,plev -v lev,plev cli_Aday_${model}_${expout}_r1i1p1_${year1}01-${year2}12.nc

# ta
cdo -r -f nc -t echam6 -setname,ta -selvar,t -sellevel,${level} ATM_dm_${expin}_${year1}_${year2} \
        ta_Aday_${model}_${expout}_r1i1p1_${year1}01-${year2}12.nc
ncatted -O -a long_name,ta,o,c,air_temperature ta_Aday_${model}_${expout}_r1i1p1_${year1}01-${year2}12.nc
ncrename -d lev,plev -v lev,plev ta_Aday_${model}_${expout}_r1i1p1_${year1}01-${year2}12.nc

# ua
cdo -r -f nc -t echam6 -setname,ua -selvar,u -sellevel,${level} ATM_dm_${expin}_${year1}_${year2} \
        ua_Aday_${model}_${expout}_r1i1p1_${year1}01-${year2}12.nc
ncatted -O -a long_name,ua,o,c,eastward_wind ua_Aday_${model}_${expout}_r1i1p1_${year1}01-${year2}12.nc
ncrename -d lev,plev -v lev,plev ua_Aday_${model}_${expout}_r1i1p1_${year1}01-${year2}12.nc

# va
cdo -r -f nc -t echam6 -setname,va -selvar,v -sellevel,${level} ATM_dm_${expin}_${year1}_${year2} \
        va_Aday_${model}_${expout}_r1i1p1_${year1}01-${year2}12.nc
ncatted -O -a long_name,va,o,c,northward_wind va_Aday_${model}_${expout}_r1i1p1_${year1}01-${year2}12.nc
ncrename -d lev,plev -v lev,plev va_Aday_${model}_${expout}_r1i1p1_${year1}01-${year2}12.nc

# hus 
cdo -r -f nc -t echam6 -setname,hus -selvar,q -sellevel,${level} ATM_dm_${expin}_${year1}_${year2} \
        hus_Aday_${model}_${expout}_r1i1p1_${year1}01-${year2}12.nc
ncatted -O -a long_name,hus,o,c,specific_humidity hus_Aday_${model}_${expout}_r1i1p1_${year1}01-${year2}12.nc
ncrename -d lev,plev -v lev,plev hus_Aday_${model}_${expout}_r1i1p1_${year1}01-${year2}12.nc

# hur 
cdo -r -f nc -t echam6 -mulc,100 -setname,hur -selvar,rhumidity -sellevel,${level} ATM_dm_${expin}_${year1}_${year2} \
        hur_Aday_${model}_${expout}_r1i1p1_${year1}01-${year2}12.nc
ncatted -O -a long_name,hur,o,c,relative_humidity hur_Aday_${model}_${expout}_r1i1p1_${year1}01-${year2}12.nc
ncrename -d lev,plev -v lev,plev hur_Aday_${model}_${expout}_r1i1p1_${year1}01-${year2}12.nc

# wap 
cdo -r -f nc -t echam6 -setname,wap -selvar,omega -sellevel,${level} ATM_dm_${expin}_${year1}_${year2} \
        wap_Aday_${model}_${expout}_r1i1p1_${year1}01-${year2}12.nc
ncatted -O -a long_name,wap,o,c,lagrangian_tendency_or_air_pressure wap_Aday_${model}_${expout}_r1i1p1_${year1}01-${year2}12.nc
ncrename -d lev,plev -v lev,plev wap_Aday_${model}_${expout}_r1i1p1_${year1}01-${year2}12.nc

# zg 
cdo -r -f nc -t echam6 -setname,zg -selvar,geopoth -sellevel,${level} ATM_dm_${expin}_${year1}_${year2} \
        zg_Aday_${model}_${expout}_r1i1p1_${year1}01-${year2}12.nc
ncatted -O -a long_name,zg,o,c,geopotential_height zg_Aday_${model}_${expout}_r1i1p1_${year1}01-${year2}12.nc
ncrename -d lev,plev -v lev,plev zg_Aday_${model}_${expout}_r1i1p1_${year1}01-${year2}12.nc