#!/bin/ksh
set -ax

#------------------------------------------------------------------------------
#
#   NAME:  ocd5_run.ksh 
#
#   LANGUAGE: K-shell
#
#   PURPOSE:  End to end processing of the SHIPS Intensity Hurricane Model. 
#
#   DESCRIPTION:  Primary script used to drive the SHIPS Intensity model
#                 processing.
#
#   INFORMATION:
#            National Hurricane Center / Tropical Prediction Center
#            Department of Commerce - NOAA / NWS
#
#   CATEGORY:  Statistical Model
#
#   CALLING SEQUENCE:  ocd5_run.ksh <basin>CCYYYY.com a<basin>CCYYYY.dat
#
#   INPUTS:
#        <basin>CCYYYY.com   <-- computes CARQs 
#        a<basin>CCYYYY.dat  <-- A-deck format, containing CLP5 and SHF5 forecasts 
#                                for the date/time in the CARQ file
#
#   OUTPUTS: ocd5.dat - OCD5 forecast in A-deck format
#
#   EXTERNAL:
#
#   SIDE EFFECTS: Unknown
#
#   RESTRICTIONS:
#
#   EXAMPLE or USAGE: ocd5_run.ksh al152013.com aal152013.dat
#
#   OTHER NOTES:  None
#
#   MODIFICATION HISTORY:
#        Written by:    M. DeMaria 15-JUN-2013 (NOAA-NESDIS)
#
#   REVISION INFORMATION (RCS Keyword):
#
#     $Id$
#
#-----------------------------------------------------------------------------
#... DEFINE PATHS ...
#--------------------
# note - user should define these paths 

#
# common path
#commonp="/nhc/save/guidance/prgms"                           ### need to move to j-job script ?
commonp="${HOMEnhcg}"  

#
# Specify name of directory containing executables
xdir=${commonp}"/exec/"

# Specify name of directory containing static data files
#ydir=${commonp}"/coefs/"
ydir=${commonp}"/parm/"

#
# Log directory
#ldir="/nhc/save/guidance/storm-data/log/"
ldir="${WORK_DIR}/storm-data/log/${stormid}/"
#ldir="/nhc/save/guidance/prgms/test/ocd5/"

mkdir -p -m 775 $ldir

#---------------------------------------
#... EXPORT ENV for COEFS/DATA files ...
#---------------------------------------

#
# Coefs and Data files ...
#export SHIPS_COEF="/nhc/save/guidance/prgms/coefs/"
export SHIPS_COEF="${HOMEnhcg}/parm/"


#------------
#... LOGS ...
#------------

#
# Set LOG equal to the storm ID
LOG="OCD5-"`echo ${1##/*/} | sed -e 's/\.com//'`
LOG=${ldir}${LOG}.`date +"%Y%m%d.%H%M"`".log"

# -----------------------------------------------------------
# BEGIN SHIPS PREP :  Track, Grib and Input Deck Creation 
# -----------------------------------------------------------

echo "-------------------------------------------------------" > $LOG
echo "             ***** OCD5 RUN STARTED ***** "             >> $LOG
echo "-------------------------------------------------------">> $LOG
echo "start time: "`date +"%Y/%m/%d %T %Z doy=(%y%j)"`        >> $LOG
echo ""                                                       >> $LOG

#
# Create the ocd5.com input file for the model run
rm -f ocd5.com
sleep 1
cabname=${1}

# diagnostic info
echo "[${cabname}] [${1}]"                                     >> $LOG
hostname                                                       >> $LOG
cp -f $cabname ocd5.com                                        >> $LOG 2>&1

echo ""                                                        >> $LOG
echo "***** COM FILE *****"                                    >> $LOG
echo "-------------------------------------------------------" >> $LOG
echo " [${cabname}]                                          " >> $LOG
cat ocd5.com                                                   >> $LOG
echo "-------------------------------------------------------" >> $LOG
echo ""                                                        >> $LOG
#
# Check to make sure ocd5.com file is not empty
if [ ! -s ocd5.com ]; then
   
   echo "*** ERROR: ocd5.com file is empty, halt ocd5 processing." \
   >> $LOG
   exit 1
fi

echo "---- getting a CLP5/SHF5 forecast file----"              >> $LOG
echo ""                                                        >> $LOG

#
# Create the ocd5.inp file containing SHF5 and CLP5 forecasts in A-deck format            
rm -f ocd5.inp
itrname=${2}
cp ${itrname} ocd5.inp
#
# Check to make sure itrack.dat file is not empty
if [ ! -s ocd5.inp ]
then
   echo "*** ERROR ***: ocd5.inp file is empty, halt OCD5 processing." \
   >> $LOG
   exit 1
fi

#
echo " "                                                        >> $LOG
echo "---- CLP5/SHF5 forecasts available in ocd5.inp file ----" >> $LOG
grep CLP5 ocd5.inp                                              >> $LOG
grep SHF5 ocd5.inp                                              >> $LOG

echo " "                                                        >> $LOG
echo "---- Running the OCD5 fortran code  ----"  >> $LOG

rm -f ocd5.dat
rm -f ocd5.log

export pgm=ocd5.x
 prep_step
 startmsg
${xdir}"ocd5.x"                >> $LOG
export err=$?;  err_chk;

#
echo ""                          >> $LOG
echo "---- OCD5.x log file ----" >> $LOG
cat ocd5.log                     >> $LOG

#
echo ""                                     >> $LOG
echo "---- OCD5.x output file for ATCF----" >> $LOG
cat ocd5.dat                                >> $LOG

# --- Perform Cleanup ---

# Clean up after ships model run
rm -f ocd5.inp
rm -f ocd5.com
rm -f ocd5.log

echo ""                                                        >> $LOG
echo "-------------------------------------------------------" >> $LOG
echo "             ***** OCD5 RUN COMPLETED *****"             >> $LOG
echo "-------------------------------------------------------" >> $LOG
echo "end time: "`date +"%Y/%m/%d %T %Z doy=(%y%j)"`           >> $LOG
echo ""                                                        >> $LOG

# compress $LOG
#gzip $LOG

# --- End of Model Run ---
#
# exit with success stat
exit 0
#

#------------------------------------------------------------------------------
#
# EXPANDED REVISION INFORMATION (RCS Keywords):
#
#   $Id$
#   $Author$
#   $Log$
#
#------------------------------------------------------------------------------