scale_atmos_phy_cp_kf_jmapplib Module Reference

module atmosphere / physics / cumulus / kf-jmapplib More...

Functions/Subroutines
subroutine, public	atmos_phy_cp_kf_jmapplib_setup (KA, KS, KE, IA, JA, dx, dt)
	ATMOS_PHY_CP_KF_JMAPPLIB_setup Setup. More...
subroutine, public	atmos_phy_cp_kf_jmapplib_finalize
subroutine, public	atmos_phy_cp_kf_jmapplib_tendency (KA, KS, KE, IA, IS, IE, JA, JS, JE, DENS, U, V, W, TEMP, POTT, PRES, EXNER, QDRY, QV, QC, QI, us, PBLH, SFLX_BUOY, CZ, FZ, DENS_t, RHOT_t, RHOQV_t, RHOQ_t, w0avg, nca, SFLX_rain, SFLX_snow, SFLX_prec, cloudtop, cloudbase)
	ATMOS_PHY_CP_KF_JMAPPLIB_tendency calculate tendency by the virtical eddy viscosity. More...

Detailed Description

module atmosphere / physics / cumulus / kf-jmapplib

Description

Cumulus parameterization Kain-Fritsch scheme implemented in JMA Physics Process library

Author

Team SCALE

NAMELIST

• PARAM_ATMOS_PHY_CP_KF_JMAPPLIB

  name	type	default value	comment
  ATMOS_PHY_CP_KF_JMAPPLIB_TYPE	character(len=6)	"KF1701"	> KF1701 for MSM; KF for LFM
  ATMOS_PHY_CP_KF_JMAPPLIB_RADIUS_TYPE	integer	3	> 3 for KF1701
  ATMOS_PHY_CP_KF_JMAPPLIB_DTLCL_FACT	real(RP)	0.5_RP	> ignored for KF1701
  ATMOS_PHY_CP_KF_JMAPPLIB_DPMIN	real(RP)	5000.0_RP	> 5000 for KF1701; 1000 for KF
  ATMOS_PHY_CP_KF_JMAPPLIB_DLIFETIME	real(RP)	600.0_RP	> 600 for KF1701; 3600 for KF
  ATMOS_PHY_CP_KF_JMAPPLIB_SLIFETIME	real(RP)	600.0_RP	> 600 for KF1701; 3600 for KF
  ATMOS_PHY_CP_KF_JMAPPLIB_DETLQ_QR_RATE	real(RP)	0.0_RP	> 0.0 for KF1701; 0.0 for KF
  ATMOS_PHY_CP_KF_JMAPPLIB_DETIC_QS_RATE	real(RP)	1.0_RP	> 1.0 for KF1701; 0.0 for KF

  
  

History Output

No history output

Function/Subroutine Documentation

atmos_phy_cp_kf_jmapplib_setup()

subroutine, public scale_atmos_phy_cp_kf_jmapplib::atmos_phy_cp_kf_jmapplib_setup	(	integer, intent(in)	KA,
		integer, intent(in)	KS,
		integer, intent(in)	KE,
		integer, intent(in)	IA,
		integer, intent(in)	JA,
		real(rp), intent(in)	dx,
		real(dp), intent(in)	dt
	)

ATMOS_PHY_CP_KF_JMAPPLIB_setup Setup.

Definition at line 62 of file scale_atmos_phy_cp_kf_jmapplib.F90.

    use scale_prc, only: &
       prc_abort
    use scale_const, only: &
       pre00 => const_pre00
#ifdef JMAPPLIB
    use conv_grid, only: &
       conv_grid_ini
    use conv_kf_const, only: &
       conv_kf_const_ini
    use conv_kf_parm, only: &
       conv_kf_parm_ini, &
       conv_kf_parm_get_stepcu
    use conv_kf_parm1701, only: &
       conv_kf_parm_ini1701, &
       conv_kf_parm_get_stepcu1701
    use conv_kf_lut, only: &
       conv_kf_lut_ini
