module urban / grid / icosahedralA More...

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

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

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

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

Detailed Description

module urban / grid / icosahedralA

Description: Grid module for icosahedral coordinate for urban

Author: Team SCALE

NAMELIST

PARAM_URBAN_GRID_ICOA

name type default value comment

URBAN_GRID_ICOA_IN_BASENAME character(len=H_LONG) ''

URBAN_GRID_ICOA_IN_AGGREGATE logical

UDZ real(RP), dimension(200)

History Output: No history output

Function/Subroutine Documentation

◆ urban_grid_icoa_setup()

subroutine, public scale_urban_grid_icoa::urban_grid_icoa_setup ( )

Setup.

Definition at line 60 of file scale_urban_grid_icoA.F90.

References scale_file::file_aggregate, scale_file::file_open(), scale_io::io_fid_conf, scale_prc::prc_abort(), scale_prc::prc_myrank, scale_urban_grid_icoa_index::uke, scale_urban_grid_icoa_index::ukmax, scale_urban_grid_icoa_index::uks, urban_grid_icoa_cdz, urban_grid_icoa_cz, and urban_grid_icoa_fz.

     use scale_prc, only: &
        prc_abort
     use scale_file, only: &
        file_aggregate
     implicit none
 
     namelist / param_urban_grid_icoa / &
        urban_grid_icoa_in_basename,  &
        urban_grid_icoa_in_aggregate, &
        udz
 
     integer :: ierr
     integer :: k
     !---------------------------------------------------------------------------
 
     log_newline
     log_info("URBAN_GRID_ICOA_setup",*) 'Setup'
 
     if ( ukmax < 1 ) then
        log_info("URBAN_GRID_ICOA_setup",*) 'Skip because UKMAX < 1'
        return
     end if
 
     udz(:) = 0.0_rp
 
     urban_grid_icoa_in_aggregate = file_aggregate
 
     !--- read namelist
     rewind(io_fid_conf)
     read(io_fid_conf,nml=param_urban_grid_icoa,iostat=ierr)
     if( ierr < 0 ) then !--- missing
        log_info("URBAN_GRID_ICOA_setup",*) 'Not found namelist. Default used.'
     elseif( ierr > 0 ) then !--- fatal error
        log_error("URBAN_GRID_ICOA_setup",*) 'Not appropriate names in namelist PARAM_URBAN_GRID_ICOA. Check!'
        call prc_abort
     endif
     log_nml(param_urban_grid_icoa)
 
     allocate( urban_grid_icoa_cz(uks  :uke) )
     allocate( urban_grid_icoa_fz(uks-1:uke) )
     allocate( urban_grid_icoa_cdz(uks  :uke) )
 
     log_newline
     log_info("URBAN_GRID_ICOA_setup",*) 'Urban grid information '
 
     if ( urban_grid_icoa_in_basename /= '' ) then
        call urban_grid_icoa_read
     else
        log_info("URBAN_GRID_ICOA_setup",*) 'Not found input grid file. Grid position is calculated.'
 
        call urban_grid_icoa_generate
     endif
 
     if ( uke == uks ) then
        log_newline
        log_info("URBAN_GRID_ICOA_setup",*) 'Single layer. LDZ = ', udz(1)
     else
        log_newline
        log_info("URBAN_GRID_ICOA_setup",'(1x,A)') 'Vertical Coordinate'
        log_info_cont('(1x,A)')                  '|   k       z      zh      dz   k |'
        log_info_cont('(1x,A)')                  '|         [m]     [m]     [m]     |'
        k = uks-1
        log_info_cont('(1x,A,F8.3,A,I4,A)')      '|            ',urban_grid_icoa_fz(k),'        ',k,' | Atmosphere interface'
        do k = uks, uke-1
        log_info_cont('(1x,A,I4,F8.3,A,F8.3,A)') '|',k,urban_grid_icoa_cz(k),'        ',urban_grid_icoa_cdz(k),'     | '
        log_info_cont('(1x,A,F8.3,A,I4,A)')      '|            ',urban_grid_icoa_fz(k),'       |',k,' | '
        enddo
        k = uke
        log_info_cont('(1x,A,I4,F8.3,A,F8.3,A)') '|',k,urban_grid_icoa_cz(k),'        ',urban_grid_icoa_cdz(k),'     | '
        log_info_cont('(1x,A,F8.3,A,I4,A)')      '|            ',urban_grid_icoa_fz(k),'        ',k,' | bedrock'
        log_info_cont('(1x,A)')                  '|=================================|'
     endif
 
     return