#!/bin/csh -f
#
# This command file will create a surface map for the specified forecast hour and model 
# on X windows.
#
# Philip Schumacher	12/96
#
# Check if a data has been input.  If not print instructions.
if ($1x == x || $2x == x || $3x == x) then
        echo "This script creates a mean sea level pressure map and boundary temps    "
	echo "for a map centered on North Dakota."
	echo "make_map."
	echo " " 
	echo "Syntax:  look_sfc YYMMDDHH FHR MDL                          "
	echo " "
	echo "  where:                                               "
	echo "         YYMMDDHH is the year, month, day and hour of the model "
	echo "         to be examined       "
	echo " "
	echo "         FHR is the forecast hour of the map to be created " 
        echo " "
        echo "         MDL is the model to look at (in small letters)"
	echo "         ngm, ruc, avn, eta                            "
        echo " " 
	echo "   Example:  look_sfc 96102912 F12 eta                    "
	echo " "
	echo "   This create a map from the 29th of the month at 00Z"
	exit

endif

# Set up the name of the ASCII model file.

set GEMDAT = ${2}
set GEMFIL = $METDAT/gempak/grids/${1}_${3}.gem 
#set GEMFIL = /usr2/phil/oct29/gempak/surface/${1}_sao.gem
echo ${GEMFIL}

# Now run sfmap to create quadrant 1

    
 gdplot << EOF
 res $GEMNTS/${3}looksfc
 GDFILE   = $GEMFIL
 GDATTIM  = $GEMDAT
 r
   
exit

EOF

echo `********************************`
echo `*** TYPE gpend WHEN FINISHED ***'
echo `********************************`

exit
