mod_urban_driver Module Reference

module URBAN driver More...

Functions/Subroutines
subroutine, public	urban_driver_setup
	Setup. More...
subroutine, public	urban_driver_calc_tendency (force)
	Calclate tendency. More...
subroutine, public	urban_driver_update
	Urban step. More...
subroutine, public	urban_surface_get
	Get surface boundary. More...
subroutine, public	urban_surface_set (countup)
	Set surface boundary to other model. More...

Detailed Description

module URBAN driver

Description

Urban module driver

Author

Team SCALE

NAMELIST

• No namelist group

History Output

name	description	unit	variable
URBAN_QC_t	tendency of URBAN_QC	kg/kg/s	URBAN_QC_t
URBAN_RAINB_t	tendency of URBAN_RAINB	kg/m2/s	URBAN_RAINB_t
URBAN_RAING_t	tendency of URBAN_RAING	kg/m2/s	URBAN_RAING_t
URBAN_RAINR_t	tendency of URBAN_RAINR	kg/m2/s	URBAN_RAINR_t
URBAN_ROFF	urban runoff water	kg/m2/s	URBAN_ROFF
URBAN_TBL_t	tendency of URBAN_TBL	K/s	URBAN_TBL_t
URBAN_TB_t	tendency of URBAN_TB	K/s	URBAN_TB_t
URBAN_TC_t	tendency of URBAN_TC	K/s	URBAN_TC_t
URBAN_TGL_t	tendency of URBAN_TGL	K/s	URBAN_TGL_t
URBAN_TG_t	tendency of URBAN_TG	K/s	URBAN_TG_t
URBAN_TRL_t	tendency of URBAN_TRL	K/s	URBAN_TRL_t
URBAN_TR_t	tendency of URBAN_TR	K/s	URBAN_TR_t
URBAN_UC_t	tendency of URBAN_UC	m/s2	URBAN_UC_t

Function/Subroutine Documentation

urban_driver_setup()

subroutine, public mod_urban_driver::urban_driver_setup

Setup.

Definition at line 54 of file mod_urban_driver.F90.

    use scale_prc, only: &
       prc_abort
    use scale_const, only: &
       undef => const_undef
    use mod_urban_admin, only: &
       urban_do, &
       urban_dyn_type, &
       urban_sfc_type
    use scale_urban_dyn_kusaka01, only: &
       urban_dyn_kusaka01_setup
    use scale_landuse, only: &
       landuse_fact_urban
    use mod_urban_vars, only: &
       urban_z0m, &
       urban_z0h, &
       urban_z0e, &
       urban_zd
    implicit none
    !---------------------------------------------------------------------------
 
    log_newline
    log_info("URBAN_driver_setup",*) 'Setup'
 
    if ( urban_do ) then
 
       allocate( ah_urb(uia,uja,1:24) )
       allocate( ahl_urb(uia,uja,1:24) )
       ah_urb(:,:,:) = undef
       ahl_urb(:,:,:) = undef
 
       select case ( urban_dyn_type )
       case ( 'KUSAKA01' )
          call urban_dyn_kusaka01_setup( uia, uis, uie, uja, ujs, uje,                   & ! [IN]
                                         landuse_fact_urban(:,:),                        & ! [IN]
                                         urban_z0m(:,:), urban_z0h(:,:), urban_z0e(:,:), & ! [OUT]
                                         urban_zd(:,:),                                  & ! [OUT]
                                         ah_urb(:,:,:), ahl_urb(:,:,:)                   ) ! [OUT]
 
          urban_sfc_type = 'KUSAKA01'
       case default
          log_error("URBAN_driver_setup",*) 'LAND_DYN_TYPE is invalid: ', trim(urban_dyn_type)
          call prc_abort
       end select
 
       select case ( urban_sfc_type )
       case ( 'KUSAKA01' )
          ! do nothing
       case default
          log_error("URBAN_driver_setup",*) 'LAND_SFC_TYPE is invalid: ', trim(urban_sfc_type)
          call prc_abort
       end select
 
    end if
 
    return

References scale_const::const_undef, scale_landuse::landuse_fact_urban, scale_prc::prc_abort(), scale_urban_grid_cartesc_index::uia, scale_urban_grid_cartesc_index::uie, scale_urban_grid_cartesc_index::uis, scale_urban_grid_cartesc_index::uja, scale_urban_grid_cartesc_index::uje, scale_urban_grid_cartesc_index::ujs, mod_urban_admin::urban_do, scale_urban_dyn_kusaka01::urban_dyn_kusaka01_setup(), mod_urban_admin::urban_dyn_type, mod_urban_admin::urban_sfc_type, mod_urban_vars::urban_z0e, mod_urban_vars::urban_z0h, mod_urban_vars::urban_z0m, and mod_urban_vars::urban_zd.

Referenced by mod_rm_driver::rm_driver().

Here is the call graph for this function:

Here is the caller graph for this function:

urban_driver_calc_tendency()

subroutine, public mod_urban_driver::urban_driver_calc_tendency

(

logical, intent(in)

force

)

Calclate tendency.

Definition at line 114 of file mod_urban_driver.F90.

    use scale_const, only: &
       d2r => const_d2r
    use scale_statistics, only: &
       statistics_checktotal, &
       statistics_total
    use scale_urban_grid_cartesc_real, only: &
       urban_grid_cartesc_real_vol,    &
       urban_grid_cartesc_real_totvol, &
       urban_grid_cartesc_real_area,   &
       urban_grid_cartesc_real_totarea
    use scale_topography, only: &
       tansl_x => topography_tansl_x, &
       tansl_y => topography_tansl_y
    use scale_file_history, only: &
       file_history_in
    use mod_atmos_admin, only: &
       atmos_sw_phy_ch
    use mod_atmos_phy_ch_driver, only: &
       atmos_phy_ch_driver_urban_flux
    use mod_urban_vars, only: &
       atmos_temp,      &
       atmos_pres,      &
       atmos_u,         &
       atmos_v,         &
       atmos_dens,      &
       atmos_qv,        &
       atmos_pbl,       &
       atmos_sfc_dens,  &
       atmos_sfc_pres,  &
       atmos_sflx_lw,   &
       atmos_sflx_sw,   &
       atmos_cossza,    &
       atmos_sflx_water, &
       atmos_sflx_engi, &
       urban_trl_t,     &
       urban_tbl_t,     &
       urban_tgl_t,     &
       urban_tr_t,      &
       urban_tb_t,      &
       urban_tg_t,      &
       urban_tc_t,      &
       urban_qc_t,      &
       urban_uc_t,      &
       urban_rainr_t,   &
       urban_rainb_t,   &
       urban_raing_t,   &
       urban_sfc_temp,    &
       urban_sfc_albedo,  &
       urban_sflx_mw,     &
       urban_sflx_mu,     &
       urban_sflx_mv,     &
       urban_sflx_sh,     &
       urban_sflx_lh,     &
       urban_sflx_shex,   &
       urban_sflx_qvex,   &
       urban_sflx_qtrc,   &
       urban_sflx_gh,     &
       urban_z0m,         &
       urban_z0h,         &
       urban_z0e,         &
       urban_zd,          &
       urban_ah,          &
       urban_ahl,         &
       urban_ustar,       &
       urban_tstar,       &
       urban_qstar,       &
       urban_wstar,       &
       urban_rlmo,        &
       urban_u10,         &
       urban_v10,         &
       urban_t2,          &
       urban_q2,          &
       urban_tr,          &
       urban_tb,          &
       urban_tg,          &
       urban_tc,          &
       urban_qc,          &
       urban_uc,          &
       urban_trl,         &
       urban_tbl,         &
       urban_tgl,         &
       urban_rainr,       &
       urban_rainb,       &
       urban_raing,       &
       urban_roff
    use scale_atmos_hydrometeor, only: &
       hydrometeor_lhv => atmos_hydrometeor_lhv, &
       atmos_hydrometeor_dry,                    &
       i_qv
    use scale_time, only: &
       dt => time_dtsec_urban, &
       nowdate => time_nowdate
    use scale_mapprojection, only: &
       base_lon => mapprojection_basepoint_lon, &
       base_lat => mapprojection_basepoint_lat
    use scale_atmos_grid_cartesc_real, only: &
       real_z1 => atmos_grid_cartesc_real_z1
    use scale_urban_grid_cartesc, only: &
       cdz => urban_grid_cartesc_cdz
    use scale_landuse, only: &
       landuse_fact_urban, &
       exists_urban => landuse_exists_urban
    use mod_urban_admin, only: &
       urban_sfc_type
    use scale_urban_dyn_kusaka01, only: &
       urban_dyn_kusaka01
    implicit none
    logical, intent(in) :: force
 
    real(RP) :: TRL(UKA,UIA,UJA), TBL(UKA,UIA,UJA), TGL(UKA,UIA,UJA)
    real(RP) :: TR(UIA,UJA), TB(UIA,UJA), TG(UIA,UJA)
    real(RP) :: TC(UIA,UJA), QC(UIA,UJA), UC(UIA,UJA)
    real(RP) :: RAINR(UIA,UJA), RAINB(UIA,UJA), RAING(UIA,UJA), ROFF(UIA,UJA)
 
    real(RP) :: LHV(UIA,UJA) ! latent heat of vaporization [J/kg]
 
    real(RP) :: URBAN_SFLX_LHEX(UIA,UJA)
 
    real(RP) :: LAT, LON ! [deg]
    integer  :: tloc     ! local time (1-24h)
    real(RP) :: dsec     ! second [s]
 
    integer :: k, i, j, iq
    !---------------------------------------------------------------------------
 
    call prof_rapstart('URB_CalcTend', 1)
 
    !########## Get Surface Boundary from coupler ##########
    call urban_surface_get
 
    !########## initialize tendency ##########
!OCL XFILL
    do j = ujs, uje
    do i = uis, uie
    do k = uks, uke
       urban_trl_t(k,i,j) = 0.0_rp
       urban_tbl_t(k,i,j) = 0.0_rp
       urban_tgl_t(k,i,j) = 0.0_rp
    end do
    end do
    end do
 
!OCL XFILL
    do j = ujs, uje
    do i = uis, uie
       urban_tr_t(i,j) = 0.0_rp
       urban_tb_t(i,j) = 0.0_rp
       urban_tg_t(i,j) = 0.0_rp
       urban_tc_t(i,j) = 0.0_rp
       urban_qc_t(i,j) = 0.0_rp
       urban_uc_t(i,j) = 0.0_rp
 
       urban_rainr_t(i,j) = 0.0_rp
       urban_rainb_t(i,j) = 0.0_rp
       urban_raing_t(i,j) = 0.0_rp
    enddo
    enddo
 
!OCL XFILL
    !$omp parallel do
    do iq = 1, qa
    do j  = ujs, uje
    do i  = uis, uie
       urban_sflx_qtrc(i,j,iq) = 0.0_rp
    enddo
    enddo
    enddo
 
    select case ( urban_sfc_type )
    case ( 'KUSAKA01' )
 
!OCL XFILL
       !$omp parallel do
       do j = ujs, uje
       do i = uis, uie
       do k = uks, uke
          trl(k,i,j) = urban_trl(k,i,j)
          tbl(k,i,j) = urban_tbl(k,i,j)
          tgl(k,i,j) = urban_tgl(k,i,j)
       end do
       end do
       end do
 
!OCL XFILL
       !$omp parallel do
       do j = ujs, uje
       do i = uis, uie
          tr(i,j) = urban_tr(i,j)
          tb(i,j) = urban_tb(i,j)
          tg(i,j) = urban_tg(i,j)
          tc(i,j) = urban_tc(i,j)
          qc(i,j) = urban_qc(i,j)
          uc(i,j) = urban_uc(i,j)
          rainr(i,j) = urban_rainr(i,j)
          rainb(i,j) = urban_rainb(i,j)
          raing(i,j) = urban_raing(i,j)
       end do
       end do
 
 
       ! local time
       lat = base_lat
       lon = base_lon
       if (lon < 0.0_rp )   lon = mod(lon, 360.0_rp) + 360.0_rp
       if (lon > 360.0_rp ) lon = mod(lon, 360.0_rp)
       tloc = mod( (nowdate(4) + int(lon/15.0_rp)),24 )
       dsec = real( nowdate(5)*60.0_rp + nowdate(6), kind=rp ) / 3600.0_rp
       if( tloc == 0 ) tloc = 24
 
       !--- Calculate AH at LST
       if ( tloc == 24 ) then
          do j = ujs, uje
          do i = uis, uie
          if ( exists_urban(i,j) ) then
             urban_ah(i,j)  = ( 1.0_rp-dsec ) * ah_urb(i,j,tloc  ) &
                            + (        dsec ) * ah_urb(i,j,1     )
             urban_ahl(i,j) = ( 1.0_rp-dsec ) * ahl_urb(i,j,tloc  ) &
                            + (        dsec ) * ahl_urb(i,j,1     )
          end if
          enddo
          enddo
       else
          do j = ujs, uje
          do i = uis, uie
          if ( exists_urban(i,j) ) then
             urban_ah(i,j)  = ( 1.0_rp-dsec ) * ah_urb(i,j,tloc  ) &
                            + (        dsec ) * ah_urb(i,j,tloc+1)
             urban_ahl(i,j) = ( 1.0_rp-dsec ) * ahl_urb(i,j,tloc  ) &
                            + (        dsec ) * ahl_urb(i,j,tloc+1)
          end if
          enddo
          enddo
       endif
 
       call hydrometeor_lhv( uia, uis, uie, uja, ujs, uje, &
                             atmos_temp(:,:), lhv(:,:) )
 
       call urban_dyn_kusaka01( uka, uks, uke, uia, uis, uie, uja, ujs, uje, &
                                atmos_temp(:,:), atmos_pres(:,:),                            & ! [IN]
                                atmos_u(:,:), atmos_v(:,:),                                  & ! [IN]
                                atmos_dens(:,:), atmos_qv(:,:), lhv(:,:),                    & ! [IN]
                                real_z1(:,:), atmos_pbl(:,:),                                & ! [IN]
                                atmos_sfc_dens(:,:), atmos_sfc_pres(:,:),                    & ! [IN]
                                atmos_sflx_lw(:,:,:), atmos_sflx_sw(:,:,:),                  & ! [IN]
                                atmos_sflx_water(:,:), atmos_sflx_engi(:,:),                 & ! [IN]
                                urban_z0m(:,:), urban_z0h(:,:), urban_z0e(:,:),              & ! [IN]
                                urban_zd(:,:),                                               & ! [IN]
                                cdz(:),                                                      & ! [IN]
                                tansl_x(:,:), tansl_y(:,:),                                  & ! [IN]
                                landuse_fact_urban(:,:),                                     & ! [IN]
                                dt,                                                          & ! [IN]
                                trl(:,:,:), tbl(:,:,:), tgl(:,:,:),                          & ! [INOUT]
                                tr(:,:), tb(:,:), tg(:,:), tc(:,:), qc(:,:), uc(:,:),        & ! [INOUT]
                                rainr(:,:), rainb(:,:), raing(:,:), urban_roff(:,:),         & ! [INOUT]
                                urban_sfc_temp(:,:),                                         & ! [OUT]
                                urban_sfc_albedo(:,:,:,:),                                   & ! [OUT]
                                urban_sflx_mw(:,:), urban_sflx_mu(:,:), urban_sflx_mv(:,:),  & ! [OUT]
                                urban_sflx_sh(:,:), urban_sflx_lh(:,:), urban_sflx_gh(:,:),  & ! [OUT]
                                urban_ustar(:,:), urban_tstar(:,:), urban_qstar(:,:),        & ! [OUT]
                                urban_wstar(:,:),                                            & ! [OUT]
                                urban_rlmo(:,:),                                             & ! [OUT]
                                urban_u10(:,:), urban_v10(:,:), urban_t2(:,:), urban_q2(:,:) ) ! [OUT]
 
       !-----------------------------------------------------------
       ! anthropogenic heat fluxes
       !-----------------------------------------------------------
       !$omp parallel do
       do j = ujs, uje
       do i = uis, uie
       if ( exists_urban(i,j) ) then
          urban_sflx_shex(i,j) = urban_ah(i,j) / landuse_fact_urban(i,j) ! Sensible heat flux [W/m2]
          urban_sflx_lhex(i,j) = urban_ahl(i,j) / landuse_fact_urban(i,j) ! Latent heat flux   [W/m2]
          urban_sflx_sh(i,j) = urban_sflx_sh(i,j) + urban_sflx_shex(i,j)
          urban_sflx_lh(i,j) = urban_sflx_lh(i,j) + urban_sflx_lhex(i,j)
       end if
       end do
       end do
 
!OCL XFILL
       !$omp parallel do
       do j = ujs, uje
       do i = uis, uie
       if ( exists_urban(i,j) ) then
          do k = uks, uke
             urban_trl_t(k,i,j) = ( trl(k,i,j) - urban_trl(k,i,j) ) / dt
             urban_tbl_t(k,i,j) = ( tbl(k,i,j) - urban_tbl(k,i,j) ) / dt
             urban_tgl_t(k,i,j) = ( tgl(k,i,j) - urban_tgl(k,i,j) ) / dt
          end do
       end if
       end do
       end do
 
!OCL XFILL
       !$omp parallel do
       do j = ujs, uje
       do i = uis, uie
       if ( exists_urban(i,j) ) then
          urban_tr_t(i,j) = ( tr(i,j) - urban_tr(i,j) ) / dt
          urban_tb_t(i,j) = ( tb(i,j) - urban_tb(i,j) ) / dt
          urban_tg_t(i,j) = ( tg(i,j) - urban_tg(i,j) ) / dt
          urban_tc_t(i,j) = ( tc(i,j) - urban_tc(i,j) ) / dt
          urban_qc_t(i,j) = ( qc(i,j) - urban_qc(i,j) ) / dt
          urban_uc_t(i,j) = ( uc(i,j) - urban_uc(i,j) ) / dt
          urban_rainr_t(i,j) = ( rainr(i,j) - urban_rainr(i,j) ) / dt
          urban_rainb_t(i,j) = ( rainb(i,j) - urban_rainb(i,j) ) / dt
          urban_raing_t(i,j) = ( raing(i,j) - urban_raing(i,j) ) / dt
       end if
       end do
       end do
 
       if ( .NOT. atmos_hydrometeor_dry ) then
          !$omp parallel do
          do j = ujs, uje
          do i = uis, uie
          if ( exists_urban(i,j) ) then
             urban_sflx_qtrc(i,j,i_qv) = urban_sflx_lh(i,j) / lhv(i,j)
             urban_sflx_qvex(i,j)      = urban_sflx_lhex(i,j) / lhv(i,j)
          end if
          enddo
          enddo
       endif
 
    end select
 
    ! Surface flux for chemical tracers
    if ( atmos_sw_phy_ch ) then
       call atmos_phy_ch_driver_urban_flux( urban_sflx_qtrc(:,:,:) ) ! [INOUT]
    endif
 
    call file_history_in( urban_tr_t(:,:), 'URBAN_TR_t', 'tendency of URBAN_TR', 'K/s',     dim_type='XY' )
    call file_history_in( urban_tb_t(:,:), 'URBAN_TB_t', 'tendency of URBAN_TB', 'K/s',     dim_type='XY' )
    call file_history_in( urban_tg_t(:,:), 'URBAN_TG_t', 'tendency of URBAN_TG', 'K/s',     dim_type='XY' )
    call file_history_in( urban_tc_t(:,:), 'URBAN_TC_t', 'tendency of URBAN_TC', 'K/s',     dim_type='XY' )
    call file_history_in( urban_qc_t(:,:), 'URBAN_QC_t', 'tendency of URBAN_QC', 'kg/kg/s', dim_type='XY' )
    call file_history_in( urban_uc_t(:,:), 'URBAN_UC_t', 'tendency of URBAN_UC', 'm/s2',    dim_type='XY' )
 
    call file_history_in( urban_trl_t(:,:,:), 'URBAN_TRL_t', 'tendency of URBAN_TRL', 'K/s', dim_type='UXY' )
    call file_history_in( urban_tbl_t(:,:,:), 'URBAN_TBL_t', 'tendency of URBAN_TBL', 'K/s', dim_type='UXY' )
    call file_history_in( urban_tgl_t(:,:,:), 'URBAN_TGL_t', 'tendency of URBAN_TGL', 'K/s', dim_type='UXY' )
 
    call file_history_in( urban_rainr_t(:,:), 'URBAN_RAINR_t', 'tendency of URBAN_RAINR', 'kg/m2/s', dim_type='XY' )
    call file_history_in( urban_rainb_t(:,:), 'URBAN_RAINB_t', 'tendency of URBAN_RAINB', 'kg/m2/s', dim_type='XY' )
    call file_history_in( urban_raing_t(:,:), 'URBAN_RAING_t', 'tendency of URBAN_RAING', 'kg/m2/s', dim_type='XY' )
    call file_history_in( urban_roff(:,:), 'URBAN_ROFF',    'urban runoff water',      'kg/m2/s', dim_type='XY' )
 
    if ( statistics_checktotal ) then
 
       call statistics_total( uka, uks, uke, uia, uis, uie, uja, ujs, uje, &
                              urban_trl_t(:,:,:), 'URBAN_TRL_t', &
                              urban_grid_cartesc_real_vol(:,:,:), &
                              urban_grid_cartesc_real_totvol      )
       call statistics_total( uka, uks, uke, uia, uis, uie, uja, ujs, uje, &
                              urban_tbl_t(:,:,:), 'URBAN_TBL_t', &
                              urban_grid_cartesc_real_vol(:,:,:), &
                              urban_grid_cartesc_real_totvol      )
       call statistics_total( uka, uks, uke, uia, uis, uie, uja, ujs, uje, &
                              urban_tgl_t(:,:,:), 'URBAN_TGL_t', &
                              urban_grid_cartesc_real_vol(:,:,:), &
                              urban_grid_cartesc_real_totvol      )
 
       call statistics_total( uia, uis, uie, uja, ujs, uje, &
                              urban_tr_t(:,:), 'URBAN_TR_t',     &
                              urban_grid_cartesc_real_area(:,:), &
                              urban_grid_cartesc_real_totarea    )
       call statistics_total( uia, uis, uie, uja, ujs, uje, &
                              urban_tb_t(:,:), 'URBAN_TB_t',     &
                              urban_grid_cartesc_real_area(:,:), &
                              urban_grid_cartesc_real_totarea    )
       call statistics_total( uia, uis, uie, uja, ujs, uje, &
                              urban_tg_t(:,:), 'URBAN_TG_t',     &
                              urban_grid_cartesc_real_area(:,:), &
                              urban_grid_cartesc_real_totarea    )
       call statistics_total( uia, uis, uie, uja, ujs, uje, &
                              urban_tc_t(:,:), 'URBAN_TC_t',     &
                              urban_grid_cartesc_real_area(:,:), &
                              urban_grid_cartesc_real_totarea    )
       call statistics_total( uia, uis, uie, uja, ujs, uje, &
                              urban_qc_t(:,:), 'URBAN_QC_t',     &
                              urban_grid_cartesc_real_area(:,:), &
                              urban_grid_cartesc_real_totarea    )
       call statistics_total( uia, uis, uie, uja, ujs, uje, &
                              urban_uc_t(:,:), 'URBAN_UC_t',     &
                              urban_grid_cartesc_real_area(:,:), &
                              urban_grid_cartesc_real_totarea    )
 
       call statistics_total( uia, uis, uie, uja, ujs, uje, &
                              urban_rainr_t(:,:), 'URBAN_RAINR_t', &
                              urban_grid_cartesc_real_area(:,:),   &
                              urban_grid_cartesc_real_totarea      )
       call statistics_total( uia, uis, uie, uja, ujs, uje, &
                              urban_rainb_t(:,:), 'URBAN_RAINB_t', &
                              urban_grid_cartesc_real_area(:,:),   &
                              urban_grid_cartesc_real_totarea      )
       call statistics_total( uia, uis, uie, uja, ujs, uje, &
                              urban_raing_t(:,:), 'URBAN_RAING_t', &
                              urban_grid_cartesc_real_area(:,:),   &
                              urban_grid_cartesc_real_totarea      )
       call statistics_total( uia, uis, uie, uja, ujs, uje, &
                              urban_roff(:,:),    'URBAN_ROFF',    &
                              urban_grid_cartesc_real_area(:,:),   &
                              urban_grid_cartesc_real_totarea      )
    endif
 
 
    call prof_rapend  ('URB_CalcTend', 1)
 
    !########## Set Surface Boundary to coupler ##########
    call urban_surface_set( countup=.true. )
 
    return

References mod_urban_vars::atmos_cossza, mod_urban_vars::atmos_dens, scale_atmos_grid_cartesc_real::atmos_grid_cartesc_real_z1, scale_atmos_hydrometeor::atmos_hydrometeor_dry, mod_urban_vars::atmos_pbl, mod_atmos_phy_ch_driver::atmos_phy_ch_driver_urban_flux(), mod_urban_vars::atmos_pres, mod_urban_vars::atmos_qv, mod_urban_vars::atmos_sfc_dens, mod_urban_vars::atmos_sfc_pres, mod_urban_vars::atmos_sflx_engi, mod_urban_vars::atmos_sflx_lw, mod_urban_vars::atmos_sflx_sw, mod_urban_vars::atmos_sflx_water, mod_atmos_admin::atmos_sw_phy_ch, mod_urban_vars::atmos_temp, mod_urban_vars::atmos_u, mod_urban_vars::atmos_v, scale_const::const_d2r, scale_atmos_hydrometeor::i_qv, scale_tracer::k, scale_landuse::landuse_exists_urban, scale_landuse::landuse_fact_urban, scale_mapprojection::mapprojection_basepoint_lat, scale_mapprojection::mapprojection_basepoint_lon, scale_prof::prof_rapend(), scale_prof::prof_rapstart(), scale_tracer::qa, scale_statistics::statistics_checktotal, scale_time::time_dtsec_urban, scale_time::time_nowdate, scale_topography::topography_tansl_x, scale_topography::topography_tansl_y, scale_urban_grid_cartesc_index::uia, scale_urban_grid_cartesc_index::uie, scale_urban_grid_cartesc_index::uis, scale_urban_grid_cartesc_index::uja, scale_urban_grid_cartesc_index::uje, scale_urban_grid_cartesc_index::ujs, scale_urban_grid_cartesc_index::uka, scale_urban_grid_cartesc_index::uke, scale_urban_grid_cartesc_index::uks, mod_urban_vars::urban_ah, mod_urban_vars::urban_ahl, scale_urban_dyn_kusaka01::urban_dyn_kusaka01(), scale_urban_grid_cartesc::urban_grid_cartesc_cdz, scale_urban_grid_cartesc_real::urban_grid_cartesc_real_area, scale_urban_grid_cartesc_real::urban_grid_cartesc_real_totarea, scale_urban_grid_cartesc_real::urban_grid_cartesc_real_totvol, scale_urban_grid_cartesc_real::urban_grid_cartesc_real_vol, mod_urban_vars::urban_q2, mod_urban_vars::urban_qc, mod_urban_vars::urban_qc_t, mod_urban_vars::urban_qstar, mod_urban_vars::urban_rainb, mod_urban_vars::urban_rainb_t, mod_urban_vars::urban_raing, mod_urban_vars::urban_raing_t, mod_urban_vars::urban_rainr, mod_urban_vars::urban_rainr_t, mod_urban_vars::urban_rlmo, mod_urban_vars::urban_roff, mod_urban_vars::urban_sfc_albedo, mod_urban_vars::urban_sfc_temp, mod_urban_admin::urban_sfc_type, mod_urban_vars::urban_sflx_gh, mod_urban_vars::urban_sflx_lh, mod_urban_vars::urban_sflx_mu, mod_urban_vars::urban_sflx_mv, mod_urban_vars::urban_sflx_mw, mod_urban_vars::urban_sflx_qtrc, mod_urban_vars::urban_sflx_qvex, mod_urban_vars::urban_sflx_sh, mod_urban_vars::urban_sflx_shex, urban_surface_get(), urban_surface_set(), mod_urban_vars::urban_t2, mod_urban_vars::urban_tb, mod_urban_vars::urban_tb_t, mod_urban_vars::urban_tbl, mod_urban_vars::urban_tbl_t, mod_urban_vars::urban_tc, mod_urban_vars::urban_tc_t, mod_urban_vars::urban_tg, mod_urban_vars::urban_tg_t, mod_urban_vars::urban_tgl, mod_urban_vars::urban_tgl_t, mod_urban_vars::urban_tr, mod_urban_vars::urban_tr_t, mod_urban_vars::urban_trl, mod_urban_vars::urban_trl_t, mod_urban_vars::urban_tstar, mod_urban_vars::urban_u10, mod_urban_vars::urban_uc, mod_urban_vars::urban_uc_t, mod_urban_vars::urban_ustar, mod_urban_vars::urban_v10, mod_urban_vars::urban_wstar, mod_urban_vars::urban_z0e, mod_urban_vars::urban_z0h, mod_urban_vars::urban_z0m, and mod_urban_vars::urban_zd.

Referenced by mod_rm_driver::restart_read(), and mod_rm_driver::rm_driver().

Here is the call graph for this function:

Here is the caller graph for this function:

urban_driver_update()

subroutine, public mod_urban_driver::urban_driver_update

Urban step.

Definition at line 529 of file mod_urban_driver.F90.

    use scale_time, only: &
       dt => time_dtsec_urban
    use mod_urban_vars, only: &
       urban_trl_t,       &
       urban_tbl_t,       &
       urban_tgl_t,       &
       urban_tr_t,        &
       urban_tb_t,        &
       urban_tg_t,        &
       urban_tc_t,        &
       urban_qc_t,        &
       urban_uc_t,        &
       urban_rainr_t,     &
       urban_rainb_t,     &
       urban_raing_t,     &
       urban_tr,          &
       urban_tb,          &
       urban_tg,          &
       urban_tc,          &
       urban_qc,          &
       urban_uc,          &
       urban_trl,         &
       urban_tbl,         &
       urban_tgl,         &
       urban_rainr,       &
       urban_rainb,       &
       urban_raing,       &
       urban_vars_check
    use mod_urban_admin, only: &
       urban_dyn_type
    use scale_landuse, only: &
       exists_urban => landuse_exists_urban
    implicit none
 
    integer :: k, i, j
    !---------------------------------------------------------------------------
 
    call prof_rapstart('URB_Update', 1)
 
    !########## Get Surface Boundary from coupler ##########
    call urban_surface_get
 
    !########## Dynamics / Update variables ##########
    select case ( urban_dyn_type )
    case ( 'KUSAKA01' )
 
!OCL XFILL
       !$omp parallel do
       do j = ujs, uje
       do i = uis, uie
       if ( exists_urban(i,j) ) then
          do k = uks, uke
             urban_trl(k,i,j) = urban_trl(k,i,j) + urban_trl_t(k,i,j) * dt
             urban_tbl(k,i,j) = urban_tbl(k,i,j) + urban_tbl_t(k,i,j) * dt
             urban_tgl(k,i,j) = urban_tgl(k,i,j) + urban_tgl_t(k,i,j) * dt
          end do
       end if
       end do
       end do
 
