module GRID (cartesian) for urban More...

Functions/Subroutines
subroutine, public	urban_grid_setup
	Setup. More...

Variables
real(rp), dimension(:), allocatable, public	grid_ucz
	center coordinate [m]: z, local=global More...

real(rp), dimension(:), allocatable, public	grid_ufz
	face coordinate [m]: z, local=global More...

real(rp), dimension(:), allocatable, public	grid_ucdz
	z-length of control volume [m] More...

Detailed Description

module GRID (cartesian) for urban

Description: Grid module for cartesian coordinate for urban

Author: Team SCALE

History

NAMELIST

PARAM_URBAN_GRID

name type default value comment

URBAN_GRID_IN_BASENAME character(len=H_LONG) ''

URBAN_GRID_OUT_BASENAME character(len=H_LONG) ''

UDZ real(RP), dimension(200)

History Output: No history output

Function/Subroutine Documentation

◆ urban_grid_setup()

subroutine, public scale_urban_grid::urban_grid_setup ( )

Setup.

Definition at line 61 of file scale_urban_grid.F90.

References grid_ucdz, grid_ucz, grid_ufz, scale_stdio::io_fid_conf, scale_stdio::io_fid_log, scale_stdio::io_fid_nml, scale_stdio::io_l, scale_stdio::io_nml, scale_process::prc_mpistop(), scale_process::prc_myrank, scale_urban_grid_index::uke, and scale_urban_grid_index::uks.

Referenced by mod_rm_driver::scalerm(), and mod_rm_prep::scalerm_prep().

     use scale_process, only: &
        prc_mpistop
     implicit none
 
     namelist / param_urban_grid / &
        urban_grid_in_basename,  &
        urban_grid_out_basename, &
        udz
 
     integer :: ierr
     integer :: k
     !---------------------------------------------------------------------------
 
     if( io_l ) write(io_fid_log,*)
     if( io_l ) write(io_fid_log,*) '++++++ Module[GRID] / Categ[URBAN GRID] / Origin[SCALElib]'
 
     udz(:) = 0.0_rp
 
     !--- read namelist
     rewind(io_fid_conf)
     read(io_fid_conf,nml=param_urban_grid,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_GRID. Check!'
        call prc_mpistop
     endif
     if( io_nml ) write(io_fid_nml,nml=param_urban_grid)
 
     allocate( grid_ucz(uks  :uke) )
     allocate( grid_ufz(uks-1:uke) )
     allocate( grid_ucdz(uks  :uke) )
 
     if( io_l ) write(io_fid_log,*)
     if( io_l ) write(io_fid_log,*) '*** Urban grid information ***'
 
     if ( urban_grid_in_basename /= '' ) then
        call urban_grid_read
     else
        if( io_l ) write(io_fid_log,*) '*** Not found input grid file. Grid position is calculated.'
 
        call urban_grid_generate
     endif
 
     if ( uke == uks ) then
        if( io_l ) write(io_fid_log,*)
        if( io_l ) write(io_fid_log,*) '*** Single layer. LDZ = ', udz(1)
     else
        if( io_l ) write(io_fid_log,*)
        if( io_l ) write(io_fid_log,'(1x,A)') &
        '|====== Vertical Coordinate ======|'
        if( io_l ) write(io_fid_log,'(1x,A)') &
        '|   k       z      zh      dz   k |'
        if( io_l ) write(io_fid_log,'(1x,A)') &
        '|         [m]     [m]     [m]     |'
        k = uks-1
        if( io_l ) write(io_fid_log,'(1x,A,F8.3,A,I4,A)') &
        '|            ',grid_ufz(k),'        ',k,' | Atmosphere interface'
        do k = uks, uke-1
        if( io_l ) write(io_fid_log,'(1x,A,I4,F8.3,A,F8.3,A)') &
        '|',k,grid_ucz(k),'        ',grid_ucdz(k),'     | '
        if( io_l ) write(io_fid_log,'(1x,A,F8.3,A,I4,A)') &
        '|            ',grid_ufz(k),'       |',k,' | '
        enddo
        k = uke
        if( io_l ) write(io_fid_log,'(1x,A,I4,F8.3,A,F8.3,A)') &
        '|',k,grid_ucz(k),'        ',grid_ucdz(k),'     | '
        if( io_l ) write(io_fid_log,'(1x,A,F8.3,A,I4,A)') &
        '|            ',grid_ufz(k),'        ',k,' | bedrock'
        if( io_l ) write(io_fid_log,'(1x,A)') &
        '|=================================|'
     endif
 
     return