#/bin/ksh
###########################################################################
echo "--------------------------------------------------------------------"
echo "exjtwc_submit.ksh.sms - Creates specific grib files for JTWC based  "
echo "                        off of GFS Master and UKMET output files"
echo "--------------------------------------------------------------------"
echo "History: Mar 03 2005 -  Entirely based off of script by Mike Fiorino"
echo "                                                      - Justin Cooke"
###########################################################################


set -x

# List and wgrib 

for tau in $ftaus
do
     
   dfile="${ddir}/$model.t${cyc}z.master.grbf${tau}"
         
   if [ $model == "ukmet" ]
      then
         dfile="${ddir}/${model}.t${cyc}z.ukmet${tau}"
   fi
         
   lfile="$DATA/${model}.p.${tau}.lst"
   $utilexec/wgrib $dfile > $lfile

done

rm $ofile 

touch $ofile

for tau in $ftaus
do
    lfile="$DATA/${model}.p.${tau}.lst"
    dfile="${ddir}/${model}.t${cyc}z.master.grbf${tau}"
    
    if [ $model == "ukmet" ]
       then
          dfile="${ddir}/${model}.t${cyc}z.ukmet${tau}"
    fi
    
    rm $tlfile $tgfile
    touch $tlfile

    for v in $vars
    do
       if [ $v == "zg" ]
       then
          kpds5=7
          kpds6=100
          kpds7="1000 500 200"
          var="HGT"
       elif [ $v == "ua" ]
       then
          kpds5=33
          kpds6=100
          kpds7="850 700 500 300 200"
          var="UGRD"
       elif [ $v == "va" ]
       then
          kpds5=34
          kpds6=100
          kpds7="850 700 500 300 200"
          var="VGRD"
       elif [ $v == "ta" ]
       then
          kpds5=11
          kpds6=100
          kpds7="850 700 500 300"
          var="TMP"
       elif [ $v == "vrta" ]
       then
          kpds5=41
          kpds6=100
          kpds7=""
          var="ABSV"
       elif [ $v == "uas" ]
       then
          kpds5=33
          kpds6=105
          kpds7=10
          var="UGRD"
       elif [ $v == "vas" ]
       then
          kpds5=34
          kpds6=105
          kpds7=10
          var="VGRD"
       elif [ $v == "tas" ]
       then
          kpds5=11
          kpds6=105
          kpds7=2
          var="TMP"
       elif [ $v == "pr" ]
       then
          kpds5=61
          kpds6=1
          kpds7=0
          var="APCP"
       elif [ $v == "hur" ]
       then
          kpds5=52
          kpds6=100
          kpds7="850 700 500 300"
          var="RH"
       elif [ $v == "wa" ]
       then
          kpds5=39
          kpds6=100
          kpds7=700
          var="VVEL"
       elif [ $v == "psl" ]
       then
          kpds5=2
          kpds6=102
          kpds7=0
          var="PRMSL"
       elif [ $v == "tasmax" ]
       then
          kpds5=15
          kpds6=105
          kpds7=2
          var="TMAX"
       elif [ $v == "tasmin" ]
       then
          kpds5=16
          kpds6=105
          kpds7=2
          var="TMIN"
       fi
       
#      If ukmet, special case

       if [ $model == "ukmet" ]
       then
          if [ $v == "uas" ]
          then
             kpds5=33
             kpds6=1
             kpds7=0
             var="UGRD"
          elif [ $v == "vas" ]
          then
             kpds5=34
             kpds6=1
             kpds7=0
             var="VGRD"
          fi
      fi
     
      for lev in $kpds7
      do
         grep "kpds5=$kpds5:kpds6=$kpds6:kpds7=$lev" $lfile >> $tlfile
      done
       
   done
   
   cat $tlfile | sort +0 -t: -n | $utilexec/wgrib $dfile -i -grib -o $tgfile

# Convert the GFS grids to 1.0 resolution
   
   if [ $model == "gfs" ]
   then
      $utilexec/copygb -g 3 -a -x $tgfile $ofile
   else
      cat $tgfile  >> $ofile
   fi
   
   rm $tgfile

done

# Convert GRIB1 to GRIB2

/nwprod/util/exec/cnvgrib -g12 -p40 $ddir/$model.$res.$PDY$cyc.jtwc_grb $ddir/$model.$res.$PDY$cyc.jtwc_grb.grib2

# Send DBNet alert

if [ $SENDDBN == "YES" ]
   then
      $DBNROOT/bin/dbn_alert MODEL JTWC_GFS jtwc_data $ddir/$model.$res.$PDY$cyc.jtwc_grb
fi

if [ $SENDDBN_GB2 == "YES" ]
   then
      $DBNROOT/bin/dbn_alert MODEL JTWC_GFS_GB2 jtwc_data $ddir/$model.$res.$PDY$cyc.jtwc_grb.grib2
fi