#! /bin/ksh 

# JHWRF_NMM_3DVAR_ANALYSIS -- generates WRF output and anayslis
# (wrfout, wrfanl) files needed by the JHWRF_RELOCATE and HWRF_MERGE
# jobs for storm relocation and initial storm intensity modification
# so that the initial storm will be close to that of the one specified
# in the message file.  This job also produces an initial time
# (actually, 1.5 minute forecast) WRF output file to be used by the
# JHWRF_GFS_TRACK0 job for GFS vortex location via the hwrf_gettrk.

# Last modified: Samuel Trahan, April 5, 2012.

# *** Operational Workflow ***

# Submitted: After the 12hr JHWRF_NMM_REAL completes
# Submits: JHWRF_GFS_TRACK0 to find the initial GFS vortex, and
#   JHWRF_NMM_3DVAR_ANALYSIS to interpolate or initialize all model 
#   fields to the larger ghost/dummy domain

# Primary scripts and parameter files:
#   scripts/exhwrf_nmm_3dvar_analysis.sh.sms -- passes control to 
#     ush/hwrf_nmm_3dvar.sh with the proper arguments to specify
#     the analysis case
#   ush/hwrf_nmm_3dvar.sh -- runs the hwrf_wrf in a special mode to
#     produce the wrfanl_d02 and wrfanl_d03 files in the ATMOS
#     directory.  Also produces a wrfout_d01_* file for the 1.5 minute
#     forecast to be used by the JHWRF_GFS_TRACK0 job
#   ush/hwrf_namelist.sh -- produces the namelist.input file needed by
#     the hwrf_wrf, using parm/hwrf_namelist.input
#   parm/hwrf_namelist.input -- input file to hwrf_namelist.sh, used
#     by all jobs that run WRF, to produce the namelist for that job
#   hwrf_stream_parse.pl -- a simple stream parser used to parse the
#     namelist file

# Primary executables:
#   hwrf_wrf -- the WRF atmospheric model
#   hwrf_swcorner_dynamic -- used to get the domain 2 start location

# *** EMC Suggested Job Card ***

#@ job_type = parallel 
#@ wall_clock_limit = 00:10:00
#@ total_tasks = 16
#@ node = 1
#@ task_affinity=core(1)
#@ node_resources = ConsumableMemory(55 GB)
#@ notification = never
#@ bulkxfer = yes
#@ network.MPI = csss,shared,us
#@ bulkxfer = yes


export PARAFLAG=${PARAFLAG:-NO}

if [ ${PARAFLAG} == NO ]; then
$SMSBIN/smsinit $LOADL_STEP_ID
fi
#----------------------------------------------
# JHWRF_NMM_3DVAR_ANALYSIS.sms  
#----------------------------------------------
set -xa
export MP_IOAGENT_CNT=all
export MP_IO_BUFFER_SIZE=8M

export PS4='$SECONDS + '
date

if [ ${PARAFLAG} = 'YES' ];
then
  export yymmdd=` echo $PDY | cut -c3-8`
  export yymmddhh=${yymmdd}${cyc}
  export jlogfile=${DATA}/jlogfile
  export SENDSMS=YES
  export SENDCOM=YES
  export SENDDBN=NO

  if [ ! -d $outdir ]; then mkdir -p $outdir; fi
  if [ ! -d $COMOUT ]; then mkdir -p $COMOUT; fi
else
  export NWPROD=${NWPROD:-/nwprod}
fi

###############################################################
# This block can be modified for different Production test
# environment. This is used for operational testings
###############################################################
if [ ${PARAFLAG} == NO -a $envir != prod ]; then
  export SENDDBN=${SENDDBN:-NO}
  export jlogfile=${jlogfile:-/com/logs/${envir}/jlogfile}
  export DBN_ALERT_TYPE=TBD_PARA
  export DBNROOT=/nwprod/spa_util/fakedbn
fi

# #############################################
# SETUP HWRF FCST PROCESSING VARIABLES
# #############################################

####################################
# obtain unique process id (pid) and make temp directory
####################################
export pid=$$
export DATA=${DATA:-/tmpnwprd/hwrf${storm_num}_${cyc}_${envir}}
cd $DATA

export cycle=t${cyc}z

####################################
# Specify NET and RUN Name and model
####################################
export NET=hur 
export RUN=hwrf

####################################
# File To Log Msgs
####################################
export jlogfile=${jlogfile:-/com/logs/jlogfile}

####################################
# Determine Job Output Name on System
####################################
export outid="LL$job"
export jobid="${outid}.o${pid}"
export pgmout="OUTPUT.${pid}"
export pgmerr=errfile

####################################
# SENDSMS  - Flag Events on SMS
# SENDCOM  - Copy Files From TMPDIR to $COMOUT
# SENDDBN  - Issue DBNet Client Calls
# RERUN    - Rerun fcst from beginning (default no)
####################################
export SENDSMS=${SENDSMS:-YES}
export SENDCOM=${SENDCOM:-YES}
export SENDDBN=${SENDDBN:-YES}

####################################
# Specify HWRF version number
####################################
export HWRF_VERSION=${HWRF_VERSION:-${model_ver}}

####################################
# Specify Execution Areas
####################################
export HOMEhwrf=${HOMEhwrf:-/nw${envir}/hwrf.$HWRF_VERSION}
export EXEChwrf=${EXEChwrf:-$HOMEhwrf/exec}
export FIXhwrf=${FIXhwrf:-$HOMEhwrf/fix}
export PARMhwrf=${PARMhwrf:-$HOMEhwrf/parm}
export STOREhwrf=${STOREhwrf:-/com/hwrf/${envir}/store}
export EXhwrf=${EXhwrf:-$HOMEhwrf/scripts}
export USHhwrf=${USHhwrf:-$HOMEhwrf/ush}

##############################
# Set up the UTILITIES
##############################
export utilscript=$NWPROD/util/ush
export utilities=$NWPROD/util/ush

##############################
# Run setup to initialize working directory and utility scripts
##############################
sh $utilscript/setup.sh
#export PDY=20120101

##############################
# Run setpdy and initialize PDY variables
##############################
sh $utilscript/setpdy.sh
. PDY

##############################################
# Define COM directories
##############################################
export COMIN=${COMIN:-/com/${NET}/${envir}/${RUN}.${PDY}${cyc}}
export COMOUT=${COMOUT:-/com/${NET}/${envir}/${RUN}.${PDY}${cyc}}
mkdir -m 775 -p $COMOUT

msg="HAS BEGUN on `hostname`"
postmsg "$jlogfile" "$msg"

#############################################################
# Execute the script
${EXhwrf}/exhwrf_nmm_3dvar_analysis.sh.sms

#if [ ${PARAFLAG} = 'YES' ]
#then 
#   ${EXhwrf}/exhwrf_nmm_3dvar_analysis.sh.sms
#else
#   /nw${envir}/scripts/exhwrf_nmm_3dvar_analysis.sh.sms
#fi 

#############################################################

msg="ENDED NORMALLY."
postmsg "$jlogfile" "$msg"

##############################
# Remove the Temporary working directory
##############################
cd /tmpnwprd
if [ ${PARAFLAG} = 'YES' ];
then
  echo " "
#else
#  rm -rf $DATA
fi

date

if [ ${PARAFLAG} == NO ]; then
$SMSBIN/endt
fi