#!/bin/sh 

#
# Metafile Script : ghm_meta_grid
#
# Log :
# D.W.Plummer/NCEP   2/97   Add log header
# D.W.Plummer/NCEP	12/97   Change latlon bound calculations from using
#				"GRID AREA" to "LL CORNER" & "UR CORNER"
# B. Gordon/NCO      4/00   Converted for production on the IBM-SP.
#

set -x
mkdir $DATA/grid
cd $DATA/grid
cp $FIXgempak/datatype.tbl datatype.tbl

device="nc | grid.nmeta"

minlat=90
maxlat=-90
minlon=360
maxlon=-360

gdfile=${COMIN}/${RUN}g_${PDY}${cyc}f000_${storm_name}

export pgm=gdinfo;. $DATA/prep_step; $DATA/startmsg
$GEMEXE/gdinfo << EOF
GDFILE  = $gdfile
LSTALL  = no
OUTPUT  = f/gdinfo_$$.dat
l
r

EOF
   gpend

lat1=`grep "LL CORNER" gdinfo_$$.dat | head -1 | awk -F: '{print $2}' | \
      awk '{print $1}' | awk -F. '{print $1}' `
lon1=`grep "LL CORNER" gdinfo_$$.dat | head -1 | awk -F: '{print $2}' | \
      awk '{print $2}' | awk -F. '{print $1}' `
lat2=`grep "UR CORNER" gdinfo_$$.dat | head -1 | awk -F: '{print $2}' | \
      awk '{print $1}' | awk -F. '{print $1}' `
lon2=`grep "UR CORNER" gdinfo_$$.dat | head -1 | awk -F: '{print $2}' | \
      awk '{print $2}' | awk -F. '{print $1}' `

let lat1=lat1-1
let lon1=lon1-1
let lat2=lat2+1
let lon2=lon2+1

if [ $lon1 -gt 180 ] ; then
  let lon1=lon1-360
fi
if [ $lon1 -lt -180 ] ; then
  let lon1=360+lon1
fi
if [ $lon2 -gt 180 ] ; then
  let lon2=lon2-360
fi
if [ $lon2 -lt -180 ] ; then
  let lon2=360+lon2
fi

garea="${lat1};${lon1};${lat2};${lon2}"

latlon="1/2//1;1/5;5"
PDY2=`echo $PDY | cut -c3-`

export pgm=gdplot2_nc;. $DATA/prep_step; $DATA/startmsg
gdplot2_nc << EOF >> $NAOPX
GDFILE	= F-GHM | ${PDY2}/${cyc}00
GDATTIM	= FALL
DEVICE	= $device
PANEL	= 0
TEXT	= 1/21//hw
CONTUR	= 1
MAP	    = 1
CLEAR	= yes
CLRBAR  = 1

GAREA	= $garea
PROJ	= ced/0;0;0
LATLON	= $latlon

restore $USHgempak/restore/pmsl_thkn.2.nts
TITLE	= 5/-2/~GHM ^  MSLP, 1000-500 MB THKN for ${storm_name}|^ MSLP,1000-500 THKN!0
l
ru

GLEVEL  = 850
GVCORD  = PRES
PANEL   = 0
SKIP    = 0/2;2            !           !               !            !/3
SCALE   = 5                !5          !5              !5           !0
GDPFUN  = (avor(wnd))//v   !v          !mul(v,-1)      !mul(v,-1)   !kntv(wnd)
TYPE    = c/f              !c          !c/f            !c           !b
CONTUR  = 1
CINT    = 2/10/99          !2/4/8      !2/10/99        !2/4/8    
LINE    = 7/5/1/2          ! 29/5/1/2  ! 7/5/1/2       ! 29/5/1/2 
FINT    = 16;20;24;28;32;36;40;44
FLINE   = 0;23-15
HILO    = 2;6/X;N/10-99;10-99!        !2;6/X;N/10-99;10-99!    
HLSYM   = 
CLRBAR  = 1
WIND    = 9/0.7/2/112
REFVEC  =
TITLE   = 1//GHM ~ @ WINDS and VORTICITY for ${storm_name}|~ @ WINDS and VORT!0
TEXT    = 1/21//hw
CLEAR   = yes
filter  = no
l
ru

GLEVEL  = 500
ru

GLEVEL  = 200
ru

exit
EOF
export err=$?;$DATA/err_chk

#####################################################
# GEMPAK DOES NOT ALWAYS HAVE A NON ZERO RETURN CODE
# WHEN IT CAN NOT PRODUCE THE DESIRED GRID.  CHECK
# FOR THIS CASE HERE.
#####################################################
ls -l grid.nmeta
export err=$?;export pgm="GEMPAK CHECK FILE";$DATA/err_chk

if [ $SENDCOM = "YES" ] ; then
  mv grid.nmeta ${COMOUT}/ghm_${PDY}_${cyc}_${storm_name}
  if [ $SENDDBN = "YES" ] ; then
    $DBNROOT/bin/dbn_alert MODEL ${DBN_ALERT_TYPE} $job \
     $COMOUT/ghm_${PDY}_${cyc}_${storm_name}
  fi
fi

#
