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?