Newer
Older
#!/bin/bash
#SBATCH --job-name=ablvint
#SBATCH --output=JOBS/%x-%j.out
#SBATCH --qos=nf
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
if [ -n "$1" ]; then date=$1; else exit 1; 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})/.