GEMPAK Manual | Programs

GPTRAJ computes trajectories for gridded data files.

INPUT PARAMETERS
 
    MAP       Map color/dash/width/filter flag
    GAREA     Graphics area
    PROJ      Map projection/angles/margins|drop flag
    SATFIL    Satellite image filename(s)
    RADFIL    Radar image filename(s)
    LUTFIL    Enhancement lookup table filename
    IMCBAR    Color/ornt/anch/x;y/ln;wd/freq|text_info
    LATLON    Line color/dash/width/freq/inc/label/format
    PANEL     Panel loc/color/dash/width/regn
    TITLE     Title color/line/title
    TEXT      Size/fnt/wdth/brdr/N-rot/just/hw flg
    CLEAR     Clear screen flag
    DEVICE    Device|name|x size;y size|color type
    GDFILE    Grid file
    GVECT     Vector grid
    GPOINT    Grid point
    GLEVEL    Grid level
    GVCORD    Grid vertical coordinate
    GDATTIM   Grid date/time
    TSTEP     Time step (minutes)
    MARKER    Marker color/type/size/width/hw
    LINE      Color/type/width/label/smth/fltr
    RTRAJ     Back/Reverse trajectory flag
 
 

PROGRAM DESCRIPTION
 
   GPTRAJ computes horizontal trajectories from user supplied vector
   functions on constant grid level/coordinates.

   Trajectories are computed using a user specified GPOINT as the
   starting or ending location (see RTRAJ). The advected GPOINT
   location is determined by interpolating the gridded vector
   components to the trajectory location using the nearest data time.
   The trajectory location is updated iteratively through the forecast
   times using TSTEP over the time range specified in GDATTIM.

   TSTEP is the time step of integration (in minutes).

   GPOINT may be specified using the traditional methods for
   denoting a single point, or may be intered as a semicolon
   separated lat/lon pair of min/max/increment.

   MARKER is used to define the marker characteristics used to locate
   the trajectory point at times matching the available grid fields.

   LINE is used to define the line type which connects the computed
   trajectory locations.

   RTRAJ is a logical flag to specify whether a reverse/back trajectory
   will be computed. If RTRAJ is NO, a forward trajectory is computed.
   If RTRAJ is YES, a backwards trajectory will be computed from GPOINT.

   The direction of the last data point is marked with a directional arrow.


EXAMPLES
 
    1.  Draw a Lambert Conformal map of the US.
        Using the most recent ETA model output, display a trajectory
        starting at DEN using the 10m wind for the time range F024-F060.

         MAP      = 1
         GAREA    = uslcc
         PROJ     = lcc
         SATFIL   =
         RADFIL   =
         LATLON   = 0
         PANEL    = 0
         TITLE    = 1/-1/Trajectory from DEN
         TEXT     = 1/22/1/hw
         CLEAR    = yes
         DEVICE   = xw
         GDFILE   = eta
         GVECT    = wnd
         GPOINT   = den
         GLEVEL   = 10
         GVCORD   = hght
         GDATTIM  = f024-f060
         MARKER   = 6/15/6/1
         LINE     = 7
         RTRAJ    = no
         TSTEP    = 1

    2.  Repeat the map above using all forecast times in
        the ETA model run for a back trajectory.

         MAP      = 1
         GAREA    = uslcc
         PROJ     = lcc
         SATFIL   =
         RADFIL   =
         LATLON   = 0
         PANEL    = 0
         TITLE    = 1/-1/Backward Trajectory for DEN
         TEXT     = 1/22/1/hw
         CLEAR    = yes
         DEVICE   = xw
         GDFILE   = eta
         GVECT    = wnd
         GPOINT   = den
         GLEVEL   = 10
         GVCORD   = hght
         GDATTIM  = fall
         MARKER   = 6/15/6/1
         LINE     = 7
         RTRAJ    = yes
         TSTEP    = 1

    3.  Display all 700mb backward trajectories for points ending at
        lat/lon pairs every 10 degrees across the domain using all
        forecast times in the ETA model run. Set the time step to
        30 minutes.

         MAP      = 1
         GAREA    = uslcc
         PROJ     = lcc
         SATFIL   =
         RADFIL   =
         LATLON   = 0
         PANEL    = 0
         TITLE    = 1/-1
         TEXT     = 1/22/1/hw
         CLEAR    = yes
         DEVICE   = xw
         GDFILE   = eta
         GVECT    = wnd
         GPOINT   = 30/50/10;-120/-70/10
         GLEVEL   = 700
         GVCORD   = pres
         GDATTIM  = fall
         MARKER   = 6/15/6/1
         LINE     = 7
         RTRAJ    = yes
         TSTEP    = 30


ERROR MESSAGES
 
    [GPTRAJ  -1]    Fatal error initializing TAE.
    [GPTRAJ  -2]    Fatal error reading TAE parameters.
    [GPTRAJ  -3]    Fatal error initializing GEMPLT.