!OCL XFILL
       !$omp parallel do
       do j = ujs, uje
       do i = uis, uie
       if ( exists_urban(i,j) ) then
          urban_tr(i,j) = urban_tr(i,j) + urban_tr_t(i,j) * dt
          urban_tb(i,j) = urban_tb(i,j) + urban_tb_t(i,j) * dt
          urban_tg(i,j) = urban_tg(i,j) + urban_tg_t(i,j) * dt
          urban_tc(i,j) = urban_tc(i,j) + urban_tc_t(i,j) * dt
          urban_qc(i,j) = max( real(urban_qc(i,j) + urban_qc_t(i,j) * dt, rp), 0.0_rp )
          urban_uc(i,j) = max( real(urban_uc(i,j) + urban_uc_t(i,j) * dt, rp), 0.0_rp )
          urban_rainr(i,j) = max( real(urban_rainr(i,j) + urban_rainr_t(i,j) * dt, rp), 0.0_rp )
          urban_rainb(i,j) = max( real(urban_rainb(i,j) + urban_rainb_t(i,j) * dt, rp), 0.0_rp )
          urban_raing(i,j) = max( real(urban_raing(i,j) + urban_raing_t(i,j) * dt, rp), 0.0_rp )
       end if
       end do
       end do
 
    end select
 
    call urban_vars_check
 
    call prof_rapend  ('URB_Update', 1)
 
    return

References scale_tracer::k, scale_landuse::landuse_exists_urban, scale_prof::prof_rapend(), scale_prof::prof_rapstart(), scale_time::time_dtsec_urban, scale_urban_grid_cartesc_index::uie, scale_urban_grid_cartesc_index::uis, scale_urban_grid_cartesc_index::uje, scale_urban_grid_cartesc_index::ujs, scale_urban_grid_cartesc_index::uke, scale_urban_grid_cartesc_index::uks, mod_urban_admin::urban_dyn_type, mod_urban_vars::urban_qc, mod_urban_vars::urban_qc_t, mod_urban_vars::urban_rainb, mod_urban_vars::urban_rainb_t, mod_urban_vars::urban_raing, mod_urban_vars::urban_raing_t, mod_urban_vars::urban_rainr, mod_urban_vars::urban_rainr_t, urban_surface_get(), mod_urban_vars::urban_tb, mod_urban_vars::urban_tb_t, mod_urban_vars::urban_tbl, mod_urban_vars::urban_tbl_t, mod_urban_vars::urban_tc, mod_urban_vars::urban_tc_t, mod_urban_vars::urban_tg, mod_urban_vars::urban_tg_t, mod_urban_vars::urban_tgl, mod_urban_vars::urban_tgl_t, mod_urban_vars::urban_tr, mod_urban_vars::urban_tr_t, mod_urban_vars::urban_trl, mod_urban_vars::urban_trl_t, mod_urban_vars::urban_uc, mod_urban_vars::urban_uc_t, and mod_urban_vars::urban_vars_check().

Referenced by mod_rm_driver::rm_driver().

Here is the call graph for this function:

Here is the caller graph for this function:

urban_surface_get()

subroutine, public mod_urban_driver::urban_surface_get

Get surface boundary.

Definition at line 619 of file mod_urban_driver.F90.

    use mod_urban_admin, only: &
       urban_do
    use mod_urban_vars, only: &
       atmos_temp,      &
       atmos_pres,      &
       atmos_w,         &
       atmos_u,         &
       atmos_v,         &
       atmos_dens,      &
       atmos_qv,        &
       atmos_pbl,       &
       atmos_sfc_dens,  &
       atmos_sfc_pres,  &
       atmos_sflx_lw,   &
       atmos_sflx_sw,   &
       atmos_cossza,    &
       atmos_sflx_water, &
       atmos_sflx_engi
    use mod_cpl_vars, only: &
       cpl_getatm_urb
    implicit none
 
    real(RP) :: ATMOS_SFLX_rad_dn(UIA,UJA,N_RAD_DIR,N_RAD_RGN)
 
    integer  :: i, j
    !---------------------------------------------------------------------------
 
    call prof_rapstart('URB_SfcExch', 2)
 
    if ( urban_do ) then
       call cpl_getatm_urb( atmos_temp(:,:),     & ! [OUT]
                            atmos_pres(:,:),     & ! [OUT]
                            atmos_w(:,:),     & ! [OUT]
                            atmos_u(:,:),     & ! [OUT]
                            atmos_v(:,:),     & ! [OUT]
                            atmos_dens(:,:),     & ! [OUT]
                            atmos_qv(:,:),     & ! [OUT]
                            atmos_pbl(:,:),     & ! [OUT]
                            atmos_sfc_dens(:,:),     & ! [OUT]
                            atmos_sfc_pres(:,:),     & ! [OUT]
                            atmos_sflx_rad_dn(:,:,:,:), & ! [OUT]
                            atmos_cossza(:,:),     & ! [OUT]
                            atmos_sflx_water(:,:),     & ! [OUT]
                            atmos_sflx_engi(:,:)      ) ! [OUT]
    endif
 
!OCL XFILL
    do j = ujs, uje
    do i = uis, uie
       atmos_sflx_lw(i,j,i_r_direct ) = atmos_sflx_rad_dn(i,j,i_r_direct ,i_r_ir)    ! IR, direct
       atmos_sflx_lw(i,j,i_r_diffuse) = atmos_sflx_rad_dn(i,j,i_r_diffuse,i_r_ir)    ! IR, diffuse
 
       atmos_sflx_sw(i,j,i_r_direct ) = atmos_sflx_rad_dn(i,j,i_r_direct ,i_r_nir) & ! NIR, direct
                                      + atmos_sflx_rad_dn(i,j,i_r_direct ,i_r_vis)   ! VIS, direct
       atmos_sflx_sw(i,j,i_r_diffuse) = atmos_sflx_rad_dn(i,j,i_r_diffuse,i_r_nir) & ! NIR, diffuse
                                      + atmos_sflx_rad_dn(i,j,i_r_diffuse,i_r_vis)   ! VIS, diffuse
    enddo
    enddo
 
    call prof_rapend  ('URB_SfcExch', 2)
 
    return

