scale_land_grid_index.F90

Go to the documentation of this file.

!-------------------------------------------------------------------------------
!-------------------------------------------------------------------------------
module scale_land_grid_index
  !-----------------------------------------------------------------------------
  !
  !++ used modules
  !
  use scale_precision
  use scale_stdio
  use scale_prof
  !-----------------------------------------------------------------------------
  implicit none
  private
  !-----------------------------------------------------------------------------
  !
  !++ Public procedure
  !
  public :: land_grid_index_setup

  !-----------------------------------------------------------------------------
  !
  !++ Public parameters & variables
  !
  integer, public :: lkmax = 1 ! # of computational cells: z for land

  integer, public :: lks       ! start point of inner domain: z for land, local
  integer, public :: lke       ! end   point of inner domain: z for land, local

  !-----------------------------------------------------------------------------
  !
  !++ Private procedure
  !
  !-----------------------------------------------------------------------------
  !
  !++ Private parameters & variables
  !
  !-----------------------------------------------------------------------------
contains
  !-----------------------------------------------------------------------------
  subroutine land_grid_index_setup
    use scale_process, only: &
       prc_mpistop
    implicit none

    namelist / param_land_index / &
       lkmax

    integer :: ierr
    !---------------------------------------------------------------------------

    if( io_l ) write(io_fid_log,*)
    if( io_l ) write(io_fid_log,*) '++++++ Module[GRID_INDEX] / Categ[LAND GRID] / Origin[SCALElib]'

    !--- read namelist
    rewind(io_fid_conf)
    read(io_fid_conf,nml=param_land_index,iostat=ierr)
    if( ierr < 0 ) then !--- missing
       if( io_l ) write(io_fid_log,*) '*** Not found namelist. Default used.'
    elseif( ierr > 0 ) then !--- fatal error
       write(*,*) 'xxx Not appropriate names in namelist PARAM_LAND_INDEX. Check!'
       call prc_mpistop
    endif
    if( io_nml ) write(io_fid_nml,nml=param_land_index)

    lks  = 1
    lke  = lkmax

    if( io_l ) write(io_fid_log,*)
    if( io_l ) write(io_fid_log,*) '*** Land grid index information ***'
    if( io_l ) write(io_fid_log,'(1x,A,I6,A,I6,A,I6)') '*** z-axis levels :', lkmax

    return
  end subroutine land_grid_index_setup

end module scale_land_grid_index