[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
20011031: 20011030: using stnplt to add a TC marker
- 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
>>
>