scale_urban_grid_index.F90

Go to the documentation of this file.

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

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

  integer, public :: uks       ! start point of inner domain: z for urban, local
  integer, public :: uke       ! end   point of inner domain: z for urban, local

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

    namelist / param_urban_index / &
       ukmax

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

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

    !--- read namelist
    rewind(io_fid_conf)
    read(io_fid_conf,nml=param_urban_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_URBAN_INDEX. Check!'
       call prc_mpistop
    endif
    if( io_nml ) write(io_fid_nml,nml=param_urban_index)

    uks  = 1
    uke  = ukmax

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

    return
  end subroutine urban_grid_index_setup

end module scale_urban_grid_index