program gemlatlon
  implicit none
  
  integer nigem, njgem
  
  integer ezqkdef, ezuvint, ezdefset, gdll, ezwdint
  
  character * 8   cle(1)
  character * 128 def(1), val(1)
  real a
  
  real, dimension(:,:), allocatable :: gemlat, gemlon
  
  integer gzkey,vvkey
  integer nigem, njgem, nkgem, gdin
  
  
  integer  ier, nrecs
  integer  fnom, fstouv, fclos, fstfrm, fstinf, fstluk, fstprm
  external fnom, fstouv, fclos, fstfrm, fstinf, fstluk, fstprm
  
  character*2 nomvar
  character*1 typvar,grtyp
  character*8 etiket 
  
  character *1 grsrc
  
  integer ip1,ip2,ip3
  integer extra1, extra2, extra3, deet, npas, nbits, datyp, dateo
  integer ig1src, ig2src, ig3src, ig4src, date
  integer swa, lng, dltf, ubc
  integer npts,nnpas,idt,ipos
  
  
  !     
  data cle /'IMFLDS.'/
  data def /'SCRAP'/
  data val /'SCRAP'/
  
  !------------------------------------------------------------------------------------------------

  call ccard(cle,def,val, 1, ipos)
  
  print *, val(1)
  ier = fnom(1, val(1),'STD+RND+OLD+R/O',0)
  nrecs = fstouv(1, 'RND')
  
  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))
  
  print *, ip1, ip2, ip3, typvar, nomvar, etiket, grsrc, ig1src, ig2src, ig3src, ig4src
  
  !------------------------------------------------------------------------------------------------
  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) 

  !------------------------------------------------------------------------------------------------
  stop
end program gemlatlon