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
     !---------------------------------------------------------------------------
  
     !$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)
     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
  
     return