scale_atmos_phy_sf_const Module Reference

module ATMOSPHERE / Physics Surface fluxes More...

Functions/Subroutines
subroutine, public	atmos_phy_sf_const_setup (ATMOS_PHY_SF_TYPE)
	Setup. More...
subroutine, public	atmos_phy_sf_const (ATM_TEMP, ATM_PRES, ATM_W, ATM_U, ATM_V, ATM_DENS, ATM_QTRC, ATM_Z1, dt, SFC_DENS, SFC_PRES, SFLX_LW_dn, SFLX_SW_dn, SFC_TEMP, SFC_albedo, SFC_beta, SFC_Z0M, SFC_Z0H, SFC_Z0E, SFLX_MW, SFLX_MU, SFLX_MV, SFLX_SH, SFLX_LH, SFLX_QTRC, U10, V10, T2, Q2)
	Constant flux. More...

Detailed Description

module ATMOSPHERE / Physics Surface fluxes

Description

Flux from/to bottom wall of atmosphere (surface) Constant flux

Author

Team SCALE

History

• 2011-12-03 (Y.Miyamoto) [new]
• 2011-12-11 (H.Yashiro) [mod] integrate to SCALE-LES ver.3
• 2012-03-23 (H.Yashiro) [mod] Explicit index parameter inclusion
• 2012-04-10 (Y.Miyamoto) [mod] introduce coefficients for interpolation
• 2012-09-11 (S.Nishizawa) [mod] bugfix based on the scale document
• 2012-09-12 (Y.Sato) [renew] constant FLUX version

NAMELIST

• PARAM_ATMOS_PHY_SF_CONST

  name	type	default value	comment
  ATMOS_PHY_SF_FLG_MOM_FLUX	integer	0	application type for momentum flux
  ATMOS_PHY_SF_U_MINM	real(RP)	0.0_RP	minimum limit of absolute velocity for momentum [m/s]
  ATMOS_PHY_SF_CM_MIN	real(RP)	1.0E-5_RP	minimum limit of bulk coefficient for momentum [NIL]
  ATMOS_PHY_SF_CONST_USTAR	real(RP)	0.25_RP	constant friction velocity [m/s]
  ATMOS_PHY_SF_CONST_CM	real(RP)	0.0011_RP	constant bulk coefficient for momentum [NIL]
  ATMOS_PHY_SF_CONST_SH	real(RP)	15.0_RP	constant surface sensible heat flux [W/m2]
  ATMOS_PHY_SF_CONST_LH	real(RP)	115.0_RP	constant surface latent heat flux [W/m2]
  ATMOS_PHY_SF_FLG_SH_DIURNAL	logical	.false.	diurnal modulation for sensible heat flux?
  ATMOS_PHY_SF_CONST_FREQ	real(RP)	24.0_RP	frequency of sensible heat flux modulation [hour]

  
  

History Output

No history output

Function/Subroutine Documentation

atmos_phy_sf_const_setup()

subroutine, public scale_atmos_phy_sf_const::atmos_phy_sf_const_setup

(

character(len=*), intent(in)

ATMOS_PHY_SF_TYPE

)

Setup.

Definition at line 74 of file scale_atmos_phy_sf_const.F90.

References scale_stdio::io_fid_conf, scale_stdio::io_fid_log, scale_stdio::io_l, scale_stdio::io_lnml, and scale_process::prc_mpistop().

