merge_lake_ocnmsk.F90

Go to the documentation of this file.

program merge_lake_ocnmsk

  implicit none

  character(len=200) :: pth1
  character(len=200) :: pth2,pth3
  character(len=10)  :: atmres,ocnres

  integer :: binary_lake
  integer :: lat,lon,tile

  real, allocatable :: lake_frac(:,:),lake_depth(:,:),land_frac(:,:),ocn_frac(:,:),slmsk(:,:),lat2d(:,:)

  print*,"- BEGIN OCEAN MERGE PROGRAM."

  call read_nml(pth1, pth2, atmres, ocnres, pth3,binary_lake)

  do tile=1,6

    call read_grid_dims(pth1, atmres, ocnres, tile, lon, lat)

    if (tile==1) then
      write(6,*) 'lat=',lat,'lon=',lon
      allocate (lake_frac(lon,lat),lake_depth(lon,lat),land_frac(lon,lat), &
                ocn_frac(lon,lat),slmsk(lon,lat),lat2d(lon,lat))
    endif

    call read_ocean_frac(pth1,atmres,ocnres,tile,lon,lat,ocn_frac)

    call read_lake_mask(pth2,atmres,tile,lon,lat,lake_frac, &
                        lake_depth,lat2d)

    call merge(lon, lat, binary_lake, lat2d, ocn_frac, lake_frac, lake_depth, land_frac, slmsk)

    call write_data(atmres,ocnres,pth3,tile,lon,lat,land_frac, &
                    lake_frac,lake_depth,slmsk)

  end do ! tile

  deallocate (lake_frac,lake_depth,land_frac,ocn_frac,slmsk,lat2d)

  print*,"- NORMAL TERMINATION."

end program merge_lake_ocnmsk