gdgaxes

FORTRAN

 
integer gdll, ier
ier = gdgaxes(gdid, ax, ay)
integer gdid
real ax(*), ay(*)

C

 
int c_gdgaxes(int gdid,  float *ax, float *ay)

Routine description

gdgaxes returns the deformation axes (the '^^' and '>>' records in a strandard file) of the grid defined by gdid. Useful mostly for 'Z', 'G' grids. Returns a linearly increasing array from 1 to ni or nj for regular grids. Saves the user of all the validation code needed to read manually these axes from a standard file.

 

gdid

a valid grid identifier (returned by ezqkdef or ezgdef) that will be used as the source (input) grid.

ax, ay

Contents of the deformations axes. The dimensions of ax are assumed to be ni, those of ay are assumed to be nj.

?

Sample FORTRAN code :

The following code returns the deformation axes on a GEM grid. Assume necessary memory space has been allocated.

 
      datev  = -1
      etiket = '        '
      typvar = ' '
      nomvar = 'GZ'
      ip1 = 500
      ip2 = 0
      ip3 = -1
      iun = 13
      
***   read fld info and data
      key = FSTINF(iun, NI, NJ, NK, datev, etiket, ip1, ip2, ip3, typvar, nomvar) 
      ier = FSTPRM(key, DATEO, DEET, NPAS, NI, NJ, NK, NBITS, DATYP, IP1,
     $     IP2, IP3, TYPVAR, NOMVAR, ETIKET, GRTYP, IG1, IG2, IG3,
     $     IG4, SWA, LNG, DLTF, UBC, EX1, EX2, EX3)
***   Define input grid         
      
      gdgem = ezqkdef(nigem,njgem ,'Z',ig1, ig2, ig3, ig4, iun)
***   Gets latlon values for gdid "gdgem"
 
      ier = gdgaxes(gdgem, ax, ay)

Go to the list of functions, or the general index.