References mod_urban_vars::atmos_cossza, mod_urban_vars::atmos_dens, mod_urban_vars::atmos_pbl, mod_urban_vars::atmos_pres, mod_urban_vars::atmos_qv, mod_urban_vars::atmos_sfc_dens, mod_urban_vars::atmos_sfc_pres, mod_urban_vars::atmos_sflx_engi, mod_urban_vars::atmos_sflx_lw, mod_urban_vars::atmos_sflx_sw, mod_urban_vars::atmos_sflx_water, mod_urban_vars::atmos_temp, mod_urban_vars::atmos_u, mod_urban_vars::atmos_v, mod_urban_vars::atmos_w, mod_cpl_vars::cpl_getatm_urb(), scale_cpl_sfc_index::i_r_diffuse, scale_cpl_sfc_index::i_r_direct, scale_cpl_sfc_index::i_r_ir, scale_cpl_sfc_index::i_r_nir, scale_cpl_sfc_index::i_r_vis, scale_prof::prof_rapend(), scale_prof::prof_rapstart(), scale_urban_grid_cartesc_index::uie, scale_urban_grid_cartesc_index::uis, scale_urban_grid_cartesc_index::uje, scale_urban_grid_cartesc_index::ujs, and mod_urban_admin::urban_do.

Referenced by urban_driver_calc_tendency(), and urban_driver_update().

Here is the call graph for this function:

Here is the caller graph for this function:

urban_surface_set()

subroutine, public mod_urban_driver::urban_surface_set

(

logical, intent(in)

countup

)

Set surface boundary to other model.

Definition at line 686 of file mod_urban_driver.F90.

    use mod_urban_admin, only: &
       urban_do
    use mod_urban_vars, only: &
       urban_sfc_temp,   &
       urban_sfc_albedo, &
       urban_sflx_mw,    &
       urban_sflx_mu,    &
       urban_sflx_mv,    &
       urban_sflx_sh,    &
       urban_sflx_lh,    &
       urban_sflx_shex,  &
       urban_sflx_qvex,  &
       urban_sflx_gh,    &
       urban_sflx_qtrc,  &
       urban_z0m,        &
       urban_z0h,        &
       urban_z0e,        &
       urban_u10,        &
       urban_v10,        &
       urban_t2,         &
       urban_q2
    use mod_cpl_vars, only: &
       cpl_puturb
    use scale_landuse, only: &
       exists_urban => landuse_exists_urban
    implicit none
 
    ! arguments
    logical, intent(in) :: countup
    !---------------------------------------------------------------------------
 
    call prof_rapstart('URB_SfcExch', 2)
 
    if ( urban_do ) then
       call cpl_puturb( urban_sfc_temp(:,:),     & ! [IN]
                        urban_sfc_albedo(:,:,:,:), & ! [IN]
                        urban_z0m(:,:),     & ! [IN]
                        urban_z0h(:,:),     & ! [IN]
                        urban_z0e(:,:),     & ! [IN]
                        urban_sflx_mw(:,:),     & ! [IN]
                        urban_sflx_mu(:,:),     & ! [IN]
                        urban_sflx_mv(:,:),     & ! [IN]
                        urban_sflx_sh(:,:),     & ! [IN]
                        urban_sflx_lh(:,:),     & ! [IN]
                        urban_sflx_shex(:,:),     & ! [IN]
                        urban_sflx_qvex(:,:),     & ! [IN]
                        urban_sflx_gh(:,:),     & ! [IN]
                        urban_sflx_qtrc(:,:,:),   & ! [IN]
                        urban_u10(:,:),     & ! [IN]
                        urban_v10(:,:),     & ! [IN]
                        urban_t2(:,:),     & ! [IN]
                        urban_q2(:,:),     & ! [IN]
                        exists_urban(:,:),     & ! [IN]
                        countup                    ) ! [IN]
    endif
 
    call prof_rapend  ('URB_SfcExch', 2)
 
    return

References mod_cpl_vars::cpl_puturb(), scale_landuse::landuse_exists_urban, scale_prof::prof_rapend(), scale_prof::prof_rapstart(), mod_urban_admin::urban_do, mod_urban_vars::urban_q2, mod_urban_vars::urban_sfc_albedo, mod_urban_vars::urban_sfc_temp, mod_urban_vars::urban_sflx_gh, mod_urban_vars::urban_sflx_lh, mod_urban_vars::urban_sflx_mu, mod_urban_vars::urban_sflx_mv, mod_urban_vars::urban_sflx_mw, mod_urban_vars::urban_sflx_qtrc, mod_urban_vars::urban_sflx_qvex, mod_urban_vars::urban_sflx_sh, mod_urban_vars::urban_sflx_shex, mod_urban_vars::urban_t2, mod_urban_vars::urban_u10, mod_urban_vars::urban_v10, mod_urban_vars::urban_z0e, mod_urban_vars::urban_z0h, and mod_urban_vars::urban_z0m.

Referenced by mod_rm_driver::restart_read(), mod_rm_prep::rm_prep(), and urban_driver_calc_tendency().

Here is the call graph for this function:

Here is the caller graph for this function: