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

20011031: 20011030: using stnplt to add a TC marker

This archive contains answers to questions sent to Unidata support through mid-2025. Note that the archive is no longer being updated. We provide the archive for reference; many of the answers presented here remain technically correct, even if somewhat outdated. For the most up-to-date information on the use of NSF Unidata software and data services, please consult the Software Documentation first.

  • Subject: 20011031: 20011030: using stnplt to add a TC marker
  • Date: Thu, 01 Nov 2001 08:47:27 -0700

Andrew,

The STNPLT parameter uses the station table, not a decoded
gempak surface file. As I mentioned yesterday, use your
station table with STNPLT. Or, use the MARK
parameter in the SFPARM variable when plotting your decoded
gempak file.

Steve Chiswell
Unidata User Support


>From: Andrew Levine <address@hidden>
>Organization: UCAR/Unidata
>Keywords: 200111010920.fA19JP115561

>
>Hi again,
>
>Thanks for your help.  I am sorry to keep asking, but I am still having
>trouble using stnplt.  I thought I set up the station table and file
>correctly, but when I try to print out the marker, it still does not
>work, I get an error message that says there are no data records in the
>table file.  This is what I did:
>I set up a station table as follows:
>
>F0000000  99999 GINGER STATION                   -- --  1020  17180     0
>
>Next, I created a data file for lat, lon and pressure:
>
>PARM = SLAT;SLON;PRES
>
>STID     DATTIM         SLAT    SLON    PRES
>F0000000 970921/1200    10.3   174.7    1007
>F0000000 970921/1800    10.3   174.4    1007
>F0000000 970922/0000    10.2   173.7    1007
>F0000000 970922/0600    10.2   173.0    1004
>F0000000 970922/1200    10.2   171.8    1004
>F0000000 970922/1800    10.2   170.3    1002
>F0000000 970923/0000    10.6   168.8    1002
>F0000000 970923/0600    11.2   167.4    1002
>F0000000 970923/1200    12.1   166.3    1000
>F0000000 970923/1800    12.9   165.3    1000
>F0000000 970924/0000    13.7   164.5     997
>F0000000 970924/0600    14.6   163.8     994
>F0000000 970924/1200    15.3   163.2     991
>F0000000 970924/1800    16.0   162.6     987
>F0000000 970925/0000    16.9   162.1     984
>F0000000 970925/0600    17.9   161.6     980
>
>Next, I used sfcfil to create the surface file:
> SFOUTF    Output surface file               /home/ua/levine/gingertrk2.gem
> SFPRMF    Surface parameter packing file    slat;slon;pres
> STNFIL    Station information file          /home/ua/levine/gingertrkgem2.txt
> SHIPFL    Ship data file flag               n
> TIMSTN    Times/additional stations         40/1
> SFFSRC    Surface file source                
> Parameters requested: SFOUTF,SFPRMF,STNFIL,SHIPFL,TIMSTN,SFFSRC.
>
>Then, I used sfedit to bring the data into the data file:
>
> SFEFIL    Surface edit file                 /home/ua/levine/gingertrkgem.txt
> SFFILE    Surface data file                 /home/ua/levine/gingertrk2.gem
>
>when I did sflist, I got a listing like this:
>
> SFFILE    Surface data file                 /home/ua/levine/gingertrk2.gem
> AREA      Data area                         dset
> DATTIM    Date/time                         all
> SFPARM    Surface parameter list            slat;slon;pres
> OUTPUT    Output device/filename            t
> IDNTYP    STNM or STID                      STID
> Parameters requested: SFFILE,AREA,DATTIM,SFPARM,OUTPUT,IDNTYP.
> GEMPAK-SFLIST>r
> PARM = SLAT;SLON;PRES                                                        
>    
>
>    STN    YYMMDD/HHMM      SLAT     SLON     PRES
> F0000000  970921/1200     10.30   174.70  1007.00
> F0000000  970921/1800     10.30   174.40  1007.00
> F0000000  970922/0000     10.20   173.70  1007.00
> F0000000  970922/0600     10.20   173.00  1004.00
> F0000000  970922/1200     10.20   171.80  1004.00
> F0000000  970922/1800     10.20   170.30  1002.00
> F0000000  970923/0000     10.60   168.80  1002.00
> F0000000  970923/0600     11.20   167.40  1002.00
> F0000000  970923/1200     12.10   166.30  1000.00
> F0000000  970923/1800     12.90   165.30  1000.00
> F0000000  970924/0000     13.70   164.50   997.00
> F0000000  970924/0600     14.60   163.80   994.00
> F0000000  970924/1200     15.30   163.20   991.00
> F0000000  970924/1800     16.00   162.60   987.00
> F0000000  970925/0000     16.90   162.10   984.00
> F0000000  970925/0600     17.90   161.60   980.00
>
>When I ran gdplot, the marker did not print out and I got the following
>error message:
>
> GDFILE    Grid file                         /home/ua/levine/gingerpvdat.gem
> GDATTIM   Grid date/time                    970922/1200
> GLEVEL    Grid level                        350
> GVCORD    Grid vertical coordinate          thta
> PANEL     Panel loc/color/dash/width/regn   0
> SKIP      Skip_cntr/skip_plt_x;skip_plt_y   /1;1
> SCALE     Scalar scale / vector scale       6/-.9999
> GFUNC     Scalar grid                       pvth
> CTYPE     Contour type:  C/F                c/f
> CONTUR    Subbox/smooth                     5
> CINT      Contour interval/min/max          1
> LINE      Color/type/width/label/smth/fltr  32/1/2
> FINT      Fill interval/min/max             50/1
> FLINE     Fill colors                       31;2
> HILO      Color/symbol/rng/rad/cnt/intp     32/H#1;L#1///5
> HLSYM     HILO txt size/posn/font/wdth/hw   //23
> CLRBAR    Color/ornt/anch/x;y/ln;wd/freq     
> GVECT     Vector grid                       wnd
> WIND      Wind symbol/siz/wdth/typ/hdsz     am8/.5///.5
> REFVEC    Mag;x;y;txtsiz/font/wdth/HW;labl  10;.80;.02
> TITLE     Title color/line/title             
> TEXT      Size/fnt/wdth/brdr/N-rot/just/hw   
> CLEAR     Clear screen flag                 n
> GAREA     Graphics area                     0;135;50;-175
> PROJ      Map projection/angles/margins|dr  ced
> MAP       Map color/dash/width/filter flag   
> LATLON    Line color/dash/width/label/inc    
> DEVICE    Device|name|x size;y size|color   xw
> STNPLT    Txtc/txt attr|marker attr|stnfil  /|32/23/2/2|/home/ua/levine/ginge
> rtrk2.gem
> Parameters requested: GDFILE,GDATTIM,GLEVEL,GVCORD,PANEL,SKIP,SCALE,GFUNC,
> CTYPE,CONTUR,CINT,LINE,FINT,FLINE,HILO,HLSYM,CLRBAR,GVECT,WIND,REFVEC,TITLE,
> TEXT,CLEAR,GAREA,PROJ,MAP,LATLON,DEVICE,STNPLT.
> GEMPAK-GDPLOT>r
>Creating process: xw for queue 603
>
> PARAMETERS FOR GDPLOT: 
>
> Grid file: /home/ua/levine/gingerpvdat.gem                                   
>                                                               
> TIMES:    970922/1200                         
>
> GRID IDENTIFIER: 
>    LEVEL VCORD PARM                            WPARM
>      350  THTA PVTH                            WND                        
>
> GAREA:  0;135;50;-175                                                        
>                                                            
> [FL -10]  No data records in table file /home/ua/levine/gingertrk2.gem.
> Parameters requested: GDFILE,GDATTIM,GLEVEL,GVCORD,PANEL,SKIP,SCALE,GFUNC,
> CTYPE,CONTUR,CINT,LINE,FINT,FLINE,HILO,HLSYM,CLRBAR,GVECT,WIND,REFVEC,TITLE,
> TEXT,CLEAR,GAREA,PROJ,MAP,LATLON,DEVICE,STNPLT.
>
>Did I set up the station table incorrectly?  
>
>Thanks again for all of your help!
>
>Andrew
>On Wed, 31 Oct 2001, Unidata Support wrote:
>
>> 
>> Andrew,
>> 
>> My suggestion was to create a "STATION TABLE" for plotting
>> with STNPLT.  A station table looks like the
>> tables found in $GEMTBL/stns. See the GEMPAK help
>> information for more information on STNPLT.
>> 
>> Since you created a surface format gempak file,
>> you can plot the observation locations using the MARK variable
>> in the SFPARM parameter using SFMAP like:
>> 
>> SFPARM = mark:23:3:6
>> 
>> where the : characters separate the marker number, size and width.
>> 
>> Use the COLOR parameter to control the marker color.
>> 
>> Steve Chiswell
>> 
>