source: NEMO/branches/2021/dev_r14318_RK3_stage1_tsplit/tests/SEICHE/EXPREF/run_job.slurm @ 14597

#!/bin/bash
#
#SBATCH -J SEICHE 
#SBATCH -o SEICHE.%j.out
#SBATCH -e SEICHE.%j.err
#SBATCH -N 1 
#SBATCH -p mono
#SBATCH --ntasks-per-node:32
#SBATCH --time=00:20:00
#SBATCH --exclusive
#SBATCH -x node-07
#
#
OCEANCORES=30
XIOCORES=2
#
###############################################################
export MPIRUN="mpirun -np $OCEANCORES"
export RUN_DIR=/data/rd_exchange/jchanut/dev_r14318_RK3_stage1_tsplit/tests/SEICHE
export TOOL_DIR=/data/rd_exchange/jchanut/dev_r14318_RK3_stage1_tsplit/tools
#export RUN_DIR="$(dirname "$(pwd)")"
export EXE_DIR=/kara/data1/jchanut/SEICHE
###############################################################
#
# change to the working directory 
#
cd ${EXE_DIR}
cp ${RUN_DIR}/EXPREF/*.xml .
cp ${RUN_DIR}/EXPREF/namelist_ref .
cp ${RUN_DIR}/EXPREF/namelist_cfg_skel_marsaleix ./namelist_cfg_skel
cp ${RUN_DIR}/BLD/bin/nemo.exe .
cp ${XIOS_BIN}/xios_server.exe .
cp ${TOOL_DIR}/REBUILD_NEMO/rebuild_nemo .
cp ${TOOL_DIR}/REBUILD_NEMO/rebuild_nemo.exe .
#
#ntest=0: test time convergence ; =1: spatial convergence ; =2: both at constant courant number
#ntest=-1: is a sanity check with constant barotropic time step and varying baroclinic time step
ntest=-1
exp_base='SEICHE_mar_rk3_dt'
rdt_base=960
nbaro_base=112
nstep_base=1080
npts_base=1600
nlev_base=30
atfp_base=0.1
btav='.false.'
btfw='.true.'
ralpha_base=0.07 # = 0.5*theta
nexp=5
for (( i=1; i<=$nexp; i++))
do
   rdt=$rdt_base
   nstep=$nstep_base
   npts=$npts_base
   nlev=$nlev_base
   nbaro=$nbaro_base
   atfp=$atfp_base
   ralpha=$ralpha_base
   exp=$exp_base'_'$i
if   [[ $ntest -eq  0 || $ntest -eq  2 || $ntest -eq -1 ]] ; then
# Temporal convergence
   for (( j=1; j<$i; j++)) do rdt=` expr  $rdt / 2 `; done
   for (( j=1; j<$i; j++)) do nstep=` expr  $nstep \\* 2 `; done
elif [[ $ntest -eq  1 ]] ; then
# Spatial convergence only:
# smallest time step:
   for (( j=1; j<$nexp; j++)) do rdt=` expr  $rdt / 2 `; done
   for (( j=1; j<$nexp; j++)) do nstep=` expr  $nstep \\* 2 `; done
fi
if   [[ $ntest -eq  1 || $ntest -eq  2 ]] ; then
# With spatial and temporal convergence:
   for (( j=1; j<$i; j++)) do npts=` expr  $npts \\* 2 `; done
fi
if   [[ $ntest -eq  -1 ]] ; then
# This is a verification test: whatever the baroclinic time step is
# number of barotropic iterations is modified accordingly
# such that barotropic time step is constant. If no time averaging in the
# barotropic mode, the solution should basically be unchanged
   for (( j=1; j<$i; j++)) do nbaro=` expr  $nbaro / 2 `; done
fi
if   [[ $ntest -eq  3 ]] ; then
# With spatial and temporal convergence:
   for (( j=1; j<$i; j++)) do nlev=` expr  $nlev \\* 2 `; done
fi
   sed -e "s/EXPNAME/'$exp'/"  \
       -e "s/RDT/${rdt}/"      \
       -e "s/NPTS/${npts}/"    \
       -e "s/NLEV/${nlev}/"    \
       -e "s/NSTEP/${nstep}/"  \
       -e "s/NBARO/$nbaro/"    \
       -e "s/RALPHA/$ralpha/"  \
       -e "s/BTAV/$btav/"      \
       -e "s/BTFW/$btfw/"      \
       -e "s/ASS/$atfp/" namelist_cfg_skel > namelist_cfg 
   time mpirun -bynode -np $OCEANCORES ./nemo.exe : -np $XIOCORES ./xios_server.exe
   mv run.stat.nc ${i}_run.stat.nc
   mv run.stat ${i}_run.stat
   mv ocean.output ${i}_ocean.output
   stepout=`printf "%08d\n" $nstep`
   ./rebuild_nemo ${exp}_${stepout}_restart $OCEANCORES 
   rm -rf ${exp}_${stepout}_restart_*
   rm nam_rebuild*
   mv ${exp}_${stepout}_restart.nc ${exp}_restart.nc
done
#
exit