## Programe daily output of RegCM
## contact TrinhTuanLong, longtt(trinhtuanlong@gmail.com)
## HUS, Sep 2013
## requirement: cdo, nco, sigm2p (for Reg 4.3.5), and the run_post_RegCM4.sh script written by Thanh Ngo-Duc
###############################################
# NOTE please select directory of output and experiment number
lk='/work/users/longtt/SEACLID/'
exp='15'
dir=${lk}output_${exp}/
ename='sea_erain_vn_'${exp}'_'
# please select to convert dd (daily), mm (monthly), both (daily and monthly)
label='both' 
# please select years and months
year='1989 1990 1991 1992 1993 1994 1995 1996 1997 1998 1999 2000 2001 2002 2003 2004 2005 2006 2007 2008'
month='01 02 03 04 05 06 07 08 09 10 11 12'
# for first test
#year='1989'
#month=01
################################################
mkdir -p tmp_$exp
mkdir -p DAILY_$exp
mkdir -p MONTHLY_$exp
for y in $year ; do
   for m in $month ; do
	# 1. sigma to p (ATM,RAD)
	./sigma2p ${dir}/${ename}ATM.${y}${m}0100.nc 
	./sigma2p ${dir}/${ename}RAD.${y}${m}0100.nc 
	mv ${ename}ATM.${y}${m}0100_pressure.nc tmp_$exp/
	mv ${ename}RAD.${y}${m}0100_pressure.nc tmp_$exp/
	#
	if [[ $label == dd || $label == both ]]
	then
		# 2. convert to daily
		cdo daymean tmp_$exp/${ename}ATM.${y}${m}0100_pressure.nc tmp_$exp/ATM.daily.${y}${m}.nc
		cdo daymean tmp_$exp/${ename}RAD.${y}${m}0100_pressure.nc tmp_$exp/RAD.daily.${y}${m}.nc
		cdo daymean ${dir}/${ename}SRF.${y}${m}0100.nc tmp_$exp/SRF.daily.${y}${m}.nc
		#  add lat,lon for grads 
		./2run_post_RegCM4.sh tmp_$exp/ATM.daily.${y}${m}.nc DAILY_$exp/${ename}dd_ATM.${y}${m}.nc
		./2run_post_RegCM4.sh tmp_$exp/RAD.daily.${y}${m}.nc DAILY_$exp/${ename}dd_RAD.${y}${m}.nc
		./2run_post_RegCM4.sh tmp_$exp/SRF.daily.${y}${m}.nc DAILY_$exp/${ename}dd_SRF.${y}${m}.nc
		./run_post_RegCM4.sh ${dir}/${ename}STS.${y}${m}0100.nc DAILY_$exp/${ename}dd_STS.${y}${m}.nc
	fi
	if [[ $label == mm || $label == both ]]
	then
		# 3. convert to monthly
		cdo monmean tmp_$exp/${ename}ATM.${y}${m}0100_pressure.nc tmp_$exp/ATM.monthly.${y}${m}.nc
		cdo monmean tmp_$exp/${ename}RAD.${y}${m}0100_pressure.nc tmp_$exp/RAD.monthly.${y}${m}.nc
		cdo monmean ${dir}/${ename}SRF.${y}${m}0100.nc tmp_$exp/SRF.monthly.${y}${m}.nc
		#  add lat,lon for grads 
		./2run_post_RegCM4.sh tmp_$exp/ATM.monthly.${y}${m}.nc  MONTHLY_$exp/${ename}mm_ATM.${y}${m}.nc
		./2run_post_RegCM4.sh tmp_$exp/RAD.monthly.${y}${m}.nc  MONTHLY_$exp/${ename}mm_RAD.${y}${m}.nc
		./2run_post_RegCM4.sh tmp_$exp/SRF.monthly.${y}${m}.nc  MONTHLY_$exp/${ename}mm_SRF.${y}${m}.nc
		./run_post_RegCM4.sh ${dir}/${ename}STS.${y}${m}0100.nc MONTHLY_$exp/${ename}mm_STS.${y}${m}.nc
	fi
  done
rm -rf tmp_$exp/*
done
remote temple file
rm -rf tmp_$exp
############################################### 
# END