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, Ustar, Tstar, Qstar, Wstar, RLmo, 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)	Ustar,
		real(rp), dimension(lia,lja), intent(out)	Tstar,
		real(rp), dimension(lia,lja), intent(out)	Qstar,
		real(rp), dimension(lia,lja), intent(out)	Wstar,
		real(rp), dimension (lia,lja), intent(out)	RLmo,
		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 59 of file scale_land_phy_snow_diagnos.F90.

     use scale_const, only: &
       eps   => const_eps,   &
       undef => const_undef, &
       pre00 => const_pre00, &
       rdry  => const_rdry,  &
       cpdry => const_cpdry, &
       rvap  => const_rvap,  &
       stb   => const_stb
     use scale_atmos_saturation, only: &
       qsat => atmos_saturation_dens2qsat_all
 !      qsat => ATMOS_SATURATION_pres2qsat_all
     use scale_bulkflux, only: &
       bulkflux, &
       bulkflux_diagnose_surface
     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) :: Ustar(LIA,LJA) ! friction velocity         [m/s]
     real(RP), intent(out) :: Tstar(LIA,LJA) ! temperature scale         [K]
     real(RP), intent(out) :: Qstar(LIA,LJA) ! moisture scale            [kg/kg]
     real(RP), intent(out) :: Wstar(LIA,LJA) ! convective velocity scale [m/s]
     real(RP), intent(out) :: RLmo (LIA,LJA) ! inversed Obukhov length   [1/m]
     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) :: 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(LIA,LJA) ! 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(LIA,LJA) ! calculation parameter for U10 [-]
     real(RP) :: FracT2 (LIA,LJA) ! calculation parameter for T2 [-]
     real(RP) :: FracQ2 (LIA,LJA) ! calculation parameter for Q2 [-]
  
     real(RP) :: MFLUX
  
     logical :: calc_flag(LIA,LJA)
  
     integer  :: i, j
     !---------------------------------------------------------------------------
  
     log_info("LAND_PHY_SNOW_DIAGS",*) 'Snow surface diagnostic'
  
     ! calculate surface flux
     !$omp parallel do &
     !$omp private(QVsat,Uabs,MFLUX)
     do j = ljs, lje
     do i = lis, lie
  
       if( snow_frac(i,j) > 0.0_rp ) then
  
          calc_flag(i,j) = .true.
  
 !        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]
         call qsat( lst1(i,j),  rhos(i,j), & ! [IN]
                    qvsat                  ) ! [OUT]
  
         qvs(i,j)  = ( 1.0_rp - qvef(i,j) ) * qva(i,j) + qvef(i,j) * qvsat
  
         uabs = sqrt( wa(i,j)**2 + ua(i,j)**2 + va(i,j)**2 )
  
         call bulkflux( tmpa(i,j), lst1(i,j),                  & ! [IN]
                        prsa(i,j), prss(i,j),                  & ! [IN]
                        qva(i,j), qvs(i,j),                  & ! [IN]
                        uabs, z1(i,j), pbl(i,j),               & ! [IN]
                        z0m(i,j), z0h(i,j), z0e(i,j),          & ! [IN]
                        ustar(i,j), tstar(i,j), qstar(i,j),    & ! [OUT]
                        wstar(i,j), rlmo(i,j), ra,             & ! [OUT]
                        fracu10(i,j), fract2(i,j), fracq2(i,j) ) ! [OUT]
  
         if ( uabs < eps ) then
            zmflx(i,j) = 0.0_rp
            xmflx(i,j) = 0.0_rp
            ymflx(i,j) = 0.0_rp
         else
            mflux = - rhos(i,j) * ustar(i,j)**2
            zmflx(i,j) = mflux * wa(i,j) / uabs
            xmflx(i,j) = mflux * ua(i,j) / uabs
            ymflx(i,j) = mflux * va(i,j) / uabs
         end if
  
       else
  
          ! not calculate surface flux
  
          calc_flag(i,j) = .false.
  
          zmflx(i,j) = undef
          xmflx(i,j) = undef
          ymflx(i,j) = undef
          ustar(i,j) = undef
          tstar(i,j) = undef
          qstar(i,j) = undef
          wstar(i,j) = undef
          rlmo(i,j) = undef
          u10(i,j) = undef
          v10(i,j) = undef
          t2(i,j) = undef
          q2(i,j) = undef
  
       end if
  
     end do
     end do
  
     call bulkflux_diagnose_surface( lia, lis, lie, lja, ljs, lje, &
                                     ua(:,:), va(:,:),                      & ! (in)
                                     tmpa(:,:), qva(:,:),                   & ! (in)
                                     lst1(:,:), qvs(:,:),                   & ! (in)
                                     z1(:,:), z0m(:,:), z0h(:,:), z0e(:,:), & ! (in)
                                     u10(:,:), v10(:,:), t2(:,:), q2(:,:),  & ! (out)
                                     mask = calc_flag(:,:),                 & ! (in)
                                     fracu10 = fracu10(:,:),                & ! (in)
                                     fract2 = fract2(:,:),                  & ! (in)
                                     fracq2 = fracq2(:,:)                   ) ! (in)
  
     return

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

Referenced by mod_land_driver::land_driver_calc_tendency().

Here is the caller graph for this function:

Functions/Subroutines

Detailed Description

Function/Subroutine Documentation

◆ land_phy_snow_diags()