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

20010201: 20010201: GEMPAK's gdplot2

Stephen,

I think I see what you are seeing.

The cmcreg data I have is 6 hourly with accumulated precipitations at
each forecast time, eg P06M^f06, P12M^f12, P18M^f18, P24M^f24....etc to P48M^f48

When I use just the 12hr_mm restore file (eg p12m, I get displays at
F012, F024 and f048 with blank screens at the other times.

You are correct that the P--M parameter allows the software to accumulate over 
time.
This is done through both PRXX grids (precipitation rate grids) or major
and subinterval P--M grids, eg: DG_PRFT(). 
However, the code only has the subinterval accumulations at :
     cdt  / '03', '06', '09', '12', '24', '48' /

That is why I only get 12 hour precip grids at 12, 24 and 48. There is no
18 hour or 30 hour accumulation interval etc.

What the code is trying to do is create a p12m grid at each forecast time
step from the accumulating value at each forecast time, eg p12M, p24m, p48m.

The eta for instance gives us P06m at F006, F018, F030, F042, F054 and
P12M at F012, F024, F036, F048, F060. The gempak code actually subtracts
the P06M from the P12M to get the P06M values at all the forecast times. 

So, I think the problem lies in whether you want a constant period of 
accumulation
(eg 6 hours) or just the total accumulation (growing with time). To accomplish
the first case, you could add these times into the source code (both major and
subinterval times).

To see the the accumulation values, its probably easiest just to rename the
grids so that the code doesn't interpret the data as [PCS]--M grids
(gempak will actually try to use P = S + C algebra to calculate total
precip from largescale and convective grids too).

For instance:
#!/bin/csh -f

@ COUNT = 6
while ($COUNT <= 48)
   if($COUNT < 10) then
      set HOUR=0${COUNT}
   else
      set HOUR=$COUNT
   endif
   echo $HOUR
   gddiag << EOF
      gdfile = $MODEL/2001020112_cmcreg.gem
      gdoutf = $MODEL/2001020112_cmcreg.gem
      glevel = 0
      gvcord = none
      GFUNC  = p${HOUR}m
      GDATTIM = f${HOUR}
      GRDNAM = ptot
      gpack  =
      r

      e
EOF

   @ COUNT = $COUNT + 6

end 



In the restore file, I just plotted GFUNC = PTOT. This gives me the growing
blobs for the total accumulation as CMC is creating.

So I guess the question is, what do you expect to see: totals, or the 
sub(PyyM, PxxM) type of plots?

Steve Chiswell
Unidata User Support




     



>From: Stephen Sinnis <address@hidden>
>Organization: UCAR/Unidata
>Keywords: 200102012237.f11MbOX28437

>Steve,
>
>I think I might of missed something in one of the configuration files which
>identifies the GEM to the NAWIPS.  Everything is ok if I use the other
>models like eta, mrf, ruc....
>
>I can use gdplot2 and see any on the GEM's precip grids if I specify
>PXXM^fXX for gdpfunc, but the same assignment does not work with NMAP2.  I
>have included the contents of the NTS file I use ok with gdplot2.  What's
>also neat is that if I specify just p03m, I get the same as if i specify
>sub(pYYm,pXXm), where YY is XX+3hr. for each time step (from gdplot2).
>
>!
>! Restore File : reg_precip_3_48
>!
>! Log:
>! S.W.Sinnis/Pelmorex    12/00  Initial creation
>!
>GDATTIM f003-f012-003
>GLEVEL  0
>
>GVCORD  none
>
>PANEL   0
>
>SKIP    0
>
>SCALE   0
>
>GDPFUN  p03m^f3!p06m^f6!p09m^f9!p12m^12
>TYPE    f
>CONTUR  1
>
>CINT    
>LINE    
>FINT    .01;.1;.25;.5;.75;1;1.25;1.5;1.75;2;2.5;3;4;5;6;7;8;9
>FLINE   0;21-30;14-20;5
>HILO    31;0/x#2/0.1-99//100;0/y
>HLSYM   1.5
>CLRBAR  1
>
>WIND
>
>REFVEC
>
>TITLE   5/-1/~ Total Cummulative Precipitation (mm) 
>TEXT    1/21//hw
>
>GAREA   grid
>CLEAR   yes
>
>STNPLT
>
>SATFIL
>
>RADFIL
>
>STREAM
>
>POSN    4
>
>COLORS  2
>
>MARKER  2
>
>GRDLBL  5
>
>LUTFIL  none
>FILTER  yes
>-----Original Message-----
>From: Unidata Support [mailto:address@hidden]
>Sent: Thursday, February 01, 2001 5:12 PM
>To: Stephen Sinnis
>Cc: address@hidden
>Subject: 20010201: GEMPAK's gdplot2 
>
>
>
>Stephen,
>
>I was able to plot the 12mm precip in NMAP2 by adding cmc entry to the
>existing
>PRECIP_12hr_MM entry in mod_res.tbl. You mentioned creating your own
>entry in $NMAP_RESTORE/prcp/- so not using the existing product?
>
>When you load, you will not have a 12 hour precip frame for F000, but
>will have data plotted for frames 3, 5 and 9 (12, 24 and 48 hours as you
>mention below with gdplot). You can select the time widget to start with
>F012
>instead of F000 (that over rides the GDATTIM in the restore file). 
>
>Since you don't mention seeing any data for any frames, I guess I should ask
>you
>to either provide me with the entry in the prcp directory you created, 
>or ask you to try using the restore file that already exists.
>
>
>Steve Chiswell
>Unidata User SUpport
>
>
>
>>From: Stephen Sinnis <address@hidden>
>>Organization: UCAR/Unidata
>>Keywords: 200102011635.f11GZwX13624
>
>>Good day.
>>
>>I am having a problem displaying the precipitation of a GRID using NMAP2,
>>(GEMPAK6a HP-UX 10.20).
>>
>>I created entries in $GEMTBL/config/datatype.tbl, $GEMTBL/mod_res.tbl and
>>$NMAP_RESTORE/prcp/ for NMAP2, and can request the loading of the precip,
>>but nothing appears on the screen.
>>
>>If I use GDPLOT2 and restore the same precip NTS file from
>>$NMAP_RESTORE/prcp/precip_12hr_mm, I get images for 12, 24 and 48.
>>
>>Any suggestions would be appreciated.
>>
>>Stephen Sinnis - The Weather Network
>>
>>DETAILS:
>>
>>I get a Canadian Regional (short-term F000-F048-003) GEM grib file, (I
>>included the gdinfo at the end of this e-mail).  I use dcgrib2, with the
>>wmogrib2 file that has all 255 components active.  The parameter 61 (Total
>>Precipitation) is available for each of the 3hr grid sets i.e. P00M, P03M,
>>P06M,...,P48M, and is defined by CMS to contain the total accumulated
>>precipitation at each successive grid.  That is the value for P48M will
>>always be greater than or equal to any other PXXM in the file.
>>
>>(GDINFO FALL ALL)
>>          
>>
>> GRID NAVIGATION: 
>>     PROJECTION:          CED                 
>>     GRID SIZE:          410 104
>>     LL CORNER:              40.00    220.02
>>     UR CORNER:              62.66    310.00
>>
>> GRID ANALYSIS BLOCK: 
>>      UNKNOWN ANALYSIS TYPE 
>>
>> Number of grids in file:   264
>>
>> Maximum number of grids in file:   5000
>>
>>  NUM       TIME1              TIME2           LEVL1 LEVL2  VCORD PARM
>>    1     010201/0000F000                      10000         ELVL TMPK
>>
>>    2     010201/0000F000                      10000         ELVL DWPK
>>
>>    3     010201/0000F000                      10000         ELVL UREL
>>
>>    4     010201/0000F000                      10000         ELVL VREL
>>
>>    5     010201/0000F000                      10000         ELVL WXTR
>>
>>    6     010201/0000F000                        850         PRES UREL
>>
>>    7     010201/0000F000                        850         PRES VREL
>>
>>    8     010201/0000F000                        700         PRES UREL
>>
>>    9     010201/0000F000                        700         PRES VREL
>>
>>   10     010201/0000F000                        500         PRES UREL
>>
>>   11     010201/0000F000                        500         PRES VREL
>>
>>   12     010201/0000F000                          0         NONE PMSL
>>
>>   13     010201/0000F000                          0         NONE PXXM
>>
>>   14     010201/0000F000                          0         NONE C00M
>>
>>   15     010201/0000F000                          0         NONE CLD
>>
>>   16     010201/0000F000                          0         NONE LFT4
>>
>>   17     010201/0000F003                      10000         ELVL TMPK
>>
>>   18     010201/0000F003                      10000         ELVL DWPK
>>
>>   19     010201/0000F003                      10000         ELVL UREL
>>
>>   20     010201/0000F003                      10000         ELVL VREL
>>
>>   21     010201/0000F003                      10000         ELVL WXTR
>>
>>   22     010201/0000F003                        850         PRES UREL
>>
>>   23     010201/0000F003                        850         PRES VREL
>>
>>   24     010201/0000F003                        700         PRES UREL
>>
>>   25     010201/0000F003                        700         PRES VREL
>>
>>   26     010201/0000F003                        500         PRES UREL
>>
>>   27     010201/0000F003                        500         PRES VREL
>>
>>   28     010201/0000F003                          0         NONE PMSL
>>
>>   29     010201/0000F003                          0         NONE P03M
>>
>>   30     010201/0000F003                          0         NONE C03M
>>
>>   31     010201/0000F003                          0         NONE CLD
>>
>>   32     010201/0000F006                      10000         ELVL TMPK
>>
>>   33     010201/0000F006                      10000         ELVL DWPK
>>
>>   34     010201/0000F006                      10000         ELVL UREL
>>
>>   35     010201/0000F006                      10000         ELVL VREL
>>
>>...
>> 234     010201/0000F045                      10000         ELVL TMPK
>>
>>  235     010201/0000F045                      10000         ELVL DWPK
>>
>  236     010201/0000F045                      10000         ELVL UREL
>>
>>  237     010201/0000F045                      10000         ELVL VREL
>>
>>  238     010201/0000F045                      10000         ELVL WXTR
>>
>>  239     010201/0000F045                        850         PRES UREL
>>
>>  240     010201/0000F045                        850         PRES VREL
>>
>>  241     010201/0000F045                        700         PRES UREL
>>
>>  242     010201/0000F045                        700         PRES VREL
>>
>>  243     010201/0000F045                        500         PRES UREL
>>
>>  244     010201/0000F045                        500         PRES VREL
>>
>>  245     010201/0000F045                          0         NONE PMSL
>>
>>  246     010201/0000F045                          0         NONE P45M
>>
>>  247     010201/0000F045                          0         NONE C45M
>>
>>  248     010201/0000F045                          0         NONE CLD
>>
>>  249     010201/0000F048                      10000         ELVL TMPK
>>
>>  250     010201/0000F048                      10000         ELVL DWPK
>>
>>  251     010201/0000F048                      10000         ELVL UREL
>>
>>  252     010201/0000F048                      10000         ELVL VREL
>>
>>  253     010201/0000F048                      10000         ELVL WXTR
>>
>>  254     010201/0000F048                        850         PRES UREL
>>
>>  255     010201/0000F048                        850         PRES VREL
>>
>>  256     010201/0000F048                        700         PRES UREL
>>
>>  257     010201/0000F048                        700         PRES VREL
>>
>>  258     010201/0000F048                        500         PRES UREL
>>
>>  259     010201/0000F048                        500         PRES VREL
>>
>>  260     010201/0000F048                          0         NONE PMSL
>>
>>  261     010201/0000F048                          0         NONE P48M
>>
>>  262     010201/0000F048                          0         NONE C48M
>>
>>  263     010201/0000F048                          0         NONE CLD
>>
>>  264     010201/0000F048                          0         NONE LFT4
>>
>
>****************************************************************************
><
>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/
><
>****************************************************************************
><
>