#!/bin/sh
set -aeux

IM0=IM1       # Dimension of each face, for cube the total number of grid points in one layer is IM0xIM0x6
LM=LM1        # Number of vertical layers, should always be 38 for this version
NSUB=NSUB1    # Subdivision number, e.g. if NSUB 2, each face are divided
NMIN=NMIN1    # Is the number of faces, should always be 6 for cube


TMP_DIR=__TMP_DIR__
PRP_DIR=__PRP_DIR__
ETA_RUN=__ETA_RUN__
FInitBC=__FInitBC__

INIM0=`cat ${PRP_DIR}/../configure | grep "IM0=" | cut -d= -f2`
INLM0=`cat ${PRP_DIR}/../configure | grep "LM=" | cut -d= -f2`
INSUB=`cat ${PRP_DIR}/../configure | grep "NSUB=" | cut -d= -f2`
INNMI=`cat ${PRP_DIR}/../configure | grep "NMIN=" | cut -d= -f2`


if [ ! -s ${PRP_DIR}/data_in/grid/topo.dat ] ; then


cd ${PRP_DIR}/data_in/grid/

${PRP_DIR}/exe/initCube.exe > ${ETA_RUN}/../err/err.initcube 2>&1

${PRP_DIR}/exe/metrics.exe > ${ETA_RUN}/../err/err.metrics 2>&1

export PBS_SERVER=eslogin13
qsub ${TMP_DIR}/run_topo.sh


LastArq="False"
while [ "${LastArq}" != "True" ] ; do

    tam=`du -b ${PRP_DIR}/data_in/grid/topo.dat | awk '{print $1}'` 

      if [ -s ${PRP_DIR}/data_in/grid/topo.dat ] ; then

         if (( ${tam} > 7795600 )) ; then

            ${TMP_DIR}/run_init.sh ${LM}
            LastArq="True" 
         
	 else
  
            sleep 200
          
	 fi
   
            else
 
          sleep 200			
      fi
	
done

else

cd ${TMP_DIR}		
./run_init.sh ${LM}

fi