scale_atmos_thermodyn Module Reference

module ATMOSPHERE / Thermodynamics More...

Functions/Subroutines
subroutine, public	atmos_thermodyn_setup
	Setup. More...
subroutine	atmos_thermodyn_qd_0d (qdry, q)
	calc dry air mass (0D) More...
subroutine, public	atmos_thermodyn_tempre (temp, pres, Ein, dens, qdry, q)
subroutine, public	atmos_thermodyn_tempre2 (temp, pres, dens, pott, qdry, q)

Variables
real(rp), public	thermodyn_emask = 0.0_RP
	=0: SIMPLE, 1: EXACT More...
real(rp), dimension(:), allocatable, public	aq_cp
	CP for each hydrometeors [J/kg/K]. More...
real(rp), dimension(:), allocatable, public	aq_cv
	CV for each hydrometeors [J/kg/K]. More...

Detailed Description

module ATMOSPHERE / Thermodynamics

Description

Thermodynamics module

Author

Team SCALE

History

• 2011-10-24 (T.Seiki) [new] Import from NICAM
• 2012-02-10 (H.Yashiro) [mod] Reconstruction
• 2012-03-23 (H.Yashiro) [mod] Explicit index parameter inclusion
• 2012-12-22 (S.Nishizawa) [mod] Use thermodyn macro set

Function/Subroutine Documentation

atmos_thermodyn_setup()

subroutine, public scale_atmos_thermodyn::atmos_thermodyn_setup

(

)

Setup.

Definition at line 159 of file scale_atmos_sub_thermodyn.F90.

References aq_cp, aq_cv, scale_const::const_ci, scale_const::const_cl, scale_const::const_cpvap, scale_const::const_cvvap, scale_const::const_thermodyn_type, scale_tracer::i_qv, scale_stdio::io_fid_log, scale_stdio::io_l, scale_process::prc_mpistop(), scale_tracer::qie, scale_tracer::qis, scale_tracer::qqa, scale_tracer::qwe, scale_tracer::qws, and thermodyn_emask.

Referenced by mod_rm_driver::scalerm(), and mod_rm_prep::scalerm_prep().

    use scale_process, only: &
       prc_mpistop
    use scale_const, only: &
       cpvap          => const_cpvap,          &
       cvvap          => const_cvvap,          &
       cl             => const_cl,             &
       ci             => const_ci,             &
       thermodyn_type => const_thermodyn_type
    implicit none

    integer :: n
    !---------------------------------------------------------------------------

    if( io_l ) write(io_fid_log,*)
    if( io_l ) write(io_fid_log,*) '++++++ Module[THERMODYN] / Categ[ATMOS SHARE] / Origin[SCALElib]'

    allocate( aq_cp(qqa) )
    allocate( aq_cv(qqa) )

    if ( thermodyn_type == 'EXACT' ) then
       thermodyn_emask = 1.0_rp

       aq_cp(i_qv) = cpvap
       aq_cv(i_qv) = cvvap

       if ( qws /= 0 ) then
          do n = qws, qwe
             aq_cp(n) = cl
             aq_cv(n) = cl
          enddo
       endif

       if ( qis /= 0 ) then
          do n = qis, qie
             aq_cp(n) = ci
             aq_cv(n) = ci
          enddo
       endif
    elseif( thermodyn_type == 'SIMPLE' ) then
       thermodyn_emask = 0.0_rp

       aq_cp(i_qv) = cpdry
       aq_cv(i_qv) = cvdry

       if ( qws /= 0 ) then
          do n = qws, qwe
             aq_cp(n) = cvdry
             aq_cv(n) = cvdry
          enddo
       endif

       if ( qis /= 0 ) then
          do n = qis, qie
             aq_cp(n) = cvdry
             aq_cv(n) = cvdry
          enddo
       endif
    endif

    return

Here is the call graph for this function:

Here is the caller graph for this function:

atmos_thermodyn_qd_0d()

subroutine scale_atmos_thermodyn::atmos_thermodyn_qd_0d	(	real(rp), intent(out)	qdry,
		real(rp), dimension(qa), intent(in)	q
	)

calc dry air mass (0D)

Parameters

[out]	qdry	dry mass concentration [kg/kg]
[in]	q	mass concentration [kg/kg]

Definition at line 226 of file scale_atmos_sub_thermodyn.F90.

References aq_cp, aq_cv, scale_tracer::i_qc, scale_tracer::i_qg, scale_tracer::i_qi, scale_tracer::i_qr, scale_tracer::i_qs, scale_tracer::i_qv, scale_grid_index::ieb, scale_grid_index::isb, scale_grid_index::jeb, scale_grid_index::jsb, scale_grid_index::ke, scale_grid_index::ks, dc_log::log(), scale_tracer::qqe, scale_tracer::qqs, and thermodyn_emask.

    implicit none

    real(RP), intent(out) :: qdry
    real(RP), intent(in)  :: q(qa)

    integer :: iqw
    !-----------------------------------------------------------------------------

    qdry = 1.0_rp
#ifndef DRY
    do iqw = qqs, qqe
       qdry = qdry - q(iqw)
    enddo
#endif

    return

Here is the call graph for this function:

atmos_thermodyn_tempre()

subroutine, public scale_atmos_thermodyn::atmos_thermodyn_tempre	(	real(rp), dimension(ka,ia,ja), intent(out)	temp,
		real(rp), dimension(ka,ia,ja), intent(out)	pres,
		real(rp), dimension (ka,ia,ja), intent(in)	Ein,
		real(rp), dimension(ka,ia,ja), intent(in)	dens,
		real(rp), dimension(ka,ia,ja), intent(in)	qdry,
		real(rp), dimension (ka,ia,ja,qa), intent(in)	q
	)

Definition at line 1282 of file scale_atmos_sub_thermodyn.F90.

References aq_cv, scale_tracer::i_qv, scale_grid_index::ieb, scale_grid_index::isb, scale_grid_index::jeb, scale_grid_index::jsb, scale_grid_index::ke, and scale_grid_index::ks.

    implicit none

    real(RP), intent(out) :: temp(ka,ia,ja)    ! temperature
    real(RP), intent(out) :: pres(ka,ia,ja)    ! pressure
    real(RP), intent(in)  :: ein (ka,ia,ja)    ! internal energy
    real(RP), intent(in)  :: dens(ka,ia,ja)    ! density
    real(RP), intent(in)  :: qdry(ka,ia,ja)    ! dry concentration
    real(RP), intent(in)  :: q   (ka,ia,ja,qa) ! water concentration

    real(RP) :: cv, rmoist

    integer :: i, j, k, iqw
    !---------------------------------------------------------------------------

    !$omp parallel do private(i,j,k,iqw,cv,Rmoist) OMP_SCHEDULE_ collapse(2)
    do j = jsb, jeb
    do i = isb, ieb
    do k = ks, ke

       calc_cv(cv, qdry(k,i,j), q, k, i, j, iqw, cvdry, aq_cv)
       calc_r(rmoist, q(k,i,j,i_qv), qdry(k,i,j), rdry, rvap)

       temp(k,i,j) = ein(k,i,j) / cv

       pres(k,i,j) = dens(k,i,j) * rmoist * temp(k,i,j)

    enddo
    enddo
    enddo

    return

atmos_thermodyn_tempre2()

subroutine, public scale_atmos_thermodyn::atmos_thermodyn_tempre2	(	real(rp), dimension(ka,ia,ja), intent(out)	temp,
		real(rp), dimension(ka,ia,ja), intent(out)	pres,
		real(rp), dimension(ka,ia,ja), intent(in)	dens,
		real(rp), dimension(ka,ia,ja), intent(in)	pott,
		real(rp), dimension(ka,ia,ja), intent(in)	qdry,
		real(rp), dimension (ka,ia,ja,qa), intent(in)	q
	)

Definition at line 1319 of file scale_atmos_sub_thermodyn.F90.

References aq_cp, scale_tracer::i_qv, scale_grid_index::ieb, scale_grid_index::isb, scale_grid_index::jeb, scale_grid_index::jsb, scale_grid_index::ke, and scale_grid_index::ks.

    implicit none

    real(RP), intent(out) :: temp(ka,ia,ja)    ! temperature
    real(RP), intent(out) :: pres(ka,ia,ja)    ! pressure
    real(RP), intent(in)  :: dens(ka,ia,ja)    ! density
    real(RP), intent(in)  :: pott(ka,ia,ja)    ! potential temperature
    real(RP), intent(in)  :: qdry(ka,ia,ja)    ! dry concentration
    real(RP), intent(in)  :: q   (ka,ia,ja,qa) ! water concentration

    real(RP) :: rmoist, cp

    integer :: i, j, k, iqw
    !---------------------------------------------------------------------------

    !$omp parallel do private(i,j,k,iqw,cp,Rmoist) OMP_SCHEDULE_ collapse(2)
    do j = jsb, jeb
    do i = isb, ieb
    do k = ks, ke

       calc_cp(cp, qdry(k,i,j), q, k, i, j, iqw, cpdry, aq_cp)
       calc_r(rmoist, q(k,i,j,i_qv), qdry(k,i,j), rdry, rvap)
       calc_pre(pres(k,i,j), dens(k,i,j), pott(k,i,j), rmoist, cp, pre00)

       temp(k,i,j) = pres(k,i,j) / ( dens(k,i,j) * rmoist )

    enddo
    enddo
    enddo

    return

Variable Documentation

thermodyn_emask

real(rp), public scale_atmos_thermodyn::thermodyn_emask = 0.0_RP

=0: SIMPLE, 1: EXACT

Definition at line 140 of file scale_atmos_sub_thermodyn.F90.

Referenced by atmos_thermodyn_qd_0d(), and atmos_thermodyn_setup().

  real(RP), public              :: thermodyn_emask = 0.0_rp 

aq_cp

real(rp), dimension(:), allocatable, public scale_atmos_thermodyn::aq_cp

CP for each hydrometeors [J/kg/K].

Definition at line 142 of file scale_atmos_sub_thermodyn.F90.

Referenced by scale_atmos_phy_ae_kajino13::atmos_phy_ae_kajino13(), scale_atmos_phy_mp_sn14::atmos_phy_mp_sn14(), atmos_thermodyn_qd_0d(), atmos_thermodyn_setup(), atmos_thermodyn_tempre2(), mod_atmos_vars::atmos_vars_history(), and scale_atmos_phy_mp_sn14::update_by_phase_change_kij().

  real(RP), public, allocatable :: aq_cp(:)

aq_cv

real(rp), dimension(:), allocatable, public scale_atmos_thermodyn::aq_cv

CV for each hydrometeors [J/kg/K].

Definition at line 143 of file scale_atmos_sub_thermodyn.F90.

Referenced by mod_atmos_dyn_driver::atmos_dyn_driver(), scale_atmos_phy_ae_kajino13::atmos_phy_ae_kajino13(), scale_atmos_phy_mp_common::atmos_phy_mp_precipitation(), scale_atmos_phy_mp_common::atmos_phy_mp_saturation_adjustment(), scale_atmos_phy_mp_sn14::atmos_phy_mp_sn14(), atmos_thermodyn_qd_0d(), atmos_thermodyn_setup(), atmos_thermodyn_tempre(), mod_atmos_vars::atmos_vars_history(), mod_atmos_vars::atmos_vars_monitor(), mod_atmos_vars::atmos_vars_total(), mod_mkinit::rect_setup(), and scale_atmos_phy_mp_sn14::update_by_phase_change_kij().

  real(RP), public, allocatable :: aq_cv(:)