scale_atmos_phy_rd_common Module Reference

module atmosphere / physics / radiation / common More...

Functions/Subroutines
subroutine, public	atmos_phy_rd_calc_heating (KA, KS, KE, IA, IS, IE, JA, JS, JE, flux_rad, DENS, TEMP, CVtot, FZ, RHOH, TEMP_t)
	Calc heating rate. More...

Variables
integer, parameter, public	i_up = 1
integer, parameter, public	i_dn = 2
integer, parameter, public	i_lw = 1
integer, parameter, public	i_sw = 2

Detailed Description

module atmosphere / physics / radiation / common

Description

Common module for Radiation

Author

Team SCALE

Function/Subroutine Documentation

atmos_phy_rd_calc_heating()

subroutine, public scale_atmos_phy_rd_common::atmos_phy_rd_calc_heating	(	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,2,2), intent(in)	flux_rad,
		real(rp), dimension (ka,ia,ja), intent(in)	DENS,
		real(rp), dimension (ka,ia,ja), intent(in)	TEMP,
		real(rp), dimension (ka,ia,ja), intent(in)	CVtot,
		real(rp), dimension (0:ka,ia,ja), intent(in)	FZ,
		real(rp), dimension(ka,ia,ja), intent(out)	RHOH,
		real(rp), dimension(ka,ia,ja,3), intent(out), optional	TEMP_t
	)

Calc heating rate.

Definition at line 60 of file scale_atmos_phy_rd_common.F90.

    implicit none
    integer, intent(in) :: KA, KS, KE
    integer, intent(in) :: IA, IS, IE
    integer, intent(in) :: JA, JS, JE
 
    real(RP), intent(in)  :: flux_rad(KA,IA,JA,2,2)
    real(RP), intent(in)  :: DENS    (KA,IA,JA)
    real(RP), intent(in)  :: TEMP    (KA,IA,JA)
    real(RP), intent(in)  :: CVtot   (KA,IA,JA)
    real(RP), intent(in)  :: FZ      (0:KA,IA,JA)
 
    real(RP), intent(out) :: RHOH(KA,IA,JA)
 
    real(RP), intent(out), optional :: TEMP_t(KA,IA,JA,3)
 
    real(RP) :: RHOH_LW, RHOH_SW
 
    integer :: k, i, j
    !---------------------------------------------------------------------------
 
    if ( present(temp_t) ) then
 
       !$omp parallel do default(none) OMP_SCHEDULE_ collapse(2) &
       !$omp private(i,j,k, &
       !$omp         RHOH_LW,RHOH_SW) &
       !$omp shared(RHOH,TEMP_t,flux_rad,DENS,CVtot,FZ, &
       !$omp        KS,KE,IS,IE,JS,JE)
       !$acc kernels
       do j = js, je
       do i = is, ie
       do k = ks, ke
 
          rhoh_lw = ( ( flux_rad(k,i,j,i_lw,i_dn) - flux_rad(k-1,i,j,i_lw,i_dn) ) &
                    - ( flux_rad(k,i,j,i_lw,i_up) - flux_rad(k-1,i,j,i_lw,i_up) ) &
                    ) / ( fz(k,i,j) - fz(k-1,i,j) )
 
          rhoh_sw = ( ( flux_rad(k,i,j,i_sw,i_dn) - flux_rad(k-1,i,j,i_sw,i_dn) ) &
                    - ( flux_rad(k,i,j,i_sw,i_up) - flux_rad(k-1,i,j,i_sw,i_up) ) &
                    ) / ( fz(k,i,j) - fz(k-1,i,j) )
 
          rhoh(k,i,j) = rhoh_lw + rhoh_sw
 
          temp_t(k,i,j,i_lw) = rhoh_lw / dens(k,i,j) / cvtot(k,i,j) * 86400.0_rp ! [K/day]
          temp_t(k,i,j,i_sw) = rhoh_sw / dens(k,i,j) / cvtot(k,i,j) * 86400.0_rp ! [K/day]
 
          temp_t(k,i,j,3)    = temp_t(k,i,j,i_lw) + temp_t(k,i,j,i_sw)
 
       enddo
       enddo
       enddo
       !$acc end kernels
 
    else
 
       !$omp parallel do default(none) OMP_SCHEDULE_ collapse(2) &
       !$omp private(i,j,k, &
       !$omp         RHOH_LW,RHOH_SW) &
       !$omp shared(RHOH,flux_rad,FZ, &
       !$omp        KS,KE,IS,IE,JS,JE)
       !$acc kernels
       do j = js, je
       do i = is, ie
       do k = ks, ke
 
          rhoh_lw = ( ( flux_rad(k,i,j,i_lw,i_dn) - flux_rad(k-1,i,j,i_lw,i_dn) ) &
                    - ( flux_rad(k,i,j,i_lw,i_up) - flux_rad(k-1,i,j,i_lw,i_up) ) &
                    ) / ( fz(k,i,j) - fz(k-1,i,j) )
 
          rhoh_sw = ( ( flux_rad(k,i,j,i_sw,i_dn) - flux_rad(k-1,i,j,i_sw,i_dn) ) &
                    - ( flux_rad(k,i,j,i_sw,i_up) - flux_rad(k-1,i,j,i_sw,i_up) ) &
                    ) / ( fz(k,i,j) - fz(k-1,i,j) )
 
          rhoh(k,i,j) = rhoh_lw + rhoh_sw
 
       enddo
       enddo
       enddo
       !$acc end kernels
 
    end if
 
    return

References i_dn, i_lw, i_sw, and i_up.

Referenced by mod_atmos_phy_rd_driver::atmos_phy_rd_driver_calc_tendency().

Here is the caller graph for this function:

Variable Documentation

i_up

integer, parameter, public scale_atmos_phy_rd_common::i_up = 1

Definition at line 34 of file scale_atmos_phy_rd_common.F90.

  integer,  public, parameter :: I_up = 1

Referenced by atmos_phy_rd_calc_heating(), mod_atmos_phy_rd_driver::atmos_phy_rd_driver_calc_tendency(), scale_atmos_phy_rd_mstrnx::atmos_phy_rd_mstrnx_flux(), and scale_atmos_phy_rd_offline::atmos_phy_rd_offline_flux().

i_dn

integer, parameter, public scale_atmos_phy_rd_common::i_dn = 2

Definition at line 35 of file scale_atmos_phy_rd_common.F90.

  integer,  public, parameter :: I_dn = 2

Referenced by atmos_phy_rd_calc_heating(), mod_atmos_phy_rd_driver::atmos_phy_rd_driver_calc_tendency(), scale_atmos_phy_rd_mstrnx::atmos_phy_rd_mstrnx_flux(), and scale_atmos_phy_rd_offline::atmos_phy_rd_offline_flux().

i_lw

integer, parameter, public scale_atmos_phy_rd_common::i_lw = 1

Definition at line 37 of file scale_atmos_phy_rd_common.F90.

  integer,  public, parameter :: I_LW = 1

Referenced by atmos_phy_rd_calc_heating(), mod_atmos_phy_rd_driver::atmos_phy_rd_driver_calc_tendency(), scale_atmos_phy_rd_mstrnx::atmos_phy_rd_mstrnx_flux(), and scale_atmos_phy_rd_offline::atmos_phy_rd_offline_flux().

i_sw

integer, parameter, public scale_atmos_phy_rd_common::i_sw = 2

Definition at line 38 of file scale_atmos_phy_rd_common.F90.

  integer,  public, parameter :: I_SW = 2

Referenced by atmos_phy_rd_calc_heating(), mod_atmos_phy_rd_driver::atmos_phy_rd_driver_calc_tendency(), scale_atmos_phy_rd_mstrnx::atmos_phy_rd_mstrnx_flux(), and scale_atmos_phy_rd_offline::atmos_phy_rd_offline_flux().