NAWIPS & GEMPAK Scripts (No changes from version 5.2 to 5.4. PB/NWS 10/96) Prepared by Peggy Bruehl based on earlier work by Ron Miller. ---------------------------------------------------------------------------- This file provides an annotated list of many GEMPAK & NAWIPS Scripts. For a quick summary of the products of these scripts, see the GEMPAK Scripts Quick Reference. These UNIX shell scripts (mostly C-shell) can be run to to plot a variety of graphics using surface and upper air observational data, model output fields, and satellite imagery. Users are encouraged to copy or modify these scripts to address their own particular requirements for data analysis and display. Some of these scripts require command line parameters. These command line parameters and defaults are not documented in this file. You can get help in specifying the command line parameters for any script by using the help option. For example, to get help on command line options for the samap script, type samap help at the UNIX prompt: % samap help ---------------------------------------------------------------------------- Table of Contents * GEMPAK Scripts for Plotting Observational Data * MORE GEMPAK Scripts for Plotting Observational Data * GEMPAK Scripts for Plotting Model Data: * Miscellaneous GEMPAK Scripts * Command Station Models The [ ] symbols indicate a command line option. Defaults may exist. Use the help option for a command line summary. You can resize the GEMPAK window at any point, however you may have to re-run script to refresh the graphic. GEMPAK Scripts for Plotting Observational Data: * samap Plots a standard surface map from surface airways observations. * gale_samap Plots a standard surface map from surface airways observations PLUS plots any obs from ship/buoys. * uamap Plots a map of standard upper air radiosonde observations. * samore Plots a map of extra surface parameters from surface airways observations. * uamore Plots a map of stability indicies and miscellaneous parameters from radiosonde observations. * sazoom and uazoom Zoom in on a surface (sazoom) or upper air (uazoom) map. The cursor will turn into a "+" . Select one corner of the area to zoom in on by positioning the cursor in the desired location and depressing the left mouse button. Hold the mouse button down and drag the cursor so the box describes the desired area and release the mouse button. This is called a rubberband box. * sfc Plots a surface analysis using the Rapid Update Cycle (RUC) model output for the requested area, product and observation hour (or most recent available). Overlays this with a plot of the current surface airways observations. sfc [product]: tmp temperature (F), dwp dew point (F), the Equivalent Potential Temperature (K), rlh Relative Humidity, qdv Moisture Convergence, qtr Moisture Transport, fnt Frontogenesis, div Divergence, vor Relative Vorticity, prf 3-hr Pressure Change, dwpc Dew Point (C), tmpc Temperature (C). * anl Takes a mandatory pressure level, product, and model name as input. Plots an analysis from the specified model initial analysis (must be available) overlayed with a plot of the current upper air observations. anl [level - product] [mdl] [level-product]: sfcpwt sfctmp sfclft sfckin 850dwp 850tmp 850all 850spd 850rlh 850qtr 850fnt 700rlh 700tmp 700fnt 700ipv 500htf 500tmp 500rlh 500abv 500spd 500ipv 400spd 400abv 400ipv 300spd 300abv 300ipv 250spd 250abv 250ipv 200spd 200abv 200ipv 150spd 150abv 100spd (bold denotes default product for each level) [mdl]: lfm eta ngm avn * skewt Plots a Skew-T-ln p diagram for a radiosonde station. More GEMPAK Scripts for Plotting Observational Data: * hodo hodod Plots a hodograph for a radiosonde station with height or pressure as the vertical coordinate, respectively. * sasat Overlays a standard surface map from surface airways observations on top of a satellite image. * sfcsat Display a Rapid Update Cycle surface analysis for the requested area, product and obervation hour (or most recent available) combined with a plot of the current surface airways observations. * uasat Plot a map of standard upper air radiosonde observations over a satellite image. * anlsat Display an analysis of the requested mandatory pressure level and product from the specified model initial analysis (must be available) combined with a plot of the current upper air observations. * skewtm Plot a Skew-T-ln p diagram from model grid point data for any station location. GEMPAK Scripts for Plotting Model Data: * crossx Plot a cross section from sounding data. You pick the cross section line. * xsjet * xsjetm Create a vertical cross section at the location specified. To choose the location of the cross section, move the cursor (+) to the desired location for one end point of the section. Depress the left mouse button and drag the cursor to the other section end point. Release the left mouse button. The cross section is of potential temperature, isotachs, mixing ratio, and ageostrophic circulation vectors. * xscsi cyscsi Create a vertical cross section at the location specified. Choose the location of the cross section once the cursor becomes a (+) (see xsjet). The cross section is of potential temperature, m-sfcs, relative humidity, equivalent potential temperature, and total ciculation vectors. * xsomega Create a vertical cross section at the location specified. Choose the location of the cross section once the cursor becomes a (+) (see xsjet). The cross section is of potential temperature, omega, relative humidity, and ageostrophic ciculation vectors. * xsmixr Create a vertical cross section at the location specified. Choose the location of the cross section once the cursor becomes a (+) (see xsjet). The cross section is of potential temperature, mixing ratio, and total wind barbs. * xsptype Create a vertical cross section at the location specified. Choose the location of the cross section once the cursor becomes a (+) (see xsjet). The cross section is of temperature, relative humidity, and total ciculation vectors. ---------------------------------------------------------------------------- Additional GEMPAK scripts: ([ ] indicate optional parameters) hard copy printout: * uaprint Print upper air maps for 925, 850, 700, 500, 300, and 250 mb and stability indices from rawinsonde observations. Uses the same plotting model as uamap. * local product metafiles: (View created metafiles as a Local Product under ntrans, Select: File; Open: Local Products) These scripts create a metafile on your local workstation that can then be viewed with ntrans. Production of the metafile will take between 1/2 minute (skewtmeta) and up to 5 minutes (qlmeta). o sameta Plot a time sequence of surface maps using the SAMAP and SAMORE station models from surface airways observations and place them in a metafile labeled SA_META. The default is the last six hours for the area specified and ending at the current time. o qlmeta Creates a metafile file containing a time-sequence of Sea-level Pressure & 100-500 mb Thickness and of 500 mb Height & Absolute Vorticity for every forecast time available for the selected numerical model. This script can be run before the entire model is complete to get a "quick look" at the model output (esp. for the AVN). The resulting metafile is labeled QL_[model]_[mdlarea]_META in the ntrans local product list. mdlarea can be: usa United States can Canada car Carribean pac Pacific atl Atlantic epa East Pac. wat West Atl. nam North America sam South America o skewtmeta Plot a Skew-T-ln p diagram for every radiosonde station listed or every radiosonde station in the listed states and place them in a metafile labeled SKEWT_META. ---------------------------------------------------------------------------- Command Station Models TMPF & C: Temperature (F & C) PMSL: Sea-level Pressure WCIF: Wind Chill Index (F) DWPF & C: Dew Point (F & C) PTND: 3-hr Pressure Tendency THTE: Equiv. Pot. Temp. ( K ) WSYM: Weather Symbols STNID: 3-letter Station ID STDZ: Geopotential Height ( dm ) samap GUST ( kt ) TMPF PMSL WSYM O PTND DWPF STNID Visibility ( .1 mi ) samore WCIF 3-hr Precip [.01"] WSYM O 24-hr Precip [.01"] THTE STNID Snow Cover ( inches) uamap TMPC STDZ O DWPC STNID uamore CAPE ( J kg-1 ) LIFT KINX STNID TOTL PWTR (.01" ) CINS ( J Kg-1 ) Command Legend area: Center map on a state (two letter ID) or station ID. A "-" sign will zoom out (eg. wv-) whereas a "+" sign will zoom in. For upper air maps you can also specify: usa (USA) can (CANADA) car (CARRIBEAN) afos (AFOS size USA) nam (N. AMERICA) hour: Hour, hh, in GMT for plot (eg. 12 = 1200 UTC) filter: Determines Data Density: 0 = plot all data, 1 = avoid all overplots of data (enter any rational number, 0 or greater) yymmdd: Date within last two weeks you want to plot. (940126 = 26 January 1994) level: Mandatory isobaric level for plot (1000, 925, 850, 700, 500, 400, 300, 250, 200, 150, 100) stnid: Radiosonde station ID for a vertical profile. Can be entered as a number (eg. 72403 = Dulles) or 3-letter station ID (eg. IAD = Dulles) < help > To get help on a command type help or h after the command name on press (eg. samap help = help on the command samap) < print > Optional: Can be placed at the end of any of the above data plotting scripts to get a B&W printout. ** Printer must be attached to the workstation !!!! **