source: trunk/LMDZ.MARS/deftank/occigen/run_month1

#!/bin/bash
#SBATCH --nodes=6                           #A
#SBATCH --ntasks-per-node=4                 #B
#SBATCH --cpus-per-task=5                   #C
#SBATCH -J clim
#SBATCH --time=13:00:00
#SBATCH --output run_month1.%j.out
##SBATCH --mem=8GB
#SBATCH --constraint=[BDW28]
#SBATCH --exclusive

## Script to run chained simulations
## (uses script "run0" and reference file "run.def.ref")
## Set values of "num_now" and "num_end" in the script below
## to set initial month # and final month # of the simulation

# A few parameters that might need be changed depending on your setup:
# Path to the arch.env to source
source ../trunk/LMDZ.COMMON/arch.env
# Environment setup for OpenMP:
export I_MPI_DOMAIN=auto
export I_MPI_PIN_RESPECT_CPUSET=0
#Make sure that OMP_NUM_THREADS = cpus-per-task * KMP_HW_SUBSET
export KMP_HW_SUBSET=1T
export OMP_NUM_THREADS=$SLURM_CPUS_PER_TASK #C
export OMP_STACKSIZE=400M
export KMP_AFFINITY=granularity=fine,compact,1,0,verbose
#
########################################################################
set -exv
ls -al
trap 'echo -e "Error at line $LINENO!"' ERR

# Set starting month and ending month below:
num_now=1
num_end=12
num_previous=$(( num_now - 1 ))

echo "$num_previous" > num_run
# next month number
num_next=$(( num_now + 1 ))
# true (i.e. modulo 12) month number
true_num=$(( num_now % 12 ))

# Check if required files exist
\rm -f  error; touch error

if [ ! -f run.def.ref ]; then
  echo "Error: file \"run.def.ref\" not found in current directory!" > error
  exit 1
fi

if [ ! -x run0 ]; then
  echo "Error: file \"run0\" not found or not executable in current directory!" > error
  exit 1
fi

if [ ! -f run_month$num_now ]; then
  echo "Error: file \"run_month${num_now}\" not found in current directory!" > error
  exit 1
fi

# Run model depending on current month
case $true_num in
  1 ) sed s/9999/61/ run.def.ref > run.def ; ./run0 >> error ;;    #1
  2 ) sed s/9999/66/ run.def.ref > run.def ; ./run0 >> error ;;    #2
  3 ) sed s/9999/66/ run.def.ref > run.def ; ./run0 >> error ;;    #3
  4 ) sed s/9999/65/ run.def.ref > run.def ; ./run0 >> error ;;    #4
  5 ) sed s/9999/60/ run.def.ref > run.def ; ./run0 >> error ;;    #5
  6 ) sed s/9999/54/ run.def.ref > run.def ; ./run0 >> error ;;    #6
  7 ) sed s/9999/50/ run.def.ref > run.def ; ./run0 >> error ;;    #7
  8 ) sed s/9999/46/ run.def.ref > run.def ; ./run0 >> error ;;    #8
  9 ) sed s/9999/47/ run.def.ref > run.def ; ./run0 >> error ;;    #9
 10 ) sed s/9999/47/ run.def.ref > run.def ; ./run0 >> error ;;    #10
 11 ) sed s/9999/51/ run.def.ref > run.def ; ./run0 >> error ;;    #11
  0 ) sed s/9999/56/ run.def.ref > run.def ; ./run0 >> error ;;    #12
  * ) echo "Error: Invalid value of true_num ($true_num)" ; exit 1 ;;
esac

# Launch job for next month
if (( num_next <= num_end )) ; then
  cp -f run_month$num_now tmp
  sed -e "s@run_month${num_now}@run_month${num_next}@" \
      -e "s@num_now=${num_now}@num_now=${num_next}@" tmp > run_month$num_next
  rm tmp
  sbatch run_month$num_next
fi