scale_land_grid_cartesC_index.F90

Go to the documentation of this file.

!-------------------------------------------------------------------------------
!-------------------------------------------------------------------------------
#include "scalelib.h"
module scale_land_grid_cartesc_index
  !-----------------------------------------------------------------------------
  !
  !++ used modules
  !
  use scale_precision
  use scale_io
  use scale_prof
  !-----------------------------------------------------------------------------
  implicit none
  private
  !-----------------------------------------------------------------------------
  !
  !++ Public procedure
  !
  public :: land_grid_cartesc_index_setup

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

  integer, public :: lka   = -1 ! # of total grids: z for land, local
  integer, public :: lia        ! # of total grids: x for land, local
  integer, public :: lja        ! # of total grids: y for land, local

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

  integer, public :: lie        ! end   point of inner domain: x for land, local
  integer, public :: ljs        ! start point of inner domain: y for land, local
  integer, public :: lje        ! end   point of inner domain: y for land, local

  !-----------------------------------------------------------------------------
  !
  !++ Private procedure
  !
  !-----------------------------------------------------------------------------
  !
  !++ Private parameters & variables
  !
  !-----------------------------------------------------------------------------
contains
  !-----------------------------------------------------------------------------
  subroutine land_grid_cartesc_index_setup
    use scale_prc, only: &
       prc_abort
    use scale_atmos_grid_cartesc_index, only: &
       imax,       &
       ia, is, ie, &
       jmax,       &
       ja, js, je
    implicit none

    namelist / param_land_grid_cartesc_index / &
       lkmax

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

    log_newline
    log_info("LAND_GRID_CARTESC_INDEX_setup",*) 'Setup'

    !--- read namelist
    rewind(io_fid_conf)
    read(io_fid_conf,nml=param_land_grid_cartesc_index,iostat=ierr)
    if( ierr < 0 ) then !--- missing
       log_info("LAND_GRID_CARTESC_INDEX_setup",*) 'Not found namelist. Default used.'
    elseif( ierr > 0 ) then !--- fatal error
       log_error("LAND_GRID_CARTESC_INDEX_setup",*) 'Not appropriate names in namelist PARAM_LAND_GRID_CARTESC_INDEX. Check!'
       call prc_abort
    endif
    log_nml(param_land_grid_cartesc_index)

    if ( lkmax < 1 ) then
       log_error("LAND_GRID_CARTESC_INDEX_setup",*) 'LKMAX must be >= 1 ', lkmax
       call prc_abort
    end if

    lks  = 1
    lke  = lkmax
    lka  = lkmax

    log_newline
    log_info("LAND_GRID_CARTESC_INDEX_setup",*) 'Land grid index information '
    log_info_cont('(1x,A,I6,A,I6,A,I6)') 'z-axis levels :', lkmax

    ! at this moment horizontal grid is same as that in atmosphere
    limax = imax
    lia = ia
    lis = is
    lie = ie

    ljmax = jmax
    lja = ja
    ljs = js
    lje = je

    return
  end subroutine land_grid_cartesc_index_setup

end module scale_land_grid_cartesc_index