#!/bin/csh

# This script shows how to use gpmap to derive GIF radar images from a WSI NOWrad composite.
# To use this script you must install GEMPAK.


# USAGE -  makeNOWrad_gif <filename>
#
# where <filename> is a WSI NOWrad filename.  File is in native WSI format (not HDF or MDV).

# OUTPUT - Generates four GIF images over specific subregions. Image filenames automatically 
# generated using prefixes 'national_', 'cent_missvly_', 'cent_plains_', and 'south_missvly_'.

# strip '.MASTER15' off the WSI MASTER15 filename and replace with '.gif' 
set outfile = `basename $1 .MASTER15`.gif



# Make sure there is no more than one argument.
if ("x$2" != "x") then
	echo "only one argument please"
	exit 1
endif



# These GAREA definitions correspond to different regions of the continental U.S.
# They may be cut-and-pasted below.  Remove the '#' comment symbol before using, and
# you might want to change the output filename too in the DEVICE line. 

# cent_missvly # GAREA    = 34.1;-95.4;42.0;-83.9 
# cent_plains # GAREA    = 34.7;-105;43.2;-92.6
# national # GAREA    = 24.4;-125;50.3;-66.5
# north_plains # GAREA    = 39.7;-105.3;49.8;-90.5 
# south_missvly # GAREA    = 28.5;-96.0;36.6;-85.0
# south_plains # GAREA    = 25.8;-107.0;37.7;-90.9
# upper_missvly # GAREA    = 40.1;-98.5;49.6;-85.3

# Initialize GEMPAK environment so you can execute gpmap.
source $NAWIPS/Gemenviron


gpmap<<END_INPUT
 MAP      = 23
 \$MAPFIL  =
 GAREA    = 24.4;-125;50.3;-66.5
 PROJ     = sat
 LATLON   = 0
 SATFIL   = $1
 LUTFIL   = wsi_mosaic.tbl
 TITLE    = 31/-2
 DEVICE   = GIF|national_$outfile|768;576
 STNPLT   = 
 TEXT     = 1/3/1/hw
 CLEAR    = YES
 LINE     = 3
 IMBAR    = 1/ / /0.0;0.1/0.8;/|0.7/3/1/hw
r

 MAP      = 23 + 1//1.5  
 \$MAPFIL  = HICNUS.NWS + MESTUS.NWS
 GAREA    = 34.1;-95.4;42.0;-83.9 
 DEVICE   = GIF|cent_missvly_$outfile|768;576
 STNPLT   = 0/0.6|1/2/0.3/|sfstns.tbl
r

 GAREA    = 34.7;-105;43.2;-92.6
 DEVICE   = GIF|cent_plains_$outfile|768;576
r

 GAREA    = 28.5;-96.0;36.6;-85.0
 DEVICE   = GIF|south_missvly_$outfile|768;576
r

exit
END_INPUT

# Gempak uses gpend to finish its graphical output.
gpend