mod_urban_driver Module Reference

module URBAN driver More...

Functions/Subroutines
subroutine, public	urban_driver_setup
	Setup. More...
subroutine, public	urban_driver_resume
	Resume. More...
subroutine, public	urban_driver
	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

Function/Subroutine Documentation

urban_driver_setup()

subroutine, public mod_urban_driver::urban_driver_setup

(

)

Setup.

Definition at line 55 of file mod_urban_driver.f90.

References scale_stdio::io_fid_log, scale_stdio::io_l, and mod_urban_phy_driver::urban_phy_driver_setup().

Referenced by mod_rm_driver::scalerm().

    use mod_urban_phy_driver, only: &
       urban_phy_driver_setup
    implicit none
    !---------------------------------------------------------------------------

    if( io_l ) write(io_fid_log,*)
    if( io_l ) write(io_fid_log,*) '++++++ Module[DRIVER] / Categ[URBAN] / Origin[SCALE-RM]'

    call urban_phy_driver_setup

    return

Here is the call graph for this function:

Here is the caller graph for this function:

urban_driver_resume()

subroutine, public mod_urban_driver::urban_driver_resume

(

)

Resume.

Definition at line 71 of file mod_urban_driver.f90.

References scale_stdio::io_fid_log, scale_stdio::io_l, scale_prof::prof_rapend(), scale_prof::prof_rapstart(), mod_urban_phy_driver::urban_phy_driver_resume(), urban_surface_get(), urban_surface_set(), mod_urban_admin::urban_sw, and mod_urban_vars::urban_vars_history().

Referenced by mod_rm_driver::resume_state().

    use mod_urban_phy_driver, only: &
       urban_phy_driver_resume
    use mod_urban_vars, only: &
       urban_vars_history
    use mod_urban_admin, only: &
       urban_sw
    implicit none
    !---------------------------------------------------------------------------

    if( io_l ) write(io_fid_log,*)
    if( io_l ) write(io_fid_log,*) '++++++ Module[DRIVER] / Categ[URBAN] / Origin[SCALE-RM]'

    !########## Get Surface Boundary from coupler ##########
    call urban_surface_get

    call urban_phy_driver_resume

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

    !########## History & Monitor ##########
    if ( urban_sw ) then
       call prof_rapstart('URB_History', 1)
       call urban_vars_history
       call prof_rapend  ('URB_History', 1)
    endif

    return

Here is the call graph for this function:

Here is the caller graph for this function:

urban_driver()

subroutine, public mod_urban_driver::urban_driver

(

)

Urban step.

Definition at line 104 of file mod_urban_driver.f90.

References scale_grid_index::ie, scale_grid_index::is, scale_grid_index::je, scale_grid_index::js, scale_prof::prof_rapend(), scale_prof::prof_rapstart(), scale_time::time_dtsec_urban, scale_urban_grid_index::uke, scale_urban_grid_index::uks, mod_urban_phy_driver::urban_phy_driver(), 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, urban_surface_get(), urban_surface_set(), mod_urban_admin::urban_sw, 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, mod_urban_vars::urban_vars_history(), and mod_urban_vars::urban_vars_total().

Referenced by mod_rm_driver::scalerm().

    use scale_time, only: &
       dt => time_dtsec_urban
    use mod_urban_admin, only: &
       urban_sw
    use mod_urban_vars, only: &
       urban_tr_t,        &
       urban_tb_t,        &
       urban_tg_t,        &
       urban_tc_t,        &
       urban_qc_t,        &
       urban_uc_t,        &
       urban_trl_t,       &
       urban_tbl_t,       &
       urban_tgl_t,       &
       urban_rainr_t,     &
       urban_rainb_t,     &
       urban_raing_t,     &
       urban_roff_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_roff,        &
       urban_vars_total,  &
       urban_vars_history
    use mod_urban_phy_driver, only: &
       urban_phy_driver
    implicit none

    integer :: k, i, j
    !---------------------------------------------------------------------------

    !########## Get Surface Boundary from coupler ##########
    call prof_rapstart('URB_SfcExch', 2)
    call urban_surface_get
    call prof_rapend  ('URB_SfcExch', 2)

    !########## Physics ##########
    if ( urban_sw ) then
       call prof_rapstart('URB_Physics', 1)
       call urban_phy_driver( update_flag = .true. )
       call prof_rapend  ('URB_Physics', 1)
    endif

    !########## Update ##########
    do j = js, je
    do i = is, ie
    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
    do j = js, je
    do i = is, ie
       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

    call urban_vars_total

    !########## Set Surface Boundary to coupler ##########
    call prof_rapstart('URB_SfcExch', 2)
    call urban_surface_set( countup=.true. )
    call prof_rapend  ('URB_SfcExch', 2)

    !########## reset tendencies ##########
