mod_atmos_phy_sf_driver Module Reference

module ATMOSPHERE / Physics Surface fluxes More...

Functions/Subroutines
subroutine, public	atmos_phy_sf_driver_setup
	Setup. More...
subroutine, public	atmos_phy_sf_driver_calc_tendency (update_flag)
	calculation tendency More...
subroutine	history_output

Detailed Description

module ATMOSPHERE / Physics Surface fluxes

Description

Flux from/to bottom boundary of atmosphere (surface)

Author

Team SCALE

NAMELIST

• No namelist group

History Output

name	description	unit	variable
SFLX_{TRACER_NAME}	surface {TRACER_NAME} flux; {TRACER_NAME} depends on the physics schemes, e.g., QV, QC, QR.	kg/m2/s	SFLX_QTRC
GHFLX	ground heat flux (downward)	W/m2	SFLX_GH
LHFLX	latent heat flux	W/m2	SFLX_LH
MSLP	mean sea-level pressure	Pa	MSLP
MUFLX	u-momentum flux	kg/m/s2	SFLX_MU
MVFLX	v-momentum flux	kg/m/s2	SFLX_MV
MWFLX	w-momentum flux	kg/m/s2	SFLX_MW
Q2	2m specific humidity	kg/kg	Q2
Q2	2m specific humidity	kg/kg	Q2
Qstar	moisuter scale	kg/kg	Qstar
RH2	2m relative humidity	%	RH2
RLmo	inverse of Obukhov length	1/m	RLmo
SFC_ALB_IR_dif	surface albedo (IR; diffuse)	1	SFC_albedo
SFC_ALB_IR_dir	surface albedo (IR; direct	1	SFC_albedo
SFC_ALB_NIR_dif	surface albedo (NIR; diffuse	1	SFC_albedo
SFC_ALB_NIR_dir	surface albedo (NIR; direct	1	SFC_albedo
SFC_ALB_VIS_dif	surface albedo (VIS; diffuse	1	SFC_albedo
SFC_ALB_VIS_dir	surface albedo (VIS; direct	1	SFC_albedo
SFC_DENS	surface atmospheric density	kg/m3	SFC_DENS
SFC_PRES	surface atmospheric pressure	Pa	SFC_PRES
SFC_TEMP	surface skin temperature	K	SFC_TEMP
SFC_Z0E	roughness length (vapor)	m	SFC_Z0E
SFC_Z0H	roughness length (heat)	m	SFC_Z0H
SFC_Z0M	roughness length (momentum)	m	SFC_Z0M
SFLX_ENGI	ground internal energy flux (merged)	W/m2	SFLX_ENGI
SHFLX	sensible heat flux	W/m2	SFLX_SH
T2	2m air temperature	K	T2
T2	2m air temperature	K	T2
Tstar	temperature scale	K	Tstar
U10	10m x-wind	m/s	U10
U10m	10m eastward wind	m/s	U10m
Uabs10	10m absolute wind	m/s	Uabs10
Ustar	friction velocity	m/s	Ustar
V10	10m y-wind	m/s	V10
V10m	10m northward wind	m/s	V10m
Wstar	convective velocity scale	m/s	Wstar

Function/Subroutine Documentation

atmos_phy_sf_driver_setup()

subroutine, public mod_atmos_phy_sf_driver::atmos_phy_sf_driver_setup

Setup.

Definition at line 53 of file mod_atmos_phy_sf_driver.F90.

    use scale_prc, only: &
       prc_abort
    use scale_atmos_phy_sf_bulk, only: &
       atmos_phy_sf_bulk_setup
    use scale_atmos_phy_sf_const, only: &
       atmos_phy_sf_const_setup
    use scale_file_history, only: &
       file_history_reg
    use mod_atmos_admin, only: &
       atmos_phy_sf_type, &
       atmos_sw_phy_sf
    use mod_atmos_phy_sf_vars, only: &
       sfc_z0m   => atmos_phy_sf_sfc_z0m,   &
       sfc_z0h   => atmos_phy_sf_sfc_z0h,   &
       sfc_z0e   => atmos_phy_sf_sfc_z0e,   &
       sflx_mw   => atmos_phy_sf_sflx_mw,   &
       sflx_mu   => atmos_phy_sf_sflx_mu,   &
       sflx_mv   => atmos_phy_sf_sflx_mv,   &
       sflx_sh   => atmos_phy_sf_sflx_sh,   &
       sflx_lh   => atmos_phy_sf_sflx_lh,   &
       sflx_shex => atmos_phy_sf_sflx_shex, &
       sflx_qvex => atmos_phy_sf_sflx_qvex, &
       sflx_qtrc => atmos_phy_sf_sflx_qtrc, &
       sflx_engi => atmos_phy_sf_sflx_engi, &
       ustar     => atmos_phy_sf_ustar,     &
       tstar     => atmos_phy_sf_tstar,     &
       qstar     => atmos_phy_sf_qstar,     &
       wstar     => atmos_phy_sf_wstar,     &
       rlmo      => atmos_phy_sf_rlmo
    use mod_cpl_admin, only: &
       cpl_sw
    implicit none
    !---------------------------------------------------------------------------
 
    log_newline
    log_info("ATMOS_PHY_SF_driver_setup",*) 'Setup'
 
    if ( atmos_sw_phy_sf ) then
 
       if ( cpl_sw ) then
          log_info("ATMOS_PHY_SF_driver_setup",*) 'Coupler is enabled.'
       else
          ! setup library component
          select case( atmos_phy_sf_type )
          case ( 'BULK' )
             call atmos_phy_sf_bulk_setup
          case ( 'CONST' )
             call atmos_phy_sf_const_setup
          case default
             log_error("ATMOS_PHY_SF_driver_setup",*) 'invalid Surface flux type(', trim(atmos_phy_sf_type), '). CHECK!'
             call prc_abort
          end select
       endif
 
    else
 
       log_info("ATMOS_PHY_SF_driver_setup",*) 'this component is never called.'
       log_info("ATMOS_PHY_SF_driver_setup",*) 'surface fluxes are set to zero.'
       !$acc kernels
       sflx_mw(:,:) = 0.0_rp
       sflx_mu(:,:) = 0.0_rp
       sflx_mv(:,:) = 0.0_rp
       sflx_sh(:,:) = 0.0_rp
       sflx_lh(:,:) = 0.0_rp
       sflx_shex(:,:) = 0.0_rp
       sflx_qvex(:,:) = 0.0_rp
       ustar(:,:) = 0.0_rp
       tstar(:,:) = 0.0_rp
       qstar(:,:) = 0.0_rp
       wstar(:,:) = 0.0_rp
       rlmo(:,:) = 0.0_rp
       !$acc end kernels
       log_info("ATMOS_PHY_SF_driver_setup",*) 'SFC_TEMP, SFC_albedo is set in ATMOS_PHY_SF_vars.'
 
    endif
 
    !$acc kernels
    sflx_qtrc(:,:,:) = 0.0_rp
    sflx_engi(:,:)   = 0.0_rp
    !$acc end kernels
 
    call file_history_reg( 'Uabs10', '10m absolute wind',         'm/s'  , hist_uabs10, ndims=2, fill_halo=.true. )
    call file_history_reg( 'U10m',   '10m eastward wind',         'm/s'  , hist_u10m,   ndims=2, fill_halo=.true. )
    call file_history_reg( 'V10m',   '10m northward wind',        'm/s'  , hist_v10m,   ndims=2, fill_halo=.true. )
    call file_history_reg( 'T2',     '2m air temperature',        'K'    , hist_t2,     ndims=2, fill_halo=.true. )
    call file_history_reg( 'Q2',     '2m specific humidity',      'kg/kg', hist_q2,     ndims=2, fill_halo=.true. )
    call file_history_reg( 'RH2',    '2m relative humidity',      '%',     hist_rh2,    ndims=2, fill_halo=.true. )
    call file_history_reg( 'MSLP',   'mean sea-level pressure',   'Pa'   , hist_mslp,   ndims=2, fill_halo=.true., standard_name='air_pressure_at_mean_sea_level' )
 
    return

References scale_atmos_phy_sf_bulk::atmos_phy_sf_bulk_setup(), scale_atmos_phy_sf_const::atmos_phy_sf_const_setup(), mod_atmos_phy_sf_vars::atmos_phy_sf_qstar, mod_atmos_phy_sf_vars::atmos_phy_sf_rlmo, mod_atmos_phy_sf_vars::atmos_phy_sf_sfc_z0e, mod_atmos_phy_sf_vars::atmos_phy_sf_sfc_z0h, mod_atmos_phy_sf_vars::atmos_phy_sf_sfc_z0m, mod_atmos_phy_sf_vars::atmos_phy_sf_sflx_engi, mod_atmos_phy_sf_vars::atmos_phy_sf_sflx_lh, mod_atmos_phy_sf_vars::atmos_phy_sf_sflx_mu, mod_atmos_phy_sf_vars::atmos_phy_sf_sflx_mv, mod_atmos_phy_sf_vars::atmos_phy_sf_sflx_mw, mod_atmos_phy_sf_vars::atmos_phy_sf_sflx_qtrc, mod_atmos_phy_sf_vars::atmos_phy_sf_sflx_qvex, mod_atmos_phy_sf_vars::atmos_phy_sf_sflx_sh, mod_atmos_phy_sf_vars::atmos_phy_sf_sflx_shex, mod_atmos_phy_sf_vars::atmos_phy_sf_tstar, mod_atmos_admin::atmos_phy_sf_type, mod_atmos_phy_sf_vars::atmos_phy_sf_ustar, mod_atmos_phy_sf_vars::atmos_phy_sf_wstar, mod_atmos_admin::atmos_sw_phy_sf, mod_cpl_admin::cpl_sw, scale_file_history::file_history_reg(), and scale_prc::prc_abort().

Referenced by mod_atmos_driver::atmos_driver_setup().

Here is the call graph for this function:

Here is the caller graph for this function:

atmos_phy_sf_driver_calc_tendency()

subroutine, public mod_atmos_phy_sf_driver::atmos_phy_sf_driver_calc_tendency

(

logical, intent(in)

update_flag

)

calculation tendency

Definition at line 148 of file mod_atmos_phy_sf_driver.F90.

    use scale_const, only: &
       undef  => const_undef, &
       pre00  => const_pre00, &
       rdry   => const_rdry,  &
       rvap   => const_rvap,  &
       cpdry  => const_cpdry, &
       cpvap  => const_cpvap, &
       epstvap => const_epstvap
    use scale_atmos_grid_cartesc_real, only: &
       cz => atmos_grid_cartesc_real_cz, &
       fz => atmos_grid_cartesc_real_fz, &
       z1 => atmos_grid_cartesc_real_z1, &
       atmos_grid_cartesc_real_area, &
       atmos_grid_cartesc_real_totarea
    use scale_topography, only: &
       tansl_x => topography_tansl_x, &
       tansl_y => topography_tansl_y
    use scale_time, only: &
       dt_sf => time_dtsec_atmos_phy_sf
    use scale_statistics, only: &
       statistics_checktotal, &
       statistics_total
    use scale_atmos_bottom, only: &
       bottom_estimate => atmos_bottom_estimate
    use scale_atmos_hydrometeor, only: &
       atmos_hydrometeor_dry, &
       lhv, &
       i_qv
    use scale_atmos_phy_sf_bulk, only: &
       atmos_phy_sf_bulk_flux
    use scale_atmos_phy_sf_const, only: &
       atmos_phy_sf_const_flux
    use scale_bulkflux, only: &
       bulkflux_diagnose_scales
    use mod_atmos_admin, only: &
       atmos_phy_sf_type
    use mod_atmos_vars, only: &
       dens   => dens_av, &
       rhot   => rhot_av, &
       pott,              &
       temp,              &
       pres,              &
       w,                 &
       u,                 &
       v,                 &
       qv,                &
       cptot,             &
       cvtot,             &
       dens_t => dens_tp, &
       momz_t => momz_tp, &
       rhou_t => rhou_tp, &
       rhov_t => rhov_tp, &
       rhoh   => rhoh_p,  &
       rhoq_t => rhoq_tp
    use mod_atmos_phy_rd_vars, only: &
       sflx_lw_dn => atmos_phy_rd_sflx_lw_dn, &
       sflx_sw_dn => atmos_phy_rd_sflx_sw_dn
    use mod_atmos_phy_bl_vars, only: &
       pbl_zi => atmos_phy_bl_zi
    use mod_atmos_phy_sf_vars, only: &
       dens_t_sf => atmos_phy_sf_dens_t,    &
       momz_t_sf => atmos_phy_sf_momz_t,    &
       rhou_t_sf => atmos_phy_sf_rhou_t,    &
       rhov_t_sf => atmos_phy_sf_rhov_t,    &
       rhoh_sf   => atmos_phy_sf_rhoh,      &
       rhoq_t_sf => atmos_phy_sf_rhoq_t,    &
       sfc_dens  => atmos_phy_sf_sfc_dens,  &
       sfc_pres  => atmos_phy_sf_sfc_pres,  &
       sfc_temp  => atmos_phy_sf_sfc_temp,  &
       sfc_z0m   => atmos_phy_sf_sfc_z0m,   &
       sfc_z0h   => atmos_phy_sf_sfc_z0h,   &
       sfc_z0e   => atmos_phy_sf_sfc_z0e,   &
       sflx_mw   => atmos_phy_sf_sflx_mw,   &
       sflx_mu   => atmos_phy_sf_sflx_mu,   &
       sflx_mv   => atmos_phy_sf_sflx_mv,   &
       sflx_sh   => atmos_phy_sf_sflx_sh,   &
       sflx_lh   => atmos_phy_sf_sflx_lh,   &
       sflx_shex => atmos_phy_sf_sflx_shex, &
       sflx_qvex => atmos_phy_sf_sflx_qvex, &
       sflx_qtrc => atmos_phy_sf_sflx_qtrc, &
       sflx_engi => atmos_phy_sf_sflx_engi, &
       ustar     => atmos_phy_sf_ustar,     &
       tstar     => atmos_phy_sf_tstar,     &
       qstar     => atmos_phy_sf_qstar,     &
       wstar     => atmos_phy_sf_wstar,     &
       rlmo      => atmos_phy_sf_rlmo,      &
       u10       => atmos_phy_sf_u10,       &
       v10       => atmos_phy_sf_v10,       &
       t2        => atmos_phy_sf_t2,        &
       q2        => atmos_phy_sf_q2
    use mod_cpl_admin, only: &
       cpl_sw
    implicit none
 
    logical, intent(in) :: update_flag
 
    real(RP) :: ATM_W   (IA,JA)
    real(RP) :: ATM_U   (IA,JA)
    real(RP) :: ATM_V   (IA,JA)
    real(RP) :: ATM_DENS(IA,JA)
    real(RP) :: ATM_TEMP(IA,JA)
    real(RP) :: ATM_PRES(IA,JA)
    real(RP) :: ATM_QV  (IA,JA)
    real(RP) :: SFC_POTV(IA,JA)
    real(RP) :: SFLX_SH2(IA,JA)
    real(RP) :: SFLX_QV (IA,JA)
    real(RP) :: CP_t, CV_t
    real(RP) :: ENGI_t
    real(RP) :: rdz
    real(RP) :: work
    real(RP) :: kappa
 
    integer  :: i, j, iq
    !---------------------------------------------------------------------------
 
    if ( update_flag ) then
 
       !$acc data create(ATM_DENS,ATM_TEMP,ATM_PRES,SFLX_QV)
 
       ! update surface density, surface pressure
       call bottom_estimate( ka, ks,  ke, ia, is, ie, ja, js, je, &
                             dens(:,:,:), pres(:,:,:), qv(:,:,:), & ! [IN]
                             sfc_temp(:,:),                       & ! [IN]
                             fz(:,:,:),                           & ! [IN]
                             sfc_dens(:,:), sfc_pres(:,:)         ) ! [OUT]
 
       !$omp parallel do
       !$acc kernels
       do j = js, je
       do i = is, ie
          atm_dens(i,j) = dens(ks,i,j)
          atm_temp(i,j) = temp(ks,i,j)
          atm_pres(i,j) = pres(ks,i,j)
       end do
       end do
       !$acc end kernels
 
       if ( cpl_sw ) then
 
          !$acc data create(SFLX_SH2)
 
          !$omp parallel do
          !$acc kernels
          do j = js, je
          do i = is, ie
             sflx_sh2(i,j) = sflx_sh(i,j) - sflx_shex(i,j)
          end do
          end do
          !$acc end kernels
 
          if ( atmos_hydrometeor_dry ) then
             !$omp parallel do
             !$acc kernels
             do j = js, je
             do i = is, ie
                sflx_qv(i,j) = 0.0_rp
                sfc_potv(i,j) = sfc_temp(i,j) * ( pre00 / sfc_pres(i,j) )**( rdry / cpdry )
             end do
             end do
             !$acc end kernels
          else
             !$omp parallel do private(kappa)
             !$acc kernels
             do j = js, je
             do i = is, ie
                sflx_qv(i,j) = sflx_qtrc(i,j,i_qv) - sflx_qvex(i,j)
                kappa = ( rdry + ( rvap - rdry ) * qv(ks,i,j) ) &
                     / ( cpdry + ( cpvap - cpdry ) * qv(ks,i,j) )
                sfc_potv(i,j) = sfc_temp(i,j) * ( pre00 / sfc_pres(i,j) )**kappa &
                              * ( 1.0_rp + epstvap * qv(ks,i,j) )
             end do
             end do
             !$acc end kernels
          end if
 
          call bulkflux_diagnose_scales( ia, is, ie, ja, js, je, &
                                         sflx_mw(:,:), sflx_mu(:,:), sflx_mv(:,:),  & ! [IN]
                                         sflx_sh2(:,:), sflx_qv(:,:),               & ! [IN]
                                         sfc_dens(:,:), sfc_potv(:,:), pbl_zi(:,:), & ! [IN]
                                         ustar(:,:), tstar(:,:), qstar(:,:),        & ! [OUT]
                                         wstar(:,:), rlmo(:,:)                      ) ! [OUT]
          !$acc end data
 
       else
 
          !$acc data create(ATM_U,ATM_V,ATM_W,ATM_QV)
 
          !$omp parallel do
          !$acc kernels
          do j = js, je
          do i = is, ie
             atm_u(i,j) = u(ks,i,j)
             atm_v(i,j) = v(ks,i,j)
             atm_w(i,j) = atm_u(i,j) * tansl_x(i,j) + atm_v(i,j) * tansl_y(i,j)
             atm_qv(i,j) = qv(ks,i,j)
          enddo
          enddo
          !$acc end kernels
 
          select case ( atmos_phy_sf_type )
          case ( 'BULK' )
 
             call atmos_phy_sf_bulk_flux( ia, is, ie, ja, js, je,                      & ! [IN]
                                          atm_w(:,:), atm_u(:,:), atm_v(:,:),          & ! [IN]
                                          atm_temp(:,:), atm_pres(:,:), atm_qv(:,:),   & ! [IN]
                                          sfc_dens(:,:), sfc_temp(:,:), sfc_pres(:,:), & ! [IN]
                                          sfc_z0m(:,:), sfc_z0h(:,:), sfc_z0e(:,:),    & ! [IN]
                                          pbl_zi(:,:), z1(:,:),                        & ! [IN]
                                          sflx_mw(:,:), sflx_mu(:,:), sflx_mv(:,:),    & ! [OUT]
                                          sflx_sh(:,:), sflx_lh(:,:), sflx_qv(:,:),    & ! [OUT]
                                          ustar(:,:), tstar(:,:), qstar(:,:),          & ! [OUT]
                                          wstar(:,:),                                  & ! [OUT]
                                          rlmo(:,:),                                   & ! [OUT]
                                          u10(:,:), v10(:,:), t2(:,:), q2(:,:)         ) ! [OUT]
 
          case ( 'CONST' )
 
             !$acc update host(ATM_W,ATM_U,ATM_V,ATM_TEMP,SFC_DENS)
             call atmos_phy_sf_const_flux( ia, is, ie, ja, js, je,                            & ! [IN]
                                           atm_w(:,:), atm_u(:,:), atm_v(:,:), sfc_temp(:,:), & ! [IN]
                                           z1(:,:), sfc_dens(:,:),                            & ! [IN]
                                           sflx_mw(:,:), sflx_mu(:,:), sflx_mv(:,:),          & ! [OUT]
                                           sflx_sh(:,:), sflx_lh(:,:), sflx_qv(:,:),          & ! [OUT]
                                           u10(:,:), v10(:,:)                                 ) ! [OUT]
             ustar(:,:) = undef
             tstar(:,:) = undef
             qstar(:,:) = undef
             rlmo(:,:) = undef
             t2(:,:) = atm_temp(:,:)
             q2(:,:) = atm_qv(:,:)
             !$acc update device(SFLX_MW,SFLX_MU,SFLX_MV,SFLX_SH,SFLX_LH,SFLX_QV,Ustar,Tstar,Qstar,Wstar,RLmo,U10,V10,T2,Q2)
 
          end select
 
          if ( .NOT. atmos_hydrometeor_dry ) then
             !$acc kernels
             sflx_qtrc(:,:,i_qv) = sflx_qv(:,:)
             sflx_engi(:,:)      = sflx_qv(:,:) * ( tracer_cv(i_qv) * sfc_temp(:,:) + lhv )
             !$acc end kernels
          endif
 
          !$acc end data
 
       endif
 
       call history_output
 
!OCL XFILL
       !$omp parallel do &
       !$omp private(rdz)
       !$acc kernels
       do j = js, je
       do i = is, ie
          rdz = 1.0_rp / ( fz(ks,i,j) - fz(ks-1,i,j) )
          momz_t_sf(i,j) = sflx_mw(i,j) / ( cz(ks+1,i,j) - cz(ks,i,j) )
          rhou_t_sf(i,j) = sflx_mu(i,j) * rdz
          rhov_t_sf(i,j) = sflx_mv(i,j) * rdz
       enddo
       enddo
       !$acc end kernels
 
       !$omp parallel do &
       !$omp private(work,rdz,CP_t,CV_t,ENGI_t)
       !$acc kernels
       do j = js, je
       do i = is, ie
          rdz = 1.0_rp / ( fz(ks,i,j) - fz(ks-1,i,j) )
          dens_t_sf(i,j) = 0.0_rp
          cp_t = 0.0_rp
          cv_t = 0.0_rp
          engi_t = sflx_engi(i,j) * rdz
          !$acc loop seq
          do iq = 1, qa
             work = sflx_qtrc(i,j,iq) * rdz
 
             rhoq_t_sf(i,j,iq) = work
             dens_t_sf(i,j)    = dens_t_sf(i,j) + work * tracer_mass(iq)
             cp_t              = cp_t + work * tracer_cp(iq)
             cv_t              = cv_t + work * tracer_cv(iq)
             engi_t            = engi_t - tracer_engi0(iq) * work
          end do
          cp_t = ( cp_t - cptot(ks,i,j) * dens_t_sf(i,j) ) / atm_dens(i,j)
          cv_t = ( cv_t - cvtot(ks,i,j) * dens_t_sf(i,j) ) / atm_dens(i,j)
 
          rhoh_sf(i,j) = sflx_sh(i,j) * rdz + engi_t &
                       - ( cp_t + log( atm_pres(i,j) / pre00 ) * ( cvtot(ks,i,j) / cptot(ks,i,j) * cp_t - cv_t ) ) * atm_dens(i,j) * atm_temp(i,j)
       enddo
       enddo
       !$acc end kernels
 
       !$acc end data
 
    endif
 
    !$omp parallel do
    !$acc kernels
    do j = js, je
    do i = is, ie
       momz_t(ks,i,j) = momz_t(ks,i,j) + momz_t_sf(i,j)
       rhou_t(ks,i,j) = rhou_t(ks,i,j) + rhou_t_sf(i,j)
       rhov_t(ks,i,j) = rhov_t(ks,i,j) + rhov_t_sf(i,j)
       rhoh(ks,i,j) = rhoh(ks,i,j) + rhoh_sf(i,j)
       dens_t(ks,i,j) = dens_t(ks,i,j) + dens_t_sf(i,j)
    enddo
    enddo
    !$acc end kernels
 
    !$omp parallel
    !$acc kernels
    do iq = 1, qa
    !$omp do
    do j = js, je
    do i = is, ie
       rhoq_t(ks,i,j,iq) = rhoq_t(ks,i,j,iq) + rhoq_t_sf(i,j,iq)
    enddo
    enddo
    !$omp end do nowait
    enddo
    !$acc end kernels
    !$omp end parallel
 
    if ( statistics_checktotal ) then
 
       if ( .NOT. atmos_hydrometeor_dry ) then
          do iq = 1, qa
             call statistics_total( ia, is, ie, ja, js, je, &
                                    sflx_qtrc(:,:,iq), 'SFLX_'//trim(tracer_name(iq)), &
                                    atmos_grid_cartesc_real_area(:,:),                 &
                                    atmos_grid_cartesc_real_totarea                    )
          enddo
       endif
 
       call statistics_total( ia, is, ie, ja, js, je, &
                              dens_t_sf(:,:), 'DENS_t_SF',       &
                              atmos_grid_cartesc_real_area(:,:), &
                              atmos_grid_cartesc_real_totarea    )
       call statistics_total( ia, is, ie, ja, js, je, &
                              momz_t_sf(:,:), 'MOMZ_t_SF',       &
                              atmos_grid_cartesc_real_area(:,:), &
                              atmos_grid_cartesc_real_totarea    )
       call statistics_total( ia, is, ie, ja, js, je, &
                              rhou_t_sf(:,:), 'RHOU_t_SF',       &
                              atmos_grid_cartesc_real_area(:,:), &
                              atmos_grid_cartesc_real_totarea    )
       call statistics_total( ia, is, ie, ja, js, je, &
                              rhov_t_sf(:,:), 'RHOV_t_SF',       &
                              atmos_grid_cartesc_real_area(:,:), &
                              atmos_grid_cartesc_real_totarea    )
       call statistics_total( ia, is, ie, ja, js, je, &
                              rhoh_sf(:,:), 'RHOH_SF',         &
                              atmos_grid_cartesc_real_area(:,:), &
                              atmos_grid_cartesc_real_totarea    )
 
       do iq = 1, qa
          call statistics_total( ia, is, ie, ja, js, je, &
                                 rhoq_t_sf(:,:,iq), trim(tracer_name(iq))//'_t_SF', &
                                 atmos_grid_cartesc_real_area(:,:),                 &
                                 atmos_grid_cartesc_real_totarea                    )
       enddo
 
    endif
 
    return

References scale_atmos_bottom::atmos_bottom_estimate(), scale_atmos_grid_cartesc_real::atmos_grid_cartesc_real_area, scale_atmos_grid_cartesc_real::atmos_grid_cartesc_real_cz, scale_atmos_grid_cartesc_real::atmos_grid_cartesc_real_fz, scale_atmos_grid_cartesc_real::atmos_grid_cartesc_real_totarea, scale_atmos_grid_cartesc_real::atmos_grid_cartesc_real_z1, scale_atmos_hydrometeor::atmos_hydrometeor_dry, mod_atmos_phy_bl_vars::atmos_phy_bl_zi, mod_atmos_phy_rd_vars::atmos_phy_rd_sflx_lw_dn, mod_atmos_phy_rd_vars::atmos_phy_rd_sflx_sw_dn, scale_atmos_phy_sf_bulk::atmos_phy_sf_bulk_flux(), scale_atmos_phy_sf_const::atmos_phy_sf_const_flux(), mod_atmos_phy_sf_vars::atmos_phy_sf_dens_t, mod_atmos_phy_sf_vars::atmos_phy_sf_momz_t, mod_atmos_phy_sf_vars::atmos_phy_sf_q2, mod_atmos_phy_sf_vars::atmos_phy_sf_qstar, mod_atmos_phy_sf_vars::atmos_phy_sf_rhoh, mod_atmos_phy_sf_vars::atmos_phy_sf_rhoq_t, mod_atmos_phy_sf_vars::atmos_phy_sf_rhou_t, mod_atmos_phy_sf_vars::atmos_phy_sf_rhov_t, mod_atmos_phy_sf_vars::atmos_phy_sf_rlmo, mod_atmos_phy_sf_vars::atmos_phy_sf_sfc_dens, mod_atmos_phy_sf_vars::atmos_phy_sf_sfc_pres, mod_atmos_phy_sf_vars::atmos_phy_sf_sfc_temp, mod_atmos_phy_sf_vars::atmos_phy_sf_sfc_z0e, mod_atmos_phy_sf_vars::atmos_phy_sf_sfc_z0h, mod_atmos_phy_sf_vars::atmos_phy_sf_sfc_z0m, mod_atmos_phy_sf_vars::atmos_phy_sf_sflx_engi, mod_atmos_phy_sf_vars::atmos_phy_sf_sflx_lh, mod_atmos_phy_sf_vars::atmos_phy_sf_sflx_mu, mod_atmos_phy_sf_vars::atmos_phy_sf_sflx_mv, mod_atmos_phy_sf_vars::atmos_phy_sf_sflx_mw, mod_atmos_phy_sf_vars::atmos_phy_sf_sflx_qtrc, mod_atmos_phy_sf_vars::atmos_phy_sf_sflx_qvex, mod_atmos_phy_sf_vars::atmos_phy_sf_sflx_sh, mod_atmos_phy_sf_vars::atmos_phy_sf_sflx_shex, mod_atmos_phy_sf_vars::atmos_phy_sf_t2, mod_atmos_phy_sf_vars::atmos_phy_sf_tstar, mod_atmos_admin::atmos_phy_sf_type, mod_atmos_phy_sf_vars::atmos_phy_sf_u10, mod_atmos_phy_sf_vars::atmos_phy_sf_ustar, mod_atmos_phy_sf_vars::atmos_phy_sf_v10, mod_atmos_phy_sf_vars::atmos_phy_sf_wstar, scale_const::const_cpdry, scale_const::const_cpvap, scale_const::const_epstvap, scale_const::const_pre00, scale_const::const_rdry, scale_const::const_rvap, scale_const::const_undef, mod_cpl_admin::cpl_sw, mod_atmos_vars::cptot, mod_atmos_vars::cvtot, mod_atmos_vars::dens, mod_atmos_vars::dens_av, mod_atmos_vars::dens_tp, history_output(), scale_atmos_hydrometeor::i_qv, scale_atmos_grid_cartesc_index::ia, scale_atmos_grid_cartesc_index::ie, scale_atmos_grid_cartesc_index::is, scale_atmos_grid_cartesc_index::ja, scale_atmos_grid_cartesc_index::je, scale_atmos_grid_cartesc_index::js, scale_atmos_grid_cartesc_index::ka, scale_atmos_grid_cartesc_index::ke, scale_atmos_grid_cartesc_index::ks, scale_atmos_hydrometeor::lhv, mod_atmos_vars::momz_tp, mod_atmos_vars::pott, mod_atmos_vars::pres, scale_tracer::qa, mod_atmos_vars::qv, mod_atmos_vars::rhoh_p, mod_atmos_vars::rhoq_tp, mod_atmos_vars::rhot, mod_atmos_vars::rhot_av, mod_atmos_vars::rhou_tp, mod_atmos_vars::rhov_tp, scale_statistics::statistics_checktotal, mod_atmos_vars::temp, scale_time::time_dtsec_atmos_phy_sf, scale_topography::topography_tansl_x, scale_topography::topography_tansl_y, scale_tracer::tracer_cp, scale_tracer::tracer_cv, scale_tracer::tracer_engi0, scale_tracer::tracer_mass, scale_tracer::tracer_name, mod_atmos_vars::u, mod_atmos_vars::v, and mod_atmos_vars::w.

Referenced by mod_atmos_driver::atmos_driver_calc_tendency(), and mod_atmos_driver::atmos_driver_calc_tendency_from_sflux().

Here is the call graph for this function:

Here is the caller graph for this function:

history_output()

subroutine mod_atmos_phy_sf_driver::history_output

Definition at line 515 of file mod_atmos_phy_sf_driver.F90.

    use scale_const, only: &
       undef => const_undef, &
       rvap  => const_rvap
    use scale_file_history, only: &
       file_history_query, &
       file_history_in, &
       file_history_put
    use scale_atmos_hydrostatic, only: &
       barometric_law_mslp => atmos_hydrostatic_barometric_law_mslp
    use scale_atmos_saturation, only: &
       atmos_saturation_psat_liq
    use scale_atmos_grid_cartesc_real, only: &
       real_cz => atmos_grid_cartesc_real_cz
    use scale_atmos_grid_cartesc_metric, only: &
       rotc => atmos_grid_cartesc_metric_rotc
    use scale_atmos_hydrometeor, only: &
       atmos_hydrometeor_dry, &
       i_qv
    use mod_atmos_vars, only: &
       temp, &
       pres, &
       qv
    use mod_atmos_phy_sf_vars, only: &
       sfc_dens   => atmos_phy_sf_sfc_dens,   &
       sfc_pres   => atmos_phy_sf_sfc_pres,   &
       sfc_temp   => atmos_phy_sf_sfc_temp,   &
       sfc_albedo => atmos_phy_sf_sfc_albedo, &
       sfc_z0m    => atmos_phy_sf_sfc_z0m,    &
       sfc_z0h    => atmos_phy_sf_sfc_z0h,    &
       sfc_z0e    => atmos_phy_sf_sfc_z0e,    &
       sflx_mw    => atmos_phy_sf_sflx_mw,    &
       sflx_mu    => atmos_phy_sf_sflx_mu,    &
       sflx_mv    => atmos_phy_sf_sflx_mv,    &
       sflx_sh    => atmos_phy_sf_sflx_sh,    &
       sflx_lh    => atmos_phy_sf_sflx_lh,    &
       sflx_gh    => atmos_phy_sf_sflx_gh,    &
       sflx_qtrc  => atmos_phy_sf_sflx_qtrc,  &
       sflx_engi  => atmos_phy_sf_sflx_engi,  &
       ustar      => atmos_phy_sf_ustar,      &
       tstar      => atmos_phy_sf_tstar,      &
       qstar      => atmos_phy_sf_qstar,      &
       wstar      => atmos_phy_sf_wstar,      &
       rlmo       => atmos_phy_sf_rlmo,       &
       u10        => atmos_phy_sf_u10,        &
       v10        => atmos_phy_sf_v10,        &
       t2         => atmos_phy_sf_t2,         &
       q2         => atmos_phy_sf_q2
    implicit none
 
    real(RP) :: work(IA,JA)
 
    logical :: do_put
    integer :: i, j, iq
    !---------------------------------------------------------------------------
 
    call file_history_in( sfc_dens(:,:),                     'SFC_DENS',        'surface atmospheric density',          'kg/m3'   )
    call file_history_in( sfc_pres(:,:),                     'SFC_PRES',        'surface atmospheric pressure',         'Pa'      )
    call file_history_in( sfc_temp(:,:),                     'SFC_TEMP',        'surface skin temperature',             'K'       )
    call file_history_in( sfc_albedo(:,:,i_r_direct ,i_r_ir ), 'SFC_ALB_IR_dir' , 'surface albedo (IR; direct',           '1'       , fill_halo=.true. )
    call file_history_in( sfc_albedo(:,:,i_r_diffuse,i_r_ir ), 'SFC_ALB_IR_dif' , 'surface albedo (IR; diffuse)',         '1'       , fill_halo=.true. )
    call file_history_in( sfc_albedo(:,:,i_r_direct ,i_r_nir), 'SFC_ALB_NIR_dir', 'surface albedo (NIR; direct',          '1'       , fill_halo=.true. )
    call file_history_in( sfc_albedo(:,:,i_r_diffuse,i_r_nir), 'SFC_ALB_NIR_dif', 'surface albedo (NIR; diffuse',         '1'       , fill_halo=.true. )
    call file_history_in( sfc_albedo(:,:,i_r_direct ,i_r_vis), 'SFC_ALB_VIS_dir', 'surface albedo (VIS; direct',          '1'       , fill_halo=.true. )
    call file_history_in( sfc_albedo(:,:,i_r_diffuse,i_r_vis), 'SFC_ALB_VIS_dif', 'surface albedo (VIS; diffuse',         '1'       , fill_halo=.true. )
    call file_history_in( sfc_z0m(:,:),                     'SFC_Z0M',         'roughness length (momentum)',           'm'       , fill_halo=.true. )
    call file_history_in( sfc_z0h(:,:),                     'SFC_Z0H',         'roughness length (heat)',               'm'       , fill_halo=.true. )
    call file_history_in( sfc_z0e(:,:),                     'SFC_Z0E',         'roughness length (vapor)',              'm'       , fill_halo=.true. )
    call file_history_in( sflx_mw(:,:),                     'MWFLX',           'w-momentum flux',                       'kg/m/s2' )
    call file_history_in( sflx_mu(:,:),                     'MUFLX',           'u-momentum flux',                       'kg/m/s2' )
    call file_history_in( sflx_mv(:,:),                     'MVFLX',           'v-momentum flux',                       'kg/m/s2' )
    call file_history_in( sflx_sh(:,:),                     'SHFLX',           'sensible heat flux',                    'W/m2'    , fill_halo=.true. )
    call file_history_in( sflx_lh(:,:),                     'LHFLX',           'latent heat flux',                      'W/m2'    , fill_halo=.true. )
    call file_history_in( sflx_gh(:,:),                     'GHFLX',           'ground heat flux (downward)',           'W/m2'    , fill_halo=.true. )
    do iq = 1, qa
       call file_history_in( sflx_qtrc(:,:,iq), 'SFLX_'//trim(tracer_name(iq)), &
                             'surface '//trim(tracer_name(iq))//' flux',        &
                             'kg/m2/s' , fill_halo=.true.                       )
    enddo
    call file_history_in( sflx_engi(:,:),                     'SFLX_ENGI',        'ground internal energy flux (merged)', 'W/m2'    , fill_halo=.true. )
 
    call file_history_in( ustar(:,:), 'Ustar',  'friction velocity',         'm/s'  , fill_halo=.true. )
    call file_history_in( tstar(:,:), 'Tstar',  'temperature scale',         'K'    , fill_halo=.true. )
    call file_history_in( qstar(:,:), 'Qstar',  'moisuter scale',            'kg/kg', fill_halo=.true. )
    call file_history_in( wstar(:,:), 'Wstar',  'convective velocity scale', 'm/s',   fill_halo=.true. )
    call file_history_in( rlmo(:,:), 'RLmo',   'inverse of Obukhov length', '1/m'  , fill_halo=.true. )
 
    call file_history_in( u10(:,:), 'U10',    '10m x-wind',                'm/s'  , fill_halo=.true. )
    call file_history_in( v10(:,:), 'V10',    '10m y-wind',                'm/s'  , fill_halo=.true. )
    call file_history_in( t2(:,:), 'T2 ',    '2m air temperature',        'K'    , fill_halo=.true. )
    call file_history_in( q2(:,:), 'Q2 ',    '2m specific humidity',      'kg/kg', fill_halo=.true. )
 
 
    !$acc data create(work)
 
    call file_history_query( hist_uabs10, do_put )
    if ( do_put ) then
       !$omp parallel do
       !$acc kernels
       do j = js, je
       do i = is, ie
          work(i,j) = sqrt( u10(i,j)**2 + v10(i,j)**2 )
       enddo
       enddo
       !$acc end kernels
       call file_history_put( hist_uabs10, work(:,:) )
    end if
 
    call file_history_query( hist_u10m, do_put )
    if ( do_put ) then
       !$omp parallel do
       !$acc kernels
       do j = js, je
       do i = is, ie
          work(i,j) = u10(i,j) * rotc(i,j,1) - v10(i,j) * rotc(i,j,2)
       enddo
       enddo
       !$acc end kernels
       call file_history_put( hist_u10m, work(:,:) )
    end if
 
    call file_history_query( hist_v10m, do_put )
    if ( do_put ) then
       !$omp parallel do
       !$acc kernels
       do j = js, je
       do i = is, ie
          work(i,j) = u10(i,j) * rotc(i,j,2) + v10(i,j) * rotc(i,j,1)
       enddo
       enddo
       !$acc end kernels
       call file_history_put( hist_v10m, work(:,:) )
    end if
 
    call file_history_query( hist_rh2, do_put )
    if ( do_put ) then
       call atmos_saturation_psat_liq( &
            ia, is, ie, ja, js, je, &
            t2(:,:), & ! (in)
            work(:,:)  ) ! (out)
       !$omp parallel do
       !$acc kernels
       do j = js, je
       do i = is, ie
          work(i,j) = sfc_dens(i,j) * q2(i,j) &
                     / work(i,j) * rvap * t2(i,j) &
                     * 100.0_rp
       enddo
       enddo
       !$acc end kernels
       call file_history_put( hist_rh2, work(:,:) )
    end if
 
 
    call file_history_query( hist_mslp, do_put )
    if ( do_put ) then
       call barometric_law_mslp( ka, ks, ke, ia, is, ie, ja, js, je,  & ! [IN]
                                 pres(:,:,:), temp(:,:,:), qv(:,:,:), & ! [IN]
                                 real_cz(:,:,:),                      & ! [IN]
                                 work(:,:)                            ) ! [OUT]
       call file_history_put( hist_mslp, work(:,:) )
    end if
 
    !$acc end data
 
    return

References scale_atmos_grid_cartesc_metric::atmos_grid_cartesc_metric_rotc, scale_atmos_grid_cartesc_real::atmos_grid_cartesc_real_cz, scale_atmos_hydrometeor::atmos_hydrometeor_dry, mod_atmos_phy_sf_vars::atmos_phy_sf_q2, mod_atmos_phy_sf_vars::atmos_phy_sf_qstar, mod_atmos_phy_sf_vars::atmos_phy_sf_rlmo, mod_atmos_phy_sf_vars::atmos_phy_sf_sfc_albedo, mod_atmos_phy_sf_vars::atmos_phy_sf_sfc_dens, mod_atmos_phy_sf_vars::atmos_phy_sf_sfc_pres, mod_atmos_phy_sf_vars::atmos_phy_sf_sfc_temp, mod_atmos_phy_sf_vars::atmos_phy_sf_sfc_z0e, mod_atmos_phy_sf_vars::atmos_phy_sf_sfc_z0h, mod_atmos_phy_sf_vars::atmos_phy_sf_sfc_z0m, mod_atmos_phy_sf_vars::atmos_phy_sf_sflx_engi, mod_atmos_phy_sf_vars::atmos_phy_sf_sflx_gh, mod_atmos_phy_sf_vars::atmos_phy_sf_sflx_lh, mod_atmos_phy_sf_vars::atmos_phy_sf_sflx_mu, mod_atmos_phy_sf_vars::atmos_phy_sf_sflx_mv, mod_atmos_phy_sf_vars::atmos_phy_sf_sflx_mw, mod_atmos_phy_sf_vars::atmos_phy_sf_sflx_qtrc, mod_atmos_phy_sf_vars::atmos_phy_sf_sflx_sh, mod_atmos_phy_sf_vars::atmos_phy_sf_t2, mod_atmos_phy_sf_vars::atmos_phy_sf_tstar, mod_atmos_phy_sf_vars::atmos_phy_sf_u10, mod_atmos_phy_sf_vars::atmos_phy_sf_ustar, mod_atmos_phy_sf_vars::atmos_phy_sf_v10, mod_atmos_phy_sf_vars::atmos_phy_sf_wstar, scale_const::const_rvap, scale_const::const_undef, scale_atmos_hydrometeor::i_qv, 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_atmos_grid_cartesc_index::ie, scale_atmos_grid_cartesc_index::is, scale_atmos_grid_cartesc_index::je, scale_atmos_grid_cartesc_index::js, scale_atmos_grid_cartesc_index::ka, scale_atmos_grid_cartesc_index::ke, scale_atmos_grid_cartesc_index::ks, mod_atmos_vars::pres, scale_tracer::qa, mod_atmos_vars::qv, mod_atmos_vars::temp, and scale_tracer::tracer_name.

Referenced by atmos_phy_sf_driver_calc_tendency().

Here is the caller graph for this function: