#!/usr/bin/env bash
# Wenchang Yang (wenchang@princeton.edu)
# Wed Jun 19 14:19:49 EDT 2024
##SBATCH --nodes=1                # node count
##SBATCH --ntasks-per-node=1      # number of tasks per node
# 
#SBATCH --ntasks=1               # total number of tasks across all nodes = nodes x ntasks-per-node
#SBATCH --cpus-per-task=1        # cpu-cores per task (>1 if multi-threaded tasks)
#SBATCH --mem-per-cpu=16G         # memory per cpu-core (4G is default)
#SBATCH --time=24:00:00          # total run time limit (HH:MM:SS)
#SBATCH --mail-type=all          # send email when job begins/ends/fails
#SBATCH --mail-user=wenchang@princeton.edu
# 
##SBATCH --array=1-100#%32        # job array with index values 1, 2, ...,; max job # is 32 if specified
##SBATCH --output=slurm-%A.%a.out # stdout file
##SBATCH --error=slurm-%A.%a.err  # stderr file
#set -ev
##env settings
#export PATH=/tigress/wenchang/miniconda3/bin:$PATH
#export PYTHONPATH=/tigress/wenchang/wython
#export PYTHONUNBUFFERED=TRUE # see https://stackoverflow.com/questions/230751/how-to-flush-output-of-print-function
#export OMP_NUM_THREADS=$SLURM_CPUS_PER_TASK #for multi-threaded job
#ii_job=$SLURM_ARRAY_TASK_ID #for job array

#make restore mask
./wydo_FLOR_restore_mask.py

#obs tau_x and tau_y mclim
./wydo_mclim_era5_tau_x.py
./wydo_mclim_era5_tau_y.py

#model tau_x and tau_y mclim
./wydo_mclim_FLOR_tau_x.py
./wydo_mclim_FLOR_tau_y.py

#model bias: model - obs
./wydo_bias_FLOR_tau_x.py
./wydo_bias_FLOR_tau_y.py

#model correction: obs - model
./wydo_correction_FLOR_tau_x.py
./wydo_correction_FLOR_tau_y.py


#flux corrections from restore values
#mclim
./wydo_mclim_FLOR_pme_restore.py
./wydo_mclim_FLOR_sfc_hflux_restore.py
#prepare for model running
./wydo_correction_FLOR_salt.py
./wydo_correction_FLOR_temp.py