job_ablvint.sbatch

#!/bin/bash
#SBATCH --job-name=ablvint
#SBATCH --output=JOBS/%x-%j.out 
#SBATCH --qos=nf
#SBATCH --time=01:00:00
#SBATCH --gres=ssdtmp:15G


module purge
module load prgenv/gnu gcc/12.2.0 hpcx-openmpi/may23 hdf5-parallel/may23 netcdf4-parallel/may23 xios/trunk/r2515_gcc-xmpi.may23
module list

set +xv

if [ -n "$1" ]; then date=$1; else exit 1; fi
if [ -n "$2" ]; then var_name=$2; fi

p2D="${SCRATCH_DIR}/NC_${GRID/N/F}/${date:0:6}/surf2D_L1_${BASE_NAME}_${date}.nc"
p3D="${SCRATCH_DIR}/NC_${GRID/N/F}/${date:0:6}/inst3D_${VLEVEL}_${BASE_NAME}_${date}.nc"
pmsk="/ec/res5/project/mercator/ATMO_FORCINGS/ERA5_ABL50_F320_v1/mask0_${BASE_NAME%_ana*}_NNint_20120101.nc"
fabl=$(basename ${p3D/inst3D/abl}) ; fabl=${fabl/${VLEVEL}/L50Z10}

time cp -av ${p2D} ${p3D} ${pmsk} $TMPDIR/.
if [[ -n "${var_name}" && ${var_name} =~ "hpg" ]] || [ -z "${var_name}" ]; then time cp -av ${p3D/inst3D/hpg} $TMPDIR/.; fi
cp -v ${DIR_NEMO}/main_vinterp.exe $TMPDIR/.

sed -e "s/fsurf2D/$(basename ${p2D})/g; s/finst3D/$(basename ${p3D})/g; s/fmask/$(basename ${pmsk})/g; s/fhpg/$(basename ${p3D/inst3D/hpg})/g; s/fabl/${fabl}/g; s~fdir~$TMPDIR~g" < ${DIR_NEMO}/namelist_abl_tools.base > $TMPDIR/namelist_abl_tools.vint.$date
sed -i "s/ln_read_mask    =   .false./ln_read_mask    =    .true./g" $TMPDIR/namelist_abl_tools.vint.$date
if [ -n "${var_name}" ]; then
  sed -i "s/    var_name  = '',/    var_name  = '${var_name}',/g" $TMPDIR/namelist_abl_tools.vint.$date
  if [[ ! ${var_name} =~ "hpg" ]]; then sed -i "s/ln_hpg_frc      =   .true. /ln_hpg_frc      =   .false./g" $TMPDIR/namelist_abl_tools.vint.$date; fi
fi

ls -l $TMPDIR
cat $TMPDIR/namelist_abl_tools.vint.$date

$TMPDIR/main_vinterp.exe $TMPDIR/namelist_abl_tools.vint.$date

ls -l $TMPDIR
time mv -v $TMPDIR/${var_name:+"${var_name}_"}${fabl} $(dirname ${p3D})/.

exit 0