#!/bin/sh -evx 

#YEAR1=1960
YEAR1=1994
YEARN=2000
YEARD=1
IBASE=/projects/NS9039K/shared/norcpm/cases/NorESM1-CMIP6/noresm1-cmip6_hindcast
IPREF=noresm1-cmip6_hindcast_
NMEM=10 
ZTOP=0000
#ZBOT=1000
ZBOT=2000
#ZBOT=10000
MATLAB='/tos-project4/NS9039K/software/matlab_R2013b_nird/bin/matlab -nodisplay -nojvm -singleCompThread'

COMP_TS2D () { 
 $MATLAB -r "comp_TS2D('"${IBASE}/${IPREF}${YEAR}1015"','"${IPREF}${YEAR}1015_mem${MEM}"',"$ZTOP","$ZBOT","`expr $YEAR + 1`","`expr $YEAR + 10`");exit"

}

for YEAR in `seq $YEAR1 $YEARD $YEARN` 
do
  for BATCH in `seq 1 1` 
  do 
    M1=`expr 1 + \( $BATCH - 1 \) \* 10` 
    M1=`echo 0$M1 | tail -3c`  
    M2=`expr $BATCH \* 10` 
    M2=`echo 0$M2 | tail -3c`  
    for MEM in `seq -w $M1 $M2`
    do
      OFILE=Data_TS2D/TS2D_${IPREF}${YEAR}_mem${MEM}_${ZTOP}-${ZBOT}_`expr ${YEAR} + 1`-`expr $YEAR + 10`.matXXX
      echo $OFILE
      if [ -e $OFILE ]
      then
        echo $OFILE exists. skipping ...
      else
        COMP_TS2D & 
      fi
    done
    wait 
  done 
done