20000221: GDVINT

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
>