module land / physics / snow / diagnostics More...

Functions/Subroutines
subroutine, public	land_phy_snow_diags (LIA, LIS, LIE, LJA, LJS, LJE, SNOW_frac, TMPA, PRSA, WA, UA, VA, RHOA, QVA, Z1, PBL, RHOS, PRSS, LST1, QVEF, Z0M, Z0H, Z0E, ZMFLX, XMFLX, YMFLX, U10, V10, T2, Q2)

Detailed Description

module land / physics / snow / diagnostics

Description: Momentum fluxes and Diagnostic variables

Author: Team SCALE

Function/Subroutine Documentation

◆ land_phy_snow_diags()

subroutine, public scale_land_phy_snow_diagnos::land_phy_snow_diags	(	integer, intent(in)	LIA,
		integer, intent(in)	LIS,
		integer, intent(in)	LIE,
		integer, intent(in)	LJA,
		integer, intent(in)	LJS,
		integer, intent(in)	LJE,
		real(rp), dimension(lia,lja), intent(in)	SNOW_frac,
		real(rp), dimension(lia,lja), intent(in)	TMPA,
		real(rp), dimension(lia,lja), intent(in)	PRSA,
		real(rp), dimension (lia,lja), intent(in)	WA,
		real(rp), dimension (lia,lja), intent(in)	UA,
		real(rp), dimension (lia,lja), intent(in)	VA,
		real(rp), dimension(lia,lja), intent(in)	RHOA,
		real(rp), dimension (lia,lja), intent(in)	QVA,
		real(rp), dimension (lia,lja), intent(in)	Z1,
		real(rp), dimension (lia,lja), intent(in)	PBL,
		real(rp), dimension(lia,lja), intent(in)	RHOS,
		real(rp), dimension(lia,lja), intent(in)	PRSS,
		real(rp), dimension (lia,lja), intent(in)	LST1,
		real(rp), dimension (lia,lja), intent(in)	QVEF,
		real(rp), dimension (lia,lja), intent(in)	Z0M,
		real(rp), dimension (lia,lja), intent(in)	Z0H,
		real(rp), dimension (lia,lja), intent(in)	Z0E,
		real(rp), dimension(lia,lja), intent(out)	ZMFLX,
		real(rp), dimension(lia,lja), intent(out)	XMFLX,
		real(rp), dimension(lia,lja), intent(out)	YMFLX,
		real(rp), dimension (lia,lja), intent(out)	U10,
		real(rp), dimension (lia,lja), intent(out)	V10,
		real(rp), dimension (lia,lja), intent(out)	T2,
		real(rp), dimension (lia,lja), intent(out)	Q2
	)

Definition at line 56 of file scale_land_phy_snow_diagnos.F90.

References scale_bulkflux::bulkflux, scale_const::const_cpdry, scale_const::const_pre00, scale_const::const_rdry, scale_const::const_rvap, and scale_const::const_stb.

