GEMPAK Comparison Guide

The popular meteorological plotting and analysis tool GEMPAK is no longer formally supported and Unidata is only producing occasional maintenance releases. To help GEMPAK users convert to the Python ecosystem we have compiled this list of GEMPAK programs, commands, parameters, etc. and their Python equivalents. These tables are continually updated as we add more functionality to MetPy as part of an National Science Foundation grant. Where issues have been filed, we have linked to them. Please feel free to tackle of these, comment on their importance in the issue tracker, or file a new issue.

Green rows represent implemented functionality, red is not implemented (with issue listed), blue is uncertain of parity, and white is unevaluated.

GEMPAK Scalar Output/Scalar Grid

GEMPAK Function Description Python/MetPy Equivalent Grid Compatible? Tested against GEMPAK? Units?
ABS(S) Absolute value numpy.absolute Yes Yes Yes
ACOS(S) Arc cosine numpy.arccos Yes Yes Yes
ASIN(S) Arc sine numpy.arcsin Yes Yes Yes
ATAN(S) Arc tangent numpy.arctan Yes Yes Yes
ATN2(S1, S2) Arc tangent numpy.arctan2 Yes Yes Yes
COS(S) Cosine numpy.cos Yes Yes Yes
EXP(S1, S2) Exponential to real S1**S2 Yes Yes Yes
EXPI(S1, S2) Exponential to integer (uses NINT) S1 ** S2.astype(int) Yes Yes Yes
INT(S) Convert to integer S.astype(int) Yes Yes Yes
LN(S) Natural logarithm numpy.log Yes Yes Yes
LOG(S) Base 10 logarithm numpy.log10 Yes Yes Yes
NINT(S) Round to nearest integer numpy.round Yes Yes Yes
SIN(S) Sine numpy.sin Yes Yes Yes
SQRT(S) Square root numpy.sqrt Yes Yes Yes
TAN(S) Tangent numpy.tan Yes Yes Yes
NCDF(S1, S2, S3) Cumulative normal distribution for value, mean, std dev scipy.stats.norm.cdf
INCD(S1, S2, S3) Inverse cumulative normal distribution value given cumulative probability, mean, std dev scipy.stats.norm.ppf(scipy.stats.norm.cdf)
BNCDF(S1, S2, S3, S4) Cumulative binormal mixture distribution given value, mode, left sigma, and right sigma scipy.stats.multivariate_normal
IBNCDF(S1, S2, S3, S4) Inverse cumulative binormal mixture distribution value given cumulative probability, mode, left sigma, and right sigma
PAUB(S1, S2, S3) Combine probabilities (S1 or S2) with dependency parameter S3
ADD(S1, S2) Addition S1 + S2 Yes Yes Yes
MUL(S1, S2) Multiplication S1 * S2 Yes Yes Yes
QUO(S1, S2) Division S1 / S2 Yes Yes Yes
SUB(S1, S2) Subtraction S1 - S2 Yes Yes Yes
BOOL(S) Boolean function
SLT(S1, S2) Less than function S1 < S2 Yes Yes Yes
SLE(S1, S2) Less than or equal to function S1 <= S2 Yes Yes Yes
SGT(S1, S2) Greater than function S1 > S2 Yes Yes Yes
SGE(S1, S2) Greater than or equal to function S1 >= S2 Yes Yes Yes
SBTW(S1, S2, S3) Between two values function S1 > S2 and S1 < S3 Yes Yes Yes
SMAX(S1, S2) Maximum of S1 and S2
SMIN(S1, S2) Minimum of S1 and S2
MASK(S1, S2) Mask S1 based on S2
MISS(S1, S2) Set missing values in S1 to S2
ADV(S, V) Advection metpy.calc.advection Yes Yes Yes
AVG(S1, S2) Average xr.concat([S1, S2], dim='new').mean(axis=0) Yes Yes Yes
AVOR(V) Absolute vorticity metpy.calc.absolute_vorticity Yes Yes* Still needs spherical correction terms Yes
BVSQ(THTA) Brunt-Vaisala frequency squared in a layer metpy.calc.brunt_vaisala_frequency_squared (wrap with np.average to obtain layer average) Yes Yes Yes
CROS(V1, V2) Vector cross product magnitude numpy.cross(V1, V2) Yes Yes No
DDEN(PRES, TMPC) Density of dry air metpy.calc.density(mixing=0) Yes Yes Yes
DDR(S) Partial derivative with respect to R metpy.calc.first_derivative No
DDT(S) Time derivative metpy.calc.first_derivative No
DDX(S) Partial derivative with respect to X metpy.calc.first_derivative Yes Yes Yes
DDY(S) Partial derivative with respect to Y metpy.calc.first_derivative Yes Yes Yes
DEF(V) Total deformation metpy.calc.total_deformation Yes Yes* Still needs spherical correction terms Yes
DIRN(V) Direction relative to north
DIRR(V) Direction relative to grid
DIV(V) Divergence metpy.calc.h_divergence Yes Yes* Still needs spherical correction terms Yes
DIVT(S, V) Divergence tendency (only for cylindrical grids)
DOT(V1, V2) Vector dot product xr.dot(xr.concat([V1u, V1v], dim='component'), xr.concat([V2u, V2v], dim='component'), dims='component') Yes Yes Yes
DTH(S) Partial derivative with respect to theta metpy.calc.first_derivative No
FCNT(S) Coriolis force at grid center (polar) metpy.calc.coriolis_parameter No
FOSB(TMPC, RELH, SPED) Fosberg index (Fire weather index frontogenesis) Issue #636
FRNT(THETA, V) Frontogenesis metpy.calc.frontogenesis Yes Yes* Still needs spherical correction terms Yes
GWFS(S, N) Filter with normal distribution of weights metpy.calc.smooth_gaussian Yes Yes Yes
RDFS(S, DX) GWFS applied for given output resolution
HIGH(S, RADIUS) Relative maxima over a grid
JCBN(S1, S2) Jacobian determinant
KNTS(S) Convert m/s to knots .metpy.convert_units('knots') Yes Yes Yes
LAP(S) Laplacian operator metpy.calc.laplacian Yes Yes - Different Answer Yes
LAV(S) Layer average (2-levels)
LDF(S) Layer difference (2-levels)
LOWS(S, RADIUS) Relative minima over a grid
MAG(V) Magnitude of a vector numpy.linalg.norm Yes Yes Yes
MASS Mass per unit volume in a layer from PRES
MDIV(V) Layer-avg. mass divergence
MIXR(DWPC, PRES) Mixing ratio metpy.calc.mixing_ratio Yes Yes Yes
MRAD(V, LAT, LON, DIR, SPD) Magnitude of radial wind Issue #653
MSDV(S, V) Layer-avg. mass-scalar flux divergence Issue #655
MSFC(V) Pseudo angular momentum (cross-sections) metpy.calc.absolute_momentum No
MTNG(V, LAT, LON, DIR, SPD) Magnitude of tangential wind Issue #653
NORM(V) Normal component (cross-sections) metpy.calc.normal_component No
NMAX(S, ROI) Neighborhood maximum value for a radius of influence
NMIN(S, ROI) Neighborhood minimum value for a radius of influence
PLAT(S) Latitude at each point (polar)
PLON(S) Longitude at each point (polar)
PLCL(PRES, TMPC, DWPC) Pressure of the lifting condensation level metpy.calc.lcl No
POIS(S1, S2) Solve Poisson equation of forcing function with boundary conditions Issue #651
POLF(S) Coriolis force at each point (polar) metpy.calc.coriolis_parameter No
PVOR(S, V) Potential vorticity in a layer metpy.calc.potential_vorticity_baroclinic
metpy.calc.potential_vorticity_barotropic
Yes Yes* Still needs spherical correction terms Yes
PVR3(S, V) 3-D potential vorticity for a level Issue #649
RELH(TMPC, DWPT) Relative humidity metpy.calc.relative_humidity_from_dewpoint Yes Yes Yes
RICH(V) Richardson stability number in a layer Issue #647
ROSS(V1, V2) Rossby number Issue #637
SAVG(S) Average over whole grid numpy.mean Yes Yes Yes
SAVS(S) Average over display subset area of grid numpy.mean Yes Yes Yes
SDIV(S, V) Flux divergence of a scalar Issue #646
SGMX(S) Maximum of S over GAREA numpy.nanmax Yes Yes Yes
SGMN(S) Minimum of S over GAREA numpy.nanmin Yes Yes Yes
SHR(V) Shearing deformation metpy.calc.shearing_deformation Yes Yes* Still needs spherical correction terms Yes
SM5S(S) 5-point smoother metpy.calc.smooth_n_point Yes Yes Yes
SM9S(S) 9-point smoother metpy.calc.smooth_n_point Yes Yes Yes
STAB(TMPC) Lapse rate over a layer in K/km
STR(V) Stretching deformation metpy.calc.stretching_deformation Yes Yes* Still needs spherical correction terms Yes
TANG(V) Tangential component (cross-sections) metpy.calc.tangential_component No
TAV(S) Time average
TDF(S) Time difference
THES(PRES, TMPC) Saturated equivalent potential temperature metpy.calc.saturation_equivalent_potential_temperature Yes Yes Yes
THTA(TMPC, PRES) Potential temperature metpy.calc.potential_temperature Yes Yes Yes
THTE(PRES, TMPC, DWPC) Equivalent potential temperature metpy.calc.equivalent_potential_temperature Yes Yes Yes
THWC(PRES, TMPC, DWPC) Wet bulb temperature metpy.calc.wet_bulb_temperature Yes for individual grid points Yes Yes
TLCL(TMPC, DWPC) Temperature of the lifting condensation level metpy.calc.lcl Yes for individual grid points Yes Yes
TMST(THTE, PRES) Parcel temperature along a moist adiabat metpy.calc.moist_lapse No
TMWK(PRES, TMPK, RMIX) Web bulb temperature in Kelvin metpy.calc.wet_bulb_temperature Yes for individual grid points Yes Yes
DSUB(V1, V2) DIRN(V1) - DIRN(V2)
UN(V) North relative u component
UR(V) Grid relative u component
VN(V) North relative v component
VOR(V) Vorticity metpy.calc.v_vorticity Yes Yes* Still needs spherical correction terms Yes
VR(V) Grid relative v component
WNDX(S1, S2, S3, S4) WINDEX (index for microburst potential) Issue #639
WSHR(V) Magnitude of the vertical wind shear in a layer
XAV(S) Average along a display subset grid row numpy.mean(a, axis=1) Yes Yes Yes
XSUM(S) Sum along a display subset grid row numpy.sum(a, axis=1) Yes Yes Yes
YAV(S) Average along a display subset grid column numpy.mean(a, axis=0) Yes Yes Yes
YSUM(S) Sum along a display subset grid column numpy.sum(a, axis=0) Yes Yes Yes

Vector Output Grid

GEMPAK Function Description Python/MetPy Equilvalent Grid Compatible? Tested against GEMPAK? Units?
AGE(S) Ageostrophic wind metpy.calc.ageostrophic_wind Yes Yes Yes
CIRC(V, S) Circulation (cross-section) Issue #662
DVDX(V) Partial x derivative of V metpy.calc.first_derivative Yes Yes Yes
DVDY(V) Partial y derivative of V metpy.calc.first_derivative Yes Yes Yes
GCIR(LAT, LON) Great circle from point to antipodal point
GCWV(LAT, LON, MASK) Great circle from point to antipodal point with land blocking algorithm MASK will usually be SEA
GEO(S) Geostrophic wind metpy.calc.geostrophic_wind Yes Yes Yes
GRAD(S) Gradient of a scalar metpy.calc.gradient Yes Yes Yes
GWFV(V, N) Filter with normal distribution of weights Issue #664
RDFV(V, DX) GWFV applied for given output resolution
INAD(V1, V2) Inertial advective wind metpy.calc.inertial_advective_wind Yes Yes, slightly different answer in some instances Yes
ISAL(S) Isallobaric wind Issue #641
KCRS(V) Curl of V
KNTV(V) Convert m/s to knots .metpy.convert_units('knots') Yes Yes Yes
LTRN(S, V) Layer-averaged transport of a scalar
NORMV(V) Vector normal wind (cross-section) metpy.calc.normal_component No
QVEC(S, V) Q-vector at a level metpy.calc.q_vector Yes Yes Yes
QVCL(THTA, V) Q-vector of a layer metpy.calc.q_vector (use with metpy.calc.mixed_layer to obtain layer average of each component)
RAD(V, LAT, LON, DIR, SPD) Storm relative radial wind
ROT(angle, V) Coordinate rotation
SMUL(S, V) Multiply a vector's components by a scalar
SM5V(V) 5-point smoother metpy.calc.smooth_n_point Yes Yes Yes
SQUO(S, V) Vector division by a scalar
TANGV(V) Vector tangential wind (cross-section) metpy.calc.tangential_component No
THRM(S) Thermal wind over a layer
TNG(V, LAT, LON, DIR,SPD) Storm relative tangential wind Issue #658
VADD(V1, V2) Add the components of two vectors
VASV(V1, V2) Vector component of V1 along V2
VAVG(V) Average over whole grid
VAVS(V) Average vector over display subset grid
VECN(S1, S2) Create vector from north relative components
VECR(S1, S2) Create vector from grid relative components
VESD(SPD, DIR) Create vector from speed and north-rel direction
VLAV(V) Layer average for a vector
VLDF(V) Layer difference for a vector
VMUL(V1, V2) Multiple the components of two vectors
VQUO(V1, V2) Divide the components of two vectors
VSUB(V1, V2) Subtract the components of two vectors
VLT(V, S) Less than function
VLE(V, S) Less than or equal to function
VGT(V, S) Greater than function
VGE(V, S) Greater than or equal to function
VBTW(V, S1, S2) Between two values function
VMSK(V, S) Mask V based on S

Logical Output Grid

GEMPAK Function Description Python/MetPy Equilvalent Grid Compatible? Tested against GEMPAK? Units?
LT(S1, S2) Returns 1 if S1 < S2
LE(S1, S2) Returns 1 if S1 <= S2
GT(S1, S2) Returns 1 if S1 > S2
GE(S1, S2) Returns 1 if S1 >= S2
GTLT(S1, S2, S3) Returns 1 if S1 > S2 and S1 < S3
GELT(S1, S2, S3) Returns 1 if S1 >= S2 and S1 < S3
GTLE(S1, S2, S3) Returns 1 if S1 > S2 and S1 <= S3
GELE(S1, S2, S3) Returns 1 if S1 >= S2 and S1 <= S3
EQ(S1, S2, S3) Returns 1 if |S1-S2| <= S3
NE(S1, S2, S3) Returns 1 if |S1-S2| > S3
AND(S1, S2, ..., Sn) Returns 1 if all Sx > 0
OR(S1, S2, ..., Sn) Returns 1 if at least one Sx > 0
EOR(S1, S2, ..., Sn) Returns 1 if exactly one Sx > 0
NOT(S) Returns 1 if S ==0

Ensemble Functions

GEMPAK Function Description Python/MetPy Equilvalent Grid Compatible? Tested against GEMPAK? Units?
ENS_SAVG(input_arg1 & input_arg2) Compute the average of a scalar diagnostic field over an ensemble. This average is an equally weighted ensemble mean unless weights have been specified in the GDFILE entry as described above or in input_arg2. Input_arg1 specifies the scalar field to average and is a dynamic argument. Input_arg2 specifies the name of a weight grid to be used in the calculation.
ENS_SSPRD(input_arg1) Compute the spread (sample standard deviation) of a scalar diagnostic field over an ensemble. Input_arg1 specifies the scalar field and is a dynamic argument. The formula used to compute the spread is SQRT { [ SUM ( Xi**2 ) - ( SUM (Xi) ** 2 ) / N ] / ( N -1 ) } where SUM is the summation over the members, Xi are values of the scalar field from individual ensemble members, and N is the number of members. Note that the average of the squared deviations from the mean is computed by dividing the sum of the squared deviations by ( N - 1 ).
ENS_SSUM(input_arg1) Compute the non-weighted sum of the values of a scalar field over an ensemble. Note that this returns the number of members in the ensemble if "input_arg1"" is a logical function that always evaluates to 1 [e.g. gt(2<1)]
ENS_SWSPRD(input_arg1 & input_arg2) Compute the spred, similar to ENS_SSPRD, of a scalar diagnostic field over an ensemble. The spread is weighted and input_arg2 specifies the name of a weight grid to be used in the calculation. NOTE: For uniform weights ens_ssprd (input_arg1) might be expected to give the same result as ens_swsprd(input_arg1 & w1) where w1 is uniform field of 1. This does not happen because of the division by (N-1) in ens_ssprd. The same is also true in comparing ens_vsprd and ens_vwsprd results.
ENS_VAVG(input_arg1 & input_arg2) Compute the average of a vector diagnostic field over an ensemble. This average is an equally weighted ensemble mean unless weights have been specified in the GDFILE entry as described above or in input_arg2. Input_arg1 specifies the vector field to average and is a dynamic argument. Input_arg2 specifies the name of a weight grid to be used in the calculation.
ENC_VSPRD(input_arg1) Compute the spread (sample standard deviation) of a vector diagnostic field over an ensemble. Input_arg1 specifies the vector field and is a dynamic argument. The formula used to compute the variance for each individual component of a vector is the same as for ENS_SSPRD. The variances of the components are added, then the square root is taken to compute the vector spread, which is a scalar.
ENS_VWSPRD(input_arg1 & input_arg2) Compute the spred, similar to ENS_VSPRD, of a vector diagnostic field over an ensemble. The spread is weighted and input_arg2 specifies the name of a weight grid to be used in the calculation. Also, see NOTE for function ENS_SWSPRD.
ENS_SMAX(input_arg1) Compute the maximum of a scalar diagnostic field over an ensemble.
ENS_SMIN(input_arg1) Compute the minimum of a scalar diagnostic field over an ensemble.
ENS_SRNG(input_arg1) Compute the range of a scalar diagnostic field over an ensemble. The range is defined as the difference between the maximum and the minimum.
ENS_PRCNTL(input_arg1 & input_arg2 & input_arg3) At each grid point determine the value of input_arg1 such that it exceeds that found in P% of the members of the ensemble, where P is given by input_arg2. The first argument is a dynamic argument, the second is a static argument. Input_arg3 specifies the name of a weight grid to be used in the calculation. The value of P must range between 0 and 100, and it may vary from grid point to grid point. If the ensemble members are weighted, it is the sum of the weights associated with the order statistics that determines the value corresponding to the percentile. If the percentile falls between two members, the value is determined by interpolation
ENS_MODE(input_arg1 & input_arg2) At each grid point, determine the mode of a weighted ensemble distribution. The mode value is the first one found by a trisection algorithm. The algorithm uses a moving, shrinking data interval in which the weight sums for the data values are made in each of three subintervals. The largest weight sum determines the next interval, with ties decided in the direction of the data mean. The algorithm terminates when only the middle subinterval has a non-zero sum of weights. The weighted average of the data values in this lone subinterval is the mode. Input_arg2 specifies the name of a weight grid to be used in the calculation.
ENS_PROB(input_arg1) Compute the probability of the multi-variate expression given in the input argument. The expression in the input argument is composed of the logical functions AND, OR, and EOR operating on the logical arithmetic comparison functions LT, LE, GT, GE, GTLT, GELT, GTLE, and GELE. ENS_PROB computes the weighted average of the 1 or 0 logical evaluations over the members of the ensemble to yield the relative frequency, or probability, of occurrence.
ENS_CPRB(input_arg1 & input_arg2 & input_arg3 & input_arg4) At each grid point, determine the univariate cumulative probability that the value of a given function is less than or equal to a given value in an ideal ensemble. The input arguments are, respectively, the function to evaluate, the given value whose cumulative probability is sought, the lower bound below which values of the function are impossible, and the upper bound above which values of the function are impossible. If a value for input_arg2 is outside of a given bound, the result is a missing value. The first argument is a dynamic argument; the others are static. The last two are optional. To omit input_arg3 but include input_arg4, specify input as ( input_arg1 & input_arg2 & & input_arg4 ).All arguments may be functions that vary from one grid point to another. For an ideal ensemble, there is a 2/(N+1) probability of values falling outside of the ensemble envelope. The probability density function is piecewise linear with triangular extrapolations beyond the ensemble, limited by any bounding values given. The following example finds the cumulative probability o
ENS_CVAL(input_arg1 & input_arg2 & input_arg3 & input_arg4) At each grid point, determine the value of a given function such that the univariate cumulative probability of values less than or equal to the computed value is the given probability, assuming an ideal ensemble. The input arguments are, respectively, the function to evaluate, the given cumulative probability ranging from 0 to 1, the lower bound below which values of the function are impossible, and the upper bound above which values of the function are impossible. A computed value will never fall outside of a bounding value. The first argument is a dynamic argument; the others are static. The last two are optional. To omit input_arg3 but include input_arg4, specify input as ( input_arg1 & input_arg2 & & input_arg4 ). All arguments may be functions that vary from one grid point to another. If the value of the second argument is not between 0 and 1, inclusive, then a missing value is assigned at that grid point. For an ideal ensemble, there is a 2/(N+1) probability of values falling outside of the ensemble envelope. The probability density function is piecewise linear with triangular extrapolations beyond the ensemble, limited by any bounding values given.

Layer Functions

GEMPAK Function Description Python/MetPy Equilvalent Grid Compatible? Tested against GEMPAK? Units?
LYR_SWTM(input_arg1 | levels) Computes the weighted mean of a scalar field. If the input argument is a single-level parameter or function, the weight is the depth of the layer extending from point halfway between the level and the next lower level and the point halfway between the level and the next upper level. If the input argument is a two-level function, then the weight is the depth of the two-level layer used in calculating the function. In isobaric coordinates, the logarithm of pressure is used to compute the depths. Input_arg1 specifies diagnostic function to average over the layers and output_arg is not required.
LYR_MXMN(argin & fldout [&argout] [|lvls] Computes the following over the specified range of levels... 1.The maximum or minimum value of a scalar quantity. 2.The value of a second output function coincident with the extrema of the input function. It might be a level, such as HGHT or PRES. It might be some other function, e.g. absolute vorticity.
LYR_FVONISFC(fvalu & fisfc & visfc & n & gradflag [|lvls] Computes the value of a function on a specified isosurface of a second function. The function will have the ability to traverse the atmosphere starting from the bottom level upwards or top level downwards and return either the 1st, 2nd, or nth level at which the isosurface threshold value exist.

GUI Programs

Program Description Python/MetPy Equivalent
nfax Displays NMX 6 bit raster fax products
nmap Displays and animates different types of meteorological data on a geographic background.
nsharp Advanced interactive sounding analysis program
ntrans Displays and animates N-AWIPS graphics metafiles in an X-Window.
nwx Displays text products from the Family of Services (FOS) data feed.

Decoders

GEMPAK Function Description Python/MetPy Equilvalent Grid Compatible? Tested against GEMPAK? Units?
dcacars NetCDF ACARS data
dcacft Raw AIREP, PIREP, RECCO, and AMDAR reports
dcairm AIRMET (AIRman's METeorological Information) reports
dccosmic
dccsig Convective signet and convective outlook reports
dcffa Flash flood watch reports
dcffg Flash flood guidance data
dcgmos GFS MOS
dcgrib GRIB
dcgrib2 GRIB2
dchrcn Forecast/advisory reports for tropical depressions, tropical storms, and hurricanes for the Atlantic and Pacific oceans.
dcigdr IGDR data in BUFR format
dcisig SIGMET reports
dclsfc Land surface synoptic reports
dcmetr Raw SAO and METAR reports
dcmsfc Raw buoy, ship, C-MAN, and Coast Guard reports
dcncon Non-convective SIGMET reports
dcncprof NetCDF format profiler and RASS reports from NOAA/FSL
dcnexr2 CRAFT IDD data stream
dcnldn NLDN lightning data reports
dcnmos NGM MOS
dcprof BUFR format profiler reports
dcrdf Regional digital forecast (RDF)
dcredbook Creates displays of Redbook graphic format products
dcscd Supplemental Climatological Data reports
dcshef Raw SHEF reports
dcshef_comet Data feed (LDM) raw SHEF reports
dcstorm Severe storm reports from SPC
dcsuomi NetCDF format SUOMINET data
dcsvrl Severe local storm reports (tornado and severe thunderstorm watch reports)
dctaf Raw terminal aerodrome forecast (TAF) reports
dctama TAMDAR data in BUFR format
dctrop Hurricane/tropical storm reprots
dcuair Upper air sounding data
dcuspln USPLN lightning data reports
dcwarn Flash flood, tornado, and severe thunderstorm warning reports
dcwatch WWUS40 format severe thunderstorm and tornado watch box bulletins
dcwcn Watch county notification reports
dcwcp Tornado and severe thunderstorm Watch Corner Points (WCP) reports
dcwou Watch outline update reports
dcwstm Winter storm reports
dcwtch Tornado and severe thunderstorm watch box reports and watch status reports
dcxmos GFSX MOS

Programs

GEMPAK Function Description Python/MetPy Equilvalent Grid Compatible? Tested against GEMPAK? Units?
acprof Draws profiles of ACARS data
bfr2gp Transfers data from a Jack Woollen BUFR file into GEMPAK sounding and surface data files
bufrenc Processes ASCII input file to produce one or more BUFR output files
cpcgstn Searches for stations located inside specified area from a vector graphics file
gd2ndfd Converts a GEMPAK grid to an NDFD GRIB2 file
gdbiint Interpolates grids from one projection to another
gdcfil Creates a GEMPAK grid file
gdcntr Draws contour lines through a scalar grid matplotlib.axes.Axes.contour No
gdcross Displays a vertical cross section of scalar and/or vector grids
gddelt Deletes grids from GEMPAK grid files
gddiag Computes a scalar/vector diagnostic grid and adds it to the grid file
gdedit Reads grids from a sequential edit file and adds them to a GEMPAK grid file
gdfrzl Generates GFA FZLVLs in VG format from a scalar grid
gdgrib2 Computes a scalar diagnostic grid and adds it to a GRIB2 file
gdgrib Computes a scalar diagnostic grid and adds it to a GRIB file
gdgsfc Computes grid data and interpolates to stations in a GEMPAK surface file
gdinfo Lists information about GEMPAK grid files
gdlist Lists data from a scalar grid
gdmap Plots data from a scalar grid
gdmod Moves grids from one GEMPAK grid file to another
gdobscnt Creates a gridded sampling of the number of surface observations within a specified radius of each grid point.
gdomeg Computes grids of vertical motion and adds them to the grid file
gdplot2 Draws contour lines through scalar grids and/or wind barbs or arrows or streamlines through vector grids. Multiple sets of contours, vectors, and/or streamlines can be generated for each frame.
gdplot3 Draws contour lines through scalar grids and/or wind barbs or arrows or streamlines through vector grids. the program also plots contents of a text file and/or objects. Multiple sets of contours, vectors, streamlines, objects and/or text files can be generated for each frame.
gdplot Draws contour lines through scalar grids and/or wind barbs or arrows through vector grids. Multiple sets of contours and vectors can be generated for each frame.
gdprof Draws profiles of a scalar grid and/or a vector grid.
gdptpdf Raw cumulative probability or probability density at selected point and time at a single level.
gdpvsf Vertically interpolates grid data to an arbitrary surface.
gdradrc Creates a gridded composite of Canadian ASCII products
gdradr Creates a gridded composite of NEXRAD level III products
gdstat Computes statistics on a time series of scalar grids
gdstream Draws streamlines through a vector grid
gdthgt Draws contours and wind barbs or arrows on a time selection at a point within a grid
gdtopo Creates a GEMPAK GRID file from a raster file of topography or land use values.
gdtser Draws a time series of a scalar at a single level
gdtspdf Draws contour or color fill of cumulative probability or probability density as a function of grid values and forecast time at a selected point and level or layer.
gdvint Performs interpolation between vertical coordinates
gdwind Displays a vector grid using wind barbs or arrows
gpanot Will allow the user to place objects at any location on the current graphic device
gpbox Draws a box around a region
gpclear Clears the current graphics device
gpclose Closes the specified graphics output window or file
gpcolor Changes the colors on a color device
gpend Terminates the GEMPLT subprocesses.
gpfax Creates a postscript, GIF, or TIFF file, or an X Windows display from a 6-bit fax file.
gpfront Version of GPMAP that plots map symbols interpreted from ASUS1 bulletins. Forecast front positions from FSUS2 bulletins can be plotted by specifying the forecast hour desired.
gpltln Draws a map, LAT/LON lines with a selected marker, and various image and graphic products
gpmap Draws a map, LAT/LON lines, and various images and graphic products.
gpnexr2 Displays NEXRAD level II products.
gpnids Plots NEXRAD level III products
gprcm Version of GPMAP that plots Radar Coded Message (RCM) bulletins
gpscat Draws scatterometer wave height and wind maps, LAT/LON lines, and various controllable image and graphic customizations.
gptcww Generates the track error watch/warn pilot for the TPC, using breakpoint information contained in a VGF file.
gptext Draws the contents of a text file to the output device.
gptpc Generates four hurricane graphics from the TPC. 1) Wind swath plot 2) Strike probability plot 3) Wind intensity graph 4) Wind speed table
gptraj Computes trajectories for gridded data files
gpvad Plots NEXRAD level III VAD wind profile
gpwarn Version of GPMAP that plots filled county/zone regions from reports which use the universal generic county/zone identified lines.
grphgd Generates a grid either from contours drawn via NMAP product generation, or from a provided VGF file
img2gd Converts an image to a grid. The grid navigation specified may be different than that of the source image.
mkelev NOT FOUND
nagrib2 Converts gridded data in GRIB2 files to GEMPAK gridded data
nagrib Converts gridded data in GRIB files to GEMPAK gridded data
namsnd Transfers model profile output in BUFR to GEMPAK sounding and surface data files.
ndfdg2 Converts NDFD gridded data in GRIB2 files to GEMPAK gridded data
nex2gini Creates as GINI format image composite of NEXRAD level III products
nex2img Creates a GIF format image composite of NEXRAD level III products
nexr2rhi Displays NEXRAD level II vertical cross sections
nexrcomp Creates a GINI format image composite of NEXRAD level III products
oabox Draws a box around an objective analysis region
oabsfc Performs a Barnes objective analysis on surface data
oabsnd Performs a Barnes objective analysis on upper air data
oagrid Creates a GEMPAK grid file which can be used in a Barnes objective analysis program
sector Subsets a GINI satellite image in area and pixel resolution.
sfcfil Creates a new GEMPAK surface file
sfchck Reads a GEMPAK surface data file and produces a table of stations and an indicator showing whether each station reported data at each time in the file
sfcntr Plots surface station data on a map and optionally contours one of the fields being plotted
sfdelt Deletes data from a surface data file
sfdsl604 Lists data from a GEMPAK surface file in a fixed format for use on the AMS DataStreme website
sfedit Adds or changes data in a surface file using a sequential text file
sfgram Draws a meteorogram for surface data
sfl604 Lists data from a GEMPAK surface file in a fixed format
sflist Surface data from a GEMPAK surface data file
sfmap Plots surface station data on a map
sfmod Moves selected surface data from an input surface file to an output surface file
sfstns Modifies the station information in a surface file
sfvgsf Adds or changes data in a surface file using the elements found in a Vector Graphics File
sncfil Creates a new GEMPAK sounding file
sncross Draws cross sections through sounding data
sndelt Deletes data from a sounding data file
sndslist Lists upper air data from a sounding file for specified vertical levels and stations in a format used for the AMS DataStreme web site.
snedit Adds data in a sequential edit file to a sounding file
snhodo Draws a hodograph of upper air data
snlist Lists upper air data from a sounding file for specified vertical levels and stations
snmap Plots sounding data on a map
snmod Moves selected sounding data from an input sounding file to an output sounding file
snprof Draws profiles of upper air data
snstns Modifies the station information in the upper air file
sntser Draws a time series at a sounding station
vg2uka Converts VG files to ASCII files, using the UKMET browsable ASCII format.

Parameters

GEMPAK Description Python/MetPy Equilvalent Grid Compatible? Tested against GEMPAK? Units?
Temperature Parameters
TMPC Temperature in Celsius
TMPF Temperature in Fahrenheit
TMPK Temperature in Kelvin
STHA Surface potential temperature in Kelvin
STHK Surface potential temperature in Kelvin
STHC Surface potential temperature in Celsius
STHE Surface equivalent potential temperature in Kelvin
STHS Surface saturation equivalent potential temperature in Kelvin
THTA Potential temperature in Kelvin
THTK Potential temperature in Kelvin
THTC Potential temperature in Celsius
THTE Equivalent potential temperature in Kelvin
THTS Saturation equivalent potential temperature in Kelvin
TVRK Virtual temperature in Kelvin
TVRC Virtual temperature in Celsius
TVRF Virtual temperature in Fahrenheit
THTV Virtual potential temperature in Kelvin
TDXC Maximum 24 hour temperature in Celsius
TDNC Minimum 24 hour temperature in Celsius
TDXF Maximum 24 hour temperature in Fahrenheit
TDNF Minimum 24 hour temperature in Fahrenheit
T6XC Maximum 6 hour temperature in Celsius
T6NC Minimum 6 hour temperature in Celsius
T6XF Maximum 6 hour temperature in Fahrenheit
T6NF Minimum 6 hour temperature in Fahrenheit
DMAX Daily weather map maximum temperature in Fahrenheit
DMIN Daily weather map minimum temperature in Fahrenheit
SSTC Sea surface temperature in Celsius
SSTF Sea surface temperature in Fahrenheit
LTMP Temperature in Celsius of surface air lifted to 500 or !x mb
Moisture Parameters
DWPC Dew point in Celsius
DWPF Dew point in Fahrenheit
DWPK Dew point in Kelvin
DPDC Dew point depression in Celsius
DPDF Dew point depress in Fahrenheit
DPDK Dew point depression in Kelvin
MIXR Mixing ratio metpy.calc.mixing_ratio
MIXS Saturated mixing ratio metpy.calc.saturation_mixing_ratio
SMXR Surface mixing ratio
SMXS Surface saturated mixing ratio
RELH Relative humidity metpy.calc.relative_humidity_from_dewpoint No
TMWK Wet bulb temperature in Kelvin metpy.calc.wet_bulb_temperature(pressures, temperatures, dewpoints).to('K') No
TMWC Wet bulb temperature in Celsius metpy.calc.wet_bulb_temperature(pressures, temperatures, dewpoints).to('degC') No
TMWF Wet bulb temperature in Fahrenheit metpy.calc.wet_bulb_temperature(pressures, temperatures, dewpoints).to('degF') No
VAPR Vapor pressure in millibars metpy.calc.vapor_pressure
VAPS Saturation vapor pressure in millibars metpy.calc.saturation_vapor_pressure
LHVP Latent heat of vaporization metpy.constants.water_heat_vaporization
PWTR Precipitable water at a given level metpy.calc.precipitable_water
Height Parameters
HGHT Height in meters
HGTM Height in meters
HGTK Height in kilometers
HGTD Height in decameters
HGFT Height in feet
HGFH Height in hundreds of feet
HGFK Height in thousands of feet
HGML Height in miles
DHGT Dry hydrostatic height in meters
MHGT Moist hydrostatic height in meters
STDZ Character standard height convention used on ua charts
RSTZ Numeric standard height convention used on ua charts
ZMSL Estimated height at a pressure level
Pressure and Altimeter Parameters
PRES Station pressure in millibars
PMSL Mean sea level pressure
PALT Surface pressure in millibars from ALTI
ALTI Altimeter setting in inches of mercury
ALTM Altimeter setting converted to millibars
SALT Abbreviated standard altimeter setting
SMSL Abbreviated mean sea level pressure in millibars
SALI Abbreviated altimeter setting in inches of mercury
RMSL First 3 digits left of decimal of PMSL * 10
RSLI First 3 digits left of decimal of ALTI * 100
RSLT First 3 digits let of decimal of ALTM * 10
PTND Pressure tendency
PTSY Graphics symbol for pressure tendency
P03C 3-h numeric pressure change
P03D Pressure tendency and change group, appp
P24C 2y-h numeric pressure change
PANY Returns PMSL if available, if not, returns ALTM
RANY Computes the 3 digit display of pressure
SANY Creates a 3 character string from integral part of PMSL or ALTM
Winds
UWND u wind in m/s
VWND V wind in m/s
UKNT U wind in knots
VKNT V wind in knows
DRCT Wind direction in degrees
SPED Wind speed m/s
SKNT Wind speed knots
SMPH Wind speed mph
PSPD Packed direction and speed m/s
PKNT Packed direction and speed in knots
GUST Wind gusts in knots
GUMS Wind gusts in m/s
PWDR Peak 5 -second wind direction in degrees
PWSP Peak 5-second wind speed in m/s
PWHR Hour of 5-second peak wind
PWMN Minutes of 5-second peak wind
WNML Wind component toward a direction 90 degrees CCW from specified direction
WCMP Wind component toward a specified direction
BARB Barb m/s
BRBM Barb m/s
BRBK Barb knots
BRBS Barb mi/hr
ARRW Arrows m/s
ARRM Arrows m/s
ARRK Arrows knots
DARR Wind direction arrows of uniform length
Lifted Condenstaion Level
TLCL Temperature in Kelvin at the LCL from the given level metpy.calc.lcl
PLCL Pressure in millibar at the LCL from the given level metpy.calc.lcl
Stability Indices
BRCH Bulk Richardson number Issue #628
BRCV BRCH computed using CAPV
BVFQ Brunt-Vaisala frequency in a layer metpy.calc.brunt_vaisala_frequency (use with metpy.calc.mixed_layer to obtain layer average)
BVPD Brunt-Vaisala period in a layer metpy.calc.brunt_vaisala_period (use with metpy.calc.mixed_layer to obtain layer average)
BVSQ Brunt-Vaisala frequency squared in a layer metpy.calc.brunt_vaisala_frequency_squared (use with metpy.calc.mixed_layer to obtain layer average)
CAPE Convective available potential energy metpy.calc.cape_cin
CAPV CAPE computed by using virtual temperature
CINS Convective inhibition metpy.calc.cape_cin
CINV CINS computed using virtual temperature
CTOT cross totals index Issue #633
EQLV equilibrium level metpy.calc.el
EQTV EQLV computed using virtual temperature
KINX K index Issue #631
LAPS temperature lapse rate in a layer Issue #638
LCLP pressure in millibars at the LCL from surface metpy.calc.lcl
LCLT temperature in Kelvin at the lcl from the surface metpy.calc.lcl
LFCT level of free convection by comparing temperature between a parcel and the environment metpy.calc.lfc
LFCV LFCT computed using the virtual temperature
LIFT lifted index Issue #632
LFTV LIFT computed using the virtual temperature Issue #632
LHAN low elevation Haines index Issue #635
MHAN middle elevation Haines index Issue #635
HHAN high elevation Haines index Issue #635
MLMR mean mixed layer MIXR
MLTH mean mixed layer THTA
PWAT Precipitable water for the entire sounding metpy.calc.precipitable_water
RICH Richardson number in a layer Issue #647
SEPA Isentropic pressure thickness in a layer
SHOW Showalter index Issue #630
SHRD wind shear direction in a layer
SHRM wind shear magnitude in a layer
STAB THTA lapse rate in a layer
STAP THTA change with pressure in a layer
SWET SWEAT index Issue #634
TOTL Total totals index Issue #633
VTOT vertical totals index Issue #633
Cloud Parameters
xCLD Character cloud coverage code
TCLD xCLD at maximum cloud coverage
xCLO Fractional cloud coverage
TCLO xCLO at maximum cloud coverage
CLCx Numeric cloud coverage
CLCT CLCx at maximum cloud coverage
CLDS combined cloud coverage short code from three levels
CMBC combined cloud coverage numeric from three levels
CLHx cloud height in hundreds of feet
CLDx combined cloud height and short code
CLDT CLDx at maximum coverage level
CLDB CLDx at the lowest ceiling/layer
COMx numeric combined cloud height and coverage combined
COMT COMx at maximum coverage level
CHCx Numeric combined cloud height and coverage combined
CHDx combined cloud height and short code
CEIL ceiling in hundreds of feet
CFRL fraction of celestial dome covered by all low and mid level clouds from WMO code 2700
CTYL low-level cloud genera from WMO code 0513
CTYM mid-level cloud genera from WMO code 0515
CTYH high-level cloud genera from WMO code 0509
CBAS cloud base height from WMO code 1600
CSYL cloud graphics symbol for CTYL
CSYM cloud graphics symbol for CTYM
CSYH cloud graphics symbol for CTYH
CSYT cloud graphics symbol for first level reporting clouds
CFRT cloud coverage number from CLCT (maximum clouds)
SKYC cloud coverage graphics symbol for CFRT
SKYM sky coverage symbol with wind barbs m/s
SKYK sky coverage symbol with wind barbs knots
XVFR categorical identification of flight rules
Weather Codes
WCOD character weather code
WNUM numeric weather code
WTMO character WMO weather code
WWMO numeric WMO weather code
WSYM graphics weather symbol corresponding to WWMO
PWTH character past weather WMO code or graphics symbol for it
PWWM numeric past weather WMO code
Station Parameters
STID character station identifier
STNM station number
SLAT station latitiude in degrees
SLON station longitude in degrees
SELV station elevation in meters
RANG range in kilometers
AZIM azimuth in kilometers
LATI latitude in degrees from range/azimuth
LONG longitude in degrees from range/azimuth
DELT delta time in seconds
Model Output Statistics
MXMN maximum or minimum temperature in Fahrenheit
TNTF night temperature fcst in Fahrenheit
TNCF night temperature climatology in Fahrenheit
TNAF night temperature anomaly in Fahrenheit
TDYF day temperature fcst in Fahrenheit
TDCF fay temperature climatology in Fahrenheit
TDAF day temperature anomaly in Fahrenheit
CL12 prevailing total sky cover fcst for a 12-hr period
SK12 maximum sustain surface wind spped fcst for 12-hr period
PP06 probability of precip fcst in 6-hr period
PP12 probability of precip fcst in 12-hr period
PP1C probability of precipitation climatology in a 12-hr period
PP1A probability of precipitation anomaly in a 12-hr period
PP24 probability of precipitation fcsr in a 24-hr period
PP2C probability of precipitation climatology in a 24-hr period
PP2A probability of precipitation anomaly in a 24-hr period
QP06 Quantitative precipitation fcst in a 6-hr period
QPX2 Maximum amount of precipitation in inches fcst in a 12-hr period. Values are same as QP12.
QP12 Quantitative precipitation fcst in a 12-hr period
QP24 Quantitative precipitation fcst in a 24-hr period
TS06 Unconditional probability of thunderstorms occurring in a 6-hr period
TS12 Unconditional probability of thunderstorms occurring in a 12-hr period
TS24 Unconditional probability of thunderstorms occurring in a 24-hr period
TC06 Unconditional probability of severe weather occurring in a 6-hr period
TC12 Unconditional probability of severe weather occurring in a 12-hr period
PCPT Categorical forecast of precipitation
POZP Conditional probability of freezing precipitation
POSN Conditional probability of snow
SN06 Categorical forecast of snow amount falling in a 6-hr period
SN12 Categorical forecast of snow amount falling in a 12-hr period
SN24 Categorical forecast of snow amount falling in a 24-hr period
PZ12 Conditional probability of freezing precipitation in a 12-hr period
PS12 Conditional probability of snow in a 12-hr period
PR12 Conditional probability of mixed liquid/frozen precipitation in a 12-hr period
PC12 Categorical forecast of precipitation type in a 12-hr period
FCIG Categorical forecast of ceiling height conditions
FVIS Categorical forecast of visibility conditions
FVSA Categorical forecast of visibility conditions (for new MOS)
OVIS Categorical forecast in plain language of obstructions to vision
WXPB Categorical weather precipitation probability or areal coverage determined by the precipitation parameter having the highest probability or areal coverage in WNUM.
TAF Parameters
TDRC Temporary/probability wind direction in degrees
TSKN Temporary/probability wind speed in knots
TGST Temporary/probability wind gusts in knots
BRGK Gust barb feathered in knots
TCHx Temporary/probability numeric combined cloud height and coverage, as for CHCx
TCEL Temporary/probability ceiling in hundreds of feet, as for CEIL
TSKC Temporary/probability cloud coverage graphics symbol, as for SKYC
TXVF Temporary/probability categorical identification of flight rules, as for XVFR
TWNM Temporary/probability numeric weather code, as for WNUM
TWSY Temporary/probability graphics weather symbol corresponding to TWNM, as for WSYM
TVSB Temporary/probability visibility in statute miles
PPRB Probability for TAF forecast change indicator
VWNM Vicinity numeric weather code, as for WNUM
VWSY Vicinity graphics weather symbol corresponding to VWNM, as for WSYM
TVWN Temporary/probability vicinity numeric weather code, as for WNUM
WSKC Worst case cloud coverage graphics symbol, as for SKYC
WXVF Worst case categorical identification of flight rules, as for XVFR
TPWN Temporary/probability/vicinity numeric weather code, as for WNUM
TPWS Temporary/probability/vicinity graphics weather symbol corresponding to TPWN, as for WSYM
AWNM Prevailing/temporary/probability/vicinity numeric weather code, as for WNUM
AWSY Prevailing/temporary/probability/vicinity graphics weather symbol corresponding to AWNM, as for WSYM
LLWS Low level wind shear forecast flag
MOTV Mountain obscuration threshold value in hundreds of feet
CMSL Ceiling converted to mean sea level in hundreds of feet
MOBS Mountain obscuration threshold met indicator
TCMS Temporary/probability ceiling converted to mean sea level in hundreds of feet
TMOB Temporary/probability mountain obscuration threshold met indicator
WCMS Worst case ceiling converted to mean sea level in hundreds of feet
WMOB Worst case mountain obscuration threshold met indicator
TCTL Temporary/probability low-level cloud genera from WMO Code 0513, as for CTYL
TCSL Temporary/probability cloud graphics symbol for TCTL, as for CSYL
Marine Parameters
WHGT Wave height in meters
WHFT Wave height in feet
WPER Wave period in seconds
HOWW Height of wind wave in meters
POWW Period of wind wave in seconds
HOSW Height of predominant swell wave in meters
POSW Period of predominant swell wave in seconds
DOSW Direction of predominant swell wave in degrees
HOS2 Height of secondary swell wave in meters
POS2 Period of secondary swell wave in seconds
DOS2 Direction of secondary swell wave in degrees
WAV2 Combined wind wave period and height in feet ("2 group")
WAV3 Combined predominant and secondary swell wave direction in tens of degrees ("3 group")
WAV4 Combined predominant swell wave period and height in feet ("4 group")
WAV5 Combined secondary swell wave period and height in feet ("5 group")
WPHM Combined wave period and height in half meters
WVSW Combined swell wave direction, period and height in half meters
SWEL Character combined swell wave direction, period and height in half meters
DAWV Swell wave direction arrows of uniform length
IDTH Thickness of ice on ship in meters
ROIA Rate of ice accretion on ship from WMO Code 3551
IGRO Rate of ice accretion on vessel in salt water in inches per three hours
DIGR Character rate of ice accretion in inches per three hours
SHPD True direction from which ship is moving (for 3 hours before obs) in degrees
SHPK Ship's average speed (for 3 hours before obs) in knots
DASH Ship's true direction arrows of uniform length
PMN1 Lowest 1-minute average pressure in previous hour in mb
PMNT Time of lowest 1-minute average pressure, as hhmm
PKWD Direction of 1-minute peak wind in previous hour in degrees
PKWK Highest 1-minute mean wind speed in previous hour in knots
PKWS Highest 1-minute mean wind speed in previous hour in m/sec
PKWT Time of highest peak 1-minute wind in previous hour, as hhmm
BRPK Wind barb (knots) for highest peak 1-minute wind
Aircraft Parameters
TURB Amount of turbulence
TBSE Base of turbulence in feet
TTOP Top of turbulence in feet
HBOT Base of turbulence in meters
HTOT Top of turbulence in meters
FQOT Frequency of turbulence
TPOT Type of turbulence
TBSY Graphics symbol for turbulence
ICNG Amount of airframe icing
IBSE Base of icing in feet
ITOP Top of icing in feet
HBOI Base of icing in meters
HTOI Top of icing in meters
TPOI Type of icing
ICSY Graphics symbol for icing
WBSE Base of weather in feet
WTOP Top of weather in feet
HBWX Base of weather in meters
HTWX Top of weather in meters
CLC1 Numeric cloud coverage 1
CBS1 Cloud base 1 in feet
CTP1 Cloud top 1 in feet
CB1M Cloud base 1 in meters
CT1M Cloud top 1 in meters
CLC2 Numeric cloud coverage 2
CBS2 Cloud base 2 in feet
CTP2 Cloud top 2 in feet
CB2M Cloud base 2 in meters
CT2M Cloud top 2 in meters
ACRT Aircraft report type
SELV Flight level in meters
FELV Flight level in hundreds of feet
ITSY Icing type symbol
TTSY Turbulence type symbol
TFSY Turbulence frequency symbol
ACTP Character aircraft type
ATP1 Numeric aircraft type
Miscellaneous Parameters
VSBY Visibility in statute miles
VSBK Visibility in kilometers
VSBN Visibility in nautical miles
VSBF Character visibility in fractions of statute miles for visibilities between 0. and 1.
VSBC Character visibility in fractions of statute miles for all visibility numbers
PnnI Precipitation over last nn hours in inches
PnnM Precipitation over last nn hours in millimeters
DPRC Character daily weather map precipitation in inches
PR24 Precipitation over last 24 hours in inches, as sum of four successive 6-hour precip amounts
SNOW Snow depth in inches
SNEW Amount of new snow in inches
SNRT Forecast snow and ice pellet accumulation to watch threshold ratio
SI12 Forecast snow and ice pellet 12-h accumulation in inches
SNIP Snow and ice pellet watch threshold in inches
FZRT Forecast freezing rain accumulation to watch threshold ratio
FZ12 Forecast Freezing rain 12-h accumulation in inches
FZRN Freezing rain watch threshold in inches
WEQS Water equivalent of snow on the ground in inches
HAIL Hail flag
HLSZ Hail size in centimeters
DDEN Density of dry air in kg/(m**3) metpy.calc.density(mixing=0)
PSYM Montgomery stream function in m**2/(100*s**2) metpy.calc.montgomery_streamfunction
HEAT Heat index in Fahrenheit metpy.calc.heat_index
HMTR Humiture (apparent temperature) in Fahrenheit
WCEQ Wind chill equivalent temperature in Fahrenheit metpy.calc.windchill
WCHT Revised wind chill temperature in Fahrenheit
MSUN Duration of sunshine in minutes
FFnn Flash flood guidance for next nn hours in inches
TOST Type of station (manned or automatic)
STIM Report hour and minutes as hhmm
TEXT Undecoded data
SPCL Undecoded special reports
MARK Markers
FOSB Fosberg Index, also called Fire Weather Index
Spacing Parameters
BLNK Plot a blank, not accounted for in FILTER
SPAC Plot a space, accounted for in FILTER
Additional GEMPAK Parameters
ADDSTN Logical variable which indicates whether stations which are in STNFIL, but not already included in the data file, should be added to the data file.
AFOSFL The name of the AFOS graphics file to be displayed.
AIRM The ending valid time for the airmet, the colors for the instrument flight rules, mountain obscuration, turbulence, icing, sustained winds and low-level wind shear, and flags for plotting symbols or airmet type, the end time, and the flight levels on the map.
ANLYSS The average station spacing and the grid extend region.
ANOTLN The line attributes for annotation.
ANOTYP Specifies the fill type for annotation.
AREA The data area. Only data within the area specified will be processed.
ASCT The ending valid time for the AScat wind data, the speed intervals and colors, the wind barb size and width and plotting flags. Skip is a value that indicates how many rows and data points to skip when plotting. The flags include High Wind Speed, Low Wind Speed, KNMI Quality Control Fail, Redundant Data.
ASDI The ending valid time for ASDI data, the mode of the display, the time increments (in minutes going back from the ending time) and the corresponding colors for the ASDI data.
ATCF The initial time for the ATCF forecast tracks, the colors for each model track, the model names, flags for plotting the time, the storm name or number, the forecast wind speeds and markers at each forecast time, and an optional specific storm identifier.
AWPSFL The name of the AWIPS graphics file to be displayed.
BND Specifies the parameters needed for processing bounds areas.
BORDER The color, line type and line width of the background and underground plot.
BOUNDS Specifies the bound area(s) to consider when performing the graph-to-grid procedure.
BOXLIN Line attributes for a box drawn around a region
BUFRFIL The BUFR output file names.
CALIMG Allows the user to select whether to use the calibration values supplied with an image, or raw pixel values.
CATMAP A string that contains "label=value" pairs.
CBTOP User estimated cloud top height in meters. The default value is 8000 m.
CENTER Allows the GD2NDFD user to specify the originating or generating center ID and sub-center ID, as well as the Generating Process/Model ID.
CINT The contour interval, minimum and maximum values and number of digits.
CLEAR A logical variable which determines whether the graphics screen is cleared before plotting.
CLRBAR Specifies the characteristics of a color bar associated with contour fill.
CNTRFL The name of the file containing contour information.
CNTR Specifies the plot attributes for Cell Centroid movement barbs (knots).
CNTRPRM The SFPARM to contour.
COLORS Specifies a list of color numbers.
COLUMN Specifies the number of columns for plotting the contents of an ASCII text file specified by TXTFIL.
COMPRESS A flag to determine whether the output will be written in compressed format.
CONTUR Sets attributes for the contour algorithms.
CPYFIL Identifies the location of the grid navigation and analysis information to be stored in a grid file, as well as an optional subarea.
CSIG The ending valid time for the convective sigmet and convective outlook, the colors for the initial hour (0-hr) CSIG, extrapolated 1-hr CSIG, extrapolated 2-hr CSIG, and outlook, and flags for plotting the 0-hr CSIG sequence number, end time, direction/speed, flight level, intensity, 1-hr CSIG sequence number, and 2-hr CSIG sequence number on the map.
CTLFLG A logical flag which indicates whether control characters are included in a raw surface data set to be decoded.
CTYPE Specifies the contouring algorithms to use.
CURVE A number corresponding to the method to be used to fit the curve.
CXSTNS Defines the x-axis for a cross-section plot.
CYCLE The cycle name tag put into the VGF files.
DATOUT The date and time which will be assigned in the output file.
DATTIM The date and time to be used by GEMPAK programs.
DELTAN The average station spacing in degrees of latitude. The Barnes objective analysis programs use this number to compute weights for data interpolation.
DELTAX The spacing between grid points in the x direction on CED grids. This value is in degrees longitude.
DELTAY The spacing between grid points in the y direction on CED grids. This value is in degrees latitude.
DELZ The user chosen average height difference between pressure levels to be used in the vertical interpolation. The default value is 500 m.
DESIRE The value of the output vertical coordinate variable to which the input variables should be interpolated.
DEVICE Specifies the graphics device.
DISCRETE A string that contains "value1-value2=value3" pairs.
DITHER Used to specify the plotting behavior of the reflectivity intensities.
DLINES Indicates whether the direction of lines (order of points) should be taken into consideration when performing a graph-to-grid contour analysis.
DTAAREA Defines the area over which station data will be input to the Barnes objective analysis.
ECHO Specifies whether to plot the grid box intesities.
EDGEOPTS An option which allows users to specify boundary conditions for the analysis.
EDR Starts with the ending valid time for EDR data followed by corresponding colors for the EDR data over the time limit specified.
ENCY The initial time for the ENCY forecast tracks, the colors for each model track, the model names, flags for plotting the time, the forecast pressures and markers at each forecast time.
EXTEND Specifies the numbers of grid points beyond the GRDAREA which define the grid extend area in the Barnes objective analysis. The first pass is computed on the extend area to reduce edge effects on the GRDAREA.
FAXFIL The name of a 6-bit FAX product.
FFA The ending valid time for the flood watches, the colors for the flash flood and areal flood watches, a flag for plotting the start and stop times, a flag for plotting the zone names for the storms on the map, a flag for plotting the immediate cause for the flooding, and a flag to outline the zone.
FHOUR The forecast hour, e.g., 18, or 24, which defines the "f" value in the BUFRFIL name.
FHR The forecast hour of freezing levels. FHR can be a single hour or a range.
FILTER A logical variable or real number which controls the filtering of data in order to eliminate plotting of overlapping data.
FILTYP The filter type.
FINT The contour fill interval, minimum and maximum values.
FLINE The color and fill type to be used for contour fill.
FXYTBL The FXY table file names.
G2DIAG Allows for detailed GRIB2 message section information, entry-by-entry, to be printed out for selected GRIB2 messages.
G2DRT Specifies the scaling and packing options when encoding the GRIB2 message.
G2IDS A list of integers that are encoded into the GRIB2 Identification Section ( Section 1 ) that identifies the source and type of data packed in the GRIB2 message.
G2IS A list of integers that are encoded into the GRIB2 Indicator Section ( Section 0 ) that identify the discipline and version of the GRIB2 message.
G2PDT Can be used to specify any or all of the values in the Product Definition Template (PDT), describing the grid in the output GRIB2 message.
G2TBLS Allows for specification of the GRIB2 decoding tables.
GAIRM Makes provision to plot the G-AIRMET snapshots in 3-hour time bins.
GAMMA The convergence parameter, is a multiplier for the weight and search radius for passes after the first pass of the Barnes analysis programs.
GAREA The graphics area.
GBDIAG Allows for detailed GRIB message section information, byte-by-byte, to be printed out for selected GRIB messages.
GBFILE The name of the file which contains gridded data in GRIB messages.
GBTBLS Allows for specification of the GRIB decoding tables.
GCENTER Sets the center latitude-longitude on the ETA model domain.
GDATTIM The date/time for the grid.
GDEFIL The name of the grid edit file which will be used to update a grid file.
GDFILE The name of the file which contains gridded data.
GDNUM Allows the user to select grids by number.
GDOUTF The output grid data file name.
GDOUTL The output vertical level in the target vertical coordinate.
GDPFUN Specifies a grid diagnostic function which yields either a scalar or vector quantity.
GFUNC Specifies a grid diagnostic function which yields a scalar quantity.
GGLIMS The parameter which controls the grid value limits and values.
GGVGF The name of the VGF file.
GLEVEL The vertical level for the grid.
GPACK The packing type and the number of bits (or data precision) to be used to pack the grid data.
GPOINT The grid location to be used for the plot.
GRDAREA Specifies the area to be covered by the grid.
GRDHDR A list of the valid grid header flags.
GRDLBL The color number to be used in plotting the grid index numbers.
GRDNAM The parameter name for the grid.
GRDTYP The type of a diagnostic grid.
GSKIP Allows the NDFDG2 user to specify a skip factor in order to reduce the resolution of the resulting GEMPAK grid.
GSPACE Sets the grid spacing of the ETA model domain.
GUESFUN Specifies one or more grid diagnostic functions that compute the first guess grid for a Barnes objective analysis.
GUESS The name of the file which contains the first guess gridded data for objective analysis programs.
GVCORD The vertical coordinate of the grid to be selected.
GVECT Specifies a grid diagnostic function which yields a vector quantity.
GVOUTC The vertical coordinate of the output grid.
HILO Contains the information for plotting relative highs and lows.
HISTGRD Toggles the writing of the graph-to-grid history grid to the GEMPAK grid file.
HLSYM Defines the characteristics for the HILO symbols.
HRCN The ending valid time for the tropical disturbance, the colors for the hurricanes, tropical storms, tropical depressions and directional arrows, the symbols for the hurricanes, tropical storms and tropical depressions, and flags for plotting the center located time, the name and minimum central pressure, the speed, the wind and sea quadrant radii, and the forecast track on the map.
IDNTYP Sets character or numeric station identifiers to be used for input or output.
IJSKIP Used to control subsetting of the internal grid by declaring bounding index values and numbers of points to skip in each index direction, I and J.
IMCBAR Specifies the characteristics of a color bar for images.
IMGFIL The name of an image file.
IMGTIM The date and time to be used.
IMJM Sets the number of grid points for workstation eta in the N-S and E-W direction.
INDXFL The name of the file which contains the GRIB message header information.
INFO The information needed to define the object to be plotted.
INTERP A logical variable which determines whether interpolation between sweeps will occur.
ISIG The ending valid time for the international SIGMET, the colors for the thunderstorms, turbulence, hurricanes, tropical storms, tropical depressions, volcanic ash clouds, marked mountain waves, tropical cyclones, squall lines, CAT, icing, hail, duststorms, sandstorms, cumulonimbus, and low level wind shear, flags for plotting symbols or storm names, the start and end times, the message id, the direction and speed, and the flight level or the central pressure and maximum wind speed associated with a tropical cyclone on the map.
KEYCOL Indicates which contour lines to process.
KEYLINE Indicates which contour lines to process based on VGTYPE and SUBTYP.
KXKY Specifies the size of a grid as two numbers.
LAT Specifies the latitude grid lines to be drawn.
LATLON Specifies the latitude and longitude grid lines to be drawn.
LEVELS Specifies the vertical levels to be extracted from the data set.
LINE The color, line type, line width, line label frequency, smoothing separated by slashes, and flag to suppress small contours.
LOCI The point(s) needed to place the object to be plotted.
LON Specifies the longitude grid lines to be drawn.
LSTALL A flag indicating whether the full contents of a file are to be listed.
LSTPRM Specifies a field to list on the side of the display. Stations are added to the list if they are filtered from the map.
LTNG The ending valid time for lightning data, the time increments (in minutes going back from the ending time) and the corresponding colors for the lightning data, and the positive and negative markers to display.
LUTFIL Specifies a lookup table file name used to enhance the colors for satellite or radar images.
LVLINCR Specifies the freezing level increment.
MAPFIL The name(s) of the map file(s) to be used for maps drawn by GEMPAK programs.
MAP The map color, line type and line width.
MARKER Specifies the marker color, type, size, line width and hardware/software flag.
MAXGRD The maximum number of grids that can be stored in the grid file being created.
MAXTOP Specifies the color and filter attributes for cell top annotations (feet).
MESO Specifies the plot symbol and filter attributes for mesocyclones.
MIXRLN Specifies the color, line type, line width, minimum, maximum, and increment for the background mixing ratio lines on thermodynamic diagrams.
MRGDAT A logical variable indicating whether sounding data is to be merged or unmerged.
MSCALE Specifies the characteristics of a scale legend associated with map projections.
MXDPTH The user estimated mixed layer depth in meters.
NCON The ending valid time for the non-convective sigmet, the colors for the icing, turbulence, duststorm and sandstorm, and volcanic ash, and flags for plotting symbols, the end time, the message id, and the flight levels on the map.
NDVAL The data value to be used where NIDS products report "ND" (none detected).
NPASS Controls the number of passes for the Barnes objective analysis.
NTRACE The number of traces to be drawn in SFGRAM.
OAATTR Contains attributes to use for objective analysis in the graph to grid function.
OABND Specifies the bounds file(s) to use for 'blocking'.
OAGUESS Contains the information to use as a first guess for objective analysis in the graph to grid function.
OLKDAY The day of the extended oulook.
OSCT The ending valid time for the OScat wind data, the speed intervals and colors, the wind barb size and width and plotting flags.
OUTFIL The name of an output satellite image file.
OUTPUT Determines the output devices.
OVERWR Allows the user to specify whether existing GEMPAK grids in a file should be overwritten or left unchanged.
PANEL Specifies the panel location, panel outline color, line type and width.
PDSEXT A logical flag which only becomes applicable when a PDS extension exists in the GRIB message.
PDSVAL Provides a way to enter explicitly PDS numbers identifying a grid by parameter, level, vertical coordinate type, and time.
PIXRES The input for how many pixels and lines to include in the new image.
PLUS Specifies the size and width of a plus sign.
PMAX Defines the maximum possible pressure from which data may be interpolated to the output vertical coordinate.
POSN The position number and the format of the text to be used to plot data in GDMAP.
PRBTYP Specifies the probability type.
PRECSN The binary or decimal precision to be preserved in the packing of the gridded data.
PROJ The map projection, projection angles, and margins.
PTYPE The type of y axis plot to be used, the height-to-width ratio of the plot, and the margins.
QCNTL The quality control threshold values.
QSCT The ending valid time for the QuikScat or ambiguity wind data, the speed intervals and colors, the wind barb size and width and plotting flags.
RADDUR The length of time (in minutes) prior to RADTIM for which data will be used in composites.
RADFIL The name of a radar image file.
RADFRQ tThe update frequency for RADAR composites.
RADINFO Specifies the color for radar site operational status annotations.
RADIUS The search radius (in meters) for which data will be considered.
RADMODE Allows the user to select whether to include radar data from sites operating in (P) precipitation/storm mode, (C) clear air mode, and/or (M) maintainence mode.
RADPARM The Radar parameter to be displayed.
RADTIM The date and time to be used for RADAR composites.
RCMFIL The RCM data file.
REFVEC Specifies the size and location on the screen of the.reference arrow.
REGION Specifies an areal location.
RESPOND A logical variable indicating whether the user will respond interactively to GEMPAK programs.
ROTATE The angle of rotation for the coordinate axes.
RTRAJ A logical variable which determines whether the trajectory will start or end at the specified GPOINT.
SATFIL The name of a satellite image file.
SCALE The scaling factor for the data.
SEARCH Controls the search radius in an objective analysis program.
SFEFIL The name of the surface edit file to be used to update a surface file.
SFFILE The name of the surface data file to be accessed.
SFFSRC The surface file source.
SFOUTF The output surface data file.
SFPARM A list of surface parameters to be used in a surface program.
SFPRMF Specifies the packing information for the surface file to be created.
SGWH The ending valid time for the significant wave height data, the height intervals and colors.
SHAPE The object that the user wishes to plot on the current graphics device.
SHIPFL A logical variable which indicates whether the surface file contains stations which are not at a fixed location, such as moving ships, aircraft, or floating buoys.
SKIP A variable which determines the contour points or plot points to skip.
SKPMIS A logical variable which indicates whether non-reporting stations will be listed.
SNBUFR The name of the BUFR model sounding file to be used as input to create GEMPAK sounding and surface data files using program NAMSND.
SNEFIL The name of the sounding edit file to be used to update a sounding data file.
SNFILE The filename for an upper air data set.
SNOUTF The output sounding data file.
SNPARM A list of upper air parameters to be used in an upper-air program.
SNPRMF Specifies the packing formation for the sounding file to be created.
SOURCE Indicates whether the data used to compute the average station spacing are to be read from a surface or sounding file.
SPLINE A logical for using splines to interpolate the data to height levels.
SQUALL The length of the squall line used for the air and moisture flux calculations.
STARTL The level in the input data file at which to begin the search for the output vertical level in gdpvsf.
STAT The issuing status of the GFA freezing levels.
STATION The station to use in SFGRAM.
STNCOL Specifies the color for the station identifier, time and the parameters specified in STNDEX.
STNDEX The list of stability indices or station parameters for upper-air data.
STNFIL The name of a file which contains station information which includes the character identifier, number, name, state, country, latitude, longitude and elevation for each station.
STNPLT Allows the user to plot station markers and station information.
STNTYP Used to select the data reporting characteristic of a station.
STOPL The level in the input data file at which to end the search for the output vertical level in gdpvsf.
STREAM Controls several parameters dealing with the overall streamline calculation and display.
STRMID Specifies the storm identifier used by programs GPTPC and GPTCWW.
SVRL The ending valid time for the SLS watches, the colors for the severe thunderstorm and tornado (SLS) watches, a flag for plotting the start and stop times, a flag for plotting the county names for the watches on the map, and a flag to outline the county.
SYSTEM The system (storm) speed (m/s) and direction.
TAG Used to identify a group of GFA FZLVL elements.
TAXIS Contains the range, increment and location for labels on a time axis.
TCMG The ending valid time for the tropical disturbance, the colors for the disturbance symbol, the arrows, and the storm danger area, and the name of the center issuing the graphic.
TEXT The size, font, text width and hardware/software flag for graphics text.
THTALN Specifies the color, line type, line width, minimum, maximum, and increment for the background dry adiabats (potential temperature lines) on thermodynamic diagrams
THTELN Specifies the color, line type, line width, minimum, maximum, and increment for the background moist adiabats (equivalent potential temperature lines) on thermodynamic diagrams
TILT The Radar beam elevation/tilt angle.
TIMSTN Contains the maximum number of times to include in a file and the number of stations to be included in addition to the stations in STNFIL.
TITLE The title color, title line, and title string.
TOPOFL The topographic input file.
TRACE Specifications for each trace on the meteogram.
TRAK The ending time for the Altimetric Satellite Ground Track Prediction data, and the color of the prediction.
TROPHT The user estimated tropopause height in meters.
TRPINT The user chosen distance above (and below) the tropopause which is used for layer calculations.
TSTEP Specifies the time step, in minutes, for the calculation of updated position of the parcel within the grid domain.
TVS Specifies the plot symbol and filter attributes for tornado vortex signatures.
TXTCOL Specifies the color number for text, NOAA or NWS logo.
TXTFIL Specifies an ASCII text file to be read and displayed to the current device driver.
TXTLOC Specifies the start location for plotting the contents of an ASCII text file specified by TXTFIL.
TXTYPE Specifies the text attributes.
TYPE Specifies the processing type for the GDPLOT2 GDPFUN parameter.
UKAFIL The intermediate input/output ASCII file.
VCOORD Specifies the vertical coordinate system of the levels to process.
VERCEN Allows the GDGRIB user to specify the contents of bytes 4, 5, 6, and 26 of the GRIB PDS.
VGFILE The name of the Vector Graphics File (VGF) to be displayed or processed.
WARN The ending valid time for the warnings, the colors for the severe thunderstorm, tornado and flash flood warnings, a flag for plotting the start and stop times, a flag for plotting the county names for the warning on the map, a flag to outline the county, and a flag to plot warning polygon.
WATCH The ending valid time for the watches, the colors for the severe thunderstorm and tornado watches and a flag for plotting the start and stop times for the watch on the map.
WAVLEN the wavelength for the gravity or lee wave in kilometers.
WAVSPD The wave speed in m/s. This is used for calculating the SCORER parameter.
WCN The ending valid time for the watch county notification(WCN), the colors for the county bounds, a flag for plotting the start and stop times, a flag for plotting the county names for the WCN on the map, a flag to outline the county or union, a flag to fill the county or union, and a union flag.
WCP The ending valid time for the watches, the colors for the severe thunderstorm and tornado watches and flags for plotting the start and stop times and watch numbers for the watches on the map.
WEIGHT The Barnes weighting parameter.
WIND Specifies the wind symbol, size, width, type, and head size.
WINPOS Specifies the position for plotting winds for vertical profile plots.
WMOHDR Allows specification of a WMO header for a GRIB message.
WOU The ending valid time for the watch outline update (WOU), the colors for the county bounds, a flag for plotting the start and stop times, a flag for plotting the county names for the WOU on the map, a flag to outline the county or union, a flag to fill the county or union, and a union flag.
WSAT The ending valid time for the WindSAT wind data, the speed intervals and colors, the wind barb size and width and plotting flags.
WSPDA The ending valid time for the Altika wind speed data, the speed intervals and colors, and time stamp intervals and colors.
WSPDALT The ending valid time for the altimeter-derived wind speed data, the speed intervals and colors, and time stamp intervals and colors.
WSTM The ending valid time for the winter storms, the colors for the storm warning, watch and advisory, a flag for plotting the start and stop times, a flag for plotting the zone names for the storms on the map, and a flag to outline the zone.
XAXIS Contains the left bound, right bound, label increment, and frequency information.
YAXIS Contains the lower bound, upper bound, label increment, and frequency information.