program data_prep
c
c *** Program to prepare ETA initialization fields from any gridded
c        data set.
c
c *** Original code obtained from U. of Athens and modified at FSL.
c
      implicit none
c
      include 'ecommons.h'
c
      integer*4 n
      integer*4 nxin,nyin,nzin
c_______________________________________________________________________________
c 
c *** Fill eta model common blocks.
c
      common /sectorsize/nxin,nyin,nzin
      print*, 'Passou 1 lyra'
      
      call eta_commons
      
      
c *** Interpolate data to eta grid from each input init gridded dataset.
c
      do n=1,ninit
c
c ****** Get the native grid dimensions of the input data.
c           (Dims are filled into named common block - sectorsize).
c
!         call get_sector_size(init_in(n))
	write(6,*) 'call get_sector_size'
	write(6,*) 'file= ', init_gdsdir
	call get_sector_size(init_in(n),init_gdsdir)
c       
         
	  print*,nxin,nyin,nzin
          print*, 'Passou 2 lyra'
	  
	 call eta_interp(n)
c
      enddo
      print*, 'Passou 3 lyra'
     
c
c *** Create model init and constants files.
c
	write(6,*) 'calling eta_const'
      call eta_const
c
c *** Create model boundary condition files.
c
      call eta_boco
c
      end