Programme gem-latlon.f90

• gzkey=fstinf(1, nigem, njgem,nkgem,-1,' ',500,0,-1, ' ','GZ')

• ier=fstprm(gzkey, dateo, deet, npas, nigem, njgem, nkgem, nbits, datyp, &

• ip1, ip2, ip3, typvar, nomvar, etiket, &

• grsrc, ig1src, ig2src, ig3src, ig4src, &

• swa, lng, dltf, ubc, extra1, extra2, extra3)

• allocate(gemlat(nigem,njgem))

• allocate(gemlon(nigem,njgem))

• gdin = ezqkdef(nigem,njgem,grsrc,ig1src,ig2src,ig3src,ig4src,1)

• ier = gdll(gdin, gemlat, gemlon)

• call statfld('LA','P',ip1,ip2,ip3,date,etiket,gemlat,nigem, njgem, nkgem)

• call statfld('LO','P',ip1,ip2,ip3,date,etiket,gemlon,nigem, njgem, nkgem)