scale_urban_grid_cartesC_real.F90

Go to the documentation of this file.

!-------------------------------------------------------------------------------
!-------------------------------------------------------------------------------
#include "scalelib.h"
module scale_urban_grid_cartesc_real
  !-----------------------------------------------------------------------------
  !
  !++ used modules
  !
  use scale_precision
  use scale_io
  use scale_prof
  use scale_atmos_grid_cartesc_index
  use scale_urban_grid_cartesc_index
  !-----------------------------------------------------------------------------
  implicit none
  private
  !-----------------------------------------------------------------------------
  !
  !++ Public procedure
  !
  public :: urban_grid_cartesc_real_setup

  !-----------------------------------------------------------------------------
  !
  !++ Public parameters & variables
  !
  real(RP), public, allocatable :: urban_grid_cartesc_real_area(:,:)
  real(RP), public              :: urban_grid_cartesc_real_totarea
  real(RP), public, allocatable :: urban_grid_cartesc_real_vol (:,:,:)
  real(RP), public              :: urban_grid_cartesc_real_totvol

  !-----------------------------------------------------------------------------
  !
  !++ Private procedure
  !
  !-----------------------------------------------------------------------------
  !
  !++ Private parameters & variables
  !
  !-----------------------------------------------------------------------------
contains
  !-----------------------------------------------------------------------------
  subroutine urban_grid_cartesc_real_setup
    use scale_atmos_grid_cartesc_real, only: &
       atmos_grid_cartesc_real_area, &
       atmos_grid_cartesc_real_totarea
    use scale_urban_grid_cartesc, only: &
       urban_grid_cartesc_cdz
    use scale_file_cartesc, only: &
       file_cartesc_set_coordinates_urban

    integer :: k, i, j

    ! at this moment, horizontal grid is identical to that of the atmosphere
    allocate( urban_grid_cartesc_real_area(    uia,uja) )
    allocate( urban_grid_cartesc_real_vol(uka,uia,uja) )
    urban_grid_cartesc_real_area(:,:) = atmos_grid_cartesc_real_area(:,:)
    urban_grid_cartesc_real_totarea   = atmos_grid_cartesc_real_totarea

    do j = 1,   uja
    do i = 1,   uia
    do k = uks, uke
       urban_grid_cartesc_real_vol(k,i,j) = urban_grid_cartesc_real_area(i,j) * urban_grid_cartesc_cdz(k)
    enddo
    enddo
    enddo

    urban_grid_cartesc_real_totvol = 0.0_rp
    do j = ujs, uje
    do i = uis, uie
    do k = uks, uke
       urban_grid_cartesc_real_totvol = urban_grid_cartesc_real_totvol + urban_grid_cartesc_real_vol(k,i,j)
    end do
    end do
    end do

    call file_cartesc_set_coordinates_urban( urban_grid_cartesc_real_vol(:,:,:) ) ! [IN]

    return
  end subroutine urban_grid_cartesc_real_setup

end module scale_urban_grid_cartesc_real