#!/bin/sh -evx 

YEAR1=1950
YEARN=2010
IBASE=/projects/NS9039K/shared/norcpm/cases/NorCPM/noresm1-cmip6_analysis_19500115
IPREF=noresm1-cmip6_analysis_19500115_mem
ZTOP=0000
#ZBOT=1000
ZBOT=2000
#ZBOT=10000
%MATLAB='/tos-project4/NS9039K/software/matlab_R2013b_nird/bin/matlab -nodisplay -nojvm -singleCompThread'
MATLAB='matlab -nodisplay -nojvm -singleCompThread'

COMP_TS2D () { 
 $MATLAB -r "comp_TS2D('"$1"','"$2$3"',"$4","$5","$6","$7");exit" < /dev/null 
}

  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_${IPREF}${MEM}_${ZTOP}-${ZBOT}_${YEAR1}-${YEARN}.nc
    echo $OFILE
    if [ -e $OFILE ]
    then
      echo $OFILE exists. skipping ...
    else
      COMP_TS2D $IBASE $IPREF $MEM $ZTOP $ZBOT $YEAR1 $YEARN & 
    fi
    done
    wait 
  done