Referenced by scale_atmos_phy_sf::atmos_phy_sf_setup().

    use scale_process, only: &
       prc_mpistop
    implicit none

    character(len=*), intent(in) :: atmos_phy_sf_type

    namelist / param_atmos_phy_sf_const / &
       atmos_phy_sf_flg_mom_flux,   &
       atmos_phy_sf_u_minm,         &
       atmos_phy_sf_cm_min,         &
       atmos_phy_sf_const_ustar,    &
       atmos_phy_sf_const_cm,       &
       atmos_phy_sf_const_sh,       &
       atmos_phy_sf_const_lh,       &
       atmos_phy_sf_flg_sh_diurnal, &
       atmos_phy_sf_const_freq

    integer :: ierr
    !---------------------------------------------------------------------------

    if( io_l ) write(io_fid_log,*)
    if( io_l ) write(io_fid_log,*) '++++++ Module[SURFACE FLUX] / Categ[ATMOS PHYSICS] / Origin[SCALElib]'
    if( io_l ) write(io_fid_log,*) '+++ Constant flux'

    if ( atmos_phy_sf_type /= 'CONST' ) then
       write(*,*) 'xxx ATMOS_PHY_SF_TYPE is not CONST. Check!'
       call prc_mpistop
    endif

    !--- read namelist
    rewind(io_fid_conf)
    read(io_fid_conf,nml=param_atmos_phy_sf_const,iostat=ierr)
    if( ierr < 0 ) then !--- missing
       if( io_l ) write(io_fid_log,*) '*** Not found namelist. Default used.'
    elseif( ierr > 0 ) then !--- fatal error
       write(*,*) 'xxx Not appropriate names in namelist PARAM_ATMOS_PHY_SF_CONST. Check!'
       call prc_mpistop
    endif
    if( io_lnml ) write(io_fid_log,nml=param_atmos_phy_sf_const)

    return

Here is the call graph for this function:

Here is the caller graph for this function:

atmos_phy_sf_const()

subroutine, public scale_atmos_phy_sf_const::atmos_phy_sf_const	(	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_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_DENS,
		real(rp), dimension (ia,ja,qa), intent(in)	ATM_QTRC,
		real(rp), dimension (ia,ja), intent(in)	ATM_Z1,
		real(dp), intent(in)	dt,
		real(rp), dimension (ia,ja), intent(in)	SFC_DENS,
		real(rp), dimension (ia,ja), intent(in)	SFC_PRES,
		real(rp), dimension(ia,ja), intent(in)	SFLX_LW_dn,
		real(rp), dimension(ia,ja), intent(in)	SFLX_SW_dn,
		real(rp), dimension (ia,ja), intent(in)	SFC_TEMP,
		real(rp), dimension(ia,ja,2), intent(in)	SFC_albedo,
		real(rp), dimension (ia,ja), intent(in)	SFC_beta,
		real(rp), dimension (ia,ja), intent(inout)	SFC_Z0M,
		real(rp), dimension (ia,ja), intent(inout)	SFC_Z0H,
		real(rp), dimension (ia,ja), intent(inout)	SFC_Z0E,
		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,qa), intent(out)	SFLX_QTRC,
		real(rp), dimension (ia,ja), intent(out)	U10,
		real(rp), dimension (ia,ja), intent(out)	V10,
		real(rp), dimension (ia,ja), intent(out)	T2,
		real(rp), dimension (ia,ja), intent(out)	Q2
	)

Constant flux.

Definition at line 131 of file scale_atmos_phy_sf_const.F90.

References scale_const::const_pi, scale_tracer::i_qv, scale_grid_index::ie, scale_stdio::io_fid_log, scale_stdio::io_l, scale_grid_index::is, scale_grid_index::je, scale_grid_index::js, and scale_time::time_nowsec.

