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, CZ, Zsfc, Z1, 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)	CZ,
		real(rp), dimension (ia,ja), intent(in)	Zsfc,
		real(rp), dimension (ia,ja), intent(in)	Z1,
		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 50 of file scale_atmos_bottom.F90.

References scale_const::const_grav, and lagrange_interp().

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().

    use scale_const, only: &
       grav  => const_grav
    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)  :: cz      (ka,ia,ja)
    real(RP), intent(in)  :: zsfc    (ia,ja)
    real(RP), intent(in)  :: z1      (ia,ja)
    real(RP), intent(out) :: sfc_dens(ia,ja)
    real(RP), intent(out) :: sfc_pres(ia,ja)

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

    ! estimate surface density (extrapolation)
    !$omp parallel do OMP_SCHEDULE_
    do j = js, je
    do i = is, ie
       sfc_dens(i,j) = lagrange_interp( zsfc(i,j),         & ! [IN]
                                        cz(ks:ks+2,i,j), & ! [IN]
                                        dens(ks:ks+2,i,j)  ) ! [IN]
    enddo
    enddo

    ! estimate surface pressure (hydrostatic balance)
    !$omp parallel do OMP_SCHEDULE_
    do j = js, je
    do i = is, ie
       sfc_pres(i,j) = pres(ks,i,j) &
                     + 0.5_rp * ( sfc_dens(i,j) + dens(ks,i,j) ) * grav * z1(i,j)
    enddo
    enddo

    return

Here is the call graph for this function:

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 92 of file scale_atmos_bottom.F90.

Referenced by atmos_bottom_estimate().

    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

Here is the caller graph for this function: