mod_urban_phy_driver Module Reference

module URBAN / Physics Urban Canopy Model (UCM) More...

Functions/Subroutines
subroutine, public	urban_phy_driver_setup
	Setup. More...
subroutine, public	urban_phy_driver_resume
	Resume. More...
subroutine, public	urban_phy_driver (update_flag)
	Driver. More...

Detailed Description

module URBAN / Physics Urban Canopy Model (UCM)

Description

Urban physics module

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_phy_driver_setup()

subroutine, public mod_urban_phy_driver::urban_phy_driver_setup

(

)

Setup.

Definition at line 52 of file mod_urban_phy_driver.f90.

References scale_stdio::io_fid_log, scale_stdio::io_l, scale_urban_phy::urban_phy_setup(), mod_urban_admin::urban_sw, mod_urban_admin::urban_type, mod_urban_vars::urban_z0e, mod_urban_vars::urban_z0h, and mod_urban_vars::urban_z0m.

Referenced by mod_urban_driver::urban_driver_setup().

    use scale_urban_phy, only: &
       urban_phy_setup
    use mod_urban_admin, only: &
       urban_type, &
       urban_sw
    use mod_urban_vars, only: &
       urban_z0m, &
       urban_z0h, &
       urban_z0e
    implicit none
    !---------------------------------------------------------------------------

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

    if ( urban_sw ) then

       ! setup library component
       call urban_phy_setup( urban_type,     & ! [IN]
                             urban_z0m(:,:), & ! [OUT]
                             urban_z0h(:,:), & ! [OUT]
                             urban_z0e(:,:)  ) ! [OUT]

    else
       if( io_l ) write(io_fid_log,*) '*** this component is never called.'
    endif

    return

Here is the call graph for this function:

Here is the caller graph for this function:

urban_phy_driver_resume()

subroutine, public mod_urban_phy_driver::urban_phy_driver_resume

(

)

Resume.

Definition at line 85 of file mod_urban_phy_driver.f90.

References scale_prof::prof_rapend(), scale_prof::prof_rapstart(), mod_admin_restart::restart_run, urban_phy_driver(), and mod_urban_admin::urban_sw.

Referenced by mod_urban_driver::urban_driver_resume().

    use mod_admin_restart, only: &
       restart_run
    use mod_urban_admin, only: &
       urban_sw
    implicit none

    if ( urban_sw ) then

       if ( .NOT. restart_run ) then ! tentative
          ! run once (only for the diagnostic value)
          call prof_rapstart('URB_Physics', 1)
          call urban_phy_driver( update_flag = .true. )
          call prof_rapend  ('URB_Physics', 1)
       end if

    end if

    return

Here is the call graph for this function:

Here is the caller graph for this function:

urban_phy_driver()

subroutine, public mod_urban_phy_driver::urban_phy_driver

(

logical, intent(in)

update_flag

)

Driver.

Definition at line 107 of file mod_urban_phy_driver.f90.

References 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, scale_grid_index::ie, scale_grid_index::is, scale_grid_index::je, scale_grid_index::js, scale_mapproj::mprj_basepoint_lat, scale_mapproj::mprj_basepoint_lon, scale_grid_real::real_z1, scale_rm_statistics::statistics_checktotal, scale_time::time_dtsec_urban, scale_time::time_nowdate, scale_urban_grid_index::uke, scale_urban_grid_index::uks, scale_urban_phy::urban_phy, 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_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_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_urban_driver::urban_driver(), and urban_phy_driver_resume().

    use scale_atmos_hydrometeor, only: &
       hydrometeor_lhv => atmos_hydrometeor_lhv
    use scale_time, only: &
       nowdate => time_nowdate,     &
       dt      => time_dtsec_urban
    use scale_rm_statistics, only: &
       statistics_checktotal, &
       stat_total
    use scale_history, only: &
       hist_in
    use scale_mapproj, only: &
       base_lon => mprj_basepoint_lon, &
       base_lat => mprj_basepoint_lat
    use scale_grid_real, only: &
       real_z1
    use scale_urban_phy, only: &
       urban_phy
    use mod_urban_vars, only: &
       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_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_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,         &
       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_sflx_prec
    implicit none

    ! arguments
    logical, intent(in) :: update_flag

    ! works
    real(RP) :: LHV(IA,JA) ! latent heat of vaporization [J/kg]
    real(RP) :: total      ! dummy

    character(len=2) :: sk

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

    if ( update_flag ) then

       call urban_phy( urban_tr_t(:,:),      & ! [OUT]
                       urban_tb_t(:,:),      & ! [OUT]
                       urban_tg_t(:,:),      & ! [OUT]
                       urban_tc_t(:,:),      & ! [OUT]
                       urban_qc_t(:,:),      & ! [OUT]
                       urban_uc_t(:,:),      & ! [OUT]
                       urban_trl_t(:,:,:),    & ! [OUT]
                       urban_tbl_t(:,:,:),    & ! [OUT]
                       urban_tgl_t(:,:,:),    & ! [OUT]
                       urban_rainr_t(:,:),      & ! [OUT]
                       urban_rainb_t(:,:),      & ! [OUT]
                       urban_raing_t(:,:),      & ! [OUT]
                       urban_roff_t(:,:),      & ! [OUT]
                       urban_sfc_temp(:,:),      & ! [OUT]
                       urban_sfc_albedo(:,:,i_lw), & ! [OUT]
                       urban_sfc_albedo(:,:,i_sw), & ! [OUT]
                       urban_sflx_mw(:,:),      & ! [OUT]
                       urban_sflx_mu(:,:),      & ! [OUT]
                       urban_sflx_mv(:,:),      & ! [OUT]
                       urban_sflx_sh(:,:),      & ! [OUT]
                       urban_sflx_lh(:,:),      & ! [OUT]
                       urban_sflx_gh(:,:),      & ! [OUT]
                       urban_z0m(:,:),      & ! [OUT]
                       urban_z0h(:,:),      & ! [OUT]
                       urban_z0e(:,:),      & ! [OUT]
                       urban_u10(:,:),      & ! [OUT]
                       urban_v10(:,:),      & ! [OUT]
                       urban_t2(:,:),      & ! [OUT]
                       urban_q2(:,:),      & ! [OUT]
                       atmos_temp(:,:),      & ! [IN]
                       atmos_pres(:,:),      & ! [IN]
                       atmos_w(:,:),      & ! [IN]
                       atmos_u(:,:),      & ! [IN]
                       atmos_v(:,:),      & ! [IN]
                       atmos_dens(:,:),      & ! [IN]
                       atmos_qv(:,:),      & ! [IN]
                       real_z1(:,:),      & ! [IN]
                       atmos_pbl(:,:),      & ! [IN]
                       atmos_sfc_pres(:,:),      & ! [IN]
                       atmos_sflx_lw(:,:,:),    & ! [IN]
                       atmos_sflx_sw(:,:,:),    & ! [IN]
                       atmos_sflx_prec(:,:),      & ! [IN]
                       urban_tr(:,:),      & ! [IN]
                       urban_tb(:,:),      & ! [IN]
                       urban_tg(:,:),      & ! [IN]
                       urban_tc(:,:),      & ! [IN]
                       urban_qc(:,:),      & ! [IN]
                       urban_uc(:,:),      & ! [IN]
                       urban_trl(:,:,:),    & ! [IN]
                       urban_tbl(:,:,:),    & ! [IN]
                       urban_tgl(:,:,:),    & ! [IN]
                       urban_rainr(:,:),      & ! [IN]
                       urban_rainb(:,:),      & ! [IN]
                       urban_raing(:,:),      & ! [IN]
                       urban_roff(:,:),      & ! [IN]
                       base_lon,                   & ! [IN]
                       base_lat,                   & ! [IN]
                       nowdate(:),        & ! [IN]
                       dt                          ) ! [IN]

       call hydrometeor_lhv( lhv(:,:), atmos_temp(:,:) )

!OCL XFILL
       do j = js, je
       do i = is, ie
          urban_sflx_evap(i,j) = urban_sflx_lh(i,j) / lhv(i,j)
       end do
       end do

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

       call hist_in( urban_trl_t(:,:,:), 'URBAN_TRL_t', 'tendency of URBAN_TRL', 'K', zdim='urban' )
       call hist_in( urban_tbl_t(:,:,:), 'URBAN_TBL_t', 'tendency of URBAN_TBL', 'K', zdim='urban' )
       call hist_in( urban_tgl_t(:,:,:), 'URBAN_TGL_t', 'tendency of URBAN_TGL', 'K', zdim='urban' )

       call hist_in( urban_rainr_t(:,:), 'URBAN_RAINR_t', 'tendency of URBAN_RAINR', 'K' )
       call hist_in( urban_rainb_t(:,:), 'URBAN_RAINB_t', 'tendency of URBAN_RAINB', 'K' )
       call hist_in( urban_raing_t(:,:), 'URBAN_RAING_t', 'tendency of URBAN_RAING', 'K' )
       call hist_in( urban_roff_t(:,:), 'URBAN_ROFF_t',  'tendency of URBAN_ROFF',  'K' )

    endif

    if ( statistics_checktotal ) then
       call stat_total( total, urban_tr_t(:,:), 'URBAN_TR_t' )
       call stat_total( total, urban_tb_t(:,:), 'URBAN_TB_t' )
       call stat_total( total, urban_tg_t(:,:), 'URBAN_TG_t' )
       call stat_total( total, urban_tc_t(:,:), 'URBAN_TC_t' )
       call stat_total( total, urban_qc_t(:,:), 'URBAN_QC_t' )
       call stat_total( total, urban_uc_t(:,:), 'URBAN_UC_t' )

       do k = uks, uke
          write(sk,'(I2.2)') k

          call stat_total( total, urban_trl_t(k,:,:), 'URBAN_TRL_t'//sk  )
          call stat_total( total, urban_tbl_t(k,:,:), 'URBAN_TBL_t'//sk  )
          call stat_total( total, urban_tgl_t(k,:,:), 'URBAN_TGL_t'//sk  )
       enddo

       call stat_total( total, urban_rainr_t(:,:), 'URBAN_RAINR_t' )
       call stat_total( total, urban_rainb_t(:,:), 'URBAN_RAINB_t' )
       call stat_total( total, urban_raing_t(:,:), 'URBAN_RAING_t' )
       call stat_total( total, urban_roff_t(:,:), 'URBAN_ROFF_t'  )
    endif

    return

Here is the caller graph for this function: