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

  • Subject: [GEMPAK #QHN-534327]: GEMAPK Question - vertical maximum
  • Date: Thu, 17 Apr 2014 13:30:18 -0600
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