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

[GEMPAK #CYW-235722]: gdcross question



In ~ldm/etc/ldmd.conf there will be one or more lines such as

EXEC    "pqact -f ANY-CRAFT-NNEXRAD-CONDUIT /local/ldm/etc/pqact.conf"

Each "EXEC" line will have a file like "pqact.conf" - some may have names like 
"pqact.gempak_decoders"

Please send me each file that is contained in an un-commented EXEC lind in 
ldmd.conf and I'll take a look.

Michael


> Within the pqact file itself, the only nam212 (nothing with NAM212) related 
> stuff I see is this:
> 
> # The '-strip' option to FILE, STDIOFILE, & PIPE removes control characters
> # (other than the '\n') from the data before output. This is useful for
> # cleaning up WMO format text bulletins.
> #
> ...skipping...
> #       FILE    /usr/local/ldm/data/gempak/grib/nam/nam212_\1-\2-\3_\4.grib
> #
> #HDS    grib/ncep/NAM_84/#215/(....)(..)(..)(..)(..)
> #       FILE    /usr/local/ldm/data/gempak/grib/nam/nam215_\1-\2-\3_\4.grib
> #
> #HDS    grib/ncep/GFS/#211/(....)(..)(..)(..)(..)
> #       FILE    /usr/local/ldm/data/gempak/grib/gfs/gfs211_\1-\2-\3_\4.grib
> #
> #
> #
> #
> CONDUIT ^data/nccf/com/nam/prod/nam.(....)(..)(..)/nam.t(..)z.awip3d0(.)
> FILE    /usr/local/ldm/data/gempak/grib/nam/nam212_\1-\2-\3_\4_F0\5.grib2
> CONDUIT ^data/nccf/com/nam/prod/nam.(....)(..)(..)/nam.t(..)z.awip3d1(.)
> FILE    /usr/local/ldm/data/gempak/grib/nam/nam212_\1-\2-\3_\4_F1\5.grib2
> CONDUIT ^data/nccf/com/nam/prod/nam.(....)(..)(..)/nam.t(..)z.awip3d2(.)
> FILE    /usr/local/ldm/data/gempak/grib/nam/nam212_\1-\2-\3_\4_F2\5.grib2
> CONDUIT ^data/nccf/com/nam/prod/nam.(....)(..)(..)/nam.t(..)z.awip3d3(.)
> FILE    /usr/local/ldm/data/gempak/grib/nam/nam212_\1-\2-\3_\4_F3\5.grib2
> CONDUIT ^data/nccf/com/nam/prod/nam.(....)(..)(..)/nam.t(..)z.awip3d4(.)
> FILE    /usr/local/ldm/data/gempak/grib/nam/nam212_\1-\2-\3_\4_F4\5.grib2
> 
> 
> Our ldm/unix guy (Bill Capehart) is out for the holiday.
> 
> 
> --
> 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:25 PM
> To: Clabo, Darren R.
> Cc: address@hidden
> Subject: [GEMPAK #CYW-235722]: gdcross question
> 
> 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
> 
> 

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