scale_land_sfc_thick_slab Module Reference

module LAND / Surface fluxes with thick-slab land model More...

Functions/Subroutines
subroutine, public	land_sfc_thick_slab_setup (LAND_TYPE)
	Setup. More...
subroutine, public	land_sfc_thick_slab (LST_t, ZMFLX, XMFLX, YMFLX, SHFLX, LHFLX, GHFLX, U10, V10, T2, Q2, TMPA, PRSA, WA, UA, VA, RHOA, QVA, Z1, PBL, PRSS, LWD, SWD, TG, LST, QVEF, ALB_LW, ALB_SW, DZG, TCS, Z0M, Z0H, Z0E, dt)

Detailed Description

module LAND / Surface fluxes with thick-slab land model

Description

Surface flux with thick-slab land model

Author

Team SCALE

NAMELIST

• PARAM_LAND_SFC_THICK_SLAB

  name	type	default value	comment
  DUMMY	logical

  
  

History Output

No history output

Function/Subroutine Documentation

land_sfc_thick_slab_setup()

subroutine, public scale_land_sfc_thick_slab::land_sfc_thick_slab_setup

(

character(len=*), intent(in)

LAND_TYPE

)

Setup.

Definition at line 46 of file scale_land_sfc_thick-slab.F90.

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

Referenced by scale_land_sfc::land_sfc_setup().

    use scale_process, only: &
       prc_mpistop
    implicit none

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

    logical :: dummy

    namelist / param_land_sfc_thick_slab / &
       dummy

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

    if( io_l ) write(io_fid_log,*)
    if( io_l ) write(io_fid_log,*) '++++++ Module[THICK-SLAB] / Categ[LAND SFC] / Origin[SCALElib]'

    !--- read namelist
    rewind(io_fid_conf)
    read(io_fid_conf,nml=param_land_sfc_thick_slab,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_LAND_SFC_THICK_SLAB. Check!'
       call prc_mpistop
    endif
    if( io_nml ) write(io_fid_nml,nml=param_land_sfc_thick_slab)

    return

Here is the call graph for this function:

Here is the caller graph for this function:

land_sfc_thick_slab()

subroutine, public scale_land_sfc_thick_slab::land_sfc_thick_slab	(	real(rp), dimension(ia,ja), intent(out)	LST_t,
		real(rp), dimension(ia,ja), intent(out)	ZMFLX,
		real(rp), dimension(ia,ja), intent(out)	XMFLX,
		real(rp), dimension(ia,ja), intent(out)	YMFLX,
		real(rp), dimension(ia,ja), intent(out)	SHFLX,
		real(rp), dimension(ia,ja), intent(out)	LHFLX,
		real(rp), dimension(ia,ja), intent(out)	GHFLX,
		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,
		real(rp), dimension(ia,ja), intent(in)	TMPA,
		real(rp), dimension(ia,ja), intent(in)	PRSA,
		real(rp), dimension (ia,ja), intent(in)	WA,
		real(rp), dimension (ia,ja), intent(in)	UA,
		real(rp), dimension (ia,ja), intent(in)	VA,
		real(rp), dimension(ia,ja), intent(in)	RHOA,
		real(rp), dimension (ia,ja), intent(in)	QVA,
		real(rp), dimension (ia,ja), intent(in)	Z1,
		real(rp), dimension (ia,ja), intent(in)	PBL,
		real(rp), dimension(ia,ja), intent(in)	PRSS,
		real(rp), dimension (ia,ja), intent(in)	LWD,
		real(rp), dimension (ia,ja), intent(in)	SWD,
		real(rp), dimension (ia,ja), intent(in)	TG,
		real(rp), dimension (ia,ja), intent(in)	LST,
		real(rp), dimension (ia,ja), intent(in)	QVEF,
		real(rp), dimension(ia,ja), intent(in)	ALB_LW,
		real(rp), dimension(ia,ja), intent(in)	ALB_SW,
		real(rp), dimension (ia,ja), intent(in)	DZG,
		real(rp), dimension (ia,ja), intent(in)	TCS,
		real(rp), dimension (ia,ja), intent(in)	Z0M,
		real(rp), dimension (ia,ja), intent(in)	Z0H,
		real(rp), dimension (ia,ja), intent(in)	Z0E,
		real(dp), intent(in)	dt
	)

Definition at line 113 of file scale_land_sfc_thick-slab.F90.

References scale_bulkflux::bulkflux, scale_const::const_cpdry, scale_const::const_pre00, scale_const::const_rdry, scale_const::const_stb, 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, scale_landuse::landuse_fact_land, and scale_process::prc_mpistop().

Referenced by scale_land_sfc::land_sfc_setup().

    use scale_process, only: &
      prc_mpistop
    use scale_const, only: &
      pre00 => const_pre00, &
      rdry  => const_rdry,  &
      cpdry => const_cpdry, &
      stb   => const_stb
    use scale_landuse, only: &
      landuse_fact_land
    use scale_atmos_hydrometeor, only: &
      hydrometeor_lhv => atmos_hydrometeor_lhv
    use scale_atmos_saturation, only: &
      qsat => atmos_saturation_pres2qsat_all
    use scale_bulkflux, only: &
      bulkflux
    implicit none

    ! arguments
    real(RP), intent(out) :: LST_t(IA,JA) ! tendency of LST
    real(RP), intent(out) :: ZMFLX(IA,JA) ! z-momentum flux at the surface [kg/m/s2]
    real(RP), intent(out) :: XMFLX(IA,JA) ! x-momentum flux at the surface [kg/m/s2]
    real(RP), intent(out) :: YMFLX(IA,JA) ! y-momentum flux at the surface [kg/m/s2]
    real(RP), intent(out) :: SHFLX(IA,JA) ! sensible heat flux at the surface [J/m2/s]
    real(RP), intent(out) :: LHFLX(IA,JA) ! latent heat flux at the surface [J/m2/s]
    real(RP), intent(out) :: GHFLX(IA,JA) ! ground heat flux at the surface [J/m2/s]
    real(RP), intent(out) :: U10  (IA,JA) ! velocity u at 10m [m/s]
    real(RP), intent(out) :: V10  (IA,JA) ! velocity v at 10m [m/s]
    real(RP), intent(out) :: T2   (IA,JA) ! temperature at 2m [K]
    real(RP), intent(out) :: Q2   (IA,JA) ! water vapor at 2m [kg/kg]

    real(RP), intent(in) :: TMPA(IA,JA) ! temperature at the lowest atmospheric layer [K]
    real(RP), intent(in) :: PRSA(IA,JA) ! pressure at the lowest atmospheric layer [Pa]
    real(RP), intent(in) :: WA  (IA,JA) ! velocity w at the lowest atmospheric layer [m/s]
    real(RP), intent(in) :: UA  (IA,JA) ! velocity u at the lowest atmospheric layer [m/s]
    real(RP), intent(in) :: VA  (IA,JA) ! velocity v at the lowest atmospheric layer [m/s]
    real(RP), intent(in) :: RHOA(IA,JA) ! density at the lowest atmospheric layer [kg/m3]
    real(RP), intent(in) :: QVA (IA,JA) ! ratio of water vapor mass to total mass at the lowest atmospheric layer [kg/kg]
    real(RP), intent(in) :: Z1  (IA,JA) ! cell center height at the lowest atmospheric layer [m]
    real(RP), intent(in) :: PBL (IA,JA) ! the top of atmospheric mixing layer [m]
    real(RP), intent(in) :: PRSS(IA,JA) ! pressure at the surface [Pa]
    real(RP), intent(in) :: LWD (IA,JA) ! downward long-wave radiation flux at the surface [J/m2/s]
    real(RP), intent(in) :: SWD (IA,JA) ! downward short-wave radiation flux at the surface [J/m2/s]

    real(RP), intent(in) :: TG    (IA,JA) ! soil temperature [K]
    real(RP), intent(in) :: LST   (IA,JA) ! land surface temperature [K]
    real(RP), intent(in) :: QVEF  (IA,JA) ! efficiency of evaporation (0-1)
    real(RP), intent(in) :: ALB_LW(IA,JA) ! surface albedo for LW (0-1)
    real(RP), intent(in) :: ALB_SW(IA,JA) ! surface albedo for SW (0-1)
    real(RP), intent(in) :: DZG   (IA,JA) ! soil depth [m]
    real(RP), intent(in) :: TCS   (IA,JA) ! thermal conductivity for soil [J/m/K/s]
    real(RP), intent(in) :: Z0M   (IA,JA) ! roughness length for momemtum [m]
    real(RP), intent(in) :: Z0H   (IA,JA) ! roughness length for heat [m]
    real(RP), intent(in) :: Z0E   (IA,JA) ! roughness length for vapor [m]
    real(DP), intent(in) :: dt            ! delta time

    ! works
    real(RP) :: LST1(IA,JA)

    real(RP) :: res    ! residual

    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) :: QVsat  ! saturation water vapor mixing ratio at surface [kg/kg]
    real(RP) :: QVS    ! water vapor mixing ratio at surface [kg/kg]

    real(RP) :: FracU10 ! calculation parameter for U10 [-]
    real(RP) :: FracT2  ! calculation parameter for T2 [-]
    real(RP) :: FracQ2  ! calculation parameter for Q2 [-]

    real(RP) :: LHV(IA,JA)    ! latent heat of vaporization [J/kg]

    integer  :: i, j
    !---------------------------------------------------------------------------

    if( io_l ) write(io_fid_log,*) '*** Land surface step: Thick-Slab'

    call hydrometeor_lhv( lhv(:,:), tmpa(:,:) )

    ! update surface temperature
    do j = js, je
    do i = is, ie
      lst1(i,j) = tg(i,j)
      lst_t(i,j) = ( lst1(i,j) - lst(i,j) ) / dt
    end do
    end do

    ! calculate surface flux
    do j = js, je
    do i = is, ie

      if( landuse_fact_land(i,j) > 0.0_rp ) then

        call qsat( qvsat,     & ! [OUT]
                   lst1(i,j), & ! [IN]
                   prss(i,j)  ) ! [IN]

        qvs  = ( 1.0_rp - qvef(i,j) ) * qva(i,j) + qvef(i,j) * qvsat

        call bulkflux( &
            ustar,     & ! [OUT]
            tstar,     & ! [OUT]
            qstar,     & ! [OUT]
            uabs,      & ! [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) = -rhoa(i,j) * ustar**2 / uabs * wa(i,j)
        xmflx(i,j) = -rhoa(i,j) * ustar**2 / uabs * ua(i,j)
        ymflx(i,j) = -rhoa(i,j) * ustar**2 / uabs * va(i,j)

        shflx(i,j) = - rhoa(i,j) * ustar * tstar &
                   * cpdry * ( prss(i,j) / pre00 )**( rdry/cpdry )
        lhflx(i,j) = - rhoa(i,j) * ustar * qstar * lhv(i,j)

        ghflx(i,j) = -2.0_rp * tcs(i,j) * ( lst1(i,j) - tg(i,j) ) / dzg(i,j)

        ! calculation for residual
        res = ( 1.0_rp - alb_sw(i,j) ) * swd(i,j) &
            + ( 1.0_rp - alb_lw(i,j) ) * ( lwd(i,j) - stb * lst1(i,j)**4 ) &
            - shflx(i,j) - lhflx(i,j) + ghflx(i,j)

        ! put residual in ground heat flux
        ghflx(i,j) = ghflx(i,j) - res

        ! diagnostic variables considering unstable/stable state
        !U10(i,j) = FracU10 * UA(i,j)
        !V10(i,j) = FracU10 * VA(i,j)
        !T2 (i,j) = ( 1.0_RP - FracT2 ) * LST1(i,j) + FracT2 * TMPA(i,j)
        !Q2 (i,j) = ( 1.0_RP - FracQ2 ) * QVS       + FracQ2 * QVA (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
        shflx(i,j) = 0.0_rp
        lhflx(i,j) = 0.0_rp
        ghflx(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 call graph for this function:

Here is the caller graph for this function: