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
  !$acc declare create(LKMAX, LIMAX, LJMAX)
 
  integer, public :: lka   = -1 ! # of total grids: z 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 :: lia        ! # of total grids: x 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 :: lja        ! # of total grids: y 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
  !$acc declare create(LKA, LIA, LJA, LKS, LKE, LIS, LIE, LJS, LJE)
 
  !-----------------------------------------------------------------------------
  !
  !++ 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
 
    !$acc update device(LKMAX, LIMAX, LJMAX)
    !$acc update device(LKA, LIA, LJA, LKS, LKE, LIS, LIE, LJS, LJE)
 
    return
  end subroutine land_grid_cartesc_index_setup
 
end module scale_land_grid_cartesc_index