#endif
    implicit none
 
    integer,  intent(in) :: KA, KS, KE
    integer,  intent(in) :: IA, JA
 
    real(RP), intent(in) :: dx
    real(DP), intent(in) :: dt
 
#ifdef JMAPPLIB
 
    integer  :: ATMOS_PHY_CP_KF_JMAPPLIB_radius_type   = 3
    real(RP) :: ATMOS_PHY_CP_KF_JMAPPLIB_dtlcl_fact    = 0.5_rp
    real(RP) :: ATMOS_PHY_CP_KF_JMAPPLIB_dpmin         = 5000.0_rp
    real(RP) :: ATMOS_PHY_CP_KF_JMAPPLIB_dlifetime     = 600.0_rp
    real(RP) :: ATMOS_PHY_CP_KF_JMAPPLIB_slifetime     = 600.0_rp
    real(RP) :: ATMOS_PHY_CP_KF_JMAPPLIB_detlq_qr_rate = 0.0_rp
    real(RP) :: ATMOS_PHY_CP_KF_JMAPPLIB_detic_qs_rate = 1.0_rp
 
    namelist / param_atmos_phy_cp_kf_jmapplib / &
         atmos_phy_cp_kf_jmapplib_type, &
         atmos_phy_cp_kf_jmapplib_radius_type, &
         atmos_phy_cp_kf_jmapplib_dtlcl_fact, &
         atmos_phy_cp_kf_jmapplib_dpmin, &
         atmos_phy_cp_kf_jmapplib_dlifetime, &
         atmos_phy_cp_kf_jmapplib_slifetime, &
         atmos_phy_cp_kf_jmapplib_detlq_qr_rate, &
         atmos_phy_cp_kf_jmapplib_detic_qs_rate
 
    integer :: KMAX
 
    integer :: ierr
    integer :: k
    !---------------------------------------------------------------------------
 
    log_newline
    log_info("ATMOS_PHY_CP_KF_JMAPPLIB_setup",*) 'Setup'
    log_info("ATMOS_PHY_CP_KF_JMAPPLIB_setup",*) 'KF scheme implemented in the JMA Physics Process Library'
 
    !--- read namelist
    rewind(io_fid_conf)
    read(io_fid_conf,nml=param_atmos_phy_cp_kf_jmapplib,iostat=ierr)
    if( ierr > 0 ) then !--- fatal error
       log_error("ATMOS_PHY_CP_KF_JMAPPLIB_setup",*) 'Not appropriate names in namelist PARAM_ATMOS_PHY_CP_KF_JMAPPLIB. Check!'
       call prc_abort
    endif
    log_nml(param_atmos_phy_cp_kf_jmapplib)
 
 
    kmax = ke - ks + 1
    call conv_grid_ini( kmax )
    call conv_kf_const_ini( dx, pre00, dt )
 
    select case ( atmos_phy_cp_kf_jmapplib_type )
    case ( "KF" )
 
       call conv_kf_parm_ini( &
            dt,                                                           &
            kfrad_var_in        = atmos_phy_cp_kf_jmapplib_radius_type,   &
            cudt_in             = -1.0_rp,                                &
            dtlcl_fct_in        = atmos_phy_cp_kf_jmapplib_dtlcl_fact,    &
            dpmin_in            = atmos_phy_cp_kf_jmapplib_dpmin,         &
            cu_lifetime_min_in  = atmos_phy_cp_kf_jmapplib_dlifetime,     &
            shallow_lifetime_in = atmos_phy_cp_kf_jmapplib_slifetime,     &
            detlq_qr_r_in       = atmos_phy_cp_kf_jmapplib_detlq_qr_rate, &
            detic_qs_r_in       = atmos_phy_cp_kf_jmapplib_detic_qs_rate  )
 
