NSF Unidata

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

[GEMPAK #QHN-534327]: GEMAPK Question - vertical maximum

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.



Hi Nathan,

I've been trying to figure out different ways to do this.  The first thought I 
can is quick but hack-y and requires some scripting.  You didn't mention if 
this is grid data or SN data so I'm using grid data for my example.

 GEMPAK-GDPROF>l
 GPOINT   = IAD
 GDATTIM  = f012
 GVCORD   = PRES
 GFUNC    = TMPC
 GVECT    = WND
 GDFILE   = nam12
 LINE     = 3
 MARKER   = 0
 BORDER   = 1
 PTYPE    = LOG
 SCALE    = 999
 XAXIS    =  
 YAXIS    =  
 WIND     = BM1
 REFVEC   =  
 WINPOS   = 1
 FILTER   = YES
 TITLE    = 1
 PANEL    = 0
 CLEAR    = YES
 TEXT     = 1
 DEVICE   = XW
 OUTPUT   = f/gdprof.txt
 THTALN   = 0
 THTELN   = 0
 MIXRLN   = 0
 GEMPAK-GDPROF>r

 GDPROF PARAMETERS

 Grid file:                 nam12
 Grid time:                 140417/1200F012
 Grid function:             TMPC
 Grid vector function:      WND
 X,Y grid point:            475.74  206.16
 Lat/lon:                    38.95  -77.45
 Number of levels:           29
 Data range:                -63.79    7.41
 Scaling factor:             10** 0
 Wind display:            BM1          


Two things to notice here: 1) OUTPUT = t/gdprof.txt directs text data to be 
dumped to file.  You can exit GDPROF and cat/grep/sed/whatever the new txt file 
to computer min/max.  2) "Data range:" shows you the range right there!   But 
it's not dumped to the text file.  But if you are scripting, and call GDPROF 
and log the output such as ./callGDProf.csh >& prof.log, then you can grep for 
"Data range" and get your min/max.

Like I said, it's hack-y but it works.  I put an example script here:

http://jsfiddle.net/9PDHP/

I don't have a solution yet which deals with an entire grid.  I was hoping that 
GDDIAG or GDSTAT could help, but I'm not sure how yet.

Michael James
Unidata





> Full Name: Nathan Marsili
> Email Address: address@hidden
> Organization: National Weather Service
> Package Version:
> Operating System:
> Hardware:
> Description of problem: From what I have researched this is probably not 
> possible, but is there any way to produce a maximum value of a parameter in 
> the vertical.  I have a parameter at 10 different levels and would like to 
> display the maximum value of this parameter?
> 
> 

Ticket Details
===================
Ticket ID: QHN-534327
Department: Support GEMPAK
Priority: Normal
Status: Open