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.
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/ >< >**************************************************************************** >< >