module ATMOSPHERE / Physics Radiation More...

Functions/Subroutines
subroutine, public	atmos_phy_rd_heating (flux_rad, DENS, RHOT, QTRC, FZ, dt, TEMP_t, RHOT_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

integer, parameter, public	i_direct = 1

integer, parameter, public	i_diffuse = 2

Detailed Description

module ATMOSPHERE / Physics Radiation

Description: Common module for Radiation Heating Rate

Author: Team SCALE

History

2013-02-06 (H.Yashiro) [new]

Function/Subroutine Documentation

◆ atmos_phy_rd_heating()

subroutine, public scale_atmos_phy_rd_common::atmos_phy_rd_heating	(	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)	RHOT,
		real(rp), dimension (ka,ia,ja,qa), intent(in)	QTRC,
		real(rp), dimension (0:ka,ia,ja), intent(in)	FZ,
		real(dp), intent(in)	dt,
		real(rp), dimension (ka,ia,ja,3), intent(out)	TEMP_t,
		real(rp), dimension (ka,ia,ja), intent(out)	RHOT_t
	)

Calc heating rate.

Definition at line 70 of file scale_atmos_phy_rd_common.F90.

References i_dn, i_lw, i_sw, i_up, scale_grid_index::ie, scale_grid_index::is, scale_grid_index::je, scale_grid_index::js, scale_grid_index::ke, scale_grid_index::ks, scale_tracer::tracer_cv, scale_tracer::tracer_mass, and scale_tracer::tracer_r.

Referenced by mod_atmos_phy_rd_driver::atmos_phy_rd_driver().

     use scale_tracer, only: &
        tracer_r, &
        tracer_cv, &
        tracer_mass
     use scale_atmos_thermodyn, only: &
        thermodyn_qd   => atmos_thermodyn_qd,   &
        thermodyn_cv   => atmos_thermodyn_cv,   &
        thermodyn_rhoe => atmos_thermodyn_rhoe, &
        thermodyn_rhot => atmos_thermodyn_rhot
     implicit none
 
     real(RP), intent(in)  :: flux_rad(KA,IA,JA,2,2)
     real(RP), intent(in)  :: DENS    (KA,IA,JA)
     real(RP), intent(in)  :: RHOT    (KA,IA,JA)
     real(RP), intent(in)  :: QTRC    (KA,IA,JA,QA)
     real(RP), intent(in)  :: FZ      (0:KA,IA,JA)
     real(DP), intent(in)  :: dt
     real(RP), intent(out) :: TEMP_t  (KA,IA,JA,3)
     real(RP), intent(out) :: RHOT_t  (KA,IA,JA)
 
     real(RP) :: RHOE  (KA,IA,JA)
     real(RP) :: RHOE_t(KA,IA,JA,2)
     real(RP) :: RHOT1 (KA,IA,JA)
     real(RP) :: QDRY  (KA,IA,JA)
     real(RP) :: CVtot (KA,IA,JA)
 
     integer :: k, i, j
     !---------------------------------------------------------------------------
 
     call thermodyn_rhoe( rhoe(:,:,:),   & ! [OUT]
                          rhot(:,:,:),   & ! [IN]
                          qtrc(:,:,:,:), & ! [IN]
                          tracer_cv(:),  & ! [IN]
                          tracer_r(:),   & ! [IN]
                          tracer_mass(:) ) ! [IN]
 
     !$omp parallel do default(none)                       &
     !$omp shared(JS,JE,IS,IE,KS,KE,RHOE_t,flux_rad,FZ,dt,RHOE) &
     !$omp private(i,j,k) OMP_SCHEDULE_ collapse(2)
     do j = js, je
     do i = is, ie
     do k = ks, ke
 
        rhoe_t(k,i,j,i_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) )
 
        rhoe_t(k,i,j,i_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) )
 
        rhoe(k,i,j) = rhoe(k,i,j) + dt * ( rhoe_t(k,i,j,i_lw) + rhoe_t(k,i,j,i_sw) )
 
     enddo
     enddo
     enddo
 
     call thermodyn_rhot( rhot1(:,:,:),  & ! [OUT]
                          rhoe(:,:,:),   & ! [IN]
                          qtrc(:,:,:,:), & ! [IN]
                          tracer_cv(:),  & ! [IN]
                          tracer_r(:),   & ! [IN]
                          tracer_mass(:) ) ! [IN]
 
     ! update rhot
     do j = js, je
     do i = is, ie
     do k = ks, ke
        rhot_t(k,i,j) = ( rhot1(k,i,j) - rhot(k,i,j) ) / dt
     enddo
     enddo
     enddo
 
     call thermodyn_qd( qdry(:,:,:),   & ! [OUT]
                        qtrc(:,:,:,:), & ! [IN]
                        tracer_mass(:) ) ! [IN]
 
     call thermodyn_cv( cvtot(:,:,:),   & ! [OUT]
                        qtrc(:,:,:,:), & ! [IN]
                        tracer_cv(:),   & ! [IN]
                        qdry(:,:,:)    ) ! [IN]
 
     do j = js, je
     do i = is, ie
     do k = ks, ke
        temp_t(k,i,j,i_lw) = rhoe_t(k,i,j,i_lw) / dens(k,i,j) / cvtot(k,i,j) * 86400.0_rp ! [K/day]
        temp_t(k,i,j,i_sw) = rhoe_t(k,i,j,i_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