gdll
|
FORTRAN integer gdll, ier ier = gdll(gdid, lat, lon) integer gdid real lat(*), lon(*) |
|
C int c_gdll (int gdid, float *lat, float *lon) |
Routine description
|
gdll returns the lat-lon positions of each grid points of the grid referred by gdid. |
Description of parameters
|
gdid |
a valid grid identifier (returned by ezqkdef or ezgdef) that will be used as the source (input) grid. |
|
lat, lon |
Stream of latitude and longitude points. The dimensions of the arrays lat() and lon() are assumed to be the same as the grid referred to gdid. |
|
x,y |
Position of the points on the grid, in grid coordinates. |
|
n |
Number of points |
?
Sample FORTRAN code :
The following code computes the lat-lon coordinates of each grid point on a GEM grid. It does essentially the same work as gdll. 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 = gdll(gdgem, lat, lon)
Go to the list of functions, or the general index.