Referenced by scale_atmos_phy_sf::atmos_phy_sf_setup().

    use scale_grid_index
    use scale_tracer
    use scale_const, only: &
       pi    => const_pi
    use scale_atmos_thermodyn, only: &
       atmos_thermodyn_templhv
    use scale_time, only: &
       time_nowsec
    implicit none

    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_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_dens  (ia,ja)    ! density     at the lowermost layer (cell center) [kg/m3]
    real(RP), intent(in)    :: atm_qtrc  (ia,ja,qa) ! tracer      at the lowermost layer (cell center) [kg/kg]
    real(RP), intent(in)    :: atm_z1    (ia,ja)    ! height of the lowermost grid from surface (cell center) [m]
    real(DP), intent(in)    :: dt                   ! delta time
    real(RP), intent(in)    :: sfc_dens  (ia,ja)    ! density     at the surface atmosphere [kg/m3]
    real(RP), intent(in)    :: sfc_pres  (ia,ja)    ! pressure    at the surface atmosphere [Pa]
    real(RP), intent(in)    :: sflx_lw_dn(ia,ja)    ! downward longwave  radiation flux at the surface [J/m2/s]
    real(RP), intent(in)    :: sflx_sw_dn(ia,ja)    ! downward shortwave radiation flux at the surface [J/m2/s]
    real(RP), intent(in)    :: sfc_temp  (ia,ja)    ! temperature at the surface skin [K]
    real(RP), intent(in)    :: sfc_albedo(ia,ja,2)  ! surface albedo (LW/SW) [0-1]
    real(RP), intent(in)    :: sfc_beta  (ia,ja)    ! evaporation efficiency [0-1]
    real(RP), intent(inout) :: sfc_z0m   (ia,ja)    ! surface roughness length (momentum) [m]
    real(RP), intent(inout) :: sfc_z0h   (ia,ja)    ! surface roughness length (heat) [m]
    real(RP), intent(inout) :: sfc_z0e   (ia,ja)    ! surface roughness length (vapor) [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_qtrc (ia,ja,qa) ! surface flux for tracer mass   (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), intent(out)   :: t2        (ia,ja)    ! temperature t     at  2m height
    real(RP), intent(out)   :: q2        (ia,ja)    ! water vapor q     at  2m height

    real(RP) :: atm_uabs(ia,ja) ! absolute velocity at z1 [m/s]

    real(RP) :: cm(ia,ja)       ! bulk coefficient for momentum
    real(RP) :: r10

    real(RP) :: modulation
    real(RP) :: uabs_lim
    real(RP) :: lhv(ia,ja)
    integer  :: i, j
    !---------------------------------------------------------------------------

    if( io_l ) write(io_fid_log,*) '*** Physics step: Surface flux(const)'

    do j = js, je
    do i = is, ie
       atm_uabs(i,j) = sqrt( atm_w(i,j)*atm_w(i,j) &
                           + atm_u(i,j)*atm_u(i,j) &
                           + atm_v(i,j)*atm_v(i,j) ) ! at cell center
    enddo
    enddo

    if   ( atmos_phy_sf_flg_mom_flux == 0 ) then ! Bulk coefficient is constant
       do j = js, je
       do i = is, ie
          cm(i,j) = atmos_phy_sf_const_cm
       enddo
       enddo
    elseif( atmos_phy_sf_flg_mom_flux == 1 ) then ! Friction velocity is constant
       do j = js, je
       do i = is, ie
          cm(i,j) = ( atmos_phy_sf_const_ustar / atm_uabs(i,j) )**2
          cm(i,j) = min( max( cm(i,j), atmos_phy_sf_cm_min ), atmos_phy_sf_cm_max )
       enddo
       enddo
    endif

    !-----< momentum >-----

    do j = js, je
    do i = is, ie
       uabs_lim = min( max( atm_uabs(i,j), atmos_phy_sf_u_minm ), atmos_phy_sf_u_maxm )

       sflx_mw(i,j) = -cm(i,j) * uabs_lim * sfc_dens(i,j) * atm_w(i,j)
       sflx_mu(i,j) = -cm(i,j) * uabs_lim * sfc_dens(i,j) * atm_u(i,j)
       sflx_mv(i,j) = -cm(i,j) * uabs_lim * sfc_dens(i,j) * atm_v(i,j)
    enddo
    enddo

    !-----< heat flux >-----

    if ( atmos_phy_sf_flg_sh_diurnal ) then
       modulation = sin( 2.0_rp * pi * time_nowsec / 3600.0_rp / atmos_phy_sf_const_freq )
    else
       modulation = 1.0_rp
    endif

    do j = js, je
    do i = is, ie
       sflx_sh(i,j) = atmos_phy_sf_const_sh * modulation
       sflx_lh(i,j) = atmos_phy_sf_const_lh
    enddo
    enddo

    !-----< mass flux >-----
   call atmos_thermodyn_templhv( lhv, atm_temp )

    sflx_qtrc(:,:,:) = 0.0_rp
    do j = js, je
    do i = is, ie
       sflx_qtrc(i,j,i_qv) = sflx_lh(i,j) / lhv(i,j)
    enddo
    enddo

    !-----< U10, T2, q2 >-----

    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

    do j = js, je
    do i = is, ie
       t2(i,j) = atm_temp(i,j)
       q2(i,j) = atm_qtrc(i,j,i_qv)
    enddo
    enddo

    return

Here is the caller graph for this function: