GEMPAK Manual |
Programs
GDPROF draws profiles of a scalar grid and/or a vector grid.
INPUT PARAMETERS
GPOINT Grid point
GDATTIM Grid date/time
GVCORD Grid vertical coordinate
GFUNC Scalar grid
GVECT Vector grid
GDFILE Grid file
LINE Color/type/width/label/smth/fltr/scflg
MARKER Marker color/type/size/width/hw
BORDER Background color/type/width
PTYPE Plot type/h:w ratio/margins
SCALE Scalar scale / vector scale
XAXIS Xstrt/xstop/xinc/lbl;gln;tck
YAXIS Ystrt/ystop/yinc/lbl;gln;tck
WIND Wind symbol/siz/wdth/typ/hdsz
REFVEC Mag;x;y;txtsiz/font/wdth/HW;labl
WINPOS Wind position
FILTER Filter data factor
TITLE Title color/line/title
PANEL Panel loc/color/dash/width/regn
CLEAR Clear screen flag
TEXT Size/fnt/wdth/brdr/N-rot/just/hw flg
DEVICE Device|name|x size;y size|color type
OUTPUT Output device/filename
THTALN THTA color/dash/width/mn/mx/inc
THTELN THTE color/dash/width/mn/mx/inc
MIXRLN MIXR color/dash/width/mn/mx/inc
PROGRAM DESCRIPTION
GDPROF draws vertical profiles at a point. Data from each
level in the grid file are interpolated to the point.
GPOINT specifies the location for the profile. It may be
entered as a latitude and longitude separated with a semicolon,
as a station character identifier or station number, or as a
set of x and y grid points separated with a semicolon and
preceded with an @. The profile location may also be selected
graphically with the CURSOR command. The location of the
profile is then selected by clicking on a horizontal map in
another GEMPAK XW window.
Profiles are plotted in an animation sequence by specifying more
than one time in GDATTIM.
The vertical coordinate, set in GVCORD, can be PRES, THTA, or
HGHT. The vertical axis type, set in PTYPE, can be LIN,
LOG, KAP, or SKEW. Only temperatures vs. pressure may be
plotted on a skew t plot. The height-to-width ratio of the
plot may be specified in PTYPE following a /. If no value is
entered, a value of 0 will be used.
GFUNC is computed at every level in the data set. No errors
are returned if the function cannot be computed at any level,
except at the surface. If data at the surface cannot be
computed, a warning message is written and plotting continues.
GVECT specifies a vector to be plotted in the right margin.
If GFUNC is a temperature or dewpoint temperature, background
lines of potential temperature, equivalent potential
temperature and mixing ratio can be displayed by specifying
THTALN, THTELN, and MIXRLN. If these background lines are
requested and GFUNC is not a temperature, an error message
is generated.
GVECT specifies a vector grid to be displayed as arrows or
barbs in one of three display columns on the right of the
profile plot. For a wind vector field, the default is north
relative direction. A "/G" indicates grid relative direction.
The number of the column is specified in WINPOS.
If M is entered in WIND, winds will be displayed in m/s unless
the KNOTV operator has been specified in GVECT, in which case the
winds will be displayed in knots. If K is entered in WIND,
the wind is displayed in knots.
EXAMPLES
1. Plot temperature in Celsius on a log P chart at IAD for
the latest time in the file. Label the Y axis from 1000 mb
to 100 mb in increments of 100 mb. Label the X axis from
-40 degrees to +40 degrees in increments of 10 degrees. Add
all the background lines. Plot wind barbs in the right margin.
GPOINT = iad
GDATTIM = last
GVCORD = pres
GFUNC = tmpc
GVECT = wnd
GDFILE = $GEMPAK/data/hrcbob.grd
LINE = 2//3
MARKER = 5//1.25
BORDER = 1
PTYPE = skewt
SCALE = 0
XAXIS = -40/40/10
YAXIS = 1000/100/100/;1
WIND = bm6//2
REFVEC = 10
WINPOS = 1
FILTER = no
TITLE = 1/3
PANEL = 0
CLEAR = yes
TEXT = 1
DEVICE = xw
OUTPUT = n
THTALN = 8/3/1
THTELN = 23/1/1
MIXRLN = 16/10/2
2. Overlay the dewpoint for IAD. Do not redraw the background
lines. Plot the wind in grid relative direction.
GPOINT = iad
GDATTIM = last
GVCORD = pres
GFUNC = dwpc
GVECT = wnd/G
GDFILE = $GEMPAK/data/hrcbob.grd
LINE = 3//3
MARKER = 5//1.25
BORDER = 1
PTYPE = skewt
SCALE = 0
XAXIS = -40/40/10
YAXIS = 1000/100/100/;1
WIND = 0
REFVEC = 10
WINPOS = 1
FILTER = no
TITLE = 0
PANEL = 0
CLEAR = no
TEXT = 1
DEVICE = xw
OUTPUT = n
THTALN = 0
THTELN = 0
MIXRLN = 0
3. Clear the screen and plot the vorticity at BWI in
color number 2 using a dashed line. Do not plot any
winds. Have the graph plot in a second window named
"profile2".
GPOINT = bwi
GDATTIM = last
GVCORD = pres
GFUNC = vor(wnd)
GVECT =
GDFILE = $GEMPAK/data/hrcbob.grd
LINE = 2/2/3
MARKER = 0
BORDER = 1
PTYPE = log
SCALE = 5
XAXIS =
YAXIS = 1000/300/100
WIND =
REFVEC = 10
WINPOS = 1
FILTER = no
TITLE = 1/3
PANEL = 0
CLEAR = yes
TEXT = 1
DEVICE = xw|profile2
OUTPUT = t
THTALN = 0
THTELN = 0
MIXRLN = 0
ERROR MESSAGES
[GDPROF +2] WARNING: ... not found. CONTINUING---
[GDPROF +1] Invalid parameter for background lines.
[GDPROF -1] Fatal error initializing TAE.
[GDPROF -2] Fatal error reading TAE parameters.
[GDPROF -3] Fatal error initializing GEMPLT.
[GDPROF -4] Input for GPOINT is invalid.
[GDPROF -5] Input for GDATTIM is invalid.
[GDPROF -6] Input for GVCORD is invalid.
[GDPROF -7] Input for PTYPE is invalid.
[GDPROF -8] Error defining graph coordinates.
[GDPROF -9] No points found for profile.
[GDPROF -10] There are no levels at this time.
[GDPROF -11] Input ... for YAXIS is invalid.
[GDPROF -12] Input ... for XAXIS is invalid.