scale_atmos_hydrometeor Module Reference

module Hydrometeor More...

Functions/Subroutines
subroutine, public	atmos_hydrometeor_setup
	Setup. More...
subroutine, public	atmos_hydrometeor_regist (Q0, NV, NL, NI, NAME, DESC, UNIT, ADVC)
	Regist tracer. More...
subroutine	atmos_hydrometeor_lhv_0d (lhv, temp)
subroutine	atmos_hydrometeor_lhv_1d (lhv, temp)
subroutine	atmos_hydrometeor_lhv_2d (lhv, temp)
subroutine	atmos_hydrometeor_lhv_3d (lhv, temp)
subroutine	atmos_hydrometeor_lhs_0d (lhs, temp)
subroutine	atmos_hydrometeor_lhs_1d (lhs, temp)
subroutine	atmos_hydrometeor_lhs_2d (lhs, temp)
subroutine	atmos_hydrometeor_lhs_3d (lhs, temp)
subroutine	atmos_hydrometeor_lhf_0d (lhf, temp)
subroutine	atmos_hydrometeor_lhf_1d (lhf, temp)
subroutine	atmos_hydrometeor_lhf_2d (lhf, temp)
subroutine	atmos_hydrometeor_lhf_3d (lhf, temp)
subroutine	atmos_hydrometeor_entr_0d (entr, temp, pres, q, Rq)
	calc temp, pres, q -> entropy (0D) More...
subroutine	atmos_hydrometeor_entr_2d (entr, temp, pres, q, Rq)
	calc temp, pres, q -> entropy (2D) More...
subroutine	atmos_hydrometeor_entr_3d (entr, temp, pres, q, Rq)
	calc temp, pres, q -> entropy (3D) More...
subroutine, public	atmos_hydrometeor_diagnose_number_concentration (QTRC)

Variables
integer, public	i_qv = -1
integer, parameter, public	n_hyd = 6
integer, parameter, public	i_hc = 1
	liquid water cloud More...
integer, parameter, public	i_hr = 2
	liquid water rain More...
integer, parameter, public	i_hi = 3
	ice water cloud More...
integer, parameter, public	i_hs = 4
	snow More...
integer, parameter, public	i_hg = 5
	graupel More...
integer, parameter, public	i_hh = 6
	hail More...
integer, public	i_qc = -1
integer, public	i_qr = -1
integer, public	i_qi = -1
integer, public	i_qs = -1
integer, public	i_qg = -1
integer, public	i_qh = -1
integer, public	i_nc = -1
integer, public	i_nr = -1
integer, public	i_ni = -1
integer, public	i_ns = -1
integer, public	i_ng = -1
integer, public	i_nh = -1
integer, public	qhs = -1
integer, public	qhe = -2
integer, public	qls = -1
integer, public	qle = -2
integer, public	qis = -1
integer, public	qie = -2
real(rp), public	lhv
	latent heat of vaporization for use [J/kg] More...
real(rp), public	lhs
	latent heat of sublimation for use [J/kg] More...
real(rp), public	lhf
	latent heat of fusion for use [J/kg] More...

Detailed Description

module Hydrometeor

Description

Hydrometeor module

Author

Team SCALE

History

• 2016-08-06 (S.Nishizawa) [new]

Function/Subroutine Documentation

atmos_hydrometeor_setup()

subroutine, public scale_atmos_hydrometeor::atmos_hydrometeor_setup

(

)

Setup.

Definition at line 129 of file scale_atmos_sub_hydrometeor.F90.

References scale_const::const_ci, scale_const::const_cl, scale_const::const_cpvap, scale_const::const_cvvap, scale_const::const_lhf0, scale_const::const_lhf00, scale_const::const_lhs0, scale_const::const_lhs00, scale_const::const_lhv0, scale_const::const_lhv00, scale_const::const_thermodyn_type, scale_stdio::io_fid_log, scale_stdio::io_l, lhf, lhs, lhv, and scale_process::prc_mpistop().

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

    use scale_const, only: &
       cpvap          => const_cpvap,          &
       cvvap          => const_cvvap,          &
       cl             => const_cl,             &
       ci             => const_ci,             &
       lhv00          => const_lhv00,          &
       lhs00          => const_lhs00,          &
       lhf00          => const_lhf00,          &
       lhv0           => const_lhv0,           &
       lhs0           => const_lhs0,           &
       lhf0           => const_lhf0,           &
       thermodyn_type => const_thermodyn_type
    use scale_process, only: &
       prc_mpistop
    implicit none
    !---------------------------------------------------------------------------

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

    if ( thermodyn_type == 'EXACT' ) then

       cv_vapor = cvvap
       cp_vapor = cpvap
       cv_water = cl
       cp_water = cv_water
       cv_ice   = ci
       cp_ice   = cv_ice

       lhv = lhv00
       lhs = lhs00
       lhf = lhf00
       thermodyn_emask = 1.0_rp

    elseif( thermodyn_type == 'SIMPLE' ) then

       cv_vapor = cvvap
       cp_vapor = cpvap
       cv_water = cvvap
       cp_water = cv_water
       cv_ice   = cvvap
       cp_ice   = cv_ice

       lhv = lhv0
       lhs = lhs0
       lhf = lhf0
       thermodyn_emask = 0.0_rp

    else
       write(*,*) 'xxx Not appropriate ATMOS_THERMODYN_ENERGY_TYPE. Check!', trim(thermodyn_type)
       call prc_mpistop
    endif

    return

Here is the call graph for this function:

Here is the caller graph for this function:

atmos_hydrometeor_regist()

subroutine, public scale_atmos_hydrometeor::atmos_hydrometeor_regist	(	integer, intent(out)	Q0,
		integer, intent(in)	NV,
		integer, intent(in)	NL,
		integer, intent(in)	NI,
		character(len=*), dimension(nv+nl+ni), intent(in)	NAME,
		character(len=*), dimension(nv+nl+ni), intent(in)	DESC,
		character(len=*), dimension(nv+nl+ni), intent(in)	UNIT,
		logical, dimension(nv+nl+ni), intent(in), optional	ADVC
	)

Regist tracer.

Parameters

[in]	nv	number of vapor
[in]	nl	number of liquid water tracers
[in]	ni	number of ice water tracers

Definition at line 196 of file scale_atmos_sub_hydrometeor.F90.

References scale_const::const_rvap, i_qv, scale_process::prc_mpistop(), qhe, qhs, qie, qis, qle, qls, and scale_tracer::tracer_regist().

Referenced by scale_atmos_phy_mp_kessler::atmos_phy_mp_kessler_config(), scale_atmos_phy_mp_sn14::atmos_phy_mp_sn14_config(), scale_atmos_phy_mp_suzuki10::atmos_phy_mp_suzuki10_config(), and scale_atmos_phy_mp_tomita08::atmos_phy_mp_tomita08_config().

    use scale_process, only: &
      prc_mpistop
    use scale_tracer, only: &
      tracer_regist
    use scale_const, only: &
      rvap => const_rvap
    implicit none

    integer,          intent(out) :: Q0
    integer,          intent(in)  :: NV
    integer,          intent(in)  :: NL
    integer,          intent(in)  :: NI
    character(len=*), intent(in)  :: NAME(NV+NL+NI)
    character(len=*), intent(in)  :: DESC(NV+NL+NI)
    character(len=*), intent(in)  :: UNIT(NV+NL+NI)

    logical,          intent(in), optional :: ADVC(NV+NL+NI)

    real(RP) :: CV   (NV+NL+NI)
    real(RP) :: CP   (NV+NL+NI)
    real(RP) :: R    (NV+NL+NI)
    logical  :: MASS (NV+NL+NI)
    logical  :: ADVC_(NV+NL+NI)

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

    if ( i_qv > 0 ) then
       write(*,*) 'xxx tracer for hydrometeor is already registerd'
       call prc_mpistop
    endif

    if ( nv /= 1 ) then
       write(*,*) 'xxx number of vapor must be 1 at this moment'
       call prc_mpistop
    endif

    nq = 0

    do n = 1, nv
       nq = nq + 1
       cv(nq) = cv_vapor
       cp(nq) = cp_vapor
       r(nq) = rvap
    end do

    do n = 1, nl
       nq = nq + 1
       cv(nq) = cv_water
       cp(nq) = cp_water
       r(nq) = 0.0_rp
    end do

    do n = 1, ni
       nq = nq + 1
       cv(nq) = cv_ice
       cp(nq) = cp_ice
       r(nq) = 0.0_rp
    end do

    ! NQ = 1 + NL + NI,   vapor + liqid + ice

    if ( present(advc) ) then
       advc_(:) = advc(:)
    else
       advc_(:) = .true.
    endif

    do n = 1, nq
       mass(n) = .true.
    end do

    call tracer_regist( q0,         & ! [OUT]
                        nq,         & ! [IN]
                        name,       & ! [IN]
                        desc,       & ! [IN]
                        unit,       & ! [IN]
                        cv, cp, r,  & ! [IN], optional
                        advc_, mass ) ! [IN], optional

    i_qv = q0

    if ( nq > 1 ) then
       qhs = i_qv + 1
       qhe = qhs + nl + ni - 1
    endif

    if ( nl > 0 ) then
       qls = i_qv + 1
       qle = qls + nl - 1
    endif

    if ( ni > 0 ) then
       qis = qle + 1
       qie = qis + ni - 1
    endif

    return

Here is the call graph for this function:

Here is the caller graph for this function:

atmos_hydrometeor_lhv_0d()

subroutine scale_atmos_hydrometeor::atmos_hydrometeor_lhv_0d	(	real(rp), intent(out)	lhv,
		real(rp), intent(in)	temp
	)

Parameters

[out]	lhv	latent heat of vaporization [J/kg]
[in]	temp	temperature [K]

Definition at line 301 of file scale_atmos_sub_hydrometeor.F90.

References scale_const::const_lhv0, and scale_const::const_tem00.

    use scale_const, only: &
       tem00 => const_tem00, &
       lhv0  => const_lhv0
    implicit none

    real(RP), intent(out) :: lhv
    real(RP), intent(in)  :: temp
    !---------------------------------------------------------------------------

    lhv = lhv0 + ( cp_vapor - cp_water ) * ( temp - tem00 ) * thermodyn_emask

    return

atmos_hydrometeor_lhv_1d()

subroutine scale_atmos_hydrometeor::atmos_hydrometeor_lhv_1d	(	real(rp), dimension (ka), intent(out)	lhv,
		real(rp), dimension(ka), intent(in)	temp
	)

Parameters

[out]	lhv	latent heat of vaporization [J/kg]
[in]	temp	temperature [K]

Definition at line 319 of file scale_atmos_sub_hydrometeor.F90.

References scale_const::const_lhv0, scale_const::const_tem00, scale_grid_index::ke, and scale_grid_index::ks.

    use scale_const, only: &
       tem00 => const_tem00, &
       lhv0  => const_lhv0
    implicit none

    real(RP), intent(out) :: lhv (KA)
    real(RP), intent(in)  :: temp(KA)

    integer :: k
    !---------------------------------------------------------------------------

    do k = ks, ke
       lhv(k) = lhv0 + ( cp_vapor - cp_water ) * ( temp(k) - tem00 ) * thermodyn_emask
    enddo

    return

atmos_hydrometeor_lhv_2d()

subroutine scale_atmos_hydrometeor::atmos_hydrometeor_lhv_2d	(	real(rp), dimension (ia,ja), intent(out)	lhv,
		real(rp), dimension(ia,ja), intent(in)	temp
	)

Parameters

[out]	lhv	latent heat of vaporization [J/kg]
[in]	temp	temperature [K]

Definition at line 341 of file scale_atmos_sub_hydrometeor.F90.

References scale_const::const_lhv0, scale_const::const_tem00, scale_grid_index::ieb, scale_grid_index::isb, scale_grid_index::jeb, and scale_grid_index::jsb.

    use scale_const, only: &
       tem00 => const_tem00, &
       lhv0  => const_lhv0
    implicit none

    real(RP), intent(out) :: lhv (IA,JA)
    real(RP), intent(in)  :: temp(IA,JA)

    integer :: i, j
    !---------------------------------------------------------------------------

    do j = jsb, jeb
    do i = isb, ieb
       lhv(i,j) = lhv0 + ( cp_vapor - cp_water ) * ( temp(i,j) - tem00 ) * thermodyn_emask
    enddo
    enddo

    return

atmos_hydrometeor_lhv_3d()

subroutine scale_atmos_hydrometeor::atmos_hydrometeor_lhv_3d	(	real(rp), dimension (ka,ia,ja), intent(out)	lhv,
		real(rp), dimension(ka,ia,ja), intent(in)	temp
	)

Parameters

[out]	lhv	latent heat of vaporization [J/kg]
[in]	temp	temperature [K]

Definition at line 365 of file scale_atmos_sub_hydrometeor.F90.

References scale_const::const_lhv0, scale_const::const_tem00, 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.

    use scale_const, only: &
       tem00 => const_tem00, &
       lhv0  => const_lhv0
    implicit none

    real(RP), intent(out) :: lhv (KA,IA,JA)
    real(RP), intent(in)  :: temp(KA,IA,JA)

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

    do j = jsb, jeb
    do i = isb, ieb
    do k = ks, ke
       lhv(k,i,j) = lhv0 + ( cp_vapor - cp_water ) * ( temp(k,i,j) - tem00 ) * thermodyn_emask
    enddo
    enddo
    enddo

    return

atmos_hydrometeor_lhs_0d()

subroutine scale_atmos_hydrometeor::atmos_hydrometeor_lhs_0d	(	real(rp), intent(out)	lhs,
		real(rp), intent(in)	temp
	)

Parameters

[out]	lhs	latent heat of sublimation [J/kg]
[in]	temp	temperature [K]

Definition at line 391 of file scale_atmos_sub_hydrometeor.F90.

References scale_const::const_lhs0, and scale_const::const_tem00.

    use scale_const, only: &
       tem00 => const_tem00, &
       lhs0  => const_lhs0
    implicit none

    real(RP), intent(out) :: lhs
    real(RP), intent(in)  :: temp
    !---------------------------------------------------------------------------

    lhs = lhs0 + ( cp_vapor - cp_ice ) * ( temp - tem00 ) * thermodyn_emask

    return

atmos_hydrometeor_lhs_1d()

subroutine scale_atmos_hydrometeor::atmos_hydrometeor_lhs_1d	(	real(rp), dimension (ka), intent(out)	lhs,
		real(rp), dimension(ka), intent(in)	temp
	)

Parameters

[out]	lhs	latent heat of sublimation [J/kg]
[in]	temp	temperature [K]