!       call conv_kf_parm_get_stepcu(stepcu) ! (out)
 
    case ( "KF1701" )
 
       call conv_kf_parm_ini1701( &
            dt,                                                           &
            kfrad_var_in        = atmos_phy_cp_kf_jmapplib_radius_type,   &
            kfrad_in            = 1000.0_rp,                              &
            cudt_in             = -1.0_rp,                                &
            dpmin_in            = atmos_phy_cp_kf_jmapplib_dpmin,         &
            cu_lifetime_min_in  = atmos_phy_cp_kf_jmapplib_dlifetime,     &
            shallow_lifetime_in = atmos_phy_cp_kf_jmapplib_slifetime,     &
            detlq_qr_r_in       = atmos_phy_cp_kf_jmapplib_detlq_qr_rate, &
            detic_qs_r_in       = atmos_phy_cp_kf_jmapplib_detic_qs_rate, &
            tkemax_in           = 1.0_rp                                  )
            
!       call conv_kf_parm_get_stepcu1701(stepcu) ! (out)
 
    case default
       log_error("ATMOS_PHY_CP_KF_JMAPPLIB_setup",*) 'ATMOS_PHY_CP_KF_TYPE must be either "KF" or "KF1701"'
       call prc_abort
    end select
 
    call conv_kf_lut_ini
 
    allocate( ishall_counter(ia,ja), ideep_counter(ia,ja) )
    ishall_counter(:,:) = 0.0_rp
    ideep_counter(:,:) = 0.0_rp
 
#else
 
    log_error("ATMOS_PHY_CP_KF_JMAPPLIB_setup",*) 'To use "KF-JMAPPLIB", compile SCALE with "SCALE_ENABLE_JMAPPLIB=T" option.'
    call prc_abort
 
#endif
 
    return

References scale_const::const_pre00, scale_io::io_fid_conf, 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_jmapplib_finalize()

subroutine, public scale_atmos_phy_cp_kf_jmapplib::atmos_phy_cp_kf_jmapplib_finalize

Definition at line 192 of file scale_atmos_phy_cp_kf_jmapplib.F90.

    implicit none
 
#ifdef JMAPPLIB
    deallocate( ishall_counter, ideep_counter )
#endif
 
    return

atmos_phy_cp_kf_jmapplib_tendency()

subroutine, public scale_atmos_phy_cp_kf_jmapplib::atmos_phy_cp_kf_jmapplib_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)	W,
		real(rp), dimension (ka,ia,ja), intent(in)	TEMP,
		real(rp), dimension (ka,ia,ja), intent(in)	POTT,
		real(rp), dimension (ka,ia,ja), intent(in)	PRES,
		real(rp), dimension(ka,ia,ja), intent(in)	EXNER,
		real(rp), dimension (ka,ia,ja), intent(in)	QDRY,
		real(rp), dimension (ka,ia,ja), intent(in)	QV,
		real(rp), dimension (ka,ia,ja), intent(in)	QC,
		real(rp), dimension (ka,ia,ja), intent(in)	QI,
		real(rp), dimension (ia,ja), intent(in)	us,
		real(rp), dimension (ia,ja), intent(in)	PBLH,
		real(rp), dimension(ia,ja), intent(in)	SFLX_BUOY,
		real(rp), dimension (ka,ia,ja), intent(in)	CZ,
		real(rp), dimension (0:ka,ia,ja), intent(in)	FZ,
		real(rp), dimension (ka,ia,ja), intent(out)	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 (ka,ia,ja), intent(inout)	w0avg,
		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_prec,
		real(rp), dimension (ia,ja), intent(inout)	cloudtop,
		real(rp), dimension(ia,ja), intent(inout)	cloudbase
	)

ATMOS_PHY_CP_KF_JMAPPLIB_tendency calculate tendency by the virtical eddy viscosity.

Definition at line 217 of file scale_atmos_phy_cp_kf_jmapplib.F90.

    use scale_atmos_hydrometeor, only: &
       n_hyd, &
       i_hc, &
       i_hr, &
       i_hi, &
       i_hs, &
       i_hg
#ifdef JMAPPLIB
    use conv_kf_main, only: &
       conv_kf_run
    use conv_kf_main1701, only: &
       conv_kf_run1701
    use pbl_diag, only: &
       pbl_diag_buoyancy_excess
