scale_urban_grid_cartesC_index.F90

Go to the documentation of this file.

!-------------------------------------------------------------------------------
!-------------------------------------------------------------------------------
#include "scalelib.h"
module scale_urban_grid_cartesc_index
  !-----------------------------------------------------------------------------
  !
  !++ used modules
  !
  use scale_precision
  use scale_io
  use scale_prof
  !-----------------------------------------------------------------------------
  implicit none
  private
  !-----------------------------------------------------------------------------
  !
  !++ Public procedure
  !
  public :: urban_grid_cartesc_index_setup
 
  !-----------------------------------------------------------------------------
  !
  !++ Public parameters & variables
  !
  integer, public :: ukmax = -1 ! # of computational cells: z for urban
  integer, public :: uimax = -1 ! # of computational cells: x for urban
  integer, public :: ujmax = -1 ! # of computational cells: y for urban
  !$acc declare create(UKMAX, UIMAX, UJMAX)
 
  integer, public :: uka   = -1 ! # of total grids: z for urban, local
  integer, public :: uks        ! start point of inner domain: z for urban, local
  integer, public :: uke        ! end   point of inner domain: z for urban, local
 
  integer, public :: uia        ! # of total grids: x for urban, local
  integer, public :: uis        ! start point of inner domain: x for urban, local
  integer, public :: uie        ! end   point of inner domain: x for urban, local
 
  integer, public :: uja        ! # of total grids: Y for urban, local
  integer, public :: ujs        ! start point of inner domain: y for urban, local
  integer, public :: uje        ! end   point of inner domain: y for urban, local
  !$acc declare create(UKA, UKS, UKE, UIA, UIS, UIE, UJA, UJS, UJE)
 
  !-----------------------------------------------------------------------------
  !
  !++ Private procedure
  !
  !-----------------------------------------------------------------------------
  !
  !++ Private parameters & variables
  !
  !-----------------------------------------------------------------------------
contains
  !-----------------------------------------------------------------------------
  subroutine urban_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_urban_grid_cartesc_index / &
       ukmax
 
    integer :: ierr
    !---------------------------------------------------------------------------
 
    log_newline
    log_info("URBAN_GRID_CARTESC_INDEX_setup",*) 'Setup'
 
    !--- read namelist
    rewind(io_fid_conf)
    read(io_fid_conf,nml=param_urban_grid_cartesc_index,iostat=ierr)
    if( ierr < 0 ) then !--- missing
       log_info("URBAN_GRID_CARTESC_INDEX_setup",*) 'Not found namelist. Default used.'
    elseif( ierr > 0 ) then !--- fatal error
       log_error("URBAN_GRID_CARTESC_INDEX_setup",*) 'Not appropriate names in namelist PARAM_URBAN_GRID_CARTESC_INDEX. Check!'
       call prc_abort
    endif
    log_nml(param_urban_grid_cartesc_index)
 
    if ( ukmax < 1 ) then
       log_error("URBAN_GRID_CARTESC_INDEX_setup",*) 'UKMAX must be >= 1 ', ukmax
       call prc_abort
    end if
 
    uka  = ukmax
    uks  = 1
    uke  = ukmax
 
    log_newline
    log_info("URBAN_GRID_CARTESC_INDEX_setup",*) 'Urban grid index information '
    log_info_cont('(1x,A,I6,A,I6,A,I6)') 'z-axis levels :', ukmax
 
    ! at this moment horizontal grid is same as that in atmosphere
    uimax = imax
    uia = ia
    uis = is
    uie = ie
 
    ujmax = jmax
    uja = ja
    ujs = js
    uje = je
 
    !$acc update device(UKMAX, UIMAX, UJMAX)
 
    !$acc update device(UKA, UKS, UKE, UIA, UIS, UIE, UJA, UJS, UJE)
 
    return
  end subroutine urban_grid_cartesc_index_setup
 
end module scale_urban_grid_cartesc_index