Definition at line 409 of file scale_atmos_sub_hydrometeor.F90.

References scale_const::const_lhs0, scale_const::const_tem00, scale_grid_index::ke, and scale_grid_index::ks.

    use scale_const, only: &
       tem00 => const_tem00, &
       lhs0  => const_lhs0
    implicit none

    real(RP), intent(out) :: lhs (KA)
    real(RP), intent(in)  :: temp(KA)

    integer :: k
    !---------------------------------------------------------------------------

    do k = ks, ke
       lhs(k) = lhs0 + ( cp_vapor - cp_ice ) * ( temp(k) - tem00 ) * thermodyn_emask
    enddo

    return

atmos_hydrometeor_lhs_2d()

subroutine scale_atmos_hydrometeor::atmos_hydrometeor_lhs_2d	(	real(rp), dimension (ia,ja), intent(out)	lhs,
		real(rp), dimension(ia,ja), intent(in)	temp
	)

Parameters

[out]	lhs	latent heat of sublimation [J/kg]
[in]	temp	temperature [K]

Definition at line 431 of file scale_atmos_sub_hydrometeor.F90.

References scale_const::const_lhs0, scale_const::const_tem00, scale_grid_index::ieb, scale_grid_index::isb, scale_grid_index::jeb, and scale_grid_index::jsb.

    use scale_const, only: &
       tem00 => const_tem00, &
       lhs0  => const_lhs0
    implicit none

    real(RP), intent(out) :: lhs (IA,JA)
    real(RP), intent(in)  :: temp(IA,JA)

    integer :: i, j
    !---------------------------------------------------------------------------

    do j = jsb, jeb
    do i = isb, ieb
       lhs(i,j) = lhs0 + ( cp_vapor - cp_ice ) * ( temp(i,j) - tem00 ) * thermodyn_emask
    enddo
    enddo

    return

atmos_hydrometeor_lhs_3d()

subroutine scale_atmos_hydrometeor::atmos_hydrometeor_lhs_3d	(	real(rp), dimension (ka,ia,ja), intent(out)	lhs,
		real(rp), dimension(ka,ia,ja), intent(in)	temp
	)

Parameters

[out]	lhs	latent heat of sublimation [J/kg]
[in]	temp	temperature [K]

Definition at line 455 of file scale_atmos_sub_hydrometeor.F90.

References scale_const::const_lhs0, scale_const::const_tem00, 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.

    use scale_const, only: &
       tem00 => const_tem00, &
       lhs0  => const_lhs0
    implicit none

    real(RP), intent(out) :: lhs (KA,IA,JA)
    real(RP), intent(in)  :: temp(KA,IA,JA)

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

    do j = jsb, jeb
    do i = isb, ieb
    do k = ks, ke
       lhs(k,i,j) = lhs0 + ( cp_vapor - cp_ice ) * ( temp(k,i,j) - tem00 ) * thermodyn_emask
    enddo
    enddo
    enddo

    return

atmos_hydrometeor_lhf_0d()

subroutine scale_atmos_hydrometeor::atmos_hydrometeor_lhf_0d	(	real(rp), intent(out)	lhf,
		real(rp), intent(in)	temp
	)

Parameters

[out]	lhf	latent heat of fusion [J/kg]
[in]	temp	temperature [K]

Definition at line 481 of file scale_atmos_sub_hydrometeor.F90.

References scale_const::const_lhf0, and scale_const::const_tem00.

    use scale_const, only: &
       tem00 => const_tem00, &
       lhf0  => const_lhf0
    implicit none

    real(RP), intent(out) :: lhf
    real(RP), intent(in)  :: temp
    !---------------------------------------------------------------------------

    lhf = lhf0 + ( cp_water - cp_ice ) * ( temp - tem00 ) * thermodyn_emask

    return

atmos_hydrometeor_lhf_1d()

subroutine scale_atmos_hydrometeor::atmos_hydrometeor_lhf_1d	(	real(rp), dimension (ka), intent(out)	lhf,
		real(rp), dimension(ka), intent(in)	temp
	)

Parameters

[out]	lhf	latent heat of fusion [J/kg]
[in]	temp	temperature [K]

Definition at line 499 of file scale_atmos_sub_hydrometeor.F90.

References scale_const::const_lhf0, scale_const::const_tem00, scale_grid_index::ke, and scale_grid_index::ks.

    use scale_const, only: &
       tem00 => const_tem00, &
       lhf0  => const_lhf0
    implicit none

    real(RP), intent(out) :: lhf (KA)
    real(RP), intent(in)  :: temp(KA)

    integer :: k
    !---------------------------------------------------------------------------

    do k = ks, ke
       lhf(k) = lhf0 + ( cp_water - cp_ice ) * ( temp(k) - tem00 ) * thermodyn_emask
    enddo

    return

atmos_hydrometeor_lhf_2d()

subroutine scale_atmos_hydrometeor::atmos_hydrometeor_lhf_2d	(	real(rp), dimension (ia,ja), intent(out)	lhf,
		real(rp), dimension(ia,ja), intent(in)	temp
	)

Parameters

[out]	lhf	latent heat of fusion [J/kg]
[in]	temp	temperature [K]

Definition at line 521 of file scale_atmos_sub_hydrometeor.F90.

References scale_const::const_lhf0, scale_const::const_tem00, scale_grid_index::ieb, scale_grid_index::isb, scale_grid_index::jeb, and scale_grid_index::jsb.

    use scale_const, only: &
       tem00 => const_tem00, &
       lhf0  => const_lhf0
    implicit none

    real(RP), intent(out) :: lhf (IA,JA)
    real(RP), intent(in)  :: temp(IA,JA)

    integer :: i, j
    !---------------------------------------------------------------------------

    do j = jsb, jeb
    do i = isb, ieb
       lhf(i,j) = lhf0 + ( cp_water - cp_ice ) * ( temp(i,j) - tem00 ) * thermodyn_emask
    enddo
    enddo

    return

atmos_hydrometeor_lhf_3d()

subroutine scale_atmos_hydrometeor::atmos_hydrometeor_lhf_3d	(	real(rp), dimension (ka,ia,ja), intent(out)	lhf,
		real(rp), dimension(ka,ia,ja), intent(in)	temp
	)

Parameters

[out]	lhf	latent heat of fusion [J/kg]
[in]	temp	temperature [K]

Definition at line 545 of file scale_atmos_sub_hydrometeor.F90.

References scale_const::const_lhf0, scale_const::const_tem00, 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.

    use scale_const, only: &
       tem00 => const_tem00, &
       lhf0  => const_lhf0
    implicit none

    real(RP), intent(out) :: lhf (KA,IA,JA)
    real(RP), intent(in)  :: temp(KA,IA,JA)

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

    do j = jsb, jeb
    do i = isb, ieb
    do k = ks, ke
       lhf(k,i,j) = lhf0 + ( cp_water - cp_ice ) * ( temp(k,i,j) - tem00 ) * thermodyn_emask
    enddo
    enddo
    enddo

    return

atmos_hydrometeor_entr_0d()

subroutine scale_atmos_hydrometeor::atmos_hydrometeor_entr_0d	(	real(rp), intent(out)	entr,
		real(rp), intent(in)	temp,
		real(rp), intent(in)	pres,
		real(rp), dimension(qa), intent(in)	q,
		real(rp), dimension(qa), intent(in)	Rq
	)

calc temp, pres, q -> entropy (0D)

Parameters

[out]	entr	entropy [J/K]
[in]	temp	temperature [K]
[in]	pres	pressure [Pa]
[in]	q	mass concentration [kg/kg]
[in]	rq	gas constantt [J/kg/K]

Definition at line 575 of file scale_atmos_sub_hydrometeor.F90.

References scale_const::const_cpdry, scale_const::const_eps, scale_const::const_lhf0, scale_const::const_lhv0, scale_const::const_pre00, scale_const::const_psat0, scale_const::const_rdry, scale_const::const_rvap, scale_const::const_tem00, i_qv, scale_tracer::qa, qie, qis, qle, and qls.

    use scale_const, only: &
       eps   => const_eps,   &
       pre00 => const_pre00, &
       tem00 => const_tem00, &
       rdry  => const_rdry,  &
       cpdry => const_cpdry, &
       rvap  => const_rvap,  &
       lhv0  => const_lhv0,  &
       lhf0  => const_lhf0,  &
       psat0 => const_psat0
    use scale_tracer, only: &
       qa
    implicit none

    real(RP), intent(out) :: entr
    real(RP), intent(in)  :: temp
    real(RP), intent(in)  :: pres
    real(RP), intent(in)  :: q(QA)
    real(RP), intent(in)  :: Rq(QA)

    real(RP) :: qdry, Rtot
    real(RP) :: logT_T0, pres_dry, pres_vap

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

    logt_t0 = log( temp / tem00 )

    qdry = 1.0_rp
    rtot = 0.0_rp
    do iqw = 1, qa
       qdry = qdry - q(iqw)
       rtot = rtot + q(iqw) * rq(iqw)
    enddo
    rtot = rtot + rdry * qdry

    ! dry air + vapor
    pres_dry = max( pres * qdry * rdry / rtot, eps )
    entr = qdry * cpdry * logt_t0 &
         - qdry * rdry  * log( pres_dry / pre00 )

    if ( i_qv > 0 ) then
       pres_vap = max( pres * q(i_qv) * rvap / rtot, eps )
       entr = entr + q(i_qv) * cp_vapor * logt_t0 &
                   - q(i_qv) * rvap  * log( pres_vap / psat0 ) &
                   + q(i_qv) * lhv0 / tem00
    endif

    ! liquid water
    if ( qls > 0 ) then
       do iqw = qls, qle
          entr = entr + q(iqw) * cp_water * logt_t0
       enddo
    endif

    ! ice
    if ( qis > 0 ) then
       do iqw = qis, qie
          entr = entr + q(iqw) * cp_ice * logt_t0 &
                      - q(iqw) * lhf0 / tem00
       enddo
    endif

    return

atmos_hydrometeor_entr_2d()

subroutine scale_atmos_hydrometeor::atmos_hydrometeor_entr_2d	(	real(rp), dimension(ia,ja), intent(out)	entr,
		real(rp), dimension(ia,ja), intent(in)	temp,
		real(rp), dimension(ia,ja), intent(in)	pres,
		real(rp), dimension (ia,ja,qa), intent(in)	q,
		real(rp), dimension (qa), intent(in)	Rq
	)

calc temp, pres, q -> entropy (2D)

Parameters

[out]	entr	entropy [J/K]
[in]	temp	temperature [K]
[in]	pres	pressure [Pa]
[in]	q	mass concentration [kg/kg]
[in]	rq	gas constant [J/kg/K]

Definition at line 649 of file scale_atmos_sub_hydrometeor.F90.

References scale_const::const_cpdry, scale_const::const_eps, scale_const::const_lhf0, scale_const::const_lhv0, scale_const::const_pre00, scale_const::const_psat0, scale_const::const_rdry, scale_const::const_rvap, scale_const::const_tem00, i_qv, scale_grid_index::ieb, scale_grid_index::isb, scale_grid_index::jeb, scale_grid_index::jsb, scale_tracer::qa, qie, qis, qle, and qls.

    use scale_const, only: &
       eps   => const_eps,   &
       pre00 => const_pre00, &
       tem00 => const_tem00, &
       cpdry => const_cpdry, &
       rdry  => const_rdry,  &
       rvap  => const_rvap,  &
       lhv0  => const_lhv0,  &
       lhf0  => const_lhf0,  &
       psat0 => const_psat0
    use scale_tracer, only: &
       qa
    implicit none

    real(RP), intent(out) :: entr(IA,JA)
    real(RP), intent(in)  :: temp(IA,JA)
    real(RP), intent(in)  :: pres(IA,JA)
    real(RP), intent(in)  :: q   (IA,JA,QA)
    real(RP), intent(in)  :: Rq  (QA)

    real(RP) :: qdry, Rtot
    real(RP) :: logT_T0, pres_dry, pres_vap

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

    ! dry air + vapor
    do j = jsb, jeb
    do i = isb, ieb

       logt_t0 = log( temp(i,j) / tem00 )

       qdry = 1.0_rp
       rtot = 0.0_rp
       do iqw = 1, qa
          qdry = qdry - q(i,j,iqw)
          rtot = rtot + q(i,j,iqw) * rq(iqw)
       enddo
       rtot = rtot + rdry * qdry

       ! dry air + vapor
       pres_dry = max( pres(i,j) * qdry    * rdry / rtot, eps )
       entr(i,j) = qdry    * cpdry * logt_t0 &
                 - qdry    * rdry  * log( pres_dry / pre00 )

       if ( i_qv > 0 ) then
          pres_vap = max( pres(i,j) * q(i,j,i_qv) * rvap / rtot, eps )
          entr(i,j) = entr(i,j) + q(i,j,i_qv) * cp_vapor * logt_t0 &
                                - q(i,j,i_qv) * rvap  * log( pres_vap / psat0 ) &
                                + q(i,j,i_qv) * lhv0 / tem00
       endif

       ! liquid water
       if ( qls > 0 ) then
          do iqw = qls, qle
             entr(i,j) = entr(i,j) + q(i,j,iqw) * cp_water * logt_t0
          enddo
       endif

       ! ice
       if ( qis > 0 ) then
          do iqw = qis, qie
             entr(i,j) = entr(i,j) + q(i,j,iqw) * cp_ice * logt_t0 &
                                   - q(i,j,iqw) * lhf0 / tem00
          enddo
       endif

    enddo
    enddo

    return

atmos_hydrometeor_entr_3d()

subroutine scale_atmos_hydrometeor::atmos_hydrometeor_entr_3d	(	real(rp), dimension(ka,ia,ja), intent(out)	entr,
		real(rp), dimension(ka,ia,ja), intent(in)	temp,
		real(rp), dimension(ka,ia,ja), intent(in)	pres,
		real(rp), dimension (ka,ia,ja,qa), intent(in)	q,
		real(rp), dimension (qa), intent(in)	Rq
	)

calc temp, pres, q -> entropy (3D)

Parameters

[out]	entr	entropy [J/K]
[in]	temp	temperature [K]
[in]	pres	pressure [Pa]
[in]	q	mass concentration [kg/kg]
[in]	rq	gas constant [J/kg/K]

Definition at line 730 of file scale_atmos_sub_hydrometeor.F90.

References scale_const::const_cpdry, scale_const::const_eps, scale_const::const_lhf0, scale_const::const_lhv0, scale_const::const_pre00, scale_const::const_psat0, scale_const::const_rdry, scale_const::const_rvap, scale_const::const_tem00, 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, scale_tracer::qa, qie, qis, qle, and qls.

    use scale_const, only: &
       eps   => const_eps,   &
       pre00 => const_pre00, &
       tem00 => const_tem00, &
       cpdry => const_cpdry, &
       rdry  => const_rdry,  &
       rvap  => const_rvap,  &
       lhv0  => const_lhv0,  &
       lhf0  => const_lhf0,  &
       psat0 => const_psat0
    use scale_tracer, only: &
       qa
    implicit none

    real(RP), intent(out) :: entr(KA,IA,JA)
    real(RP), intent(in)  :: temp(KA,IA,JA)
    real(RP), intent(in)  :: pres(KA,IA,JA)
    real(RP), intent(in)  :: q   (KA,IA,JA,QA)
    real(RP), intent(in)  :: Rq  (QA)

    real(RP) :: qdry, Rtot
    real(RP) :: logT_T0, pres_dry, pres_vap

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

    ! dry air + vapor
    do j = jsb, jeb
    do i = isb, ieb
    do k = ks, ke

       logt_t0 = log( temp(k,i,j) / tem00 )

       qdry = 1.0_rp
       rtot = 0.0_rp
       do iqw = 1, qa
          qdry = qdry - q(k,i,j,iqw)
          rtot = rtot + q(k,i,j,iqw) * rq(iqw)
       enddo
       rtot = rtot + rdry * qdry

       ! dry air + vapor
       pres_dry = max( pres(k,i,j) * qdry    * rdry / rtot, eps )
       entr(k,i,j) = qdry    * cpdry * logt_t0 &
                   - qdry    * rdry  * log( pres_dry / pre00 )

       if ( i_qv > 0 ) then
          pres_vap = max( pres(k,i,j) * q(k,i,j,i_qv) * rvap / rtot, eps )
          entr(k,i,j) = entr(k,i,j) + q(k,i,j,i_qv) * cp_vapor * logt_t0 &
                                    - q(k,i,j,i_qv) * rvap  * log( pres_vap / psat0 ) &
                                    + q(k,i,j,i_qv) * lhv0 / tem00
       endif

       ! liquid water
       if ( qls > 0 ) then
          do iqw = qls, qle
             entr(k,i,j) = entr(k,i,j) + q(k,i,j,iqw) * cp_water * logt_t0
          enddo
       endif

       ! ice
       if ( qis > 0 ) then
          do iqw = qis, qie
             entr(k,i,j) = entr(k,i,j) + q(k,i,j,iqw) * cp_ice * logt_t0 &
                                       - q(k,i,j,iqw) * lhf0 / tem00
          enddo
       endif

    enddo
    enddo
    enddo

    return

atmos_hydrometeor_diagnose_number_concentration()

subroutine, public scale_atmos_hydrometeor::atmos_hydrometeor_diagnose_number_concentration

(

real(rp), dimension(:,:,:,:), intent(inout)

QTRC

)

Definition at line 808 of file scale_atmos_sub_hydrometeor.F90.

References scale_const::const_pi, i_nc, i_ng, i_ni, i_nr, i_ns, i_qc, i_qg, i_qi, i_qr, and i_qs.

Referenced by scale_atmos_boundary::atmos_boundary_update(), mod_atmos_phy_cp_driver::atmos_phy_cp_driver(), and mod_realinput::parentatomsetup().

    use scale_const, only: &
       pi => const_pi
    implicit none

    real(RP), intent(inout) :: QTRC(:,:,:,:)

    real(RP), parameter :: Dc   =  20.e-6_rp ! typical particle diameter for cloud  [m]
    real(RP), parameter :: Dr   = 200.e-6_rp ! typical particle diameter for rain   [m]
    real(RP), parameter :: Di   =  80.e-6_rp ! typical particle diameter for ice    [m]
    real(RP), parameter :: Ds   =  80.e-6_rp ! typical particle diameter for snow   [m]
    real(RP), parameter :: Dg   = 200.e-6_rp ! typical particle diameter for grapel [m]
    real(RP), parameter :: RHOw =  1000.0_rp ! typical density for warm particles   [kg/m3]
    real(RP), parameter :: RHOf =   100.0_rp ! typical density for frozen particles [kg/m3]
    real(RP), parameter :: RHOg =   400.0_rp ! typical density for grapel particles [kg/m3]
    real(RP), parameter :: b    =     3.0_rp ! assume spherical form

    real(RP) :: piov6
    !---------------------------------------------------------------------------

    piov6 = pi / 6.0_rp

    if ( i_nc > 0 ) qtrc(:,:,:,i_nc) = qtrc(:,:,:,i_qc) / ( (piov6*rhow) * dc**b )
    if ( i_nr > 0 ) qtrc(:,:,:,i_nr) = qtrc(:,:,:,i_qr) / ( (piov6*rhow) * dr**b )
    if ( i_ni > 0 ) qtrc(:,:,:,i_ni) = qtrc(:,:,:,i_qi) / ( (piov6*rhof) * di**b )
    if ( i_ns > 0 ) qtrc(:,:,:,i_ns) = qtrc(:,:,:,i_qs) / ( (piov6*rhof) * ds**b )
    if ( i_ng > 0 ) qtrc(:,:,:,i_ng) = qtrc(:,:,:,i_qg) / ( (piov6*rhog) * dg**b )

    return

Here is the caller graph for this function:

Variable Documentation

i_qv

integer, public scale_atmos_hydrometeor::i_qv = -1

Definition at line 68 of file scale_atmos_sub_hydrometeor.F90.

Referenced by scale_atmos_adiabat::atmos_adiabat_cape(), scale_atmos_adiabat::atmos_adiabat_liftparcel(), scale_atmos_dyn_common::atmos_dyn_numfilter_coef_q(), atmos_hydrometeor_entr_0d(), atmos_hydrometeor_entr_2d(), atmos_hydrometeor_entr_3d(), atmos_hydrometeor_regist(), scale_atmos_hydrostatic::atmos_hydrostatic_setup(), scale_atmos_phy_ae_kajino13::atmos_phy_ae_kajino13(), scale_atmos_phy_ae_kajino13::atmos_phy_ae_kajino13_mkinit(), scale_atmos_phy_cp_kf::atmos_phy_cp_kf_setup(), scale_atmos_phy_mp_convert::atmos_phy_mp_bulk2bin(), scale_atmos_phy_mp_kessler::atmos_phy_mp_kessler(), scale_atmos_phy_mp_sn14::atmos_phy_mp_sn14(), scale_atmos_phy_mp_tomita08::atmos_phy_mp_tomita08(), mod_atmos_phy_rd_driver::atmos_phy_rd_driver(), scale_atmos_phy_rd_mstrnx::atmos_phy_rd_mstrnx(), scale_atmos_phy_sf_bulk::atmos_phy_sf_bulk(), scale_atmos_phy_sf_const::atmos_phy_sf_const(), mod_atmos_phy_sf_driver::atmos_phy_sf_driver(), scale_atmos_phy_tb_mynn::atmos_phy_tb_mynn(), scale_atmos_refstate::atmos_refstate_update(), mod_atmos_vars::atmos_vars_history(), mod_atmos_vars::atmos_vars_monitor(), mod_atmos_vars::atmos_vars_total(), scale_atmos_phy_cp_kf::cp_kf_main(), mod_cpl_vars::cpl_getsfc_atm(), mod_cpl_vars::cpl_putatm(), mod_cpl_vars::cpl_vars_setup(), mod_mkinit::interporation_fact(), scale_atmos_phy_mp_sn14::mp_negativefilter(), scale_grid_nest::nest_setup(), mod_realinput_grads::parentatominputgrads(), mod_realinput_nicam::parentatominputnicam(), mod_realinput_wrfarw::parentatominputwrfarw(), mod_realinput::parentatomsetup(), and mod_mkinit::read_sounding().

  integer, public            :: I_QV = -1

n_hyd

integer, parameter, public scale_atmos_hydrometeor::n_hyd = 6

Definition at line 70 of file scale_atmos_sub_hydrometeor.F90.

Referenced by mod_atmos_phy_mp_driver::atmos_phy_mp_driver(), scale_atmos_phy_mp_dry::atmos_phy_mp_dry_effectiveradius(), scale_atmos_phy_mp_dry::atmos_phy_mp_dry_mixingratio(), scale_atmos_phy_mp_kessler::atmos_phy_mp_kessler_effectiveradius(), scale_atmos_phy_mp_kessler::atmos_phy_mp_kessler_mixingratio(), scale_atmos_phy_mp_sdm::atmos_phy_mp_sdm_effectiveradius(), scale_atmos_phy_mp_sdm::atmos_phy_mp_sdm_mixingratio(), scale_atmos_phy_mp_sn14::atmos_phy_mp_sn14_effectiveradius(), scale_atmos_phy_mp_sn14::atmos_phy_mp_sn14_mixingratio(), scale_atmos_phy_mp_suzuki10::atmos_phy_mp_suzuki10_effectiveradius(), scale_atmos_phy_mp_suzuki10::atmos_phy_mp_suzuki10_mixingratio(), scale_atmos_phy_mp_tomita08::atmos_phy_mp_tomita08_effectiveradius(), scale_atmos_phy_mp_tomita08::atmos_phy_mp_tomita08_mixingratio(), scale_atmos_phy_rd_mstrnx::atmos_phy_rd_mstrnx(), and scale_atmos_phy_rd_mstrnx::atmos_phy_rd_mstrnx_setup().

  integer, public, parameter :: N_HYD = 6

i_hc

integer, parameter, public scale_atmos_hydrometeor::i_hc = 1

liquid water cloud

Definition at line 72 of file scale_atmos_sub_hydrometeor.F90.

Referenced by mod_atmos_phy_mp_driver::atmos_phy_mp_driver(), scale_atmos_phy_mp_kessler::atmos_phy_mp_kessler_setup(), scale_atmos_phy_mp_sn14::atmos_phy_mp_sn14_setup(), scale_atmos_phy_mp_suzuki10::atmos_phy_mp_suzuki10_effectiveradius(), scale_atmos_phy_mp_suzuki10::atmos_phy_mp_suzuki10_mixingratio(), scale_atmos_phy_mp_suzuki10::atmos_phy_mp_suzuki10_setup(), scale_atmos_phy_mp_tomita08::atmos_phy_mp_tomita08_setup(), and scale_atmos_phy_rd_mstrnx::atmos_phy_rd_mstrnx().

  integer, public, parameter :: I_HC =  1

i_hr

integer, parameter, public scale_atmos_hydrometeor::i_hr = 2

liquid water rain

Definition at line 73 of file scale_atmos_sub_hydrometeor.F90.

Referenced by mod_atmos_phy_mp_driver::atmos_phy_mp_driver(), scale_atmos_phy_mp_kessler::atmos_phy_mp_kessler_setup(), scale_atmos_phy_mp_sn14::atmos_phy_mp_sn14_setup(), scale_atmos_phy_mp_suzuki10::atmos_phy_mp_suzuki10_effectiveradius(), scale_atmos_phy_mp_suzuki10::atmos_phy_mp_suzuki10_mixingratio(), scale_atmos_phy_mp_suzuki10::atmos_phy_mp_suzuki10_setup(), and scale_atmos_phy_mp_tomita08::atmos_phy_mp_tomita08_setup().

  integer, public, parameter :: I_HR =  2

i_hi

integer, parameter, public scale_atmos_hydrometeor::i_hi = 3

ice water cloud

Definition at line 74 of file scale_atmos_sub_hydrometeor.F90.

Referenced by mod_atmos_phy_mp_driver::atmos_phy_mp_driver(), scale_atmos_phy_mp_sn14::atmos_phy_mp_sn14_setup(), scale_atmos_phy_mp_suzuki10::atmos_phy_mp_suzuki10_effectiveradius(), scale_atmos_phy_mp_suzuki10::atmos_phy_mp_suzuki10_mixingratio(), scale_atmos_phy_mp_suzuki10::atmos_phy_mp_suzuki10_setup(), scale_atmos_phy_mp_tomita08::atmos_phy_mp_tomita08_setup(), and scale_atmos_phy_rd_mstrnx::atmos_phy_rd_mstrnx().

  integer, public, parameter :: I_HI =  3

i_hs

integer, parameter, public scale_atmos_hydrometeor::i_hs = 4

snow

Definition at line 75 of file scale_atmos_sub_hydrometeor.F90.

Referenced by mod_atmos_phy_mp_driver::atmos_phy_mp_driver(), scale_atmos_phy_mp_sn14::atmos_phy_mp_sn14_setup(), scale_atmos_phy_mp_suzuki10::atmos_phy_mp_suzuki10_effectiveradius(), scale_atmos_phy_mp_suzuki10::atmos_phy_mp_suzuki10_mixingratio(), scale_atmos_phy_mp_suzuki10::atmos_phy_mp_suzuki10_setup(), and scale_atmos_phy_mp_tomita08::atmos_phy_mp_tomita08_setup().

  integer, public, parameter :: I_HS =  4

i_hg

integer, parameter, public scale_atmos_hydrometeor::i_hg = 5

graupel

Definition at line 76 of file scale_atmos_sub_hydrometeor.F90.

Referenced by mod_atmos_phy_mp_driver::atmos_phy_mp_driver(), scale_atmos_phy_mp_sn14::atmos_phy_mp_sn14_setup(), scale_atmos_phy_mp_suzuki10::atmos_phy_mp_suzuki10_effectiveradius(), scale_atmos_phy_mp_suzuki10::atmos_phy_mp_suzuki10_mixingratio(), scale_atmos_phy_mp_suzuki10::atmos_phy_mp_suzuki10_setup(), and scale_atmos_phy_mp_tomita08::atmos_phy_mp_tomita08_setup().

  integer, public, parameter :: I_HG =  5

i_hh

integer, parameter, public scale_atmos_hydrometeor::i_hh = 6

hail

Definition at line 77 of file scale_atmos_sub_hydrometeor.F90.

Referenced by mod_atmos_phy_mp_driver::atmos_phy_mp_driver(), scale_atmos_phy_mp_suzuki10::atmos_phy_mp_suzuki10_effectiveradius(), scale_atmos_phy_mp_suzuki10::atmos_phy_mp_suzuki10_mixingratio(), and scale_atmos_phy_mp_suzuki10::atmos_phy_mp_suzuki10_setup().

  integer, public, parameter :: I_HH =  6

i_qc

integer, public scale_atmos_hydrometeor::i_qc = -1

Definition at line 79 of file scale_atmos_sub_hydrometeor.F90.

Referenced by scale_atmos_adiabat::atmos_adiabat_cape(), scale_atmos_adiabat::atmos_adiabat_liftparcel(), atmos_hydrometeor_diagnose_number_concentration(), scale_atmos_hydrostatic::atmos_hydrostatic_setup(), scale_atmos_phy_cp_kf::atmos_phy_cp_kf_setup(), scale_atmos_phy_mp_kessler::atmos_phy_mp_kessler(), scale_atmos_phy_mp_sn14::atmos_phy_mp_sn14(), scale_atmos_phy_mp_sn14::atmos_phy_mp_sn14_cloudfraction(), scale_atmos_phy_mp_tomita08::atmos_phy_mp_tomita08(), mod_atmos_phy_rd_driver::atmos_phy_rd_driver(), scale_atmos_phy_tb_mynn::atmos_phy_tb_mynn(), mod_atmos_vars::atmos_vars_history(), scale_atmos_phy_mp_sn14::aut_acc_slc_brk_kij(), scale_atmos_phy_cp_kf::cp_kf_main(), scale_atmos_phy_mp_sn14::freezing_water_kij(), scale_atmos_phy_mp_sn14::ice_multiplication_kij(), scale_atmos_phy_mp_sn14::mixed_phase_collection_kij(), scale_atmos_phy_mp_sn14::mp_negativefilter(), mod_realinput_grads::parentatominputgrads(), mod_realinput_wrfarw::parentatominputwrfarw(), mod_realinput::parentatomsetup(), mod_mkinit::read_sounding(), and scale_atmos_phy_mp_sn14::update_by_phase_change_kij().

  integer, public            :: I_QC = -1

i_qr

integer, public scale_atmos_hydrometeor::i_qr = -1

Definition at line 80 of file scale_atmos_sub_hydrometeor.F90.

Referenced by atmos_hydrometeor_diagnose_number_concentration(), scale_atmos_phy_cp_kf::atmos_phy_cp_kf_setup(), scale_atmos_phy_mp_kessler::atmos_phy_mp_kessler(), scale_atmos_phy_mp_sn14::atmos_phy_mp_sn14(), scale_atmos_phy_mp_tomita08::atmos_phy_mp_tomita08(), mod_atmos_phy_rd_driver::atmos_phy_rd_driver(), scale_atmos_phy_mp_sn14::aut_acc_slc_brk_kij(), scale_atmos_phy_cp_kf::cp_kf_main(), scale_atmos_phy_mp_sn14::freezing_water_kij(), scale_atmos_phy_mp_sn14::mixed_phase_collection_kij(), scale_atmos_phy_mp_sn14::mp_negativefilter(), mod_realinput_grads::parentatominputgrads(), mod_realinput_wrfarw::parentatominputwrfarw(), and scale_atmos_phy_mp_sn14::update_by_phase_change_kij().

  integer, public            :: I_QR = -1

i_qi

integer, public scale_atmos_hydrometeor::i_qi = -1

Definition at line 81 of file scale_atmos_sub_hydrometeor.F90.

Referenced by atmos_hydrometeor_diagnose_number_concentration(), scale_atmos_phy_cp_kf::atmos_phy_cp_kf_setup(), scale_atmos_phy_mp_sn14::atmos_phy_mp_sn14(), scale_atmos_phy_mp_tomita08::atmos_phy_mp_tomita08(), mod_atmos_phy_rd_driver::atmos_phy_rd_driver(), scale_atmos_phy_tb_mynn::atmos_phy_tb_mynn(), scale_atmos_phy_cp_kf::cp_kf_main(), scale_atmos_phy_mp_sn14::dep_vapor_melt_ice_kij(), scale_atmos_phy_mp_sn14::freezing_water_kij(), scale_atmos_phy_mp_sn14::mixed_phase_collection_kij(), scale_atmos_phy_mp_sn14::mp_negativefilter(), mod_realinput_grads::parentatominputgrads(), mod_realinput_wrfarw::parentatominputwrfarw(), and scale_atmos_phy_mp_sn14::update_by_phase_change_kij().

  integer, public            :: I_QI = -1

i_qs

integer, public scale_atmos_hydrometeor::i_qs = -1

Definition at line 82 of file scale_atmos_sub_hydrometeor.F90.

Referenced by atmos_hydrometeor_diagnose_number_concentration(), scale_atmos_phy_cp_kf::atmos_phy_cp_kf_setup(), scale_atmos_phy_mp_sn14::atmos_phy_mp_sn14(), scale_atmos_phy_mp_tomita08::atmos_phy_mp_tomita08(), mod_atmos_phy_rd_driver::atmos_phy_rd_driver(), scale_atmos_phy_cp_kf::cp_kf_main(), scale_atmos_phy_mp_sn14::dep_vapor_melt_ice_kij(), scale_atmos_phy_mp_sn14::freezing_water_kij(), scale_atmos_phy_mp_sn14::mixed_phase_collection_kij(), scale_atmos_phy_mp_sn14::mp_negativefilter(), mod_realinput_grads::parentatominputgrads(), mod_realinput_wrfarw::parentatominputwrfarw(), and scale_atmos_phy_mp_sn14::update_by_phase_change_kij().

  integer, public            :: I_QS = -1

i_qg

integer, public scale_atmos_hydrometeor::i_qg = -1

Definition at line 83 of file scale_atmos_sub_hydrometeor.F90.

Referenced by atmos_hydrometeor_diagnose_number_concentration(), scale_atmos_phy_mp_sn14::atmos_phy_mp_sn14(), scale_atmos_phy_mp_sn14::atmos_phy_mp_sn14_cloudfraction(), scale_atmos_phy_mp_tomita08::atmos_phy_mp_tomita08(), mod_atmos_phy_rd_driver::atmos_phy_rd_driver(), scale_atmos_phy_mp_sn14::dep_vapor_melt_ice_kij(), scale_atmos_phy_mp_sn14::freezing_water_kij(), scale_atmos_phy_mp_sn14::mixed_phase_collection_kij(), scale_atmos_phy_mp_sn14::mp_negativefilter(), mod_realinput_grads::parentatominputgrads(), mod_realinput_wrfarw::parentatominputwrfarw(), and scale_atmos_phy_mp_sn14::update_by_phase_change_kij().

  integer, public            :: I_QG = -1

i_qh

integer, public scale_atmos_hydrometeor::i_qh = -1

Definition at line 84 of file scale_atmos_sub_hydrometeor.F90.

  integer, public            :: I_QH = -1

i_nc

integer, public scale_atmos_hydrometeor::i_nc = -1

Definition at line 86 of file scale_atmos_sub_hydrometeor.F90.

Referenced by atmos_hydrometeor_diagnose_number_concentration(), scale_atmos_phy_mp_sn14::atmos_phy_mp_sn14(), scale_atmos_phy_mp_sn14::aut_acc_slc_brk_kij(), scale_atmos_phy_mp_sn14::dep_vapor_melt_ice_kij(), scale_atmos_phy_mp_sn14::freezing_water_kij(), scale_atmos_phy_mp_sn14::ice_multiplication_kij(), scale_atmos_phy_mp_sn14::mixed_phase_collection_kij(), scale_atmos_phy_mp_sn14::mp_negativefilter(), scale_atmos_phy_mp_sn14::nucleation_kij(), mod_realinput_wrfarw::parentatominputwrfarw(), mod_mkinit::read_sounding(), and scale_atmos_phy_mp_sn14::update_by_phase_change_kij().

  integer, public            :: I_NC = -1

i_nr

integer, public scale_atmos_hydrometeor::i_nr = -1

Definition at line 87 of file scale_atmos_sub_hydrometeor.F90.

Referenced by atmos_hydrometeor_diagnose_number_concentration(), scale_atmos_phy_mp_sn14::atmos_phy_mp_sn14(), scale_atmos_phy_mp_sn14::aut_acc_slc_brk_kij(), scale_atmos_phy_mp_sn14::dep_vapor_melt_ice_kij(), scale_atmos_phy_mp_sn14::freezing_water_kij(), scale_atmos_phy_mp_sn14::mixed_phase_collection_kij(), scale_atmos_phy_mp_sn14::mp_negativefilter(), mod_realinput_wrfarw::parentatominputwrfarw(), and scale_atmos_phy_mp_sn14::update_by_phase_change_kij().

  integer, public            :: I_NR = -1

i_ni

integer, public scale_atmos_hydrometeor::i_ni = -1

Definition at line 88 of file scale_atmos_sub_hydrometeor.F90.

Referenced by atmos_hydrometeor_diagnose_number_concentration(), scale_atmos_phy_mp_sn14::atmos_phy_mp_sn14(), scale_atmos_phy_mp_sn14::dep_vapor_melt_ice_kij(), scale_atmos_phy_mp_sn14::freezing_water_kij(), scale_atmos_phy_mp_sn14::mixed_phase_collection_kij(), scale_atmos_phy_mp_sn14::mp_negativefilter(), scale_atmos_phy_mp_sn14::nucleation_kij(), mod_realinput_wrfarw::parentatominputwrfarw(), and scale_atmos_phy_mp_sn14::update_by_phase_change_kij().

  integer, public            :: I_NI = -1

i_ns

integer, public scale_atmos_hydrometeor::i_ns = -1

Definition at line 89 of file scale_atmos_sub_hydrometeor.F90.

Referenced by atmos_hydrometeor_diagnose_number_concentration(), scale_atmos_phy_mp_sn14::atmos_phy_mp_sn14(), scale_atmos_phy_mp_sn14::dep_vapor_melt_ice_kij(), scale_atmos_phy_mp_sn14::freezing_water_kij(), scale_atmos_phy_mp_sn14::mixed_phase_collection_kij(), scale_atmos_phy_mp_sn14::mp_negativefilter(), mod_realinput_wrfarw::parentatominputwrfarw(), and scale_atmos_phy_mp_sn14::update_by_phase_change_kij().

  integer, public            :: I_NS = -1

i_ng

integer, public scale_atmos_hydrometeor::i_ng = -1

Definition at line 90 of file scale_atmos_sub_hydrometeor.F90.

Referenced by atmos_hydrometeor_diagnose_number_concentration(), scale_atmos_phy_mp_sn14::atmos_phy_mp_sn14(), scale_atmos_phy_mp_sn14::mp_negativefilter(), and mod_realinput_wrfarw::parentatominputwrfarw().

  integer, public            :: I_NG = -1

i_nh

integer, public scale_atmos_hydrometeor::i_nh = -1

Definition at line 91 of file scale_atmos_sub_hydrometeor.F90.

  integer, public            :: I_NH = -1

qhs

integer, public scale_atmos_hydrometeor::qhs = -1

Definition at line 94 of file scale_atmos_sub_hydrometeor.F90.

Referenced by atmos_hydrometeor_regist(), scale_atmos_phy_mp_common::atmos_phy_mp_negative_fixer(), mod_atmos_vars::atmos_vars_history(), and mod_mkinit::rect_setup().

  integer, public            :: QHS = -1

qhe

integer, public scale_atmos_hydrometeor::qhe = -2

Definition at line 95 of file scale_atmos_sub_hydrometeor.F90.

Referenced by atmos_hydrometeor_regist(), scale_atmos_phy_mp_common::atmos_phy_mp_negative_fixer(), mod_atmos_vars::atmos_vars_history(), mod_mkinit::read_sounding(), and mod_mkinit::rect_setup().

  integer, public            :: QHE = -2

qls

integer, public scale_atmos_hydrometeor::qls = -1

Definition at line 97 of file scale_atmos_sub_hydrometeor.F90.

Referenced by atmos_hydrometeor_entr_0d(), atmos_hydrometeor_entr_2d(), atmos_hydrometeor_entr_3d(), atmos_hydrometeor_regist(), mod_atmos_vars::atmos_vars_history(), and mod_realinput::parentatomsetup().

  integer, public            :: QLS = -1

qle

integer, public scale_atmos_hydrometeor::qle = -2

Definition at line 98 of file scale_atmos_sub_hydrometeor.F90.

Referenced by atmos_hydrometeor_entr_0d(), atmos_hydrometeor_entr_2d(), atmos_hydrometeor_entr_3d(), atmos_hydrometeor_regist(), mod_atmos_vars::atmos_vars_history(), and mod_realinput::parentatomsetup().

  integer, public            :: QLE = -2

qis

integer, public scale_atmos_hydrometeor::qis = -1

Definition at line 100 of file scale_atmos_sub_hydrometeor.F90.

Referenced by atmos_hydrometeor_entr_0d(), atmos_hydrometeor_entr_2d(), atmos_hydrometeor_entr_3d(), atmos_hydrometeor_regist(), mod_atmos_vars::atmos_vars_history(), mod_atmos_vars::atmos_vars_monitor(), and mod_atmos_vars::atmos_vars_total().

  integer, public            :: QIS = -1

qie

integer, public scale_atmos_hydrometeor::qie = -2

Definition at line 101 of file scale_atmos_sub_hydrometeor.F90.

Referenced by atmos_hydrometeor_entr_0d(), atmos_hydrometeor_entr_2d(), atmos_hydrometeor_entr_3d(), atmos_hydrometeor_regist(), mod_atmos_vars::atmos_vars_history(), mod_atmos_vars::atmos_vars_monitor(), and mod_atmos_vars::atmos_vars_total().

  integer, public            :: QIE = -2

lhv

real(rp), public scale_atmos_hydrometeor::lhv

latent heat of vaporization for use [J/kg]

Definition at line 103 of file scale_atmos_sub_hydrometeor.F90.

Referenced by atmos_hydrometeor_setup(), scale_atmos_phy_mp_kessler::atmos_phy_mp_kessler(), scale_atmos_phy_mp_common::atmos_phy_mp_saturation_adjustment(), scale_atmos_phy_mp_tomita08::atmos_phy_mp_tomita08(), mod_atmos_vars::atmos_vars_history(), mod_atmos_vars::atmos_vars_monitor(), and mod_atmos_vars::atmos_vars_total().

  real(RP), public           :: LHV

lhs

real(rp), public scale_atmos_hydrometeor::lhs

latent heat of sublimation for use [J/kg]

Definition at line 104 of file scale_atmos_sub_hydrometeor.F90.

Referenced by atmos_hydrometeor_setup().

  real(RP), public           :: LHS

lhf

real(rp), public scale_atmos_hydrometeor::lhf

latent heat of fusion for use [J/kg]

Definition at line 105 of file scale_atmos_sub_hydrometeor.F90.

Referenced by atmos_hydrometeor_setup(), scale_atmos_phy_mp_common::atmos_phy_mp_saturation_adjustment(), scale_atmos_phy_mp_tomita08::atmos_phy_mp_tomita08(), mod_atmos_vars::atmos_vars_history(), mod_atmos_vars::atmos_vars_monitor(), and mod_atmos_vars::atmos_vars_total().

  real(RP), public           :: LHF