#!/bin/sh

#####################################################################
# Script tpc_psurge_nhctrk.sh generates the hypothetical tracks for the
# Probability Storm Surge Model
# usage: tpc_psurge_nhctrk.sh <ScriptDir> <stormAdvDir> <storm> <advTime>
#####################################################################
set -x

if [[ -z $4 ]] || [[ ! -z $5 ]]
then
  echo "usage: $0 <ScriptDir> <stormAdvDir> <storm> <advTime>"
  exit 1
fi

ScriptDir=${1:-$DATA}
stormAdvDir=${2:-$DATA}
storm=${3:-$stormid}
advTime=${4:-$AdvTime}
if [ $RUN_ENVIR = "NCO" ]
then
   EXECpsurge=${EXECpsurge:-/nwprod/exec}
   SloshbsnDir=$PARMpsurge/tpc_sloshbsn
else
   SloshbsnDir=${ScriptDir}/../parm/tpc_sloshbsn
   EXECpsurge=${ScriptDir}/../exec
fi


# following awk script ignores first line, and first row.
bsnList=`cat ${SloshbsnDir}/ocean.txt | awk '{if (NR >= 2) {$1=""; print}}' | sed -e s#,#' '#g`

runDir=${stormAdvDir}/run
ansDir=${stormAdvDir}/ans

#######################################################
# Create all possible basin directories for this storm.
#######################################################
if [ ! -d ${runDir} ]
then
  mkdir ${runDir}
  mkdir ${ansDir}
fi

cd $runDir

for bsn in  $bsnList
do
  mkdir -p ${runDir}/${bsn}
done

/nwprod/util/ush/setup.sh


#########################################################################
# Create the storm tracks in the useful basin directories for this storm.
#########################################################################
histFile=${stormAdvDir}/b${stormid}.dat
fcstFile=${stormAdvDir}/${stormid}.fst
tideFile=${stormAdvDir}/${stormid}_${advTime}_wlevel.dat

export pgm=tpc_psurge_nhctrk
. prep_step
$EXECpsurge/tpc_psurge_nhctrk ${SloshbsnDir} ${histFile} ${fcstFile} ${tideFile} ${runDir} -mode3 >>nhctrk.out
export err=$?; err_chk

#####
# Copy weight files to the ans directory, or remove the empty basin directories 
# for this storm.
#####
for bsn in $bsnList
do
  if [ -e ${runDir}/${bsn}/weight.txt ]
  then
    if [ ! -e ${ansDir}/${bsn} ]
    then
      mkdir ${ansDir}/${bsn}
    fi
    cp ${runDir}/${bsn}/weight.txt ${ansDir}/${bsn}/weight.txt
  else
    rm -rf ${runDir}/${bsn}
  fi
done

#####
# Signal that we are finished generating tracks.
#####
echo "Done with tpc_psurge_nhctrk.sh" > ${runDir}/nhctrkDone.txt

exit 0