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

[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