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	URBAN_QC_t
URBAN_RAINB_t	tendency of URBAN_RAINB	K	URBAN_RAINB_t
URBAN_RAING_t	tendency of URBAN_RAING	K	URBAN_RAING_t
URBAN_RAINR_t	tendency of URBAN_RAINR	K	URBAN_RAINR_t
URBAN_ROFF_t	tendency of URBAN_ROFF	K	URBAN_ROFF_t
URBAN_TBL_t	tendency of URBAN_TBL	K	URBAN_TBL_t
URBAN_TB_t	tendency of URBAN_TB	K	URBAN_TB_t
URBAN_TC_t	tendency of URBAN_TC	K	URBAN_TC_t
URBAN_TGL_t	tendency of URBAN_TGL	K	URBAN_TGL_t
URBAN_TG_t	tendency of URBAN_TG	K	URBAN_TG_t
URBAN_TRL_t	tendency of URBAN_TRL	K	URBAN_TRL_t
URBAN_TR_t	tendency of URBAN_TR	K	URBAN_TR_t
URBAN_UC_t	tendency of URBAN_UC	m/s	URBAN_UC_t

Function/Subroutine Documentation

urban_driver_setup()

subroutine, public mod_urban_driver::urban_driver_setup

(

)

Setup.

Definition at line 52 of file mod_urban_driver.F90.

References 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, and mod_urban_vars::urban_z0m.

Referenced by mod_rm_driver::rm_driver().

    use scale_prc, only: &
       prc_abort
    use mod_urban_admin, only: &
       urban_do, &
       urban_dyn_type, &
       urban_sfc_type
    use scale_urban_dyn_kusaka01, only: &
       urban_dyn_kusaka01_setup
    use mod_urban_vars, only: &
       urban_z0m, &
       urban_z0h, &
       urban_z0e
    implicit none
    !---------------------------------------------------------------------------

    log_newline
    log_info("URBAN_driver_setup",*) 'Setup'

    if ( urban_do ) then

       select case ( urban_dyn_type )
       case ( 'KUSAKA01' )
          call urban_dyn_kusaka01_setup( uia, uis, uie, uja, ujs, uje, &
                                         urban_z0m(:,:), urban_z0h(:,:), urban_z0e(:,:) ) ! [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

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 99 of file mod_urban_driver.F90.

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_lw, mod_urban_vars::atmos_sflx_rain, mod_urban_vars::atmos_sflx_snow, mod_urban_vars::atmos_sflx_sw, mod_atmos_admin::atmos_sw_phy_ch, mod_urban_vars::atmos_temp, mod_urban_vars::atmos_u, mod_urban_vars::atmos_v, mod_urban_vars::atmos_w, scale_const::const_d2r, scale_atmos_hydrometeor::i_qv, 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_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, 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_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_roff, mod_urban_vars::urban_roff_t, 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_sh, 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_u10, mod_urban_vars::urban_uc, mod_urban_vars::urban_uc_t, 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(), and mod_rm_driver::rm_driver().

    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_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_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_rain, &
       atmos_sflx_snow, &
       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_roff_t,    &
       urban_sfc_temp,    &
       urban_sfc_albedo,  &
       urban_sflx_mw,     &
       urban_sflx_mu,     &
       urban_sflx_mv,     &
       urban_sflx_sh,     &
       urban_sflx_lh,     &
       urban_sflx_qtrc,   &
       urban_sflx_gh,     &
       urban_z0m,         &
       urban_z0h,         &
       urban_z0e,         &
       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
    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) :: 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
       urban_roff_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)
          roff(i,j) = urban_roff(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

       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_w(:,:), 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_rain(:,:), atmos_sflx_snow(:,:),                  & ! [IN]
                                cdz(:),                                                      & ! [IN]
                                landuse_fact_urban(:,:),                                     & ! [IN]
                                tloc, dsec, dt,                                              & ! [IN]
                                trl(:,:,:), tbl(:,:,:), tgl(:,:,:),                          & ! [INOUT]
                                tr(:,:), tb(:,:), tg(:,:), tc(:,:), qc(:,:), uc(:,:),        & ! [INOUT]
                                rainr(:,:), rainb(:,:), raing(:,:), 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_z0m(:,:), urban_z0h(:,:), urban_z0e(:,:),              & ! [OUT]
                                urban_u10(:,:), urban_v10(:,:), urban_t2(:,:), urban_q2(:,:) ) ! [OUT]

!OCL XFILL
       !$omp parallel do
       do j = ujs, uje
       do i = uis, uie
       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 do
       end do

