scale_atmos_phy_mp_kessler Module Reference

module atmosphere / physics / microphysics / Kessler More...

Functions/Subroutines
subroutine, public	atmos_phy_mp_kessler_setup
	ATMOS_PHY_MP_kessler_setup Setup. More...
subroutine, public	atmos_phy_mp_kessler_adjustment (KA, KS, KE, IA, IS, IE, JA, JS, JE, DENS, PRES, dt, TEMP, QTRC, CPtot, CVtot, RHOE_t, EVAPORATE)
	ATMOS_PHY_MP_kessler_adjustment calculate state after saturation process. More...
subroutine, public	atmos_phy_mp_kessler_terminal_velocity (KA, KS, KE, DENS0, RHOQ0, REFSTATE_dens_profile, vterm)
	Kessler-type warm rain microphysics (terminal velocity) More...
subroutine, public	atmos_phy_mp_kessler_cloud_fraction (KA, KS, KE, IA, IS, IE, JA, JS, JE, QTRC, mask_criterion, cldfrac)
	Calculate Cloud Fraction. More...
subroutine, public	atmos_phy_mp_kessler_effective_radius (KA, KS, KE, IA, IS, IE, JA, JS, JE, DENS0, TEMP0, QTRC0, Re)
	Calculate Effective Radius. More...
subroutine, public	atmos_phy_mp_kessler_qtrc2qhyd (KA, KS, KE, IA, IS, IE, JA, JS, JE, QTRC, Qe)
	Calculate mass ratio of each category. More...
subroutine, public	atmos_phy_mp_kessler_qhyd2qtrc (KA, KS, KE, IA, IS, IE, JA, JS, JE, Qe, QTRC)
	Calculate mass ratio of each category. More...

Variables
integer, parameter, public	atmos_phy_mp_kessler_ntracers = QA_MP
integer, parameter, public	atmos_phy_mp_kessler_nwaters = 2
integer, parameter, public	atmos_phy_mp_kessler_nices = 0
character(len=h_short), dimension(qa_mp), parameter, public	atmos_phy_mp_kessler_tracer_names = (/ 'QV', 'QC', 'QR' /)
character(len=h_mid), dimension(qa_mp), parameter, public	atmos_phy_mp_kessler_tracer_descriptions = (/ 'Ratio of Water Vapor mass to total mass (Specific humidity)', 'Ratio of Cloud Water mass to total mass ', 'Ratio of Rain Water mass to total mass '/)
character(len=h_short), dimension(qa_mp), parameter, public	atmos_phy_mp_kessler_tracer_units = (/ 'kg/kg', 'kg/kg', 'kg/kg' /)

Detailed Description

module atmosphere / physics / microphysics / Kessler

Description

Cloud Microphysics by Kessler-type parametarization Reference: Kessler(1969) Klemp and Wilhelmson(1978)

Author

Team SCALE

NAMELIST

• No namelist group

History Output

name	description	unit	variable
Pcsat	QC production term by satadjust	kg/kg/s	QC_t_sat

Function/Subroutine Documentation

atmos_phy_mp_kessler_setup()

subroutine, public scale_atmos_phy_mp_kessler::atmos_phy_mp_kessler_setup

ATMOS_PHY_MP_kessler_setup Setup.

Definition at line 90 of file scale_atmos_phy_mp_kessler.F90.

    use scale_prc, only: &
       prc_abort
    implicit none
 
    !---------------------------------------------------------------------------
 
    log_newline
    log_info("ATMOS_PHY_MP_kessler_setup",*) 'Setup'
    log_info("ATMOS_PHY_MP_kessler_setup",*) 'KESSLER-type 1-moment bulk 3 category'
 
    if( couple_aerosol ) then
       log_error("ATMOS_PHY_MP_kessler_setup",*) 'MP_aerosol_couple should be .false. for KESSLER type MP!'
       call prc_abort
    endif
 
    return

References scale_prc::prc_abort().

Referenced by mod_atmos_phy_mp_driver::atmos_phy_mp_driver_setup().

Here is the call graph for this function:

Here is the caller graph for this function:

atmos_phy_mp_kessler_adjustment()

subroutine, public scale_atmos_phy_mp_kessler::atmos_phy_mp_kessler_adjustment	(	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), intent(in)	DENS,
		real(rp), dimension (ka,ia,ja), intent(in)	PRES,
		real(dp), intent(in)	dt,
		real(rp), dimension (ka,ia,ja), intent(inout)	TEMP,
		real(rp), dimension (ka,ia,ja,qa_mp), intent(inout)	QTRC,
		real(rp), dimension (ka,ia,ja), intent(inout)	CPtot,
		real(rp), dimension (ka,ia,ja), intent(inout)	CVtot,
		real(rp), dimension (ka,ia,ja), intent(out)	RHOE_t,
		real(rp), dimension(ka,ia,ja), intent(out)	EVAPORATE
	)

ATMOS_PHY_MP_kessler_adjustment calculate state after saturation process.

Definition at line 118 of file scale_atmos_phy_mp_kessler.F90.

    use scale_const, only: &
       dwatr => const_dwatr, &
       pi    => const_pi
    use scale_file_history, only: &
       file_history_in
    use scale_atmos_phy_mp_common, only: &
       mp_saturation_adjustment => atmos_phy_mp_saturation_adjustment
    implicit none
 
    integer,  intent(in)    :: KA, KS, KE
    integer,  intent(in)    :: IA, IS, IE
    integer,  intent(in)    :: JA, JS, JE
    real(RP), intent(in)    :: DENS     (KA,IA,JA)       ! density [kg/m3]
    real(RP), intent(in)    :: PRES     (KA,IA,JA)       ! pressure [Pa]
    real(DP), intent(in)    :: dt                        ! time interval of microphysics [s]
    real(RP), intent(inout) :: TEMP     (KA,IA,JA)       ! temperature [K]
    real(RP), intent(inout) :: QTRC     (KA,IA,JA,QA_MP) ! tracer mass concentration [kg/kg]
    real(RP), intent(inout) :: CPtot    (KA,IA,JA)       ! total specific heat capacity at constant pressure [J/kg/K]
    real(RP), intent(inout) :: CVtot    (KA,IA,JA)       ! total specific heat capacity at constant volume [J/kg/K]
    real(RP), intent(out)   :: RHOE_t   (KA,IA,JA)       ! tendency of rhoe [J/m3/s]
    real(RP), intent(out)   :: EVAPORATE(KA,IA,JA)       ! number of evaporated cloud [#/m3/s]
 
    real(RP) :: QTRC_dummy(KA,IA,JA)
    real(RP) :: RHOE_d_sat(KA,IA,JA)
    real(RP) :: QC_t_sat  (KA,IA,JA)
 
    integer  :: k, i, j
    !---------------------------------------------------------------------------
 
    log_progress(*) 'atmosphere / physics / microphysics / Kessler'
 
    !##### MP Main #####
    call mp_kessler( &
         ka, ks, ke, ia, is, ie, ja, js, je, &
         dens(:,:,:), pres(:,:,:),   & ! [IN]
         dt,                         & ! [IN]
         temp(:,:,:), qtrc(:,:,:,:), & ! [INOUT]
         cptot(:,:,:), cvtot(:,:,:), & ! [INOUT]
         rhoe_t(:,:,:)               ) ! [OUT]
 
    ! save value before saturation adjustment
    do j = js, je
    do i = is, ie
    do k = ks, ke
       qc_t_sat(k,i,j) = qtrc(k,i,j,i_qc)
    enddo
    enddo
    enddo
 
!OCL XFILL
    qtrc_dummy(:,:,:) = -1.0_rp
 
    call mp_saturation_adjustment( &
         ka, ks, ke, ia, is, ie, ja, js, je,  &
         dens(:,:,:),                         & ! [IN]
         flag_liquid,                         & ! [IN]
         temp(:,:,:),                         & ! [INOUT]
         qtrc(:,:,:,i_qv),                    & ! [INOUT]
         qtrc(:,:,:,i_qc), qtrc_dummy(:,:,:), & ! [INOUT]
         cptot(:,:,:), cvtot(:,:,:),          & ! [INOUT]
         rhoe_d_sat(:,:,:)                    ) ! [OUT]
 
    do j = js, je
    do i = is, ie
    do k = ks, ke
       rhoe_t(k,i,j) = rhoe_t(k,i,j) + rhoe_d_sat(k,i,j) / dt
    enddo
    enddo
    enddo
    do j = js, je
    do i = is, ie
    do k = ks, ke
       qc_t_sat(k,i,j) = ( qtrc(k,i,j,i_qc) - qc_t_sat(k,i,j) ) / dt
    enddo
    enddo
    enddo
 
    call file_history_in( qc_t_sat(:,:,:), 'Pcsat', 'QC production term by satadjust', 'kg/kg/s' )
 
    do j = js, je
    do i = is, ie
    do k = ks, ke
       evaporate(k,i,j) = max( -qc_t_sat(k,i,j), 0.0_rp ) & ! if negative, condensation
                        * dens(k,i,j) / (4.0_rp/3.0_rp*pi*dwatr*re_qc**3) ! mass -> number (assuming constant particle radius as re_qc)
    enddo
    enddo
    enddo
 
    !##### END MP Main #####
 
    return

References scale_const::const_dwatr, scale_const::const_eps, scale_const::const_pi, scale_atmos_hydrometeor::cp_vapor, scale_atmos_hydrometeor::cp_water, scale_atmos_hydrometeor::cv_vapor, scale_atmos_hydrometeor::cv_water, scale_atmos_hydrometeor::lhv, scale_prof::prof_rapend(), and scale_prof::prof_rapstart().

Referenced by mod_atmos_phy_mp_driver::atmos_phy_mp_driver_calc_tendency().

Here is the call graph for this function:

Here is the caller graph for this function:

atmos_phy_mp_kessler_terminal_velocity()

subroutine, public scale_atmos_phy_mp_kessler::atmos_phy_mp_kessler_terminal_velocity	(	integer, intent(in)	KA,
		integer, intent(in)	KS,
		integer, intent(in)	KE,
		real(rp), dimension (ka), intent(in)	DENS0,
		real(rp), dimension (ka,qa_mp-1), intent(in)	RHOQ0,
		real(rp), dimension(ka), intent(in)	REFSTATE_dens_profile,
		real(rp), dimension (ka,qa_mp-1), intent(out)	vterm
	)

Kessler-type warm rain microphysics (terminal velocity)

Definition at line 356 of file scale_atmos_phy_mp_kessler.F90.

    implicit none
 
    integer,  intent(in)  :: KA, KS, KE
    real(RP), intent(in)  :: DENS0                (KA)         ! density [kg/m3]
    real(RP), intent(in)  :: RHOQ0                (KA,QA_MP-1) ! density of each hydrometeor [kg/m3]
    real(RP), intent(in)  :: REFSTATE_dens_profile(KA)         ! reference state density profile [kg/m3]
    real(RP), intent(out) :: vterm                (KA,QA_MP-1) ! terminal velocity of each hydrometeor [m/s]
 
    real(RP) :: qr
    real(RP) :: zerosw
 
    integer  :: k
    !---------------------------------------------------------------------------
 
    do k = ks, ke
       vterm(k,i_hyd_qc) = 0.0_rp
    enddo
 
    do k = ks, ke
       qr     = rhoq0(k,i_hyd_qr) / dens0(k)
       zerosw = 0.5_rp - sign(0.5_rp, qr - 1.e-12_rp )
 
       vterm(k,i_hyd_qr) = - 36.34_rp * ( dens0(k) * ( qr + zerosw ) )**0.1364_rp &
                         * refstate_dens_profile(ks) / refstate_dens_profile(k) * ( 1.0_rp - zerosw )
    enddo
 
    return

Referenced by mod_atmos_phy_mp_driver::atmos_phy_mp_driver_calc_tendency().

Here is the caller graph for this function:

atmos_phy_mp_kessler_cloud_fraction()

subroutine, public scale_atmos_phy_mp_kessler::atmos_phy_mp_kessler_cloud_fraction	(	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,qa_mp-1), intent(in)	QTRC,
		real(rp), intent(in)	mask_criterion,
		real(rp), dimension(ka,ia,ja), intent(out)	cldfrac
	)

Calculate Cloud Fraction.

Definition at line 392 of file scale_atmos_phy_mp_kessler.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)  :: QTRC(KA,IA,JA,QA_MP-1) ! hydrometeor mass concentration [kg/kg]
    real(RP), intent(in)  :: mask_criterion         ! criterion of hydrometeor [kg/kg]
 
    real(RP), intent(out) :: cldfrac(KA,IA,JA)      ! cloud fraction (0 or 1)
 
    real(RP) :: qhydro
    integer  :: k, i, j
    !---------------------------------------------------------------------------
 
    do j  = js, je
    do i  = is, ie
    do k  = ks, ke
       qhydro = qtrc(k,i,j,i_hyd_qc) + qtrc(k,i,j,i_hyd_qr)
       cldfrac(k,i,j) = 0.5_rp + sign(0.5_rp,qhydro-mask_criterion)
    enddo
    enddo
    enddo
 
    return

Referenced by mod_atmos_phy_mp_vars::atmos_phy_mp_vars_get_diagnostic().

Here is the caller graph for this function:

atmos_phy_mp_kessler_effective_radius()

subroutine, public scale_atmos_phy_mp_kessler::atmos_phy_mp_kessler_effective_radius	(	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), intent(in)	DENS0,
		real(rp), dimension(ka,ia,ja), intent(in)	TEMP0,
		real(rp), dimension(ka,ia,ja,qa_mp-1), intent(in)	QTRC0,
		real(rp), dimension (ka,ia,ja,n_hyd), intent(out)	Re
	)

Calculate Effective Radius.

Definition at line 424 of file scale_atmos_phy_mp_kessler.F90.

    use scale_atmos_hydrometeor, only: &
       n_hyd, &
       i_hc, &
       i_hr
    implicit none
    integer, intent(in) :: KA, KS, KE
    integer, intent(in) :: IA, IS, IE
    integer, intent(in) :: JA, JS, JE
 
    real(RP), intent(in)  :: DENS0(KA,IA,JA)         ! density [kg/m3]
    real(RP), intent(in)  :: TEMP0(KA,IA,JA)         ! temperature [K]
    real(RP), intent(in)  :: QTRC0(KA,IA,JA,QA_MP-1) ! hydrometeor mass concentration [kg/kg]
 
    real(RP), intent(out) :: Re   (KA,IA,JA,N_HYD)   ! effective radius [cm]
 
    real(RP), parameter :: um2cm = 100.0_rp
    !---------------------------------------------------------------------------
 
!OCL XFILL
    re(:,:,:,i_hc) =   8.e-6_rp * um2cm
!OCL XFILL
    re(:,:,:,i_hr) = 100.e-6_rp * um2cm
!OCL XFILL
    re(:,:,:,i_hr+1:) = 0.0_rp
 
    return

References scale_atmos_hydrometeor::i_hc, scale_atmos_hydrometeor::i_hr, and scale_atmos_hydrometeor::n_hyd.

Referenced by mod_atmos_phy_mp_vars::atmos_phy_mp_vars_get_diagnostic().

Here is the caller graph for this function:

atmos_phy_mp_kessler_qtrc2qhyd()

subroutine, public scale_atmos_phy_mp_kessler::atmos_phy_mp_kessler_qtrc2qhyd	(	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,qa_mp-1), intent(in)	QTRC,
		real(rp), dimension(ka,ia,ja,n_hyd), intent(out)	Qe
	)

Calculate mass ratio of each category.

Definition at line 458 of file scale_atmos_phy_mp_kessler.F90.

    use scale_atmos_hydrometeor, only: &
       n_hyd, &
       i_hc, &
       i_hr
    implicit none
    integer, intent(in) :: KA, KS, KE
    integer, intent(in) :: IA, IS, IE
    integer, intent(in) :: JA, JS, JE
 
    real(RP), intent(in)  :: QTRC(KA,IA,JA,QA_MP-1) ! hydrometeor mass concentration [kg/kg]
 
    real(RP), intent(out) :: Qe(KA,IA,JA,N_HYD)   ! mass ratio of each hydrometeor [kg/kg]
    !---------------------------------------------------------------------------
 
!OCL XFILL
    qe(:,:,:,i_hc) = qtrc(:,:,:,i_hyd_qc)
!OCL XFILL
    qe(:,:,:,i_hr) = qtrc(:,:,:,i_hyd_qr)
!OCL XFILL
    qe(:,:,:,i_hr+1:) = 0.0_rp
 
    return

References scale_atmos_hydrometeor::i_hc, scale_atmos_hydrometeor::i_hr, and scale_atmos_hydrometeor::n_hyd.

Referenced by mod_atmos_phy_mp_vars::atmos_phy_mp_vars_get_diagnostic().

Here is the caller graph for this function:

atmos_phy_mp_kessler_qhyd2qtrc()

subroutine, public scale_atmos_phy_mp_kessler::atmos_phy_mp_kessler_qhyd2qtrc	(	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,n_hyd), intent(in)	Qe,
		real(rp), dimension(ka,ia,ja,qa_mp-1), intent(out)	QTRC
	)

Calculate mass ratio of each category.

Definition at line 488 of file scale_atmos_phy_mp_kessler.F90.

    use scale_atmos_hydrometeor, only: &
       n_hyd, &
       i_hc, &
       i_hr
    implicit none
    integer, intent(in) :: KA, KS, KE
    integer, intent(in) :: IA, IS, IE
    integer, intent(in) :: JA, JS, JE
 
    real(RP), intent(in) :: Qe(KA,IA,JA,N_HYD)   ! mass ratio of each hydrometeor [kg/kg]
 
    real(RP), intent(out) :: QTRC(KA,IA,JA,QA_MP-1) ! hydrometeor mass concentration [kg/kg]
 
    !---------------------------------------------------------------------------
 
!OCL XFILL
    qtrc(:,:,:,i_hyd_qc) = qe(:,:,:,i_hc)
!OCL XFILL
    qtrc(:,:,:,i_hyd_qr) = qe(:,:,:,i_hr)
 
    ! ignore ice water
 
    return

References scale_atmos_hydrometeor::i_hc, scale_atmos_hydrometeor::i_hr, and scale_atmos_hydrometeor::n_hyd.

Referenced by mod_atmos_phy_mp_driver::atmos_phy_mp_driver_qhyd2qtrc().

Here is the caller graph for this function:

Variable Documentation

atmos_phy_mp_kessler_ntracers

integer, parameter, public scale_atmos_phy_mp_kessler::atmos_phy_mp_kessler_ntracers = QA_MP

Definition at line 44 of file scale_atmos_phy_mp_kessler.F90.

  integer,                parameter, public :: ATMOS_PHY_MP_kessler_ntracers = qa_mp

Referenced by mod_atmos_phy_mp_driver::atmos_phy_mp_driver_tracer_setup().

atmos_phy_mp_kessler_nwaters

integer, parameter, public scale_atmos_phy_mp_kessler::atmos_phy_mp_kessler_nwaters = 2

Definition at line 45 of file scale_atmos_phy_mp_kessler.F90.

  integer,                parameter, public :: ATMOS_PHY_MP_kessler_nwaters = 2

Referenced by mod_atmos_phy_mp_driver::atmos_phy_mp_driver_tracer_setup().

atmos_phy_mp_kessler_nices

integer, parameter, public scale_atmos_phy_mp_kessler::atmos_phy_mp_kessler_nices = 0

Definition at line 46 of file scale_atmos_phy_mp_kessler.F90.

  integer,                parameter, public :: ATMOS_PHY_MP_kessler_nices = 0

Referenced by mod_atmos_phy_mp_driver::atmos_phy_mp_driver_tracer_setup().

atmos_phy_mp_kessler_tracer_names

character(len=h_short), dimension(qa_mp), parameter, public scale_atmos_phy_mp_kessler::atmos_phy_mp_kessler_tracer_names = (/ 'QV', 'QC', 'QR' /)

Definition at line 47 of file scale_atmos_phy_mp_kessler.F90.

  character(len=H_SHORT), parameter, public :: ATMOS_PHY_MP_kessler_tracer_names(QA_MP) = (/ &
       'QV', &
       'QC', &
       'QR'  /)

Referenced by mod_atmos_phy_mp_driver::atmos_phy_mp_driver_tracer_setup().

atmos_phy_mp_kessler_tracer_descriptions

character(len=h_mid), dimension(qa_mp), parameter, public scale_atmos_phy_mp_kessler::atmos_phy_mp_kessler_tracer_descriptions = (/ 'Ratio of Water Vapor mass to total mass (Specific humidity)', 'Ratio of Cloud Water mass to total mass ', 'Ratio of Rain Water mass to total mass '/)

Definition at line 51 of file scale_atmos_phy_mp_kessler.F90.

  character(len=H_MID)  , parameter, public :: ATMOS_PHY_MP_kessler_tracer_descriptions(QA_MP) = (/ &
       'Ratio of Water Vapor mass to total mass (Specific humidity)', &
       'Ratio of Cloud Water mass to total mass                    ', &
       'Ratio of Rain Water mass to total mass                     '/)

Referenced by mod_atmos_phy_mp_driver::atmos_phy_mp_driver_tracer_setup().

atmos_phy_mp_kessler_tracer_units

character(len=h_short), dimension(qa_mp), parameter, public scale_atmos_phy_mp_kessler::atmos_phy_mp_kessler_tracer_units = (/ 'kg/kg', 'kg/kg', 'kg/kg' /)

Definition at line 55 of file scale_atmos_phy_mp_kessler.F90.

  character(len=H_SHORT), parameter, public :: ATMOS_PHY_MP_kessler_tracer_units(QA_MP) = (/ &
       'kg/kg',  &
       'kg/kg',  &
       'kg/kg'   /)

Referenced by mod_atmos_phy_mp_driver::atmos_phy_mp_driver_tracer_setup().