module atmosphere / physics / surface / simple More...

Functions/Subroutines
subroutine, public	atmos_phy_sf_simple_setup
	Setup.
subroutine, public	atmos_phy_sf_simple_flux (ia, is, ie, ja, js, je, atm_w, atm_u, atm_v, atm_temp, atm_pres, atm_qv, sfc_dens, sfc_temp, sfc_pres, atm_z1, sflx_mw, sflx_mu, sflx_mv, sflx_sh, sflx_lh, sflx_qv, u10, v10)
	Calculate surface fluxes with constant bulk coefficients.

Detailed Description

module atmosphere / physics / surface / simple

Description: Flux from/to bottom wall of atmosphere (surface) Constant bulk coefficient

Author: Team SCALE

NAMELIST

PARAM_ATMOS_PHY_SF_SIMPLE

name	type	default value	comment
ATMOS_PHY_SF_U_MINM	real(RP)	0.0_RP	> minimum limit of absolute velocity for momentum [m/s]
ATMOS_PHY_SF_CONST_CM	real(RP)	0.0011_RP	> constant bulk coefficient for momentum [NIL]
ATMOS_PHY_SF_CONST_CH	real(RP)	0.0044_RP	> constant bulk coefficient for heat [NIL]
ATMOS_PHY_SF_CONST_CE	real(RP)	0.0044_RP	> constant bulk coefficient for evaporation [NIL]
ATMOS_PHY_SF_BULK_BETA	real(RP)	1.0_RP	> evaporation efficiency (0-1)

History Output: No history output

Function/Subroutine Documentation

◆ atmos_phy_sf_simple_setup()

subroutine, public scale_atm_phy_sf_bulk_simple::atmos_phy_sf_simple_setup

Setup.

Definition at line 64 of file scale_atm_phy_sf_bulk_simple.F90.

    use scale_prc, only: &
       prc_abort
    implicit none
 
    namelist / param_atmos_phy_sf_simple / &
       atmos_phy_sf_u_minm,         &
       atmos_phy_sf_const_cm,       &
       atmos_phy_sf_const_ch,       &
       atmos_phy_sf_const_ce,       &
       atmos_phy_sf_bulk_beta
 
    integer :: ierr
    !---------------------------------------------------------------------------
 
    log_newline
    log_info("ATMOS_PHY_SF_simple_setup",*) 'Setup'
    log_info("ATMOS_PHY_SF_simple_setup",*) 'Simple flux'
 
    !--- read namelist
    rewind(io_fid_conf)
    read(io_fid_conf,nml=param_atmos_phy_sf_simple,iostat=ierr)
    if( ierr < 0 ) then !--- missing
       log_info("ATMOS_PHY_SF_simple_setup",*) 'Not found namelist. Default used.'
    elseif( ierr > 0 ) then !--- fatal error
       log_error("ATMOS_PHY_SF_simple_setup",*) 'Not appropriate names in namelist PARAM_ATMOS_PHY_SF_SIMPLE. Check!'
       call prc_abort
    endif
    log_nml(param_atmos_phy_sf_simple)
 
    return

◆ atmos_phy_sf_simple_flux()

subroutine, public scale_atm_phy_sf_bulk_simple::atmos_phy_sf_simple_flux	(	integer, intent(in)	ia,
		integer, intent(in)	is,
		integer, intent(in)	ie,
		integer, intent(in)	ja,
		integer, intent(in)	js,
		integer, intent(in)	je,
		real(rp), dimension (ia,ja), intent(in)	atm_w,
		real(rp), dimension (ia,ja), intent(in)	atm_u,
		real(rp), dimension (ia,ja), intent(in)	atm_v,
		real(rp), dimension(ia,ja), intent(in)	atm_temp,
		real(rp), dimension(ia,ja), intent(in)	atm_pres,
		real(rp), dimension (ia,ja), intent(in)	atm_qv,
		real(rp), dimension(ia,ja), intent(in)	sfc_dens,
		real(rp), dimension(ia,ja), intent(in)	sfc_temp,
		real(rp), dimension(ia,ja), intent(in)	sfc_pres,
		real(rp), dimension (ia,ja), intent(in)	atm_z1,
		real(rp), dimension(ia,ja), intent(out)	sflx_mw,
		real(rp), dimension(ia,ja), intent(out)	sflx_mu,
		real(rp), dimension(ia,ja), intent(out)	sflx_mv,
		real(rp), dimension(ia,ja), intent(out)	sflx_sh,
		real(rp), dimension(ia,ja), intent(out)	sflx_lh,
		real(rp), dimension(ia,ja), intent(out)	sflx_qv,
		real(rp), dimension (ia,ja), intent(out)	u10,
		real(rp), dimension (ia,ja), intent(out)	v10 )

Calculate surface fluxes with constant bulk coefficients.

Parameters

[in]	atm_w	velocity w at the lowermost layer (cell center) [m/s]
[in]	atm_u	velocity u at the lowermost layer (cell center) [m/s]
[in]	atm_v	velocity v at the lowermost layer (cell center) [m/s]
[in]	atm_temp	temperature at the lowermost layer (cell center) [K]
[in]	atm_pres	pressure at the lowermost layer (cell center) [Pa]
[in]	atm_qv	qv at the lowermost layer (cell center) [kg/kg]
[in]	sfc_dens	density at the surface atmosphere [kg/m3]
[in]	sfc_temp	temperature at the surface atmosphere [K]
[in]	sfc_pres	pressure at the surface atmosphere [Pa]
[in]	atm_z1	height of the lowermost grid from surface (cell center) [m]
[out]	sflx_mw	surface flux for z-momentum (area center) [m/s*kg/m2/s]
[out]	sflx_mu	surface flux for x-momentum (area center) [m/s*kg/m2/s]
[out]	sflx_mv	surface flux for y-momentum (area center) [m/s*kg/m2/s]
[out]	sflx_sh	surface flux for sensible heat (area center) [J/m2/s]
[out]	sflx_lh	surface flux for latent heat (area center) [J/m2/s]
[out]	sflx_qv	surface flux for qv (area center) [kg/m2/s]
[out]	u10	velocity u at 10m height
[out]	v10	velocity v at 10m height

Definition at line 100 of file scale_atm_phy_sf_bulk_simple.F90.

    use scale_const, only: &
       pi    => const_pi
    use scale_atmos_hydrometeor, only: &
       hydrometeor_lhv => atmos_hydrometeor_lhv
    use scale_atmos_saturation, only: &
       saturation_psat_all => atmos_saturation_psat_all
    use scale_time, only: &
       time_nowsec
    implicit none
    integer, intent(in) :: IA, IS, IE
    integer, intent(in) :: JA, JS, JE
 
    real(RP), intent(in) :: ATM_W   (IA,JA) !< velocity w  at the lowermost layer (cell center) [m/s]
    real(RP), intent(in) :: ATM_U   (IA,JA) !< velocity u  at the lowermost layer (cell center) [m/s]
    real(RP), intent(in) :: ATM_V   (IA,JA) !< velocity v  at the lowermost layer (cell center) [m/s]
    real(RP), intent(in) :: ATM_TEMP(IA,JA) !< temperature at the lowermost layer (cell center) [K]
    real(RP), intent(in) :: ATM_PRES(IA,JA) !< pressure    at the lowermost layer (cell center) [Pa]
    real(RP), intent(in) :: ATM_QV  (IA,JA) !< qv          at the lowermost layer (cell center) [kg/kg]
    real(RP), intent(in) :: SFC_DENS(IA,JA) !< density     at the surface atmosphere [kg/m3]
    real(RP), intent(in) :: SFC_TEMP(IA,JA) !< temperature  at the surface atmosphere [K]
    real(RP), intent(in) :: SFC_PRES(IA,JA) !< pressure    at the surface atmosphere [Pa]
    real(RP), intent(in) :: ATM_Z1  (IA,JA) !< height of the lowermost grid from surface (cell center) [m]
 
    real(RP), intent(out) :: SFLX_MW(IA,JA) !< surface flux for z-momentum    (area center)   [m/s*kg/m2/s]
    real(RP), intent(out) :: SFLX_MU(IA,JA) !< surface flux for x-momentum    (area center)   [m/s*kg/m2/s]
    real(RP), intent(out) :: SFLX_MV(IA,JA) !< surface flux for y-momentum    (area center)   [m/s*kg/m2/s]
    real(RP), intent(out) :: SFLX_SH(IA,JA) !< surface flux for sensible heat (area center)   [J/m2/s]
    real(RP), intent(out) :: SFLX_LH(IA,JA) !< surface flux for latent   heat (area center)   [J/m2/s]
    real(RP), intent(out) :: SFLX_QV(IA,JA) !< surface flux for qv            (area center)   [kg/m2/s]
    real(RP), intent(out) :: U10    (IA,JA) !< velocity u        at 10m height
    real(RP), intent(out) :: V10    (IA,JA) !< velocity v        at 10m height
 
    real(RP) :: ATM_Uabs(IA,JA) ! absolute velocity at z1 [m/s]
    real(RP) :: R10
 
    real(RP) :: SFC_PSAT (IA,JA) ! saturated water vapor pressure [Pa]
    real(RP) :: LHV(IA,JA)
 
    real(RP) :: SFC_QSAT      ! saturated water vapor mixing ratio [kg/kg]
    real(RP) :: SFC_QV(IA,JA) ! water vapor mixing ratio [kg/kg]
 
    integer  :: i, j
    !---------------------------------------------------------------------------
 
    log_progress(*) 'atmosphere / physics / surface flux / simple'
 
    !$omp parallel do
    do j = js, je
    do i = is, ie
       atm_uabs(i,j) = min( atmos_phy_sf_u_maxm, max( atmos_phy_sf_u_minm, &
            sqrt( atm_w(i,j)**2 + atm_u(i,j)**2 + atm_v(i,j)**2 ) ) ) ! at cell center
    enddo
    enddo
 
    !-----< momentum >-----
 
    !$omp parallel do
    do j = js, je
    do i = is, ie
       sflx_mw(i,j) = - atmos_phy_sf_const_cm * atm_uabs(i,j) * sfc_dens(i,j) * atm_w(i,j)
       sflx_mu(i,j) = - atmos_phy_sf_const_cm * atm_uabs(i,j) * sfc_dens(i,j) * atm_u(i,j)
       sflx_mv(i,j) = - atmos_phy_sf_const_cm * atm_uabs(i,j) * sfc_dens(i,j) * atm_v(i,j)
    enddo
    enddo
 
    !-----< heat & mass flux >-----
    call saturation_psat_all( ia, is, ie, ja, js, je, &
                              sfc_temp(:,:), & ! [IN]
                              sfc_psat(:,:)  ) ! [OUT]
 
    call hydrometeor_lhv( &
         ia, is, ie, ja, js, je, &
         sfc_temp(:,:), & ! [IN]
         lhv(:,:)       ) ! [OUT]
 
    !$omp parallel do private( SFC_QSAT )
    do j = js, je
    do i = is, ie
      sfc_qsat = epsvap * sfc_psat(i,j) / ( sfc_pres(i,j) - ( 1.0_rp-epsvap ) * sfc_psat(i,j) )
      sfc_qv(i,j) = ( 1.0_rp - atmos_phy_sf_bulk_beta ) * atm_qv(i,j) + atmos_phy_sf_bulk_beta * sfc_qsat
 
      sflx_sh(i,j) = atmos_phy_sf_const_ch * atm_uabs(i,j) * sfc_dens(i,j) * cpdry * ( sfc_temp(i,j) - atm_temp(i,j) )
      sflx_lh(i,j) = atmos_phy_sf_const_ce * atm_uabs(i,j) * sfc_dens(i,j) * lhv(i,j) * ( sfc_qv(i,j) - atm_qv(i,j) )
 
      sflx_qv(i,j) = sflx_lh(i,j) / lhv(i,j)
    enddo
    enddo
 
    !-----< U10, V10 >-----
 
    !$omp parallel do &
    !$omp private(R10)
    do j = js, je
    do i = is, ie
       r10 = 10.0_rp / atm_z1(i,j)
 
       u10(i,j) = r10 * atm_u(i,j)
       v10(i,j) = r10 * atm_v(i,j)
    enddo
    enddo
 
    return

Functions/Subroutines

Detailed Description

Function/Subroutine Documentation

◆ atmos_phy_sf_simple_setup()

◆ atmos_phy_sf_simple_flux()