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