#endif
    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) :: W    (KA,IA,JA)
    real(RP), intent(in) :: TEMP (KA,IA,JA)
    real(RP), intent(in) :: POTT (KA,IA,JA)
    real(RP), intent(in) :: PRES (KA,IA,JA)
    real(RP), intent(in) :: EXNER(KA,IA,JA)
    real(RP), intent(in) :: QDRY (KA,IA,JA)
    real(RP), intent(in) :: QV   (KA,IA,JA)
    real(RP), intent(in) :: QC   (KA,IA,JA)
    real(RP), intent(in) :: QI   (KA,IA,JA)
    real(RP), intent(in) :: us       (IA,JA)
    real(RP), intent(in) :: PBLH     (IA,JA)
    real(RP), intent(in) :: SFLX_BUOY(IA,JA)
    real(RP), intent(in) :: CZ   (KA,IA,JA)
    real(RP), intent(in) :: FZ (0:KA,IA,JA)
 
    real(RP), intent(out)   :: 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) :: w0avg  (KA,IA,JA)
    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_prec(IA,JA)
    real(RP), intent(inout) :: cloudtop (IA,JA)
    real(RP), intent(inout) :: cloudbase(IA,JA)
 
#ifdef JMAPPLIB
 
    integer, parameter :: iconv_run = 1
 
    real(RP) :: tend_pt(KS:KE)
    real(RP) :: tend_qv(KS:KE)
    real(RP) :: tend_qc(KS:KE)
    real(RP) :: tend_qr(KS:KE)
    real(RP) :: tend_qi(KS:KE)
    real(RP) :: tend_qs(KS:KE)
    real(RP) :: tend_qg(KS:KE)
    real(RP) :: dens_lc(KS:KE)
    real(RP) :: qv_lc  (KS:KE)
    real(RP) :: qc_lc
    real(RP) :: qi_lc
    real(RP) :: z_f    (KS:KE)
    real(RP) :: dz_f   (KS:KE)
 
    real(RP) :: lcl_temp(IA,JA)
    real(RP) :: abe_out
    real(RP) :: tv_ex
    real(RP) :: ptv_ex
 
    integer :: k, i, j
    !---------------------------------------------------------------------------
 
    log_progress(*) "atmosphere / physics / cumulus / KF-JMAPPLIB"
 
    !$omp parallel do &
    !$omp private(tend_pt,tend_qv,tend_qc,tend_qr,tend_qi,tend_qs,tend_qg, &
    !$omp         dens_lc,qv_lc,qc_lc,qi_lc,z_f,dz_f, &
    !$omp         ptv_ex,tv_ex,abe_out)
    do j = js, je
    do i = is, ie
 
       do k = ks, ke
          dens_lc(k) = dens(k,i,j) * qdry(k,i,j)
          qv_lc(k) = qv(k,i,j) / qdry(k,i,j)
          tend_pt(k) = rhot_t(k,i,j) / dens(k,i,j)
          tend_qv(k) = rhoqv_t(k,i,j) / dens_lc(k)
          tend_qc(k) = rhoq_t(k,i,j,i_hc) / dens_lc(k)
          tend_qr(k) = rhoq_t(k,i,j,i_hr) / dens_lc(k)
          tend_qi(k) = rhoq_t(k,i,j,i_hi) / dens_lc(k)
          tend_qs(k) = rhoq_t(k,i,j,i_hs) / dens_lc(k)
          tend_qg(k) = rhoq_t(k,i,j,i_hg) / dens_lc(k)
 
          z_f(k) = cz(k,i,j) - fz(ks-1,i,j)
          dz_f(k) = fz(k,i,j) - fz(k-1,i,j)
       end do
 
       select case ( atmos_phy_cp_kf_jmapplib_type )
       case ( "KF" )
 
          call conv_kf_run( &
               iconv_run,                                                           & ! (in)
               u(ks:ke,i,j), v(ks:ke,i,j), w(ks:ke,i,j), temp(ks:ke,i,j), qv_lc(:), & ! (in)
               pres(ks:ke,i,j), exner(ks:ke,i,j), dens_lc(:), z_f(:), dz_f(:),      & ! (in)
               nca(i,j),                                                            & ! (inout)
               sflx_prec(i,j), sflx_rain(i,j), sflx_snow(i,j),                      & ! (inout)
               cloudtop(i,j), cloudbase(i,j), abe_out,                              & ! (inout)
               ishall_counter(i,j), ideep_counter(i,j),                             & ! (inout)
               lcl_temp(i,j),                                                       & ! (inout)
               w0avg(ks:ke,i,j),                                                    & ! (inout)
               tend_pt(:), tend_qv(:), tend_qc(:),                                  & ! (inout)
               tend_qi(:), tend_qr(:), tend_qs(:)                                   ) ! (inout)
 
          tend_qg(ks:ke) = 0.0_rp
 
       case ( "KF1701" )
 
          qc_lc = qc(ks,i,j) / qdry(ks,i,j)
          qi_lc = qi(ks,i,j) / qdry(ks,i,j)
          call pbl_diag_buoyancy_excess( &
               pott(ks,i,j), qv_lc(ks),qc_lc, qi_lc,              & ! (in)
               pblh(i,j), us(i,j), sflx_buoy(i,j), exner(ks,i,j), & ! (in)
               ptv_ex, tv_ex                                      ) ! (out)
 
          call conv_kf_run1701( &
               iconv_run,                                                           & ! (in)
               u(ks:ke,i,j), v(ks:ke,i,j), w(ks:ke,i,j), temp(ks:ke,i,j), qv_lc(:), & ! (in)
               pres(ks:ke,i,j), exner(ks:ke,i,j), dens_lc(:), z_f(:), dz_f(:),      & ! (in)
               pblh(i,j), tv_ex,                                                    & ! (in)
               nca(i,j),                                                            & ! (inout)
               sflx_prec(i,j), sflx_rain(i,j), sflx_snow(i,j),                      & ! (inout)
               cloudtop(i,j), cloudbase(i,j), abe_out,                              & ! (inout)
               ishall_counter(i,j), ideep_counter(i,j),                             & ! (inout)
               lcl_temp(i,j),                                                       & ! (inout)
               w0avg(ks:ke,i,j),                                                    & ! (inout)
               tend_pt(:), tend_qv(:), tend_qc(:),                                  & ! (inout)
               tend_qi(:), tend_qr(:),                                              & ! (inout)
               tend_qs(:), tend_qg(:)                                               ) ! (inout)
 
       end select
 
       if ( nca(i,j) < -1 ) nca(i,j) = -1
 
       do k = ks, ke
          rhot_t(k,i,j)      = tend_pt(k) * dens(k,i,j)
          rhoqv_t(k,i,j)      = tend_qv(k) * dens_lc(k)
          rhoq_t(k,i,j,i_hc) = tend_qc(k) * dens_lc(k)
          rhoq_t(k,i,j,i_hr) = tend_qr(k) * dens_lc(k)
          rhoq_t(k,i,j,i_hi) = tend_qi(k) * dens_lc(k)
          rhoq_t(k,i,j,i_hs) = tend_qs(k) * dens_lc(k)
          rhoq_t(k,i,j,i_hg) = tend_qg(k) * dens_lc(k)
          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) + rhoq_t(k,i,j,i_hg)
       end do
 
    end do
    end do
 
#endif
 
    return

References scale_atmos_hydrometeor::i_hc, scale_atmos_hydrometeor::i_hg, scale_atmos_hydrometeor::i_hi, scale_atmos_hydrometeor::i_hr, scale_atmos_hydrometeor::i_hs, and scale_atmos_hydrometeor::n_hyd.

Referenced by mod_atmos_phy_cp_driver::atmos_phy_cp_driver_calc_tendency().

Here is the caller graph for this function: