20050304: Storm Report Problems

Robert,

The NWUS2[02] files have report times in minutes. If you are
plotting a single time from the files created by dcstorm, you would only see the reports at the time you specify, that could be 1, or dozens
of reports depending on the data file. If you specify a time range, each
frame would have the reports for the valid time. If you are overlaying
radar images in NMAP2, then the storm report time closest to the image
would be used when the dominant time is the image time.

Using SFLIST with SFFILE=SELS will show all the available report times,
for instance yesterday:
SFFILE   = sels|050303/0000
AREA     = uslcc
DATTIM   = all
SFPARM   = dset
OUTPUT   = t
IDNTYP   = stid
GEMPAK-SFLIST>r
PARM
EVNT;STAT;DATE;HOUR
   STN    YYMMDD/HHMM      EVNT     STAT     DATE     HOUR
A1624000  050303/1624      3.00     0.00     3.00  1624.00
A1652000  050303/1652      3.00     0.00     3.00  1652.00
A1652001  050303/1652      3.00     0.00     3.00  1652.00
A1652002  050303/1652      3.00     0.00     3.00  1652.00
A1652003  050303/1652      3.00     0.00     3.00  1652.00
A1652004  050303/1652      3.00     0.00     3.00  1652.00
A1652005  050303/1652      3.00     0.00     3.00  1652.00
....
many more reports
....
A1922006  050303/1922      3.00     0.00     3.00  1922.00
A1922007  050303/1922      3.00     0.00     3.00  1922.00
A1922008  050303/1922      3.00     0.00     3.00  1922.00
A1922009  050303/1922      3.00     0.00     3.00  1922.00
A1922010  050303/1922      3.00     0.00     3.00  1922.00
A1922011  050303/1922      3.00     0.00     3.00  1922.00
A2015000  050303/2015      3.00     0.00     3.00  2015.00
A2050000  050303/2050      3.00     0.00     3.00  2050.00
A2155000  050303/2155      3.00     0.00     3.00  2155.00
A2158000  050303/2158      3.00     0.00     3.00  2158.00
A2200000  050303/2200      3.00     0.00     3.00  2200.00
A2320000  050303/2320      3.00     0.00     3.00  2320.00


If you wanted to bin these to a single date/time, you could
use sfcfil to create a new file, and sfmod to re-assign those
times with DATOUT, eg:

% sfcfil
SFOUTF    Output surface file               newsels.gem
SFPRMF    Surface parameter packing file    sels.pack
STNFIL Station information file SHIPFL Ship data file flag no
TIMSTN    Times/additional stations         24/20000
SFFSRC Surface file source Parameters requested: SFOUTF,SFPRMF,STNFIL,SHIPFL,TIMSTN,SFFSRC.
GEMPAK-SFCFIL>r

SFCFIL PARAMETERS: New surface file: newsels.gem Parameter file: sels.pack Station file: Number of stations in STNFIL: 0
Number of additional stations:  20000
Total number of stations:       20000
Total number of times:             24

This file will be a standard surface file.

Enter <cr> to accept parameters or type EXIT:
Parameters requested: SFOUTF,SFPRMF,STNFIL,SHIPFL,TIMSTN,SFFSRC.
GEMPAK-SFCFIL>e

% sfmod
SFFILE    Surface data file                 sels
SFOUTF    Output surface file               newsels.gem
DATTIM    Date/time                         0000-2359
DATOUT    Output date/time                  0000
AREA      Data area                         dset
SFPARM    Surface parameter list            dset
Parameters requested: SFFILE,SFOUTF,DATTIM,DATOUT,AREA,SFPARM.
GEMPAK-SFMOD>r

SFMOD PARAMETERS: Input file: $GEMDATA/storm/sels/20050303_sels.gem Output file: newsels.gem Area: dset Output parameters: dset Input date/time: Output date/time: 050303/1624 050303/0000 050303/1652 050303/0000 050303/1735 050303/0000 050303/1750 050303/0000 050303/1823 050303/0000 050303/1835 050303/0000 050303/1836 050303/0000 050303/1852 050303/0000 050303/1900 050303/0000 050303/1909 050303/0000 050303/1920 050303/0000 050303/1922 050303/0000 050303/2015 050303/0000 050303/2050 050303/0000 050303/2155 050303/0000 050303/2158 050303/0000 050303/2200 050303/0000 050303/2320 050303/0000 Enter <cr> to accept parameters or type EXIT:
Parameters requested: SFFILE,SFOUTF,DATTIM,DATOUT,AREA,SFPARM.
GEMPAK-SFMOD>e

Note, you could bin hourly, daily, 15 minutes, etc.

As I binned yesterday's reports as a single time, I can plot
those obs with in sfmap such as:
GEMPAK-SFMAP>l
AREA     = dset
GAREA    = uslcc
SATFIL = RADFIL = IMCBAR = SFPARM = mark:3
DATTIM   = all
SFFILE   = newsels.gem
COLORS   = 2
MAP      = 31/1/2
MSCALE = LATLON = 0
TITLE    = 5/-1
CLEAR    = y
PANEL    = 0
DEVICE   = xw
PROJ     = def
FILTER   = n
TEXT     = 1/21//hw
LUTFIL   = default
STNPLT = CLRBAR = 1
GEMPAK-SFMAP>r

SFMAP PARAMETERS:

Data area: dset Graphics area name: uslcc Valid parameters: MARK Parameter colors: 2 Time: 050303/0000 File: newsels.gem Map: 31/1/2 Title: 5/-1 Device: xw Projection: def Clear screen: YES Filter: NO Filter factor: 0.00 Panel: 0 Enter <cr> to accept parameters or type EXIT: Parameters requested:
AREA,GAREA,SATFIL,RADFIL,IMCBAR,SFPARM,DATTIM,SFFILE,

COLORS,MAP,MSCALE,LATLON,TITLE,CLEAR,PANEL,DEVICE,PROJ,FILTER,TEXT,LUTFIL,
STNPLT,CLRBAR.
GEMPAK-SFMAP>


The redbook product is that created by SPC. The legend in the product
gives the count of reports, and is reporting zero at this time. The
rendering of this product is just what is in the graphics drawing
commands of the rebook.

Steve Chiswell
Unidata User Support




On Thu, 2005-03-03 at 18:01, Robert Dewey wrote:
I am resending this message, now that storm season is
getting near. I sent it back in August, but never
received a response:

_______________________________

Hello,

I am getting and decoding the NWUS20 and NWUS22 storm
reports from the
LDM (and decoding them), but I never get any "real"
plots.
What I mean by real plots, is that I will get a single
hail report here
and there, but thats it... When I inspect the actual
NWUS20/22 files,
there are many more reports than just one hail report.
Also, is there a
problem with the storm reports in the Redbook
graphics? Your redbook
example gifs never show any plots, and my vgf files do
not show plots as
well...

Thanks
Robert Dewey

_________________________________

Does anyone have an answer?




  • 2005 messages navigation, sorted by:
    1. Thread
    2. Subject
    3. Author
    4. Date
    5. ↑ Table Of Contents
  • Search the gembud archives: