scale_atmos_phy_cp_kf Module Reference

module atmosphere / physics / cumulus / Kain-Fritsch More...

Functions/Subroutines
subroutine, public	atmos_phy_cp_kf_setup (KA, KS, KE, IA, IS, IE, JA, JS, JE, CZ, AREA, WARMRAIN_in)
	Setup initial setup for Kain-Fritsch Cumulus Parameterization. More...
subroutine, public	atmos_phy_cp_kf_tendency (KA, KS, KE, IA, IS, IE, JA, JS, JE, DENS, U, V, RHOT, TEMP, PRES, QDRY, QV_in, w0avg, FZ, KF_DTSECD, DENS_t, RHOT_t, RHOQV_t, RHOQ_t, nca, SFLX_rain, SFLX_snow, SFLX_engi, cloudtop, cloudbase, cldfrac_dp, cldfrac_sh)
	ATMOS_PHY_CP_kf calculate Kain-Fritsch Cumulus Parameterization. More...

Detailed Description

module atmosphere / physics / cumulus / Kain-Fritsch

Description

Main module of Kain-Fritsch Cumulus Convection Parameterization

Author

Team SCALE

References

• Kain J. S. and J. M. Fritsch, 1990: A one-dimensional entraining/detraining plume model and its application in convective parameterization. J. Atmos. Sci. 47:2784–2802.
• Kain J. S., 2004: The Kain–Fritsch Convective Parameterization: An Update, J. Appl. Meteor., 43, 170-181.
• Narita, M. and S. Ohmori, 2007: Improving Precipitation Forecasts by the Operational Nonhydrostatic Mesoscale Model with the Kain-Fritsch Convective Parameterization and Cloud Microphysics. 12th Conference on Mesoscale Processes, 6-9 August 2007, Waterville Valley conference & event center, Waterville Valley, NH, USA.
• Ogura, Y., and H. R. Cho, 1973: Diagnostic determination of cumulus cloud populations from observed large-scale variables. J. Atmos. Sci., 30, 1276-1286.

This file was originally copied from WRF. The original file was published with the following notice.

WRF was developed at the National Center for Atmospheric Research (NCAR) which is operated by the University Corporation for Atmospheric Research (UCAR). NCAR and UCAR make no proprietary claims, either statutory or otherwise, to this version and release of WRF and consider WRF to be in the public domain for use by any person or entity for any purpose without any fee or charge. UCAR requests that any WRF user include this notice on any partial or full copies of WRF. WRF is provided on an "AS IS" basis and any warranties, either express or implied, including but not limited to implied warranties of non-infringement, originality, merchantability and fitness for a particular purpose, are disclaimed. In no event shall UCAR be liable for any damages, whatsoever, whether direct, indirect, consequential or special, that arise out of or in connection with the access, use or performance of WRF, including infringement actions.

WRF® is a registered trademark of the University Corporation for Atmospheric Research (UCAR).

NAMELIST

• PARAM_ATMOS_PHY_CP_KF

  name	type	default value	comment
  ATMOS_PHY_CP_KF_RATE	real(RP)	0.03_RP	ratio of cloud water and precipitation (Ogura and Cho 1973)
  ATMOS_PHY_CP_KF_TRIGGER_TYPE	integer	1	trigger function type 1:KF2004 3:NO2007
  ATMOS_PHY_CP_KF_DLCAPE	real(RP)	0.1_RP	cape decleace rate
  ATMOS_PHY_CP_KF_DLIFETIME	real(RP)	1800.0_RP	minimum lifetime scale of deep convection [sec]
  ATMOS_PHY_CP_KF_SLIFETIME	real(RP)	2400.0_RP	lifetime of shallow convection [sec]
  ATMOS_PHY_CP_KF_DEPTH_USL	real(RP)	300.0_RP	depth of updraft source layer [hPa]
  ATMOS_PHY_CP_KF_PREC_TYPE	integer	1	precipitation type 1:Ogura-Cho(1973) 2:Kessler
  ATMOS_PHY_CP_KF_DO_PREC	logical	.true.	whether enable precipitation
  ATMOS_PHY_CP_KF_THRES	real(RP)	1.E-3_RP	autoconversion rate for Kessler
  ATMOS_PHY_CP_KF_LOG	logical	.false.	output log?

  
  

History Output

name	description	unit	variable
KF_CONVFLAG	CONVECTION FLAG		R_convflag
KF_LIFETIME	lifetime of KF scheme	s	lifetime
QC_t_KF_conv_DP	QC tendency due to flux convergence (deep convection) in KF	kg/m3/s	QC_t_KF_conv_DP
QC_t_KF_conv_SH	QC tendency due to flux convergence (shallow convection) in KF	kg/m3/s	QC_t_KF_conv_SH
QC_t_KF_updet_DP	QC tendency due to detrainment of updraft (deep convection) in KF	kg/m3/s	QC_t_KF_updet_DP
QC_t_KF_updet_SH	QC tendency due to detrainment of updraft (shallow convection) in KF	kg/m3/s	QC_t_KF_updet_SH
QI_t_KF_conv_DP	QI tendency due to flux convergence (deep convection) in KF	kg/m3/s	QI_t_KF_conv_DP
QI_t_KF_conv_SH	QI tendency due to flux convergence (shallow convection) in KF	kg/m3/s	QI_t_KF_conv_SH
QI_t_KF_updet_DP	QI tendency due to detrainment of updraft (deep convection) in KF	kg/m3/s	QI_t_KF_updet_DP
QI_t_KF_updet_SH	QI tendency due to detrainment of updraft (shallow convection) in KF	kg/m3/s	QI_t_KF_updet_SH
QR_t_KF_conv_DP	QR tendency due to flux convergence (deep convection) in KF	kg/m3/s	QR_t_KF_conv_DP
QR_t_KF_conv_SH	QR tendency due to flux convergence (shallow convection) in KF	kg/m3/s	QR_t_KF_conv_SH
QR_t_KF_rainfall_DP	QR tendency due to rain fall (deep convection) in KF	kg/m3/s	QR_t_KF_rainfall_DP
QR_t_KF_rainfall_SH	QR tendency due to rain fall (shallow convection) in KF	kg/m3/s	QR_t_KF_rainfall_SH
QS_t_KF_conv_DP	QS tendency due to flux convergence (deep convection) in KF	kg/m3/s	QS_t_KF_conv_DP
QS_t_KF_conv_SH	QS tendency due to flux convergence (shallow convection) in KF	kg/m3/s	QS_t_KF_conv_SH
QS_t_KF_snowfall_DP	QS tendency due to snow fall (deep convection) in KF	kg/m3/s	QS_t_KF_snowfall_DP
QS_t_KF_snowfall_SH	QS tendency due to snow fall (shallow convection) in KF	kg/m3/s	QS_t_KF_snowfall_SH
QV_t_KF_conv_DP	QV tendency due to flux convergence (deep convection) in KF	kg/m3/s	QV_t_KF_conv_DP
QV_t_KF_conv_SH	QV tendency due to flux convergence (shallow convection) in KF	kg/m3/s	QV_t_KF_conv_SH
QV_t_KF_downdet_DP	QV tendency due to detrainment of downdraft (deep convection) in KF	kg/m3/s	QV_t_KF_downdet_DP
QV_t_KF_downent_DP	QV tendency due to entrainment of downdraft (deep convection) in KF	kg/m3/s	QV_t_KF_downent_DP
QV_t_KF_negfix_DP	QV tendency due to negative fixer (deep convection) in KF	kg/m3/s	QV_t_KF_negfix_DP
QV_t_KF_negfix_SH	QV tendency due to negative fixer (shallow convection) in KF	kg/m3/s	QV_t_KF_negfix_SH
QV_t_KF_updet_DP	QV tendency due to detrainment of updraft (deep convection) in KF	kg/m3/s	QV_t_KF_updet_DP
QV_t_KF_updet_SH	QV tendency due to detrainment of updraft (shallow convection) in KF	kg/m3/s	QV_t_KF_updet_SH
QV_t_KF_upent_DP	QV tendency due to entrainment of updraft (deep convection) in KF	kg/m3/s	QV_t_KF_upent_DP
QV_t_KF_upent_SH	QV tendency due to entrainment of updraft (shallow convection) in KF	kg/m3/s	QV_t_KF_upent_SH

Function/Subroutine Documentation

atmos_phy_cp_kf_setup()

subroutine, public scale_atmos_phy_cp_kf::atmos_phy_cp_kf_setup	(	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)	CZ,
		real(rp), dimension(ia,ja), intent(in)	AREA,
		logical, intent(in)	WARMRAIN_in
	)

Setup initial setup for Kain-Fritsch Cumulus Parameterization.

Parameters

[in]

warmrain_in

tunning parameters, original parameter set is from KF2004 and NO2007

Definition at line 179 of file scale_atmos_phy_cp_kf.F90.

    use scale_prc, only: &
       prc_abort
    use scale_file_history, only: &
       file_history_reg
    implicit none
    integer, intent(in) :: KA, KS, KE
    integer, intent(in) :: IA, IS, IE
    integer, intent(in) :: JA, JS, JE
 
    real(RP),         intent(in) :: CZ(KA,IA,JA)
    real(RP),         intent(in) :: AREA(IA,JA)
    logical,          intent(in) :: WARMRAIN_in
 
    
    integer  :: ATMOS_PHY_CP_kf_trigger_type = 1
    integer  :: ATMOS_PHY_CP_kf_prec_type    = 1
    real(RP) :: ATMOS_PHY_CP_kf_rate         =   0.03_rp 
    logical  :: ATMOS_PHY_CP_kf_do_prec      = .true.    
    real(RP) :: ATMOS_PHY_CP_kf_dlcape       =    0.1_rp 
    real(RP) :: ATMOS_PHY_CP_kf_dlifetime    = 1800.0_rp 
    real(RP) :: ATMOS_PHY_CP_kf_slifetime    = 2400.0_rp 
    real(RP) :: ATMOS_PHY_CP_kf_DEPTH_USL    =  300.0_rp 
    real(RP) :: ATMOS_PHY_CP_kf_thres        = 1.e-3_rp  
    logical  :: ATMOS_PHY_CP_kf_LOG          = .false.   
 
    namelist / param_atmos_phy_cp_kf / &
       atmos_phy_cp_kf_rate,         &
       atmos_phy_cp_kf_trigger_type, &
       atmos_phy_cp_kf_dlcape,       &
       atmos_phy_cp_kf_dlifetime,    &
       atmos_phy_cp_kf_slifetime,    &
       atmos_phy_cp_kf_depth_usl,    &
       atmos_phy_cp_kf_prec_type,    &
       atmos_phy_cp_kf_do_prec,      &
       atmos_phy_cp_kf_thres,        &
       atmos_phy_cp_kf_log
 
    integer :: k, i, j
    integer :: n
    integer :: ierr
    !---------------------------------------------------------------------------
 
    log_newline
    log_info("ATMOS_PHY_CP_kf_setup",*) 'Setup'
    log_info("ATMOS_PHY_CP_kf_setup",*) 'Kain-Fritsch scheme'
 
    warmrain = warmrain_in
 
    !--- read namelist
    rewind(io_fid_conf)
    read(io_fid_conf,nml=param_atmos_phy_cp_kf,iostat=ierr)
    if( ierr < 0 ) then !--- missing
       log_info("ATMOS_PHY_CP_kf_setup",*) 'Not found namelist. Default used.'
    elseif( ierr > 0 ) then !--- fatal error
       log_error("ATMOS_PHY_CP_kf_setup",*) 'Not appropriate names in namelist PARAM_ATMOS_PHY_CP_KF. Check!'
       call prc_abort
    endif
    log_nml(param_atmos_phy_cp_kf)
 
    call cp_kf_lutab ! set up of KF look-up table
 
    call cp_kf_param( atmos_phy_cp_kf_prec_type,   & ! [IN]
                      atmos_phy_cp_kf_rate,        & ! [IN]
                      atmos_phy_cp_kf_do_prec,     & ! [IN]
                      atmos_phy_cp_kf_dlcape,      & ! [IN]
                      atmos_phy_cp_kf_dlifetime,   & ! [IN]
                      atmos_phy_cp_kf_slifetime,   & ! [IN]
                      atmos_phy_cp_kf_depth_usl,   & ! [IN]
                      atmos_phy_cp_kf_thres,       & ! [IN]
                      atmos_phy_cp_kf_log,         & ! [IN]
                      atmos_phy_cp_kf_trigger_type ) ! [IN]
 
    ! output parameter lists
    log_newline
    log_info("ATMOS_PHY_CP_kf_setup",*) "Ogura-Cho condense material convert rate        : ", atmos_phy_cp_kf_rate
    log_info("ATMOS_PHY_CP_kf_setup",*) "Trigger function type, 1:KF2004 3:NO2007        : ", atmos_phy_cp_kf_trigger_type
    log_info("ATMOS_PHY_CP_kf_setup",*) "CAPE decrease rate                              : ", atmos_phy_cp_kf_dlcape
    log_info("ATMOS_PHY_CP_kf_setup",*) "Minimum lifetime scale of deep convection [sec] : ", atmos_phy_cp_kf_dlifetime
    log_info("ATMOS_PHY_CP_kf_setup",*) "Lifetime of shallow convection            [sec] : ", atmos_phy_cp_kf_slifetime
    log_info("ATMOS_PHY_CP_kf_setup",*) "Updraft souce layer depth                 [hPa] : ", atmos_phy_cp_kf_depth_usl
    log_info("ATMOS_PHY_CP_kf_setup",*) "Precipitation type 1:Ogura-Cho(1973) 2:Kessler  : ", atmos_phy_cp_kf_prec_type
    log_info("ATMOS_PHY_CP_kf_setup",*) "Kessler type precipitation's threshold          : ", atmos_phy_cp_kf_thres
    log_info("ATMOS_PHY_CP_kf_setup",*) "Enable sedimentation (precipitation) ?          : ", atmos_phy_cp_kf_do_prec
    log_info("ATMOS_PHY_CP_kf_setup",*) "Warm rain?                                      : ", warmrain
    log_info("ATMOS_PHY_CP_kf_setup",*) "Output log?                                     : ", atmos_phy_cp_kf_log
 
    ! output variables
    allocate( lifetime(ia,ja) )
    allocate( i_convflag(ia,ja) )
    lifetime(:,:) = 0.0_rp
    i_convflag(:,:) = 2
 
    allocate( z(ka,ia,ja) )
    z(:,:,:) = cz(:,:,:) ! because scale_atmos_phy_cp interface ,not use scale_grid
 
    allocate( deltaz(ka,ia,ja) )
    ! deltaz is the interval of between model full levels(scalar point )
    deltaz(:,:,:) = 0.0_rp
    do j = js, je
    do i = is, ie
    do k = ks, ke
       deltaz(k,i,j) = cz(k+1,i,j) - cz(k,i,j)
    enddo
    enddo
    enddo
    deltaz(ke,:,:) = 0.0_rp
 
    allocate( deltax(ia,ja) )
    do j = js, je
    do i = is, ie
       deltax(i,j) = sqrt( area(i,j) )
    end do
    end do
 
 
    ! history
    hist_id(:,:) = -1
    call file_history_reg( 'QV_t_KF_conv_DP',                                                &
                           'QV tendency due to flux convergence (deep convection) in KF',    &
                           'kg/m3/s',                                                        &
                           hist_id(i_hist_qv_fc,i_hist_d)                                    )
    call file_history_reg( 'QV_t_KF_conv_SH',                                                &
                           'QV tendency due to flux convergence (shallow convection) in KF', &
                           'kg/m3/s',                                                        &
                           hist_id(i_hist_qv_fc,i_hist_s)                                    )
    call file_history_reg( 'QV_t_KF_upent_DP',                                                     &
                           'QV tendency due to entrainment of updraft (deep convection) in KF',    &
                           'kg/m3/s',                                                              &
                           hist_id(i_hist_qv_ue,i_hist_d)                                          )
    call file_history_reg( 'QV_t_KF_upent_SH',                                                     &
                           'QV tendency due to entrainment of updraft (shallow convection) in KF', &
                           'kg/m3/s',                                                              &
                           hist_id(i_hist_qv_ue,i_hist_s)                                          )
    call file_history_reg( 'QV_t_KF_downent_DP',                                                  &
                           'QV tendency due to entrainment of downdraft (deep convection) in KF', &
                           'kg/m3/s',                                                             &
                           hist_id(i_hist_qv_de,i_hist_d)                                         )
    call file_history_reg( 'QV_t_KF_updet_DP',                                                     &
                           'QV tendency due to detrainment of updraft (deep convection) in KF',    &
                           'kg/m3/s',                                                              &
                           hist_id(i_hist_qv_ud,i_hist_d)                                          )
    call file_history_reg( 'QV_t_KF_updet_SH',                                                     &
                           'QV tendency due to detrainment of updraft (shallow convection) in KF', &
                           'kg/m3/s',                                                              &
                           hist_id(i_hist_qv_ud,i_hist_s)                                          )
    call file_history_reg( 'QV_t_KF_downdet_DP',                                                  &
                           'QV tendency due to detrainment of downdraft (deep convection) in KF', &
                           'kg/m3/s',                                                             &
                           hist_id(i_hist_qv_dd,i_hist_d)                                         )
    call file_history_reg( 'QV_t_KF_negfix_DP',                                            &
                           'QV tendency due to negative fixer (deep convection) in KF',    &
                           'kg/m3/s',                                                      &
                           hist_id(i_hist_qv_nf,i_hist_d)                                  )
    call file_history_reg( 'QV_t_KF_negfix_SH',                                            &
                           'QV tendency due to negative fixer (shallow convection) in KF', &
                           'kg/m3/s',                                                      &
                           hist_id(i_hist_qv_nf,i_hist_s)                                  )
 
    call file_history_reg( 'QC_t_KF_conv_DP',                                                &
                           'QC tendency due to flux convergence (deep convection) in KF',    &
                           'kg/m3/s',                                                        &
                           hist_id(i_hist_qc_fc,i_hist_d)                                    )
    call file_history_reg( 'QC_t_KF_conv_SH',                                                &
                           'QC tendency due to flux convergence (shallow convection) in KF', &
                           'kg/m3/s',                                                        &
                           hist_id(i_hist_qc_fc,i_hist_s)                                    )
    call file_history_reg( 'QC_t_KF_updet_DP',                                                     &
                           'QC tendency due to detrainment of updraft (deep convection) in KF',    &
                           'kg/m3/s',                                                              &
                           hist_id(i_hist_qc_ud,i_hist_d)                                          )
    call file_history_reg( 'QC_t_KF_updet_SH',                                                     &
                           'QC tendency due to detrainment of updraft (shallow convection) in KF', &
                           'kg/m3/s',                                                              &
                           hist_id(i_hist_qc_ud,i_hist_s)                                          )
 
    call file_history_reg( 'QI_t_KF_conv_DP',                                                &
                           'QI tendency due to flux convergence (deep convection) in KF',    &
                           'kg/m3/s',                                                        &
                           hist_id(i_hist_qi_fc,i_hist_d)                                    )
    call file_history_reg( 'QI_t_KF_conv_SH',                                                &
                           'QI tendency due to flux convergence (shallow convection) in KF', &
                           'kg/m3/s',                                                        &
                           hist_id(i_hist_qi_fc,i_hist_s)                                    )
    call file_history_reg( 'QI_t_KF_updet_DP',                                                     &
                           'QI tendency due to detrainment of updraft (deep convection) in KF',    &
                           'kg/m3/s',                                                              &
                           hist_id(i_hist_qi_ud,i_hist_d)                                          )
    call file_history_reg( 'QI_t_KF_updet_SH',                                                     &
                           'QI tendency due to detrainment of updraft (shallow convection) in KF', &
                           'kg/m3/s',                                                              &
                           hist_id(i_hist_qi_ud,i_hist_s)                                          )
 
    call file_history_reg( 'QR_t_KF_conv_DP',                                                &
                           'QR tendency due to flux convergence (deep convection) in KF',    &
                           'kg/m3/s',                                                        &
                           hist_id(i_hist_qr_fc,i_hist_d)                                    )
    call file_history_reg( 'QR_t_KF_conv_SH',                                                &
                           'QR tendency due to flux convergence (shallow convection) in KF', &
                           'kg/m3/s',                                                        &
                           hist_id(i_hist_qr_fc,i_hist_s)                                    )
    call file_history_reg( 'QR_t_KF_rainfall_DP',                                     &
                           'QR tendency due to rain fall (deep convection) in KF',    &
                           'kg/m3/s',                                                 &
                           hist_id(i_hist_qr_rf,i_hist_d)                             )
    call file_history_reg( 'QR_t_KF_rainfall_SH',                                     &
                           'QR tendency due to rain fall (shallow convection) in KF', &
                           'kg/m3/s',                                                 &
                           hist_id(i_hist_qr_rf,i_hist_s)                             )
 
    call file_history_reg( 'QS_t_KF_conv_DP',                                                &
                           'QS tendency due to flux convergence (deep convection) in KF',    &
                           'kg/m3/s',                                                        &
                           hist_id(i_hist_qs_fc,i_hist_d)                                    )
    call file_history_reg( 'QS_t_KF_conv_SH',                                                &
                           'QS tendency due to flux convergence (shallow convection) in KF', &
                           'kg/m3/s',                                                        &
                           hist_id(i_hist_qs_fc,i_hist_s)                                    )
    call file_history_reg( 'QS_t_KF_snowfall_DP',                                     &
                           'QS tendency due to snow fall (deep convection) in KF',    &
                           'kg/m3/s',                                                 &
                           hist_id(i_hist_qs_sf,i_hist_d)                             )
    call file_history_reg( 'QS_t_KF_snowfall_SH',                                     &
                           'QS tendency due to snow fall (shallow convection) in KF', &
                           'kg/m3/s',                                                 &
                           hist_id(i_hist_qs_sf,i_hist_s)                             )
 
    do n = 1, hist_vmax
       if ( hist_id(n,0) > 0 .or. hist_id(n,1) > 0 ) then
          allocate( hist_work(ka,ia,ja,hist_vmax,0:1) )
          hist_work(:,:,:,:,:) = 0.0_rp
          exit
       end if
    end do
 
    return

References scale_file_history::file_history_reg(), scale_atmos_grid_cartesc_index::ia, scale_atmos_grid_cartesc_index::ie, scale_io::io_fid_conf, scale_atmos_grid_cartesc_index::is, scale_atmos_grid_cartesc_index::ja, scale_atmos_grid_cartesc_index::je, scale_atmos_grid_cartesc_index::js, scale_tracer::k, scale_atmos_grid_cartesc_index::ka, scale_atmos_grid_cartesc_index::ke, scale_atmos_grid_cartesc_index::ks, and scale_prc::prc_abort().

Referenced by mod_atmos_phy_cp_driver::atmos_phy_cp_driver_setup().

Here is the call graph for this function:

Here is the caller graph for this function:

atmos_phy_cp_kf_tendency()

subroutine, public scale_atmos_phy_cp_kf::atmos_phy_cp_kf_tendency	(	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)	U,
		real(rp), dimension (ka,ia,ja), intent(in)	V,
		real(rp), dimension (ka,ia,ja), intent(in)	RHOT,
		real(rp), dimension (ka,ia,ja), intent(in)	TEMP,
		real(rp), dimension (ka,ia,ja), intent(in)	PRES,
		real(rp), dimension (ka,ia,ja), intent(in)	QDRY,
		real(rp), dimension(ka,ia,ja), intent(in)	QV_in,
		real(rp), dimension(ka,ia,ja), intent(in)	w0avg,
		real(rp), dimension (0:ka,ia,ja), intent(in)	FZ,
		real(dp), intent(in)	KF_DTSECD,
		real(rp), dimension (ka,ia,ja), intent(inout)	DENS_t,
		real(rp), dimension (ka,ia,ja), intent(inout)	RHOT_t,
		real(rp), dimension(ka,ia,ja), intent(inout)	RHOQV_t,
		real(rp), dimension (ka,ia,ja,n_hyd), intent(inout)	RHOQ_t,
		real(rp), dimension (ia,ja), intent(inout)	nca,
		real(rp), dimension(ia,ja), intent(inout)	SFLX_rain,
		real(rp), dimension(ia,ja), intent(inout)	SFLX_snow,
		real(rp), dimension(ia,ja), intent(inout)	SFLX_engi,
		real(rp), dimension (ia,ja), intent(inout)	cloudtop,
		real(rp), dimension (ia,ja), intent(inout)	cloudbase,
		real(rp), dimension (ka,ia,ja), intent(inout)	cldfrac_dp,
		real(rp), dimension (ka,ia,ja), intent(inout)	cldfrac_sh
	)

ATMOS_PHY_CP_kf calculate Kain-Fritsch Cumulus Parameterization.

Parameters

[in]	dens	Density [kg/m3]
[in]	u	velocity u [m/s]
[in]	v	velocity v [m/s]
[in]	rhot	DENS * POTT [k*kg/m3]
[in]	temp	temperature [K]
[in]	pres	pressure of dry air [Pa]
[in]	qdry	dry air [1]
[in]	qv_in	specific humidity [kg/kg]
[in]	w0avg	running mean of vertical velocity [m/s]
[in,out]	dens_t	tendency DENS [kg/m3/s]
[in,out]	rhot_t	tendency RHOT [K*kg/m3/s]
[in,out]	rhoqv_t	tendency rho*QV [kg/kg/s]
[in,out]	rhoq_t	tendency rho*QTRC [kg/kg/s]
[in,out]	nca	convection active time [sec]
[in,out]	sflx_rain	convective rain rate [kg/m2/s]
[in,out]	sflx_snow	convective snow rate [kg/m2/s]
[in,out]	sflx_engi	precipitation internal energy [J/m2/s]
[in,out]	cloudtop	cloud top height [m]
[in,out]	cloudbase	cloud base height [m]
[in,out]	cldfrac_dp	cloud fraction (deep convection)
[in,out]	cldfrac_sh	cloud fraction (shallow convection)

Definition at line 503 of file scale_atmos_phy_cp_kf.F90.

    use scale_file_history, only: &
       file_history_query, &
       file_history_put, &
       file_history_in
    use scale_const, only: &
       grav   => const_grav,  &
       pre00  => const_pre00, &
       epsvap => const_epsvap
    use scale_atmos_hydrometeor, only: &
       cp_vapor, &
       cp_water, &
       cp_ice,   &
       cv_vapor, &
       cv_water, &
       cv_ice,   &
       n_hyd, &
       i_hc, &
       i_hr, &
       i_hi, &
       i_hs
    use scale_atmos_saturation ,only :&
       saturation_psat_liq => atmos_saturation_psat_liq
    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)
    real(RP), intent(in)    :: U    (KA,IA,JA)
    real(RP), intent(in)    :: V    (KA,IA,JA)
    real(RP), intent(in)    :: RHOT (KA,IA,JA)
    real(RP), intent(in)    :: TEMP (KA,IA,JA)
    real(RP), intent(in)    :: PRES (KA,IA,JA)
    real(RP), intent(in)    :: QDRY (KA,IA,JA)
    real(RP), intent(in)    :: QV_in(KA,IA,JA)
    real(RP), intent(in)    :: w0avg(KA,IA,JA)
    real(RP), intent(in)    :: FZ   (0:KA,IA,JA)
    real(DP), intent(in)    :: KF_DTSECD
 
    real(RP), intent(inout) :: DENS_t (KA,IA,JA)
    real(RP), intent(inout) :: RHOT_t (KA,IA,JA)
    real(RP), intent(inout) :: RHOQV_t(KA,IA,JA)
    real(RP), intent(inout) :: RHOQ_t (KA,IA,JA,N_HYD)
    real(RP), intent(inout) :: nca    (IA,JA)
 
    real(RP), intent(inout) :: SFLX_rain(IA,JA)
    real(RP), intent(inout) :: SFLX_snow(IA,JA)
    real(RP), intent(inout) :: SFLX_engi(IA,JA)
    real(RP), intent(inout) :: cloudtop     (IA,JA)
    real(RP), intent(inout) :: cloudbase    (IA,JA)
    real(RP), intent(inout) :: cldfrac_dp   (KA,IA,JA)
    real(RP), intent(inout) :: cldfrac_sh   (KA,IA,JA)
 
    integer  :: k, i, j, iq
    integer  :: nic
    integer  :: k_lcl
    integer  :: k_top
    integer  :: k_ml
    integer  :: k_lc,k_let,k_pbl
    integer  :: k_lfs
 
    real(RP) :: qv    (KA)
    real(RP) :: PSAT
    real(RP) :: QSAT  (KA)
    real(RP) :: rh    (KA)
    real(RP) :: deltap(KA)
 
    real(RP) :: cldfrac_KF(KA,2)
    real(RP) :: dens_nw(KA)
    real(RP) :: time_advec
    real(RP) :: umf(KA)
    real(RP) :: umflcl
    real(RP) :: upent(KA)
    real(RP) :: updet(KA)
    ! updraft value
    real(RP) :: temp_u(KA)
    real(RP) :: tempv(KA)
    real(RP) :: qv_u(KA)
    real(RP) :: cape
    ! water
    real(RP) :: qc(KA)
    real(RP) :: qi(KA)
    real(RP) :: qvdet(KA)
    real(RP) :: qcdet(KA)
    real(RP) :: qidet(KA)
    real(RP) :: flux_qr(KA)
    real(RP) :: flux_qs(KA)
    ! upward theta
    real(RP) :: theta_eu(KA)
    real(RP) :: theta_ee(KA)
    real(RP) :: zmix
    real(RP) :: presmix
    real(RP) :: umfnewdold(KA)
    real(RP) :: dpthmx
    real(RP) :: ems(KA)
    real(RP) :: emsd(KA)
    ! output from downdraft
    real(RP) :: wspd(3)
    real(RP) :: dmf(KA)
    real(RP) :: downent(KA)
    real(RP) :: downdet(KA)
    real(RP) :: theta_d(KA)
    real(RP) :: qv_d(KA)
    real(RP) :: rain_flux
    real(RP) :: snow_flux
    real(RP) :: prec_engi
 
    ! update variables
    real(RP) :: theta_nw(KA)
    real(RP) :: qv_nw(KA)
    real(RP) :: qc_nw(KA)
    real(RP) :: qi_nw(KA)
    real(RP) :: qr_nw(KA)
    real(RP) :: qs_nw(KA)
 
    real(RP) :: RHOD(KA)
 
    real(RP) :: dQV, dQC, dQI, dQR, dQS
 
    real(RP) :: R_convflag(IA,JA)
 
    real(RP) :: KF_DTSEC
 
    logical :: flag
    integer :: n, m
    ! ------
 
    log_progress(*) 'atmosphere / physics / cumulus / KF'
    log_info("ATMOS_PHY_CP_kf_tendency",*) 'KF Convection Check '
 
    kf_dtsec = kf_dtsecd ! DP to RP
 
    call prof_rapstart('CP_kf', 3)
 
    do n = 1, hist_vmax
       call file_history_query( hist_id(n,0), hist_flag )
       if ( hist_flag ) exit
       call file_history_query( hist_id(n,1), hist_flag )
       if ( hist_flag ) exit
    end do
 
    !$omp parallel do default(none) &
    !$omp private(RHOD,tempv,qv_d,qc,qi,PSAT,QSAT,QV,rh, &
    !$omp         dens_nw,theta_nw,qv_nw,qc_nw,qi_nw,qr_nw,qs_nw, &
    !$omp         flux_qr,flux_qs,theta_eu,theta_ee,theta_d,umfnewdold,qvdet,qcdet,qidet, &
    !$omp         umf,umflcl,cape,presmix,upent,updet,temp_u,qv_u, &
    !$omp         dpthmx,zmix,rain_flux,snow_flux,prec_engi,time_advec, &
    !$omp         k_lcl,k_lc,k_lfs,k_pbl,k_top,k_let,k_ml, &
    !$omp         nic,deltap,cldfrac_KF,ems,emsd,wspd,dmf,downent,downdet, &
    !$omp         dQv,dQC,dQR,dQI,dQS) &
    !$omp shared(DENS,RHOT,U,V,QDRY,TEMP,PRES,QV_in,FZ,Z,w0avg, &
    !$omp        DENS_t,RHOT_t,RHOQV_t,RHOQ_t, &
    !$omp        GRAV,PRE00,CP_VAPOR,CP_WATER,CP_ICE,CV_VAPOR,CV_WATER,CV_ICE,EPSvap, &
    !$omp        WARMRAIN,KF_DTSEC,SHALLOWLIFETIME, &
    !$omp        KA,KS,KE,IS,IE,JS,JE, &
    !$omp        nca,cloudtop,cloudbase,SFLX_rain,SFLX_snow,SFLX_engi, &
    !$omp        lifetime,deltaz,deltax,I_convflag, &
    !$omp        cldfrac_sh,cldfrac_dp,hist_flag,hist_work)
    do j = js, je
    do i = is, ie
 
       nca(i,j) = nca(i,j) - kf_dtsec
 
       ! check convection
       if ( nca(i,j) .ge. 0.5_rp * kf_dtsec ) cycle
 
       do k = ks, ke
          ! preparing a NON Hydriometeor condition to fit assumption in KF scheme
          rhod(k) = dens(k,i,j) * qdry(k,i,j)
       enddo
 
 
       do k = ks, ke
          ! temporary: WRF TYPE equations are used to maintain consistency with kf_main
          !call SATURATION_psat_liq( PSAT, TEMP(k,i,j) )
          !QSAT(k) = 0.622_RP * PSAT / ( PRES(k,i,j) - ( 1.0_RP-0.622_RP ) * PSAT )
          psat    = aliq*exp((bliq*temp(k,i,j)-cliq)/(temp(k,i,j)-dliq))
          qsat(k) = epsvap * psat / ( pres(k,i,j) - psat )
 
          ! calculate water vaper and relative humidity
          qv(k) = max( kf_eps, min( qsat(k), qv_in(k,i,j) / qdry(k,i,j) ) ) ! compare QSAT and QV, guess lower limit
          rh(k) = qv(k) / qsat(k)
       enddo
 
       ! calculate delta P by hydrostatic balance
       ! deltap is the pressure interval between half levels(face levels) @ SCALE
       do k = ks, ke
          deltap(k) = rhod(k) * grav * ( fz(k,i,j) - fz(k-1,i,j) ) ! rho*g*dz
       enddo
 
       call cp_kf_trigger ( &
            ka, ks, ke,                   & ! [IN]
            deltaz(:,i,j), z(:,i,j),      & ! [IN]
            qv(:), qsat(:),               & ! [IN]
            pres(:,i,j),                  & ! [IN]
            deltap(:), deltax(i,j),       & ! [IN]
            temp(:,i,j),                  & ! [IN]
            w0avg(:,i,j),                 & ! [IN]
            i_convflag(i,j),              & ! [OUT]
            cloudtop(i,j),                & ! [OUT]
            temp_u(:), tempv(:),          & ! [OUT]
            qv_u(:), qc(:), qi(:),        & ! [OUT]
            qvdet(:), qcdet(:), qidet(:), & ! [OUT]
            flux_qr(:), flux_qs(:),       & ! [OUT]
            theta_eu(:), theta_ee(:),     & ! [OUT]
            cape,                         & ! [OUT]
            umf(:), umflcl,               & ! [OUT]
            upent(:), updet(:),           & ! [OUT]
            k_lcl, k_lc,  k_pbl,          & ! [OUT]
            k_top, k_let, k_ml,           & ! [OUT]
            presmix,                      & ! [OUT]
            dpthmx,                       & ! [OUT]
            cloudbase(i,j), zmix,         & ! [OUT]
            umfnewdold(:)                 ) ! [OUT]
 
       if (i_convflag(i,j) /= 2) then ! convection allowed I_convflag=0 or 1
 
          ! calc ems(box weight[kg])
          ems(k_top+1:ke) = 0._rp
          emsd(k_top+1:ke) = 0._rp
          do k = ks, k_top
             ems(k) = deltap(k) * deltax(i,j)**2 / grav
             emsd(k) = 1._rp/ems(k)
             umfnewdold(k) = 1._rp/umfnewdold(k)
          end do
 
          call cp_kf_downdraft ( &
               ka, ks, ke,                                    & ! [IN]
               i_convflag(i,j),                               & ! [IN]
               k_lcl, k_ml, k_top, k_pbl, k_let, k_lc,        & ! [IN]
               z(:,i,j), cloudbase(i,j),                      & ! [IN]
               u(:,i,j), v(:,i,j), rh(:), qv(:), pres(:,i,j), & ! [IN]
               deltap(:), deltax(i,j),                        & ! [IN]
               ems(:),                                        & ! [IN]
               theta_ee(:),                                   & ! [IN]
               umf(:), temp_u(:),                             & ! [IN]
               flux_qr(:), flux_qs(:), tempv(:),              & ! [IN]
               wspd(:), dmf(:), downent(:), downdet(:),       & ! [OUT]
               theta_d(:), qv_d(:),                           & ! [OUT]
               rain_flux, snow_flux, prec_engi,               & ! [OUT]
               k_lfs                                          ) ! [OUT]
 
          call cp_kf_compensational ( &
               ka, ks, ke,                                                   & ! [IN]
               k_top, k_lc, k_pbl, k_ml, k_lfs,                              & ! [IN]
               deltaz(:,i,j), z(:,i,j), pres(:,i,j), deltap(:), deltax(i,j), & ! [IN]
               temp(:,i,j), qv(:), ems(:), emsd(:),                          & ! [IN]
               presmix, zmix, dpthmx,                                        & ! [IN]
               cape,                                                         & ! [IN]
               temp_u(:), qvdet(:), umflcl,                                  & ! [IN]
               qc(:), qi(:), flux_qr(:), flux_qs(:),                         & ! [IN]
               umfnewdold(:),                                                & ! [IN]
               wspd(:),                                                      & ! [IN]
               qv_d(:), theta_d(:),                                          & ! [IN]
               prec_engi,                                                    & ! [IN]
               kf_dtsec,                                                     & ! [IN]
               rhod(:), i, j,                                                & ! [IN]
               i_convflag(i,j), k_lcl,                                       & ! [INOUT]
               umf(:), upent(:), updet(:),                                   & ! [INOUT]
               qcdet(:), qidet(:), dmf(:), downent(:), downdet(:),           & ! [INOUT]
               rain_flux, snow_flux,                                         & ! [INOUT]
               nic,                                                          & ! [OUT]
               theta_nw(:),                                                  & ! [OUT]
               qv_nw(:), qc_nw(:), qi_nw(:), qr_nw(:), qs_nw(:),             & ! [OUT]
               sflx_rain(i,j), sflx_snow(i,j), sflx_engi(i,j),               & ! [OUT]
               cldfrac_kf, lifetime(i,j), time_advec                         ) ! [OUT]
 
       end if
 
       if (i_convflag(i,j) == 2) then ! no convection
 
          sflx_rain(i,j)       = 0.0_rp
          sflx_snow(i,j)       = 0.0_rp
          sflx_engi(i,j)       = 0.0_rp
          cldfrac_kf(ks:ke,:)  = 0.0_rp
          lifetime(i,j)        = 0.0_rp
          cloudtop(i,j)        = 0.0_rp
          cloudbase(i,j)       = 0.0_rp
          nca(i,j)             = -100.0_rp
 
          dens_t(ks:ke,i,j)  = 0.0_rp
          rhot_t(ks:ke,i,j)  = 0.0_rp
          rhoqv_t(ks:ke,i,j) = 0.0_rp
          do iq = 1, n_hyd
             rhoq_t(ks:ke,i,j,iq) = 0.0_rp
          end do
 
          if ( hist_flag ) then
             do m = 0, 1
             do n = 1, hist_vmax
                hist_work(:,i,j,n,m) = 0.0_rp
             end do
             end do
          end if
 
       else
 
          ! compute tendencys
 
          ! check:
          ! FEEDBACK TO RESOLVABLE SCALE TENDENCIES.
          ! IF THE ADVECTIVE TIME PERIOD (time_advec) IS LESS THAN SPECIFIED MINIMUM
          ! lifetime, ALLOW FEEDBACK TO OCCUR ONLY DURING TADVEC.
          !
          if (i_convflag(i,j) == 0) then ! deep
             if (time_advec < lifetime(i,j)) nic=nint(time_advec/kf_dtsec)
             nca(i,j) = nic * kf_dtsec ! convection feed back act this time span
          elseif (i_convflag(i,j) == 1) then ! shallow
             lifetime(i,j) = max(shallowlifetime, kf_dtsec)
             nca(i,j) = kf_dtsec ! convection feed back act this time span
          end if
 
          do k = ks, k_top
             ! vapor
             dqv = rhod(k) * ( qv_nw(k) - qv(k) )
             rhoqv_t(k,i,j) = dqv / lifetime(i,j)
             ! liquid water
             dqc = qc_nw(k) * rhod(k)
             dqr = qr_nw(k) * rhod(k)
             rhoq_t(k,i,j,i_hc) = dqc / lifetime(i,j)
             rhoq_t(k,i,j,i_hr) = dqr / lifetime(i,j)
             ! ice water
             if ( .not. warmrain ) then
                dqi = qi_nw(k) * rhod(k)
                dqs = qs_nw(k) * rhod(k)
             else
                dqi = 0.0_rp
                dqs = 0.0_rp
             end if
             rhoq_t(k,i,j,i_hi) = dqi / lifetime(i,j)
             rhoq_t(k,i,j,i_hs) = dqs / lifetime(i,j)
 
             dens_t(k,i,j) = rhoqv_t(k,i,j) &
                           + rhoq_t(k,i,j,i_hc) + rhoq_t(k,i,j,i_hr) &
                           + rhoq_t(k,i,j,i_hi) + rhoq_t(k,i,j,i_hs)
             dens_nw(k) = dens(k,i,j) + dqv + dqc + dqr + dqi + dqs
 
          end do
          do iq = i_hs+1, n_hyd
          do k = ks, k_top
             rhoq_t(k,i,j,iq) = 0.0_rp
          end do
          end do
 
          ! internal energy
          do k = ks, k_top
             rhot_t(k,i,j) = ( dens_nw(k) * theta_nw(k) - rhot(k,i,j) ) / lifetime(i,j)
          end do
 
          do k=k_top+1, ke
             dens_t(k,i,j) = 0.0_rp
             rhot_t(k,i,j) = 0.0_rp
             rhoqv_t(k,i,j) = 0.0_rp
             do iq = 1, n_hyd
                rhoq_t(k,i,j,iq) = 0.0_rp
             end do
          end do
 
          ! to keep conservation
          ! if noconvection then nca is same value before call. nca only modifyed convectioned
       end if
 
       cldfrac_sh(ks:ke,i,j) = cldfrac_kf(ks:ke,1)
       cldfrac_dp(ks:ke,i,j) = cldfrac_kf(ks:ke,2)
    end do
    end do
 
    call prof_rapend('CP_kf', 3)
 
    ! history
    call file_history_in( lifetime(:,:),   'KF_LIFETIME', 'lifetime of KF scheme', 's' )
    r_convflag(:,:) = real(i_convflag(:,:),rp)
    call file_history_in( r_convflag(:,:), 'KF_CONVFLAG', 'CONVECTION FLAG',       ''  )
    do m = 0, 1
    do n = 1, hist_vmax
       call file_history_query( hist_id(n,m), flag )
       if ( flag ) then
          call file_history_put( hist_id(n,m), hist_work(:,:,:,n,m) )
       end if
    end do
    end do
 
    return

References scale_const::const_cpdry, scale_const::const_cpvap, scale_const::const_emelt, scale_const::const_epstvap, scale_const::const_epsvap, scale_const::const_grav, scale_const::const_pre00, scale_const::const_rdry, scale_const::const_rvap, scale_const::const_tem00, scale_atmos_hydrometeor::cp_ice, scale_atmos_hydrometeor::cp_vapor, scale_atmos_hydrometeor::cp_water, scale_atmos_hydrometeor::cv_ice, scale_atmos_hydrometeor::cv_vapor, scale_atmos_hydrometeor::cv_water, scale_atmos_hydrometeor::i_hc, scale_atmos_hydrometeor::i_hi, scale_atmos_hydrometeor::i_hr, scale_atmos_hydrometeor::i_hs, scale_atmos_grid_cartesc_index::ie, scale_atmos_grid_cartesc_index::is, scale_atmos_grid_cartesc_index::je, scale_atmos_grid_cartesc_index::js, scale_tracer::k, scale_atmos_grid_cartesc_index::ka, scale_atmos_grid_cartesc_index::ke, scale_atmos_grid_cartesc_index::ks, scale_atmos_hydrometeor::lhf, scale_atmos_hydrometeor::lhs, scale_atmos_hydrometeor::n_hyd, scale_prc::prc_abort(), scale_prof::prof_rapend(), scale_prof::prof_rapstart(), and scale_precision::rp.

Referenced by mod_atmos_phy_cp_driver::atmos_phy_cp_driver_calc_tendency().

Here is the call graph for this function:

Here is the caller graph for this function: