[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

20010102: 20001228: GEMPAK trajectories



Pat,

Assuming you have created the ascii files from gdlist for urel, vrel, lat and 
lon
at each time that the program will need to interpolate the frajectory points,
I'm not sure where in the process things are getting fouled up.

I think that it might be better to try another approach here using a
single program that can plot trajectories from a grid file. This eliminates
the need to produce a surface format file of data points calculated for
the trajectory to be plotted using sfmap.

I have created a tarfile in ~gbuddy/nawips-5.6/contrib called gptraj_sol.tar.Z.
I know you are still running GEMPAK 5.4, but the gptraj_gf program in the 
tarfile
is linked directly to the gf driver so that it won't have to deal with version
differences. Hopefully you have the current runtime libraries necessary.

Create a clean directory, such as $NAWIPS/gptraj and download the tarfile above.
Unpack the contents of the tarfile in the clean directory with:

zcat gptraj_sol.tar.Z | tar xvf -

I have included all the necessary files from GEMPAK 5.6 in the directory
so you should be able to run the programunder Solaris using your
existing $GEMMAPS. You should execute the program from the same directory
so that you will find the parameter, help, pdf and err files.

The program will use GPOINT to specify the trajectory point (either
a station id or lat;lon etc). Specify GDATTIM as the model run time
to be used, eg 990313/1200 or simply f000. 

As an example:

% gptraj_gf
 GEMPAK-GPTRAJ>l
 MAP      = 1
 GAREA    = uslcc
 PROJ     = lcc
 SATFIL   =  
 RADFIL   =  
 LATLON   = 0
 PANEL    = 0
 TITLE    = 1/-1/Forward Trajectories
 TEXT     = .9/22/1/hw
 CLEAR    = yes
 DEVICE   = gf|traj.gif
 GDFILE   = modelgridfile.gem
 GVECT    = wnd
 GLEVEL   = 10
 GVCORD   = hght
 GDATTIM  = f000
 MARKER   = 6/15/6/1
 LINE     = 7
 RTRAJ    = no
 GPOINT   = 40;-105
 GEMPAK-GPTRAJ> run
Enter <cr> to accept parameters or type EXIT:
 start/stop/nmin 010103/1200V000     010103/1800V006              360
 start/stop/nmin 010103/1800V006     010104/0000V012              360
 start/stop/nmin 010104/0000V012     010104/0600V018              360
 start/stop/nmin 010104/0600V018     010104/1200V024              360
 start/stop/nmin 010104/1200V024     010104/1800V030              360
 start/stop/nmin 010104/1800V030     010105/0000V036              360
 start/stop/nmin 010105/0000V036     010105/0600V042              360
 start/stop/nmin 010105/0600V042     010105/1200V048              360
 start/stop/nmin 010105/1200V048     010105/1800V054              360
 start/stop/nmin 010105/1800V054     010106/0000V060              360
 Ending location 36.45406;-94.17933  
 Parameters requested: MAP,GAREA,PROJ,SATFIL,RADFIL,LATLON,PANEL,TITLE,TEXT,
 CLEAR,DEVICE,GDFILE,GVECT,GPOINT,GLEVEL,GVCORD,GDATTIM,MARKER,LINE,RTRAJ.
 GEMPAK-GPTRAJ>

 GEMPAK-GPTRAJ> MAP = 0
 GEMPAK-GPTRAJ> GPOINT = STL
 GEMPAK-GPTRAJ> run

 <snip>

 GEMPAK-GPTRAJ> exit


The program will find the available forecast times for the model run,
and produce a trajectory (or back trajectory if RTRAJ=YES) using GPOINT.
The program is directly linked to the GF device, so upon exit, you will have 
a gif format file (called traj.gif in the above example). There is no
running of gpend, since there is no message queue associated- so you will
need to run for all of your points before exiting the program.

You can specify any GVECT that can be computed from your data (eg WND, GEO, AGE 
etc),
using the GLEVEL and GVCORD as appropriate for the field. phelp gptraj is 
available
as well.

If you don't have the current runtime libraries, the program won't run. In
that case, type "ldd gptraj_gf" to find out what libraries you need.

Steve Chiswell
Unidata user Support


>From: "Patrick S. Market" <address@hidden>
>Organization: UCAR/Unidata
>Keywords: 200101022020.f02KKco12593

>Steve:
>
>I have included output from my latest go at it.  I have checked the
>parameters that you suggested, and added the following line to my 
>create_forw.gem file:
>
>  echo $MAX $LAT $LON $TITLE $GEMTBL $MODEL
>
>with the following output:
>
> 24 35 -105 990313/1200 /asp/met/nawips/gempak5.4/tables mass
>  [FL 1018] forw_mass.gem
>  [SF -2]  File forw_mass.gem could not be opened.
> trajectory_gem  OPEN_GEM_SHIP: Created new ship file forw_mass.gem
> 24 40 -105 990313/1200 /asp/met/nawips/gempak5.4/tables mass
> 24 35 -95 990313/1200 /asp/met/nawips/gempak5.4/tables mass
> 24 40 -95 990313/1200 /asp/met/nawips/gempak5.4/tables mass
>
>The model run only goes out to 24 hours, and I requested 4 lat/lon
>points.  There is a traj.pack file in /asp/met/nawips/gempak5.4/tables
>and a forw_mass.gem does get created (210432 bytes) in the $TRAJECTORY
>directory, in spite of the error messages.
>
>The same error messages then ensue from plot_forward.gem.
>
>I believe that nothing is being written to the forw_mass.gem file.  I
>removed the 'rm forw_($MODEL).gem' line and reran create_forw.gem.
>The file was not update, even after I changed the lat/lon
>requests.  I am looking at the trajectory_gem.c code after I send this.
>
>Thanks.
>
>Pat
>
>> 
>> 
>> Patrick,
>> 
>> The sfmap output below says that there is no data in the file.
>> 
>> It could be that create_forward.gem is not suceeding if file naming 
>> isn't correct, or if the packing file isn't found etc.
>> You will see that the create_forw.gem runs:
>> trajectory_gem $MAX $LAT $LON $TITLE $GEMTBL/traj.pack forw_${MODEL}.gem $MO
> DEL
>> 
>> You might want to echo out those values and make sure that the program
>> is finding the input data as well as the packing file etc.
>> 
>> Steve Chiswell
>> Unidata User Support
>> 
>> 
>> >From: "Patrick S. Market" <address@hidden>
>> >Organization: UCAR/Unidata
>> >Keywords: 200012282236.eBSMaRo24648
>> 
>> >Hello.
>> >
>> >I have run into a snag in running the GEMPAK trajectory program.  I have
>> >worked my way through the difficulties in the create_forw.gem program,
>> >but am having difficulty with the plot_forward.gem program.  Typical output
>> >goes like so...
>> >
>> >
>> >
>> >
>> >Creating process: gplt for queue 1550
>> > AREA      Data area                         dset
>> > GAREA     Graphics area                     20;-132;55;-68
>> > SATFIL    Satellite image filename(s)       
>> > RADFIL    Radar image filename(s)           
>> > SFPARM    Surface parameter list            arrw;hour=24
>> > DATTIM    Date/time                         all
>> > SFFILE    Surface data file                 forw_mass.gem
>> > COLORS    Color list                        2;0
>> > MAP       Map color/dash/width              0
>> > LATLON    Line color/dash/width/label/inc   0
>> > MARKER    Marker color/type/size/width/hw   5/15/1/3
>> > TITLE     Title color/line/title            0
>> > CLEAR     Clear screen flag                 n
>> > PANEL     Panel loc/color/dash/width/regn   0
>> > DEVICE    Device|name|x size;y size|color   gf|forw_mass.gif
>> > PROJ      Map projection/angles/margins|dr  ced
>> > FILTER    Filter data factor                0
>> > TEXT      Text size/font/width/hw flag      0
>> > SKPMIS    Skip missing data flag            NO
>> > LUTFIL    Enhancement lookup table filenam  
>> > STNPLT    Txc|mkc/mkt/mks/mkwd/hwfl|stnfil  
>> > Parameters requested: AREA,GAREA,SATFIL,RADFIL,SFPARM,DATTIM,SFFILE,COLORS
> ,
>> > MAP,LATLON,MARKER,TITLE,CLEAR,PANEL,DEVICE,PROJ,FILTER,TEXT,SKPMIS,LUTFIL,
>> > STNPLT.
>> >Creating process: gf for queue 1551
>> > GEMPAK-SFMAP> GEMPAK-SFMAP> GEMPAK-SFMAP> GEMPAK-SFMAP> GEMPAK-SFMAP> GEMP
> AK-
>> > SFMAP> GEMPAK-SFMAP> GEMPAK-SFMAP> GEMPAK-SFMAP> GEMPAK-SFMAP> GEMPAK-SFMA
> P> 
>> > GEMPAK-SFMAP> GEMPAK-SFMAP> GEMPAK-SFMAP> GEMPAK-SFMAP> GEMPAK-SFMAP> GEMP
> AK-
>> > SFMAP> GEMPAK-SFMAP> GEMPAK-SFMAP> GEMPAK-SFMAP> [TI -5]  There are no tim
> es 
>> > in the data set.
>> > Parameters requested: AREA,GAREA,SATFIL,RADFIL,SFPARM,DATTIM,SFFILE,COLORS
> ,
>> > MAP,LATLON,MARKER,TITLE,CLEAR,PANEL,DEVICE,PROJ,FILTER,TEXT,SKPMIS,LUTFIL,
>> > STNPLT.
>> > GEMPAK-SFMAP> GEMPAK-SFMAP> GEMPAK-SFMAP> GEMPAK-SFMAP> GEMPAK-SFMAP> GEMP
> AK-
>> > SFMAP> GEMPAK-SFMAP> GEMPAK-SFMAP> GEMPAK-SFMAP> [TI -5]  There are no tim
> es 
>> > in the data set.
>> >
>> >This then repeats 6-7 times (the number of times to calculate over?)
>> >
>> >Any ideas?
>> >
>> >Thanks.
>> >
>> >Pat
>> >
>> >-- 
>> >========================================================================
>> >Dr. Patrick S. Market              
>> >Dept. of Soil & Atmospheric Sciences 
>> >University of Missouri-Columbia                 Phone:  (573) 882 - 1496  
>> >203 Gentry Hall                                 Fax:    (573) 884 - 5133 
>> >Columbia, MO 65211 USA          E-mail: address@hidden
>> >========================================================================
>> >
>> 
>> ****************************************************************************
>> Unidata User Support                                    UCAR Unidata Program
>> (303)497-8644                                                  P.O. Box 3000
>> address@hidden                                   Boulder, CO 80307
>> ----------------------------------------------------------------------------
>> Unidata WWW Service                        http://www.unidata.ucar.edu/     
>> ****************************************************************************
>> 
>
>
>-- 
>========================================================================
>Dr. Patrick S. Market              
>Dept. of Soil & Atmospheric Sciences 
>University of Missouri-Columbia                 Phone:  (573) 882 - 1496  
>203 Gentry Hall                                 Fax:    (573) 884 - 5133 
>Columbia, MO 65211 USA          E-mail: address@hidden
>========================================================================
>