!OCL XFILL
       !$omp parallel do
       do j = ujs, uje
       do i = uis, uie
          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
          urban_roff_t(i,j) = ( roff(i,j) - urban_roff(i,j) ) / dt
       end do
       end do

       if ( .NOT. atmos_hydrometeor_dry ) then
          !$omp parallel do
          do j = ujs, uje
          do i = uis, uie
             urban_sflx_qtrc(i,j,i_qv) = urban_sflx_lh(i,j) / lhv(i,j)
          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

    !########## Set Surface Boundary to coupler ##########
    call urban_surface_set( countup=.true. )


    call file_history_in( urban_tr_t(:,:), 'URBAN_TR_t', 'tendency of URBAN_TR', 'K',     dim_type='XY' )
    call file_history_in( urban_tb_t(:,:), 'URBAN_TB_t', 'tendency of URBAN_TB', 'K',     dim_type='XY' )
    call file_history_in( urban_tg_t(:,:), 'URBAN_TG_t', 'tendency of URBAN_TG', 'K',     dim_type='XY' )
    call file_history_in( urban_tc_t(:,:), 'URBAN_TC_t', 'tendency of URBAN_TC', 'K',     dim_type='XY' )
    call file_history_in( urban_qc_t(:,:), 'URBAN_QC_t', 'tendency of URBAN_QC', 'kg/kg', dim_type='XY' )
    call file_history_in( urban_uc_t(:,:), 'URBAN_UC_t', 'tendency of URBAN_UC', 'm/s',   dim_type='XY' )

    call file_history_in( urban_trl_t(:,:,:), 'URBAN_TRL_t', 'tendency of URBAN_TRL', 'K', dim_type='UXY' )
    call file_history_in( urban_tbl_t(:,:,:), 'URBAN_TBL_t', 'tendency of URBAN_TBL', 'K', dim_type='UXY' )
    call file_history_in( urban_tgl_t(:,:,:), 'URBAN_TGL_t', 'tendency of URBAN_TGL', 'K', dim_type='UXY' )

    call file_history_in( urban_rainr_t(:,:), 'URBAN_RAINR_t', 'tendency of URBAN_RAINR', 'K', dim_type='XY' )
    call file_history_in( urban_rainb_t(:,:), 'URBAN_RAINB_t', 'tendency of URBAN_RAINB', 'K', dim_type='XY' )
    call file_history_in( urban_raing_t(:,:), 'URBAN_RAING_t', 'tendency of URBAN_RAING', 'K', dim_type='XY' )
    call file_history_in( urban_roff_t(:,:), 'URBAN_ROFF_t',  'tendency of URBAN_ROFF',  'K', 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_t(:,:), 'URBAN_ROFF_t',  &
                              urban_grid_cartesc_real_area(:,:),   &
                              urban_grid_cartesc_real_totarea      )
    endif


    call prof_rapend  ('URB_CalcTend', 1)

    return

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 452 of file mod_urban_driver.F90.

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_lw, mod_urban_vars::atmos_sflx_rain, mod_urban_vars::atmos_sflx_snow, mod_urban_vars::atmos_sflx_sw, mod_urban_vars::atmos_temp, mod_urban_vars::atmos_u, mod_urban_vars::atmos_v, mod_urban_vars::atmos_w, scale_landuse::landuse_fact_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_q2, 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, mod_urban_vars::urban_roff, mod_urban_vars::urban_roff_t, 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_sh, urban_surface_get(), 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_u10, mod_urban_vars::urban_uc, mod_urban_vars::urban_uc_t, mod_urban_vars::urban_v10, mod_urban_vars::urban_vars_history(), mod_urban_vars::urban_vars_total(), mod_urban_vars::urban_z0e, mod_urban_vars::urban_z0h, and mod_urban_vars::urban_z0m.

Referenced by mod_rm_driver::rm_driver().

    use scale_time, only: &
       dt => time_dtsec_urban
    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_rain, &
       atmos_sflx_snow, &
       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_roff_t,      &
       urban_sflx_mw,     &
       urban_sflx_mu,     &
       urban_sflx_mv,     &
       urban_sflx_sh,     &
       urban_sflx_lh,     &
       urban_sflx_qtrc,   &
       urban_sflx_gh,     &
       urban_z0m,         &
       urban_z0h,         &
       urban_z0e,         &
       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,        &
       urban_sfc_temp,    &
       urban_vars_total,  &
       urban_vars_history
    use scale_landuse, only: &
       landuse_fact_urban
    use mod_urban_admin, only: &
       urban_dyn_type
    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
       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 do
       end do

!OCL XFILL
       !$omp parallel do
       do j = ujs, uje
       do i = uis, uie
          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) = urban_qc(i,j) + urban_qc_t(i,j) * dt
          urban_uc(i,j) = urban_uc(i,j) + urban_uc_t(i,j) * dt
          urban_rainr(i,j) = urban_rainr(i,j) + urban_rainr_t(i,j) * dt
          urban_rainb(i,j) = urban_rainb(i,j) + urban_rainb_t(i,j) * dt
          urban_raing(i,j) = urban_raing(i,j) + urban_raing_t(i,j) * dt
          urban_roff(i,j) = urban_roff(i,j) + urban_roff_t(i,j) * dt
       end do
       end do

    end select

    call urban_vars_total

    !########## History & Monitor ##########
    call urban_vars_history


    call prof_rapend  ('URB_Update', 1)

    return

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 576 of file mod_urban_driver.F90.

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_lw, mod_urban_vars::atmos_sflx_rain, mod_urban_vars::atmos_sflx_snow, mod_urban_vars::atmos_sflx_sw, 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().

    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_rain, &
       atmos_sflx_snow
    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_rain(:,:),     & ! [OUT]
                            atmos_sflx_snow(:,:)      ) ! [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

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 643 of file mod_urban_driver.F90.

References mod_cpl_vars::cpl_puturb(), 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_sh, 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().

    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_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
    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_gh(:,:),     & ! [IN]
                        urban_sflx_qtrc(:,:,:),   & ! [IN]
                        urban_u10(:,:),     & ! [IN]
                        urban_v10(:,:),     & ! [IN]
                        urban_t2(:,:),     & ! [IN]
                        urban_q2(:,:),     & ! [IN]
                        countup                    ) ! [IN]
    endif

    call prof_rapend  ('URB_SfcExch', 2)

    return

Here is the call graph for this function:

Here is the caller graph for this function: