scale_atmos_phy_rd_common Module Reference

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 69 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, and scale_grid_index::ks.

Referenced by mod_atmos_phy_rd_driver::atmos_phy_rd_driver().

    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]

    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]

    ! 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]

    call thermodyn_cv( cvtot(:,:,:),   & ! [OUT]
                       qtrc(:,:,:,:), & ! [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

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 39 of file scale_atmos_phy_rd_common.F90.

Referenced by mod_atmos_phy_rd_driver::atmos_phy_rd_driver(), atmos_phy_rd_heating(), and scale_atmos_phy_rd_mstrnx::atmos_phy_rd_mstrnx().

  integer,  public, parameter :: i_up = 1

i_dn

integer, parameter, public scale_atmos_phy_rd_common::i_dn = 2

Definition at line 40 of file scale_atmos_phy_rd_common.F90.

Referenced by mod_atmos_phy_rd_driver::atmos_phy_rd_driver(), atmos_phy_rd_heating(), and scale_atmos_phy_rd_mstrnx::atmos_phy_rd_mstrnx().

  integer,  public, parameter :: i_dn = 2

i_lw

integer, parameter, public scale_atmos_phy_rd_common::i_lw = 1

Definition at line 42 of file scale_atmos_phy_rd_common.F90.

Referenced by mod_atmos_phy_rd_driver::atmos_phy_rd_driver(), atmos_phy_rd_heating(), and scale_atmos_phy_rd_mstrnx::atmos_phy_rd_mstrnx().

  integer,  public, parameter :: i_lw = 1

i_sw

integer, parameter, public scale_atmos_phy_rd_common::i_sw = 2

Definition at line 43 of file scale_atmos_phy_rd_common.F90.

Referenced by mod_atmos_phy_rd_driver::atmos_phy_rd_driver(), atmos_phy_rd_heating(), and scale_atmos_phy_rd_mstrnx::atmos_phy_rd_mstrnx().

  integer,  public, parameter :: i_sw = 2

i_direct

integer, parameter, public scale_atmos_phy_rd_common::i_direct = 1

Definition at line 45 of file scale_atmos_phy_rd_common.F90.

Referenced by mod_atmos_phy_rd_driver::atmos_phy_rd_driver().

  integer,  public, parameter :: i_direct  = 1

i_diffuse

integer, parameter, public scale_atmos_phy_rd_common::i_diffuse = 2

Definition at line 46 of file scale_atmos_phy_rd_common.F90.

Referenced by mod_atmos_phy_rd_driver::atmos_phy_rd_driver().

  integer,  public, parameter :: i_diffuse = 2