set -xa
date
echo "------------------------------------------------"
echo "Executing XXXXXXXX exmodels_realtime_gempak.sh XXXXXXXX"
echo "This script reads the gempak surface file"
echo "for making multi-panel time series data counts plots (previous"
echo "9 days from cyc hour to the current day cyc hour) of current cyc"
echo "and its corresponding monthly mean dump counts for posting"
echo "to the WEB using the gempak sfgram program. This script"
echo "does not use poe features but uses only the serial job"
echo "Author: Krishna Kumar"
echo "NCO/Data Management & Quality Assessment Branch September 2004"
echo "Modified to add several new data types  June 2006"
echo "NCO/System Integration Branch July 2008"
echo "Modified to automate plotting data types in the bufr dump list"
echo "Modified to extend the time series to 30 days for all models May 2010"
echo "Modified to apply scaling factor (divide by 100 for large counts) to"
echo "new data types 'sevcsr' and 'atms' March 2012"
echo "Modified to apply scaling factor to data types 'esamub' 'eshrs3'"
echo "and 'esmhs' August 2012"
echo "Modified to apply scaling factor to data type 'cris' October 2012"
echo "------------------------------------------------"

echo "Making graphics for $NETMODEL model"

### Specify data types for single-panel plots
### by reading it from ${NET}_dataplotlist_t${cyc}z

. /nwprod/gempak/.gempak

echo $*
ntask=$1
shift
echo $*
echo "running ntask=$ntask"
echo "Display variable = $DISPLAY"

yymm=`echo $PDY | cut -c3-6`
yy=`echo $PDY | cut -c3-4`
hhmm=$cyc'00'
webadd='WWW.NCEP.NOAA.GOV'
echo $webadd > file.webadd

ipn=0
for dtp in `cat file.dtype_$ntask`; do #### all data types in one go!!!
   ipn=`expr $ipn + 1`
   plotn=`cat plot.dtype_$ntask | awk -v pltn=$ipn '{ print $pltn }'`
   if [ $NET = "nam" ]; then 
      if [ $dtp = '1BAMUB' -o $dtp = '1BMHS'  -o $dtp = 'RADW30' -o \
           $dtp = 'RADW25' -o $dtp = 'RD2WND' -o $dtp = 'GEOST1' -o \
           $dtp = 'ESAMUB' -o $dtp = 'ESHRS3' -o $dtp = 'ESMHS'  -o \
           $dtp = 'MSOAWS' -o $dtp = 'MSOAPR' ]; then
         REALTEXT='1/10/CCNT (CURR COUNT) MCNT (MON COUNT) (*1/100)'
      elif [ $dtp = 'RD2RFT' ]; then
         REALTEXT='1/10/CCNT (CURR COUNT) MCNT (MON COUNT) (*1/1000)'
      else
         REALTEXT='1/10/CCNT (CURRENT COUNT) MCNT (MONTHLY COUNT)'
      fi
   elif [ $NET = "gfs" ]; then
      if [ $dtp = 'TRMM'   -o $dtp = '1BAMUA' -o $dtp = '1BAMUB' -o \
           $dtp = '1BMHS'  -o $dtp = '1BHRS3' -o $dtp = '1BHRS4' -o \
           $dtp = 'AMSRE'  -o $dtp = 'GEOST1' -o $dtp = 'VISEU'  -o \
           $dtp = 'AIRSEV' -o $dtp = 'H20EU'  -o $dtp = 'WDSAT'  -o \
           $dtp = 'OMI'    -o $dtp = 'ESAMUB' -o $dtp = 'ASCAT'  -o \
           $dtp = 'ESHRS3' -o $dtp = 'SSMISU' -o $dtp = 'MTIASI' -o \
           $dtp = 'AVCSAM' -o $dtp = 'AVCSPM' -o $dtp = 'ESAMUA' -o \
           $dtp = 'ATMS'   -o $dtp = 'SEVCSR' -o $dtp = 'ESMHS'  -o \
           $dtp = 'CRIS' ]; then
         REALTEXT='1/10/CCNT (CURR COUNT) MCNT (MON COUNT) (*1/100)'
      else
         REALTEXT='1/10/CCNT (CURRENT COUNT) MCNT (MONTHLY COUNT)'
      fi
   fi

###############################################
# Plot gempak sfgram files of various data types
# and create gif images for NAM, GFS, GDAS models
################################################

### The file names are for example gfs model
###                    gfs_alldatatype_filem0_t12z
###                    gfs_alldatatype_filem1_t12z
###                    etc ..... ,,   _filem7_t12z

   begfile=`sh /nwprod/util/ush/finddate.sh $PDY s-30`
   byea=`echo $begfile | awk '{ print $30 }' | cut -c3-4`
   bday=`echo  $begfile | awk '{ print $30 }'| cut -c7-8`
   bmon=`echo $begfile | awk '{ print $30 }' | cut -c5-6`
   eyea=`echo $PDY | cut -c3-4`
   eday=`echo $PDY | cut -c7-8`
   emon=`echo $PDY | cut -c5-6`

sfgram << EOF > sfgram.$NETMOD
SFFILE  = ${NETMOD}_alldatatype_alldays_${cycle}.sfc
DATTIM  = all
STATION = $dtp
TRACE1  = ccnt:1:4;mcnt:2:4/2;3/ /TIC
NTRACE  = 1
TAXIS   = ${byea}${bmon}${bday}/${hhmm}-${eyea}${emon}${eday}/${hhmm}-24;3;3;1
BORDER  = 1//2
MARKER  = 2/15/5//sw
TITLE   = 4/+16/$dtp
CLEAR   = no
DEVICE  = gif|${NETMOD}.${plotn}.gif|896;716
PANEL   = 0.01;0;1;1
TEXT    = small/2/hw
r

EOF

sfmap << EOF > sfmap.$NETMOD
AREA   = -90;-180;90;180
SFPARM =
DATTIM = all
SFFILE = ${NETMOD}_alldatatype_alldays_${cycle}.sfc
CLEAR  = no
MAP    = 0
PANEL  = 0/1/1/2/view
DEVICE = gif|${NETMOD}.${plotn}.gif|896;716
LATLON = 0
COLORS = 0;0;0
TEXT   = medium/2/hw
TITLE  = 1/4/${NETMODEL} MODEL $dtp DATA COUNTS for ${cyc}Z 
r

MAP    = 0
LATLON = 0
COLORS = 0;0;0
TITLE  = 1/6/TIME SERIES ENDING IN ${PDY} ${cyc}Z
r

MAP    = 0
LATLON = 0
COLORS = 0;0;0
TEXT   = medium/2/hw
TITLE  = $REALTEXT
r

MAP    = 0
LATLON = 0
COLORS = 0;0;0
TEXT   = medium/2/hw
TITLE  = 4/+48/DOC/NOAA/NWS/NCEP/NCO
r

MAP    = 0
LATLON = 0
COLORS = 0;0;0
TEXT   = medium/2/hw
TITLE  = 2/+43/DATE 20${byea}${bmon}${bday}/${hhmm} Z-20${eyea}${emon}${eday}/${hhmm} Z
r

EOF

gptext << EOF > gptext.$NETMOD
PANEL  = 0/1/1/2/view
COLORS = 4
TEXT   = medium/2////c/hw
CLEAR  = no
DEVICE = gif|${NETMOD}.${plotn}.gif|896;716
TXTFIL = file.webadd
TXTLOC = 0.5;.08
r

TXTFIL = NOAA | 2.65 | C
TXTLOC = 0.15;0.1
r

TXTFIL = NWS | 2.5 | C
TXTLOC = 0.85;0.1
r

EOF

gpend

done

#####################################################################
# GOOD RUN
set +x
echo "************** exmodels_realtime_gempak.sh COMPLETED NORMALLY ON THE IBM"
echo "************** exmodels_realtime_gempak.sh COMPLETED NORMALLY ON THE IBM"
echo "************** exmodels_realtime_gempak.sh COMPLETED NORMALLY ON THE IBM"
set -x
#####################################################################

msg="HAS COMPLETED NORMALLY!"
echo $msg
$DATAO/postmsg "$jlogfile" "$msg"

############## END OF SCRIPT #######################