!OCL XFILL
    do j = js, je
    do i = is, ie
    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 = js, je
    do i = is, ie
       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

    !########## History & Monitor ##########
    call prof_rapstart('URB_History', 1)
    call urban_vars_history
    call prof_rapend  ('URB_History', 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 229 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_pres, mod_urban_vars::atmos_sflx_lw, mod_urban_vars::atmos_sflx_prec, 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_grid_index::ie, scale_grid_index::is, scale_grid_index::je, scale_grid_index::js, and mod_urban_admin::urban_sw.

Referenced by urban_driver(), and urban_driver_resume().

    use mod_urban_admin, only: &
       urban_sw
    use mod_urban_vars, only: &
       atmos_temp,     &
       atmos_pres,     &
       atmos_w,        &
       atmos_u,        &
       atmos_v,        &
       atmos_dens,     &
       atmos_qv,       &
       atmos_pbl,      &
       atmos_sfc_pres, &
       atmos_sflx_lw,  &
       atmos_sflx_sw,  &
       atmos_cossza,   &
       atmos_sflx_prec
    use mod_cpl_vars, only: &
       cpl_getatm_urb
    implicit none

    real(RP) :: ATMOS_SFLX_rad_dn(IA,JA,2,2)
    real(RP) :: ATMOS_SFLX_rain  (IA,JA)
    real(RP) :: ATMOS_SFLX_snow  (IA,JA)

    integer  :: i, j
    !---------------------------------------------------------------------------

    if ( urban_sw ) 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_pres(:,:),     & ! [OUT]
                            atmos_sflx_rad_dn(:,:,:,:), & ! [OUT]
                            atmos_cossza(:,:),     & ! [OUT]
                            atmos_sflx_rain(:,:),     & ! [OUT]
                            atmos_sflx_snow(:,:)      ) ! [OUT]
    endif

!OCL XFILL
    do j = js, je
    do i = is, ie
       atmos_sflx_sw(i,j,1) = atmos_sflx_rad_dn(i,j,i_sw,1) ! direct
       atmos_sflx_lw(i,j,1) = atmos_sflx_rad_dn(i,j,i_lw,1) ! direct
       atmos_sflx_sw(i,j,2) = atmos_sflx_rad_dn(i,j,i_sw,2) ! diffuse
       atmos_sflx_lw(i,j,2) = atmos_sflx_rad_dn(i,j,i_lw,2) ! diffuse

       atmos_sflx_prec(i,j)   = atmos_sflx_rain(i,j) + atmos_sflx_snow(i,j) ! liquid+ice
    enddo
    enddo

    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 290 of file mod_urban_driver.f90.

References mod_cpl_vars::cpl_puturb(), mod_urban_vars::urban_q2, mod_urban_vars::urban_sfc_albedo, mod_urban_vars::urban_sfc_temp, mod_urban_vars::urban_sflx_evap, 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_sh, mod_urban_admin::urban_sw, 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_mkinit::mkinit(), mod_rm_driver::resume_state(), urban_driver(), and urban_driver_resume().

    use mod_urban_admin, only: &
       urban_sw
    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_evap,  &
       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
    !---------------------------------------------------------------------------

    if ( urban_sw ) 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_evap(:,:),   & ! [IN]
                        urban_u10(:,:),   & ! [IN]
                        urban_v10(:,:),   & ! [IN]
                        urban_t2(:,:),   & ! [IN]
                        urban_q2(:,:),   & ! [IN]
                        countup                  ) ! [IN]
    endif

    return

Here is the call graph for this function:

Here is the caller graph for this function: