set -xa
date 

echo "------------------------------------------------"
echo "Executing exhrly_realtime_gempak.sh.sms XXXX"
echo "This script reads the gempak surface file"
echo "for making time series plots (previous 24 hour cyc" 
echo "to the current cyc) of hourly (current) and its"
echo "corresponding monthly mean dump counts for posting" 
echo "to the WEB using the gempak sfgram program"
echo "This script uses poe features with 1 node (16 tasks)" 
echo "Author: Krishna Kumar"
echo "NCO/System Integration Branch July 2008"
echo "Modified to automate all data types July 2008"
echo "Modified to extend the time series to 3 days for hrly/ruc May 2010"
echo "Modified to replace legacy RUC with Rapid Refresh (RAP) Nov 2011"
echo "Modified to apply scaling factor (divide by 100 for large counts)"
echo "for several data types under Hourly and RAP August 2012"
echo "------------------------------------------------"

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

ipn=0
for dtp in `cat file.dtype_$ntask`; do #### all data types in a single task
   ipn=`expr $ipn + 1`
   echo $ipn
   plotn=`cat plot.dtype_$ntask | awk -v pltn=$ipn '{ print $pltn }'`
   if [ $NET = "hourly" ]; then
      if [ $dtp = 'TRMM'   -o $dtp = 'GEOSTH' -o $dtp = '1BAMUB' -o \
           $dtp = 'RADW25' -o $dtp = 'RADW30' -o $dtp = 'ATMS'   -o \
           $dtp = 'RD2WND' -o $dtp = 'RD2RFT' -o $dtp = 'AIRSEV' -o \
           $dtp = 'AVCSAM' -o $dtp = 'AVCSPM' -o $dtp = '1BMHS'  -o \
           $dtp = 'WDSAT'  -o $dtp = 'ESAMUB' -o $dtp = 'ESHRS3' -o \
           $dtp = 'ESMHS'  -o $dtp = '1BHRS4' -o $dtp = 'SSMISU' -o \
           $dtp = 'OMI'    -o $dtp = 'LGYCLD' ]; then
         REALTEXT='1/10/CCNT (CURR COUNT) MCNT (MON COUNT) (*1/100)'
      else
         REALTEXT='1/10/CCNT (CURRENT COUNT) MCNT (MONTHLY COUNT)'
      fi
      NETR=HRLY
      netr=hly
   else
      NETR=RAP
      netr=rap
      if [ $dtp = 'GEOSTH' -o $dtp = 'GEOST1' -o $dtp = 'RADW25' -o \
           $dtp = '1BAMUB' -o $dtp = '1BMHS'  -o $dtp = 'ESAMUB' -o \
           $dtp = 'ESHRS3' -o $dtp = 'ESMHS'  -o $dtp = 'LGYCLD' -o \
           $dtp = 'RD2WND' -o $dtp = 'RADW30' ]; then
         REALTEXT='1/10/CCNT (CURR COUNT) MCNT (MON COUNT) (*1/100)'
      else
         REALTEXT='1/10/CCNT (CURRENT COUNT) MCNT (MONTHLY COUNT)'
      fi
   fi

   webadd='WWW.NCEP.NOAA.GOV'
   echo $webadd > file.webadd

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

sfgram << EOF > sfgram.$ntask
SFFILE  = ${NETH}_alldatatype.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-6;1;4;1
BORDER  = 1//2
MARKER  = 2/15/5//sw
TITLE   = 4/+16/$dtp
CLEAR   = no
DEVICE  = gif|${netr}.${plotn}.gif|896;716
PANEL   = 0.01;0;1;1
TEXT    = small/2////l/hw
r

EOF

sfmap << EOF > sfmap.$ntask
AREA   = -90;-180;90;180
SFPARM =
DATTIM = all
SFFILE = ${NETH}_alldatatype.sfc
CLEAR  = no
MAP    = 0
PANEL  = 0/1/1/2/view
DEVICE = gif|${netr}.${plotn}.gif|896;716
LATLON = 0
COLORS = 0;0;0
TEXT   = medium/2/hw
TITLE  = 1/4/${NETR} 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.$ntask
PANEL  = 0/1/1/2/view
COLORS = 4
TEXT   = medium/2////c/hw
CLEAR  = no
DEVICE = gif|${netr}.${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

done

gpend