NSF Unidata

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

[GEMPAK #CYW-235722]: gdcross question

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.



The file to view is pqact.conf

But what these FILE entries do is simply same the grib2 messages to the 
directory /usr/local/ldm/data/gempak/grib/nam/.  How are these being decoded 
into GEMPAK grid files?  You may need to talk to the person who set up the LDM 
over there.

I was hoping you could find the entry which has a PIPE command, such as:

HDS|CONDUIT ((/mNAM|/mNMM).*#212|prod/nam.*awip3d)
   PIPE dcgrib2 -d data/gempak/logs/dcgrib2_NAM212.log

You may need to search for NAM capitalized.

-mj


> Alright, so here is what I got the search-
> 
> [dclabo@funnel:/ldm/etc]% grep nam212 *
> pqact.conf:#    FILE    
> /usr/local/ldm/data/gempak/grib/nam/nam212_\1-\2-\3_\4.grib
> pqact.conf:     FILE    
> /usr/local/ldm/data/gempak/grib/nam/nam212_\1-\2-\3_\4_F0\5.grib2
> pqact.conf:     FILE    
> /usr/local/ldm/data/gempak/grib/nam/nam212_\1-\2-\3_\4_F1\5.grib2
> pqact.conf:     FILE    
> /usr/local/ldm/data/gempak/grib/nam/nam212_\1-\2-\3_\4_F2\5.grib2
> pqact.conf:     FILE    
> /usr/local/ldm/data/gempak/grib/nam/nam212_\1-\2-\3_\4_F3\5.grib2
> pqact.conf:     FILE    
> /usr/local/ldm/data/gempak/grib/nam/nam212_\1-\2-\3_\4_F4\5.grib2
> pqact.conf~:#   FILE    
> /usr/local/ldm/data/gempak/grib/nam/nam212_\1-\2-\3_\4.grib
> pqact.conf~:    FILE    
> /usr/local/ldm/data/gempak/grib/nam/nam212_\1-\2-\3_\4_F0\5.grib2
> pqact.conf~:    FILE    
> /usr/local/ldm/data/gempak/grib/nam/nam212_\1-\2-\3_\4_F1\5.grib2
> pqact.conf~:    FILE    
> /usr/local/ldm/data/gempak/grib/nam/nam212_\1-\2-\3_\4_F2\5.grib2
> pqact.conf~:    FILE    
> /usr/local/ldm/data/gempak/grib/nam/nam212_\1-\2-\3_\4_F3\5.grib2
> pqact.conf~:    FILE    
> /usr/local/ldm/data/gempak/grib/nam/nam212_\1-\2-\3_\4_F4\5.grib2
> 
> Do I view one of these files?
> 
> --
> Darren R. Clabo
> South Dakota State Fire Meteorologist
> O: (605) 394-1996  M: (605) 381-9253
> address@hidden
> http://www.ias.sdsmt.edu/clabo/index.htm
> 
> 
> -----Original Message-----
> From: Unidata GEMPAK Support [mailto:address@hidden]
> Sent: Tuesday, November 20, 2012 2:07 PM
> To: Clabo, Darren R.
> Cc: address@hidden
> Subject: [GEMPAK #CYW-235722]: gdcross question
> 
> Whoops, the grep command to run in ~ldm/etc/ is "grep NAM212 *"
> 
> > Sure, no problem.
> >
> > The contents of any grid file are determined by a line in the LDM config 
> > file "pqact.conf" (or similar), located in the LDM user's ~/etc/ directory.
> >
> > Each pattern action matches a regular expression in the product header, or 
> > a file name, and send the message (in this case a grib2 file) to the 
> > "dcgrib2" decoder, which decodes from grib2 format to GEMPAK grid format.  
> > The result are the grid files used in GDFILE, such as 
> > "2012112012_nam212.gem"
> >
> > So, the content of the nam212 GEMPAK grid file is determined by what 
> > product feed is being requested (such as CONDUIT, NGRID or HDS).
> >
> > I'm guessing that the pattern action for the nam212 on your end is 
> > different than what we have here, based on the number of grids in the 
> > nam212 file being significantly less than what I see on my end.
> >
> > In ~ldm/etc/, run the command "grep NAM212" to list files which contain the 
> > text "NAM212" - then view that file with "vi" or "less" and search for the 
> > NAM212 entry, copy and paste it back to me so I can see what is different.
> >
> > Michael
> >
> >
> > > You've lost me on this one, Michael. I am not familiar with the pattern 
> > > action (and frankly not familiar with our LDM feed in general). Can you 
> > > walk me through this one?
> > >
> > > --
> > > Darren R. Clabo
> > > South Dakota State Fire Meteorologist
> > > O: (605) 394-1996  M: (605) 381-9253 address@hidden
> > > http://www.ias.sdsmt.edu/clabo/index.htm
> > >
> > >
> > > -----Original Message-----
> > > From: Unidata GEMPAK Support
> > > [mailto:address@hidden]
> > > Sent: Tuesday, November 20, 2012 12:34 PM
> > > To: Clabo, Darren R.
> > > Cc: address@hidden
> > > Subject: [GEMPAK #CYW-235722]: gdcross question
> > >
> > > Well, I'm not sure why the elevation isn't showing up for F00 since the 
> > > grid exists, but perhaps you can plot elevation as point values in 
> > > GDPLOT2 to check that coverage is complete for the US.
> > >
> > > That said, the number of grids in your nam212 file is significantly lower 
> > > than in the ones we decode to at Unidata, about 20% of what I see in our 
> > > grid files.
> > >
> > > What does your pattern action look like for the nam212?  Our is pulling 
> > > both the HDS and CONDUIT feeds:
> > >
> > > # NAM #212 40km grid CONUS
> > > HDS|CONDUIT     ((/mNAM|/mNMM).*#212|prod/nam.*awip3d)
> > > PIPE    dcgrib2 -d data/gempak/logs/dcgrib2_NAM212.log
> > > -e GEMTBL=/home/gempak/NAWIPS/gempak/tables
> > >
> > >
> > > The GDINFO grid count for the 212:
> > >
> > > GRID FILE: /data/ldm/gempak/model/nam/2012112012_nam212.gem
> > >
> > > GRID NAVIGATION:
> > > PROJECTION:          LCC
> > > ANGLES:                25.0   -95.0    25.0
> > > GRID SIZE:          185 129
> > > LL CORNER:              12.19   -133.46
> > > UR CORNER:              57.29    -49.39
> > >
> > > Number of grids in file: 18975
> > >
> > >
> > >
> > > -Michael
> > >
> > >
> > > > Here is my output for gdinfo. It appears that only the f00 forecast 
> > > > hour has a hght parameter at the surface. None of the rest of the 
> > > > forecast hours had a glevel=0 at gvcord=none. Interestingly enough, I 
> > > > was looking at the f00 run for the cross section images I just sent to 
> > > > you.
> > > >
> > > > GEMPAK-GDINFO>l
> > > > GDFILE   = /squall1/gempak-data/model/nam/2012111212_nam212.gem
> > > > LSTALL   = YES
> > > > OUTPUT   =
> > > > GDATTIM  = all
> > > > GLEVEL   = 0
> > > > GVCORD   = none
> > > > GFUNC    = hght
> > > > GEMPAK-GDINFO>r
> > > >
> > > > GRID FILE: /squall1/gempak-data/model/nam/2012111212_nam212.gem
> > > >
> > > > GRID NAVIGATION:
> > > > PROJECTION:          LCC
> > > > ANGLES:                25.0   -95.0    25.0
> > > > GRID SIZE:          185 129
> > > > LL CORNER:              12.19   -133.46
> > > > UR CORNER:              57.29    -49.39
> > > >
> > > > GRID ANALYSIS BLOCK:
> > > > UNKNOWN ANALYSIS TYPE
> > > >
> > > > Number of grids in file:  3517
> > > >
> > > > Maximum number of grids in file:  24000
> > > >
> > > > NUM       TIME1              TIME2           LEVL1 LEVL2  VCORD PARM
> > > > 170     121112/1200F000                          0         NONE HGHT
> > > >
> > > > --
> > > > Darren R. Clabo
> > > > South Dakota State Fire Meteorologist
> > > > O: (605) 394-1996  M: (605) 381-9253 address@hidden
> > > > http://www.ias.sdsmt.edu/clabo/index.htm
> > > >
> > > >
> > > > -----Original Message-----
> > > > From: Unidata GEMPAK Support
> > > > [mailto:address@hidden]
> > > > Sent: Tuesday, November 20, 2012 11:19 AM
> > > > To: Clabo, Darren R.
> > > > Cc: address@hidden
> > > > Subject: [GEMPAK #CYW-235722]: gdcross question
> > > >
> > > > Ok, I believe I know what the issue is.
> > > >
> > > > GDCROSS uses the surface elevation grid contained in GDFILE.  For the 
> > > > nam212, using GDINFO, I see the following:
> > > >
> > > >
> > > > GEMPAK-GDINFO>l
> > > > GDFILE   = $MODEL/nam/2012112012_nam212.gem
> > > > LSTALL   = YES
> > > > OUTPUT   = T
> > > > GDATTIM  = all
> > > > GLEVEL   = 0
> > > > GVCORD   = none
> > > > GFUNC    = hght
> > > > GEMPAK-GDINFO>r
> > > >
> > > > NUM       TIME1              TIME2           LEVL1 LEVL2  VCORD PARM
> > > > 565     121120/1200F000                          0         NONE HGHT
> > > > 1278     121120/1200F003                          0         NONE HGHT
> > > > 1933     121120/1200F006                          0         NONE HGHT
> > > > 2567     121120/1200F009                          0         NONE HGHT
> > > > 3243     121120/1200F012                          0         NONE HGHT
> > > > 3862     121120/1200F015                          0         NONE HGHT
> > > > 4490     121120/1200F018                          0         NONE HGHT
> > > > 5136     121120/1200F021                          0         NONE HGHT
> > > > 5821     121120/1200F024                          0         NONE HGHT
> > > > 6498     121120/1200F027                          0         NONE HGHT
> > > > 7153     121120/1200F030                          0         NONE HGHT
> > > > 7791     121120/1200F033                          0         NONE HGHT
> > > > 8411     121120/1200F036                          0         NONE HGHT
> > > > 9084     121120/1200F039                          0         NONE HGHT
> > > > 9717     121120/1200F042                          0         NONE HGHT
> > > > 10425     121120/1200F045                          0         NONE HGHT
> > > > 11048     121120/1200F048                          0         NONE HGHT
> > > > 11682     121120/1200F051                          0         NONE HGHT
> > > > 12366     121120/1200F054                          0         NONE HGHT
> > > > 13008     121120/1200F057                          0         NONE HGHT
> > > > 13680     121120/1200F060                          0         NONE HGHT
> > > > 14337     121120/1200F063                          0         NONE HGHT
> > > > 14984     121120/1200F066                          0         NONE HGHT
> > > > 15647     121120/1200F069                          0         NONE HGHT
> > > > 16318     121120/1200F072                          0         NONE HGHT
> > > > 16934     121120/1200F075                          0         NONE HGHT
> > > > 17589     121120/1200F078                          0         NONE HGHT
> > > > 18226     121120/1200F081                          0         NONE HGHT
> > > > 18889     121120/1200F084                          0         NONE HGHT
> > > >
> > > >
> > > > Can you check that these grids exist in the file you're using for 
> > > > GDFILE?  My guess is that they are not there.  If they are, something 
> > > > else is amiss.
> > > >
> > > > Michael
> > > >
> > > >
> > > > > I just ran that same script for LAX to DEN. Still no surface.
> > > > > See attached ago_test.gif-
> > > > >
> > > > > Additionally, when I run gdprof on the same nam212 file, it is 
> > > > > plotting data down to 1000mb regardless of station surface elevation. 
> > > > > See attached DNR soundings.
> > > > >
> > > > > --
> > > > > Darren R. Clabo
> > > > > South Dakota State Fire Meteorologist
> > > > > O: (605) 394-1996  M: (605) 381-9253 address@hidden
> > > > > http://www.ias.sdsmt.edu/clabo/index.htm
> > > > >
> > > > >
> > > > > -----Original Message-----
> > > > > From: Unidata GEMPAK Support
> > > > > [mailto:address@hidden]
> > > > > Sent: Tuesday, November 20, 2012 10:41 AM
> > > > > To: Clabo, Darren R.
> > > > > Cc: address@hidden
> > > > > Subject: [GEMPAK #CYW-235722]: gdcross question
> > > > >
> > > > > Ah, I see, thanks.
> > > > >
> > > > > Using your GDCROSS input I am seeing elevation between CWPL (391 m) 
> > > > > and CYMO (9 m) - see attached.
> > > > >
> > > > > This being a cross-section across Ontario, there's not much elevation 
> > > > > to show.  The filled contours are shown below the surface because of 
> > > > > the cross-section interpolation.  If I expand CXSTNS out to LAX>BWI, 
> > > > > there are grid points below the surface level which are not 
> > > > > interpolated to, though there is some overlap (second attachment).
> > > > >
> > > > > Are you using a different CXSTNS definition which does not show the 
> > > > > surface?
> > > > >
> > > > > -Michael
> > > > >
> > > > > > Michael-
> > > > > > I am simply using a nam212 gempak file within gdcross to plot
> > > > > > potential vorticity and theta. Parameter list is (note gvcord)
> > > > > >
> > > > > > CXSTNS   = CWPL>CYMO
> > > > > > GDATTIM  = f00
> > > > > > GVCORD   = pres
> > > > > > GFUNC    = pvor(thta,wnd)
> > > > > > GVECT    =
> > > > > > GDFILE   = /squall1/gempak-data/model/nam/2012111212_nam212.gem
> > > > > > WIND     =
> > > > > > REFVEC   = 10
> > > > > > PTYPE    = log
> > > > > > YAXIS    =
> > > > > > IJSKIP   = 0
> > > > > > CINT     =
> > > > > > SCALE    = 6
> > > > > > LINE     = 32/1/3
> > > > > > BORDER   = 1
> > > > > > TITLE    = 1/1
> > > > > > CLEAR    = yes
> > > > > > DEVICE   = gif|ago.gif|1200;900
> > > > > > TEXT     = 1
> > > > > > PANEL    = 0
> > > > > > CONTUR   = 3
> > > > > > FINT     = 1
> > > > > > FLINE    = 13-30
> > > > > > CTYPE    = f
> > > > > >
> > > > > > Any other thoughts?
> > > > > >
> > > > > >
> > > > > > --
> > > > > > Darren R. Clabo
> > > > > > South Dakota State Fire Meteorologist
> > > > > > O: (605) 394-1996  M: (605) 381-9253 address@hidden
> > > > > > http://www.ias.sdsmt.edu/clabo/index.htm
> > > > > >
> > > > > > -----Original Message-----
> > > > > > From: Unidata GEMPAK Support
> > > > > > [mailto:address@hidden]
> > > > > > Sent: Tuesday, November 20, 2012 8:37 AM
> > > > > > To: Clabo, Darren R.
> > > > > > Cc: address@hidden
> > > > > > Subject: [GEMPAK #CYW-235722]: gdcross question
> > > > > >
> > > > > > I can confirm that the GDCROSS code is written to only draw surface 
> > > > > > elevation when GVCORD is PRES.  I'm trying to make it work for HGHT 
> > > > > > and THTA but so far it's not an easy task.  Not the news you wanted 
> > > > > > to hear, sorry.
> > > > > >
> > > > > > Michael
> > > > > >
> > > > > > > Hi Darren,
> > > > > > >
> > > > > > > Without seeing the GDCROSS input I'm making a guess about how you 
> > > > > > > went about this.  Did you run GDVINT to interpolate to theta 
> > > > > > > levels and then run this new grid file through GDCROSS?
> > > > > > >
> > > > > > > I see what you describe when I plot a linear cross-section of 
> > > > > > > temperatures by the HGHT coordinate, though there are only 6 HGHT 
> > > > > > > levels for TMPK in the GFS, so the scope is limited, but it does 
> > > > > > > not show a surface elevation cross-section at the bottom, like 
> > > > > > > you say.
> > > > > > >
> > > > > > > I'm guessing that surface elevation is only designed to show for 
> > > > > > > GVCORD=PRES.  I'm looking at the source code and will see if I 
> > > > > > > can get it to work with height.
> > > > > > >
> > > > > > > -Michael
> > > > > > >
> > > > > > >
> > > > > > >
> > > > > > > > Good day-
> > > > > > > > I am working on putting together an isentropic analysis in 
> > > > > > > > GDCROSS. However, when I am plotting a cross section from 
> > > > > > > > across the CONUS, all of the stations are showing data down to 
> > > > > > > > 1000 mb. Considering the elevation of the Intermountain West, I 
> > > > > > > > was hoping I could get a plot that resembles what is found on 
> > > > > > > > the GEMPAK tutorial page 
> > > > > > > > (http://www.unidata.ucar.edu/software/gempak/tutorial/vertical_plotting.html#12.2)
> > > > > > > >  where the surface elevation is easily seen.
> > > > > > > >
> > > > > > > > Any thoughts on what I can do?
> > > > > > > >
> > > > > > > > --
> > > > > > > > Darren R. Clabo
> > > > > > > > South Dakota State Fire Meteorologist
> > > > > > > > O: (605) 394-1996  M: (605) 381-9253
> > > > > > > > address@hidden
> > > > > > > > http://www.ias.sdsmt.edu/clabo/index.htm
> > > > > > > >
> > > > > > > >
> > > > > > > >
> > > > > > >
> > > > > >
> > > > > >
> > > > > > Ticket Details
> > > > > > ===================
> > > > > > Ticket ID: CYW-235722
> > > > > > Department: Support GEMPAK
> > > > > > Priority: Normal
> > > > > > Status: Open
> > > > > >
> > > > > >
> > > > >
> > > > >
> > > > > Ticket Details
> > > > > ===================
> > > > > Ticket ID: CYW-235722
> > > > > Department: Support GEMPAK
> > > > > Priority: Normal
> > > > > Status: Open
> > > > >
> > > >
> > > > Ticket Details
> > > > ===================
> > > > Ticket ID: CYW-235722
> > > > Department: Support GEMPAK
> > > > Priority: Normal
> > > > Status: Open
> > > >
> > > >
> > >
> > >
> > > Ticket Details
> > > ===================
> > > Ticket ID: CYW-235722
> > > Department: Support GEMPAK
> > > Priority: Normal
> > > Status: Open
> > >
> > >
> 
> 
> Ticket Details
> ===================
> Ticket ID: CYW-235722
> Department: Support GEMPAK
> Priority: Normal
> Status: Open
> 
> 

Ticket Details
===================
Ticket ID: CYW-235722
Department: Support GEMPAK
Priority: Normal
Status: Open