Referenced by mod_land_driver::land_driver_calc_tendency().

     use scale_const, only: &
       pre00 => const_pre00, &
       rdry  => const_rdry,  &
       cpdry => const_cpdry, &
       rvap  => const_rvap,  &
       stb   => const_stb
     use scale_atmos_saturation, only: &
       qsat => atmos_saturation_pres2qsat_all
     use scale_bulkflux, only: &
       bulkflux
     implicit none
     integer, intent(in) :: lia, lis, lie
     integer, intent(in) :: lja, ljs, lje
 
     real(RP), intent(in) :: snow_frac(lia,lja) ! snow fraction
     real(RP), intent(in) :: tmpa(lia,lja) ! temperature at the lowest atmospheric layer [K]
     real(RP), intent(in) :: prsa(lia,lja) ! pressure at the lowest atmospheric layer [Pa]
     real(RP), intent(in) :: wa  (lia,lja) ! velocity w at the lowest atmospheric layer [m/s]
     real(RP), intent(in) :: ua  (lia,lja) ! velocity u at the lowest atmospheric layer [m/s]
     real(RP), intent(in) :: va  (lia,lja) ! velocity v at the lowest atmospheric layer [m/s]
     real(RP), intent(in) :: rhoa(lia,lja) ! density at the lowest atmospheric layer [kg/m3]
     real(RP), intent(in) :: qva (lia,lja) ! ratio of water vapor mass to total mass at the lowest atmospheric layer [kg/kg]
     real(RP), intent(in) :: z1  (lia,lja) ! cell center height at the lowest atmospheric layer [m]
     real(RP), intent(in) :: pbl (lia,lja) ! the top of atmospheric mixing layer [m]
     real(RP), intent(in) :: rhos(lia,lja) ! density  at the surface [kg/m3]
     real(RP), intent(in) :: prss(lia,lja) ! pressure at the surface [Pa]
 
     real(RP), intent(in) :: lst1   (lia,lja) ! land surface temperature [K]
     real(RP), intent(in) :: qvef  (lia,lja) ! efficiency of evaporation (0-1)
 
     real(RP), intent(in) :: z0m   (lia,lja) ! roughness length for momemtum [m]
     real(RP), intent(in) :: z0h   (lia,lja) ! roughness length for heat [m]
     real(RP), intent(in) :: z0e   (lia,lja) ! roughness length for vapor [m]
 
     real(RP), intent(out) :: zmflx(lia,lja) ! z-momentum flux at the surface [kg/m/s2]
     real(RP), intent(out) :: xmflx(lia,lja) ! x-momentum flux at the surface [kg/m/s2]
     real(RP), intent(out) :: ymflx(lia,lja) ! y-momentum flux at the surface [kg/m/s2]
     real(RP), intent(out) :: u10  (lia,lja) ! velocity u at 10m [m/s]
     real(RP), intent(out) :: v10  (lia,lja) ! velocity v at 10m [m/s]
     real(RP), intent(out) :: t2   (lia,lja) ! temperature at 2m [K]
     real(RP), intent(out) :: q2   (lia,lja) ! water vapor at 2m [kg/kg]
 
 
     ! works
     real(RP) :: ustar  ! friction velocity [m]
     real(RP) :: tstar  ! friction potential temperature [K]
     real(RP) :: qstar  ! friction water vapor mass ratio [kg/kg]
     real(RP) :: uabs   ! modified absolute velocity [m/s]
     real(RP) :: ra     ! Aerodynamic resistance (=1/Ce) [1/s]
 
     real(RP) :: qvsat  ! saturation water vapor mixing ratio at surface [kg/kg]
     real(RP) :: qvs    ! water vapor mixing ratio at surface [kg/kg]
     real(RP) :: rtot   ! total gas constant
     real(RP) :: qdry   ! dry air mass ratio [kg/kg]
 
     real(RP) :: fracu10 ! calculation parameter for U10 [-]
     real(RP) :: fract2  ! calculation parameter for T2 [-]
     real(RP) :: fracq2  ! calculation parameter for Q2 [-]
 
     integer  :: i, j
     !---------------------------------------------------------------------------
 
     log_info("LAND_PHY_SNOW_DIAGS",*) 'Snow surface diagnostic'
 
     ! calculate surface flux
     do j = ljs, lje
     do i = lis, lie
 
       if( snow_frac(i,j) > 0.0_rp ) then
 
         qdry = 1.0_rp - qva(i,j)
         rtot = qdry * rdry + qva(i,j) * rvap
 
         call qsat( lst1(i,j),  prss(i,j), qdry, & ! [IN]
                    qvsat                        ) ! [OUT]
 
         qvs  = ( 1.0_rp - qvef(i,j) ) * qva(i,j) + qvef(i,j) * qvsat
 
         call bulkflux( &
             ustar,     & ! [OUT]
             tstar,     & ! [OUT]
             qstar,     & ! [OUT]
             uabs,      & ! [OUT]
             ra,        & ! [OUT]
             fracu10,   & ! [OUT]
             fract2,    & ! [OUT]
             fracq2,    & ! [OUT]
             tmpa(i,j), & ! [IN]
             lst1(i,j), & ! [IN]
             prsa(i,j), & ! [IN]
             prss(i,j), & ! [IN]
             qva(i,j), & ! [IN]
             qvs,       & ! [IN]
             ua(i,j), & ! [IN]
             va(i,j), & ! [IN]
             z1(i,j), & ! [IN]
             pbl(i,j), & ! [IN]
             z0m(i,j), & ! [IN]
             z0h(i,j), & ! [IN]
             z0e(i,j)  ) ! [IN]
 
         zmflx(i,j) = -rhos(i,j) * ustar * ustar / uabs * wa(i,j)
         xmflx(i,j) = -rhos(i,j) * ustar * ustar / uabs * ua(i,j)
         ymflx(i,j) = -rhos(i,j) * ustar * ustar / uabs * va(i,j)
 
         ! diagnostic variables for neutral state
         u10(i,j) = ua(i,j) * log( 10.0_rp / z0m(i,j) ) / log( z1(i,j) / z0m(i,j) )
         v10(i,j) = va(i,j) * log( 10.0_rp / z0m(i,j) ) / log( z1(i,j) / z0m(i,j) )
         t2(i,j) = lst1(i,j) + ( tmpa(i,j) - lst1(i,j) ) * ( log(  2.0_rp / z0m(i,j) ) * log(  2.0_rp / z0h(i,j) ) ) &
                                                          / ( log( z1(i,j) / z0m(i,j) ) * log( z1(i,j) / z0h(i,j) ) )
         q2(i,j) = qvs       + (  qva(i,j) - qvs       ) * ( log(  2.0_rp / z0m(i,j) ) * log(  2.0_rp / z0e(i,j) ) ) &
                                                          / ( log( z1(i,j) / z0m(i,j) ) * log( z1(i,j) / z0e(i,j) ) )
 
       else
 
         ! not calculate surface flux
         zmflx(i,j) = 0.0_rp
         xmflx(i,j) = 0.0_rp
         ymflx(i,j) = 0.0_rp
         u10(i,j) = 0.0_rp
         v10(i,j) = 0.0_rp
         t2(i,j) = 0.0_rp
         q2(i,j) = 0.0_rp
 
       end if
 
     end do
     end do
 
     return

Here is the caller graph for this function:

Functions/Subroutines

Detailed Description

Function/Subroutine Documentation

◆ land_phy_snow_diags()