#!/bin/bash
#
#SBATCH --job-name=ReBuildTS          # Job Name
#SBATCH --output=ReBuildTS.out        # standard output
#SBATCH --error=ReBuildTS.out         # error output
#SBATCH --ntasks=1                  # Number of MPI tasks
#SBATCH --time=06:00:00             # time (D-HH:MM)
#SBATCH --mem=20G
#

cd ${SLURM_SUBMIT_DIR}
module load python/meso-3.11

PERIOD="19890101_20131231"

MoveData(){
    pp="../${INST}/${EXP}/TS_DA"
    if [ ! -d ${pp} ] ; then
	/bin/mkdir -p ${pp}
	for v in ${var[@]} ; do
	    echo "Directory ${pp}, Variable ${v}"
	    ncrcat -v $v ./${EXP}.Daily_means.*.nc ${pp}/${MOD}_${EXP}_${PERIOD}_1D_${v}.nc
	    for rv in ${!renames[@]} ; do
    		if [ $v == $rv ] ; then
		    nv=${renames[$rv]}
		    echo "Renaming ${v} into ${nv}"
		    ncrename --variable ${v},${nv} ${pp}/${MOD}_${EXP}_${PERIOD}_1D_${v}.nc ${pp}/${MOD}_${EXP}_${PERIOD}_1D_${nv}.nc
		    /bin/rm -f ${pp}/${MOD}_${EXP}_${PERIOD}_1D_${v}.nc
		fi
	    done
	done
    fi
}



var=(rainf snowf Evap Soil_evap TVeg ECanop Qs Qsb SWnet LWnet Qle Qh RadT RadTmin RadTmax SoilWet LAI SnowDepth SubSnow snow_ice snow_liq PotEvap)

declare -A renames=( ["snowf"]="Snowf" ["rainf"]="Rainf")

INST="UKMO"
MOD="JULES"

EXP="ETHZ_Avg"
##/bin/rm -rf ${EXP}*2014*nc
##/bin/rm -rf ../UKMO/${EXP}/TS_DA
##/bin/rm -rf ../UKMO/${EXP}/TS_Catchment
##MoveData

EXP="WFDE5_CRU_GPCC"
##/bin/rm -rf ${EXP}*2014*nc
##/bin/rm -rf ../UKMO/${EXP}/TS_DA
##/bin/rm -rf ../UKMO/${EXP}/TS_Catchment
##MoveData

EXP="IPSL_Alt"
##/bin/rm -rf ${EXP}*2014*nc
##/bin/rm -rf ../UKMO/${EXP}/TS_DA
##/bin/rm -rf ../UKMO/${EXP}/TS_Catchment
MoveData

EXP="IPSL_Avg"
/bin/rm -rf ${EXP}*2014*nc
/bin/rm -rf ../UKMO/${EXP}/TS_DA
/bin/rm -rf ../UKMO/${EXP}/TS_Catchment
MoveData