20020521: GEMPAK Help

Rebecca,

In GEMPAK programs, comment lines within programs will begin with "!".
When you use "$" within a shell script, you will escape it, 
(eg \$ ) unless you want the shell to interpret some environmental variable.

Below, you also had a typo:
nbd=bg/24+world_bnds/8+lakes/24
should be
bnd=bg/24+world_bnds/8+lakes/24

Also, in the title string, the characters "-" and "@" are only
special characters in the title in grid programs and not sfmap.

Here is a short example that shows using GPMAP to draw the
filled boundaries using BND, and then overlays METAR data on
the map backgrouund using the SFMAP program.

Steve Chiswell
Unidata User Support

#!/bin/csh -f
#--------------------- example script --------------------
#
# source the Gemenviron script for the GEMPAK configurations
#
source /home/gempak/Gemenviron

#
# set the DEVICE environmental variable which
# will be used in the programs below
#
set DEVICE="xw"

# Note: $respond and $mapfil are GEMPAK variables, escape the $ so that
# csh will not try to interpret them as environmental variables
gpmap << EOF
 !
 !
 \$respond = yes
 \$mapfil = hipowo.cia
 MAP      = 32
 GAREA    = 22.0;-122.0;47.0;-65.0
 PROJ     = lcc/25;-103;60
 SATFIL   =
 RADFIL   =
 LATLON   = 0
 PANEL    = 0
 TITLE    =
 TEXT     =
 CLEAR    = y
 DEVICE   = ${DEVICE}
 LUTFIL   =
 STNPLT   =
 VGFILE   =
 AFOSFL   =
 AWPSFL   =
 LINE     =
 WATCH    =
 WARN     =
 HRCN     =
 ISIG     =
 LTNG     =
 ATCF     =
 AIRM     =
 NCON     =
 SVRL     =
 BND      = bg/24+world_bnds/8+lakes/24
 ATCO     =
 TCMG     =
 QSCT     =
 IMBAR    =
 run

 exit
EOF

sfmap << EOF2
 !
 ! This example will plot the standard METAR variables
 ! using the parameters and colors defined in $GEMTBL/config/prmlist.tbl.
 !
 AREA     = dset
 GAREA    = 22.0;-122.0;47.0;-65.0
 SATFIL   =
 RADFIL   =
 SFPARM   = metar
 COLORS   =
 DATTIM   = 020521/2000
 SFFILE   = metar
 !
 ! set clear to "no" to overlay on the map drawn in GPMAP
 CLEAR    = n
 !
 ! do not redraw the map here, just overlay the data on the GPMAP output
 MAP      = 0
 LATLON   = 0
 TITLE    = 1/0
 PANEL    = 0
 DEVICE   = ${DEVICE}
 PROJ     = lcc/25;-103;60
 FILTER   = .8
 TEXT     = 1
 LUTFIL   =
 STNPLT   =
 CLRBAR   =
 IMBAR    =
 run

 exit
EOF2

#
# close the gif file if necessary
if ( $DEVICE == "gf" ) then
   gpend
endif


#--------------------- end example script --------------------


>From: Rebecca Ratz <address@hidden>
>Organization: UCAR/Unidata
>Keywords: 200205212024.g4LKOBa04515

>Steve,
>
>Let me begin by saying "thank you." I appreciate your patience with me. I
>just picked up the manual last week, so I'm still learning the basics.
>
>I tried making the map that you suggested at:
>
>http://motherlode.ucar.edu/unidata/images/nids/radar_mosaic.gif
>
>This is exactly what I'm looking for, only I can't get it to work. The
>file I'm using looks like this:
>
>gpmap <<endfile
>map=32
>garea=22.0;-122.0;47.0;-65.0
>%area=41.;-91.;48.;-78.
>%garea=41.;-91.;48.;-78.
>satfil=$GARPHOME/map/areas/topo_na.area
>radfil=
>sfparm=skyc:.75;tmpf;wsym:.75;smsl;;dwpf;;;;brbk:1:2
>colors=7;7;7;20;22;18;24;6
>nbd=bg/24+world_bnds/8+lakes/24
>dattim=last
>sffile=metar
>panel=0
>title=1/-2/- @ Surface observations for $ymdh
>text=tiny/10/sw
>clear=yes
>proj=lcc/25;-103;60/3;3;0;0
>$mapfil=hipowo.cia
>filter=1
>lutfil=
>stnplt=
>%map=4/1/1
>latlon=
>device=xw
>run
>
>endfile
>
>The image turns out with a white background...no filler, but world
>boundaries are brown and lakes outlined in blue. I'm not
>advanced enough to understand "In this image, I first use GPMAP to color
>the map areas, then display the radar data over top in gdplot2." Do you
>have access to a file that produces such graphs? Especially one that
>has detailed descriptions of commands? I wouldn't mind using something
>like that to learn the basic commands.
>
>Thank you,
>
>Rebecca Ratz
>
>
>
>
>
>
>On Mon, 20 May 2002, Unidata Support wrote:
>
>> >From: Rebecca Ratz <address@hidden>
>> >Organization: UCAR/Unidata
>> >Keywords: 200205202056.g4KKu1a04652
>>
>> >
>> >I am currently conducting a summer research session on using GEMPAK. I
>> >am hoping someone might be able to talk me through some of the problems
>> >that I am encountering.
>> >
>> >My goal is to write a file that contains temperature, precipitation,
>> >and lightning data on one map. I'm having trouble tweeking the
>> >background color, contouring, etc. Everytime I try using the variable
>> >BND=bg/31,to produce a white background, the program pops back staying
>> >that "BND" is an unrecognized parameter. Some examples of the contouring
>> >and coloring that I am looking for can be seen at...
>> >
>> >http://twister.sbs.ohio-state.edu/surface/mw.gif
>> >http://www.comet.ucar.edu/resources/cases/mapindex.htm
>> >
>> >The first website contains the background and state boundary looks that I
>> >was hoping for. As for the contouring, shading and darkening of colors
>> >in second web address...can you tell what function is used to produce
>> >such a picture?
>> >
>>
>> BND is a parameter in GPMAP used for filling geographically defined
>> boundaries as show below. It is not a parameter in SFMAP for example,
>> which is the plot you mention in the first example above.
>>
>> To set the background color to white, you have 2 easy options.
>>
>> 1) if you are running interactively on the desktop, pop up
>> the ncolor program from the ntl bar. Then, click on the leftmost color
>> tile "0" and set the color to white.
>>
>> 2) An alternative to the above, and generally the easiest way if you are
>> running scripts to create maps is to copy the $GEMTBL/colors/coltbl.xwp
>> to the directory where you are running your GEMPAK programs and
>> change the first color line in the file from:
>> BLACK           BLA      0      0      0  black
>> to
>> WHITE           WHI    255    255    255  white
>>
>> By default, GEMPAK will use the coltbl.xwp file in your current working
>> directory if it exists, and if not, look for the file in $GEMTBL/colors.
>>
>>
>> Secondly, the topographic images are stored in AREA files wich
>> can be displayed using the SATFIL parameter. The topography files
>> can be found in $GARPHOME/map/areas. To use these, set PROJ=sat
>> in your display program, and then SATFIL=$GARPHOME/map/areas/topo_na.area
>> for example.
>>
>> Or, in GPMAP, you can use the BND parameter to create a filled map like:
>> http://motherlode.ucar.edu/unidata/images/nids/radar_mosaic.gif
>> In this image, I first use GPMAP to color the map areas, then display the
>> radar data over top in gdplot2.
>> In GPMAP I used:
>>  MAP      = 32
>>  GAREA    = 22.0;-122.0;47.0;-65.0
>>  PROJ     = lcc/25;-103;60/3;3;0;0
>>  $MAPFIL  = hipowo.cia
>>  BND      = bg/24+world_bnds/8+lakes/24
>>
>> then, in your overlaying program, set CLEAR=NO.
>>
>> Steve Chiswell
>> Unidata User Support
>>
>>
>>
>>
>> >
>> >Thanks for your time, it is greatly appreciated?
>> >
>> >
>> >Rebecca Ratz
>> >
>> >address@hidden
>> >University of Michigan
>> >Undergraduate in Mechanical Engineering
>> >
>> >
>> >
>> >
>>
>> ****************************************************************************
>> Unidata User Support                                    UCAR Unidata Program
>> (303)497-8643                                                  P.O. Box 3000
>> address@hidden                                   Boulder, CO 80307
>> ----------------------------------------------------------------------------
>> Unidata WWW Service                        http://www.unidata.ucar.edu/     
>> ****************************************************************************
>>
>
>
>