NSF Unidata

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

20000221: GDVINT

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.



Clint,

My interpretation of VCOORD is that it will copy over surface values
on the specified coordinates, so if you gave gvcord=pres/thta
and vcoord=none;hght, then aside from the theta coordinate grids,
you will copy over those grids on level=0 of the specified
vcoord. This won't copy over every pressure grid from the original
file to my reading. Instead, you would use gdmod as you mentioned.

Steve Chiswell
Unidata User Support

>From: Clint Rowe <address@hidden>
>Organization: .
>Keywords: 200002212216.PAA04348

>Chiz,
>
>I thought specifying VCOORD=whatever in GDVINT was supposed to copy grids
>from the specified vertical coordinates into the output file along with
>the newly interpolated grids.  This doesn't seem to be happening.  I followed
>the examples in the manual and tutorial; here's the results of my attempt:
>
>levanto% gdvint
> GDFILE    Grid file                         ../thtagrids/testeta211.gem
> GDOUTF    Output grid file                  ../thtagrids/testeta211.gem
> GDATTIM   Grid date/time                    first
> GVCORD    Grid vertical coordinate          pres/thta
> GLEVEL    Grid level                        280-350-5
> MAXGRD    Maximum number of grids           3000
> GAREA     Graphics area                     @1;1;93;65
> VCOORD    Vertical coordinate type          pres;none;sgma;pagl;hagl;eatm
> Parameters requested: GDFILE,GDOUTF,GDATTIM,GVCORD,GLEVEL,MAXGRD,GAREA,
> VCOORD.
> GEMPAK-GDVINT>glev=300
> GEMPAK-GDVINT>gdfile=/data/gempak/hds/2000022112_eta211.gem
> GEMPAK-GDVINT>l
> GDFILE   = /data/gempak/hds/2000022112_eta211.gem
> GDOUTF   = ../thtagrids/testeta211.gem
> GDATTIM  = first
> GVCORD   = pres/thta
> GLEVEL   = 300
> MAXGRD   = 3000
> GAREA    = @1;1;93;65
> VCOORD   = pres;none;sgma;pagl;hagl;eatm
> GEMPAK-GDVINT>r
> PRES         found. 
> TMPK         found. 
> HGHT         found. 
> OMEG         found. 
> UREL         found. 
> VREL         found. 
> RELH         found. 
> AVOR         found. 
> look at ip   2
> look at ip   3
> look outlev     300.000
> look at ip   4
> look outlev     300.000
> look at ip   5
> look outlev     300.000
> look at ip   6
> look outlev     300.000
> look at ip   7
> look outlev     300.000
> look at ip   8
> look outlev     300.000
>  Interpolation is done.  No errors. 
> Parameters requested: GDFILE,GDOUTF,GDATTIM,GVCORD,GLEVEL,MAXGRD,GAREA,
> VCOORD.
> GEMPAK-GDVINT>exit
> levanto% gdinfo
> GDFILE    Grid file                         
>/data/gempak/hds/2000022112_eta211.gem
> LSTALL    Full list flag                    YES
> OUTPUT    Output device/filename            T
> Parameters requested: GDFILE,LSTALL,OUTPUT.
> GEMPAK-GDINFO>gdfile=../thtagrids/testeta211.gem
> GEMPAK-GDINFO>r
>
> GRID FILE: ../thtagrids/testeta211.gem                                       
>       
>
> GRID NAVIGATION: 
>     PROJECTION:          LCC                 
>     ANGLES:                25.0   -95.0    25.0
>     GRID SIZE:           93  65
>     LL CORNER:              12.19   -133.46
>     UR CORNER:              57.29    -49.38
>
> GRID ANALYSIS BLOCK: 
>      UNKNOWN ANALYSIS TYPE 
>
> Number of grids in file:    16
>
> Maximum number of grids in file:   3000
>
>  NUM       TIME1              TIME2           LEVL1 LEVL2  VCORD PARM
>    1     000221/1200F000                          0         THTA PRES        
>    2     000221/1200F000                          0         THTA TMPK        
>    3     000221/1200F000                          0         THTA HGHT        
>    4     000221/1200F000                          0         THTA OMEG        
>    5     000221/1200F000                          0         THTA UREL        
>    6     000221/1200F000                          0         THTA VREL        
>    7     000221/1200F000                          0         THTA SPFH        
>    8     000221/1200F000                          0         THTA AVOR        
>    9     000221/1200F000                          0         THTA THTA        
>   10     000221/1200F000                        300         THTA PRES        
>   11     000221/1200F000                        300         THTA PSYM        
>   12     000221/1200F000                        300         THTA OMEG        
>   13     000221/1200F000                        300         THTA UREL        
>   14     000221/1200F000                        300         THTA VREL        
>   15     000221/1200F000                        300         THTA SPFH        
>   16     000221/1200F000                        300         THTA AVOR        
> Parameters requested: GDFILE,LSTALL,OUTPUT.
> GEMPAK-GDINFO>exit
> levanto%
> 
>Do I have to go the route of GDCFIL to create a file, GDMOD to copy the grids
>from the original eta file into a new file with enough room for the thta
>surfaces, then run GDVINT?
>
>Clint
>
>====================================================================
>Clinton M. Rowe
>Associate Professor 
>Meteorology/Climatology Program                  phone:(402)472-1946
>Department of Geosciences                          fax:(402)472-4917
>University of Nebraska-Lincoln                        address@hidden
>