#!/bin/ksh
#=============================================================
#=============================================================
# Script to repeat ARIANE preliminary REF experiments 
# from Rühs et al. (2019)
#=============================================================
#=============================================================



# to be adjusted to computing environment (here: HLRN)
#=============================================================
#PBS -N ariane_EXPe1              
#PBS -j oe                        
#PBS -q specialm2_longq
#PBS -l nodes=1:ppn=24            
#PBS -l walltime=96:00:00         
#PBS -l feature=mpp2
#PBS -M sruehs@geomar.de
#PBS -m e                        
#=============================================================


# to be adjusted to the experiment of interest
#
# EXP: experiment name relating to release year
# LMINVAR: first time step of particle release period
# LMAXVAR: last time step of particle release period
#
# note: LMINVAR and LMAXVAR have to be specified with respect 
# to the whole set of time steps in the data specified for use
# in ARIANE via ariane_namelist_REFpreliminary.GENERIC
# note2: make sure to use an ARIANE version with active
# mod_criter1 restricting the integration length
# 
# example: for using the whole hindcast (1958-2009) 5-day mean 
# data, which is provided via annual files with subsequent 
# numbering (01-52) at /gfs1/work/shkifmsr/ARIANE/DATA, the
# particle release for the backward calculations starting in 
# 2004 would be at timesteps 3357-3430 
#=============================================================
PROJ_DIR=/gfs1/work/shkifmsr/ARIANE/ColdVsWarm
EXP=EXP2004
LMINVAR=3357
LMAXVAR=3430
export PATH=/home/h/shkifmsr/bin/ariane-2.2.6_01_ColdWarmRoute_40yrint/bin:${PATH}
#=============================================================


# main part
#=============================================================

cd ${PROJ_DIR}

mkdir -p ${EXP}
cd ${EXP}
LOG_FILE=${EXP}.log

date > ${LOG_FILE}

sed -e "s/START/$LMINVAR/" -e "s/STOP/$LMAXVAR/" ../ariane_namelist_REFpreliminary.GENERIC > namelist
cp -f -v ../ariane_sections_modified.txt ./sections.txt >> ${LOG_FILE}
cp -f -v ../ariane_region_limits . >> ${LOG_FILE}
cp -f -v ../ariane_segrid_copy . >> ${LOG_FILE}

which ariane >> ${LOG_FILE}

ariane 2>&1 >> ${LOG_FILE}

date >> ${LOG_FILE}

#=============================================================

exit