[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
[GEMPAK #CYW-235722]: gdcross question
- Subject: [GEMPAK #CYW-235722]: gdcross question
- Date: Tue, 20 Nov 2012 12:33:41 -0700
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