scale_atmos_bottom Module Reference

module atmosphere / bottom boundary extrapolation More...

Functions/Subroutines
subroutine, public	atmos_bottom_estimate (KA, KS, KE, IA, IS, IE, JA, JS, JE, DENS, PRES, QV, SFC_TEMP, FZ, SFC_DENS, SFC_PRES)
	Calc bottom boundary of atmosphere (just above surface) More...
real(rp) function	lagrange_interp (p, x, y)

Detailed Description

module atmosphere / bottom boundary extrapolation

Description

Bottom boundary treatment of model domain Extrapolation of DENS, PRES, TEMP

Author

Team SCALE

Function/Subroutine Documentation

atmos_bottom_estimate()

subroutine, public scale_atmos_bottom::atmos_bottom_estimate	(	integer, intent(in)	KA,
		integer, intent(in)	KS,
		integer, intent(in)	KE,
		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 (ka,ia,ja), intent(in)	DENS,
		real(rp), dimension (ka,ia,ja), intent(in)	PRES,
		real(rp), dimension (ka,ia,ja), intent(in)	QV,
		real(rp), dimension(ia,ja), intent(in)	SFC_TEMP,
		real(rp), dimension (0:ka,ia,ja), intent(in)	FZ,
		real(rp), dimension(ia,ja), intent(out)	SFC_DENS,
		real(rp), dimension(ia,ja), intent(out)	SFC_PRES
	)

Calc bottom boundary of atmosphere (just above surface)

Definition at line 51 of file scale_atmos_bottom.F90.

    use scale_const, only: &
       grav    => const_grav, &
       rdry    => const_rdry, &
       epstvap => const_epstvap
    implicit none
    integer, intent(in) :: KA, KS, KE
    integer, intent(in) :: IA, IS, IE
    integer, intent(in) :: JA, JS, JE
 
    real(RP), intent(in)  :: DENS    (KA,IA,JA)
    real(RP), intent(in)  :: PRES    (KA,IA,JA)
    real(RP), intent(in)  :: QV      (KA,IA,JA)
    real(RP), intent(in)  :: SFC_TEMP(IA,JA)
    real(RP), intent(in)  :: FZ      (0:KA,IA,JA)
    real(RP), intent(out) :: SFC_DENS(IA,JA)
    real(RP), intent(out) :: SFC_PRES(IA,JA)
 
    real(RP) :: Rtot
    real(RP) :: dz1, dz2
    real(RP) :: F2H1, F2H2
    real(RP) :: LogP1, LogP2
    real(RP) :: PRESH
 
    integer :: i, j
    !---------------------------------------------------------------------------
 
    !$omp parallel do OMP_SCHEDULE_ &
    !$omp private (Rtot,dz1,dz2,F2H1,F2H2,LogP1,LogP2,PRESH)
    do j = js, je
    do i = is, ie
       ! interpolate half-level pressure
       dz1 = fz(ks+1,i,j) - fz(ks,i,j)
       dz2 = fz(ks,i,j) - fz(ks-1,i,j)
 
       f2h1 = dz2 / ( dz1 + dz2 )
       f2h2 = dz1 / ( dz1 + dz2 )
 
       logp1 = log( pres(ks+1,i,j) )
       logp2 = log( pres(ks,i,j) )
 
       presh = exp( f2h1 * logp1 + f2h2 * logp2 )
 
       rtot = rdry * ( 1.0_rp + epstvap * qv(ks,i,j) )
       ! ( PRESH - SFC_PRES ) / dz2 = - GRAV * DENS(KS)
       sfc_pres(i,j) = presh + grav * dens(ks,i,j) * dz2
       sfc_dens(i,j) = sfc_pres(i,j) / ( rtot * sfc_temp(i,j) )
    enddo
    enddo
 
    return

References scale_const::const_epstvap, scale_const::const_grav, and scale_const::const_rdry.

Referenced by mod_atmos_phy_sf_driver::atmos_phy_sf_driver_calc_tendency(), mod_atmos_driver::atmos_surface_set(), and mod_atmos_vars::atmos_vars_history_setpres().

Here is the caller graph for this function:

lagrange_interp()

real(rp) function scale_atmos_bottom::lagrange_interp	(	real(rp), intent(in)	p,
		real(rp), dimension(3), intent(in)	x,
		real(rp), dimension(3), intent(in)	y
	)

Definition at line 105 of file scale_atmos_bottom.F90.

    implicit none
 
    real(RP), intent(in) :: p
    real(RP), intent(in) :: x(3), y(3)
    real(RP)             :: q
    !---------------------------------------------------------------------------
 
    q = ( (p-x(2)) * (p-x(3)) ) / ( (x(1)-x(2)) * (x(1)-x(3)) ) * y(1) &
      + ( (p-x(1)) * (p-x(3)) ) / ( (x(2)-x(1)) * (x(2)-x(3)) ) * y(2) &
      + ( (p-x(1)) * (p-x(2)) ) / ( (x(3)-x(1)) * (x(3)-x(2)) ) * y(3)
 
    return