scale_atmos_phy_cp_common Module Reference

module atmosphere / physics / cumulus / Common More...

Functions/Subroutines
subroutine, public	atmos_phy_cp_common_setup
	Setup. More...
subroutine, public	atmos_phy_cp_common_wmean (KA, KS, KE, IA, IS, IE, JA, JS, JE, W, TIME_DTSEC, CP_DTSEC, W0_mean)
	ATMOS_PHY_CP_wmean running mean vertical wind velocity comment for W0 imported from WRF. More...

Detailed Description

module atmosphere / physics / cumulus / Common

Description

Common module for Cumulus convection parameterization Running mean of vertical wind velocity

Author

Team SCALE

NAMELIST

• PARAM_ATMOS_PHY_CP_COMMON

  name	type	default value	comment
  PARAM_ATMOS_PHY_CP_WADAPT	logical	.true.
  PARAM_ATMOS_PHY_CP_W_TIME	integer	16

  
  

History Output

No history output

Function/Subroutine Documentation

atmos_phy_cp_common_setup()

subroutine, public scale_atmos_phy_cp_common::atmos_phy_cp_common_setup

Setup.

Definition at line 52 of file scale_atmos_phy_cp_common.F90.

    use scale_prc, only: &
       prc_abort
    implicit none
 
    namelist / param_atmos_phy_cp_common / &
       param_atmos_phy_cp_wadapt, &
       param_atmos_phy_cp_w_time
 
    integer :: ierr
    !---------------------------------------------------------------------------
 
    log_newline
    log_info("ATMOS_PHY_CP_common_setup",*) 'Setup'
    log_info("ATMOS_PHY_CP_common_setup",*) 'CP-COMMON'
 
    !--- read namelist
    rewind(io_fid_conf)
    read(io_fid_conf,nml=param_atmos_phy_cp_common,iostat=ierr)
    if( ierr < 0 ) then !--- missing
       log_info("ATMOS_PHY_CP_common_setup",*) 'Not found namelist. Default used.'
    elseif( ierr > 0 ) then !--- fatal error
       log_error("ATMOS_PHY_CP_common_setup",*) 'Not appropriate names in namelist PARAM_ATMOS_PHY_CP_COMMON. Check!'
       call prc_abort
    endif
    log_nml(param_atmos_phy_cp_common)
 
    ! output parameter lists
    log_newline
    log_info("ATMOS_PHY_CP_common_setup",*) "Use running mean of w in adaptive timestep?     : ", param_atmos_phy_cp_wadapt
    log_info("ATMOS_PHY_CP_common_setup",*) "Fixed time scale for running mean of w          : ", param_atmos_phy_cp_w_time
 
    return

References 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_common_wmean()

subroutine, public scale_atmos_phy_cp_common::atmos_phy_cp_common_wmean	(	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)	W,
		real(dp), intent(in)	TIME_DTSEC,
		real(dp), intent(in)	CP_DTSEC,
		real(rp), dimension(ka,ia,ja), intent(inout)	W0_mean
	)

ATMOS_PHY_CP_wmean running mean vertical wind velocity comment for W0 imported from WRF.

TST IS THE NUMBER OF TIME STEPS IN 10 MINUTES...W0AVG IS CLOSE TO A RUNNING MEAN VERTICAL VELOCITY...NOTE THAT IF YOU CHANGE TST, IT WIL CHANGE THE FREQUENCY OF THE CONVECTIVE INTITIATION CHECK (SEE BELOW) NOTE THAT THE ORDERING OF VERTICAL LAYERS MUST BE REVERSED FOR W0AVG

Definition at line 102 of file scale_atmos_phy_cp_common.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)    :: W(KA,IA,JA)
    real(DP), intent(in)    :: TIME_DTSEC
    real(DP), intent(in)    :: CP_DTSEC
 
    real(RP), intent(inout) :: W0_mean(KA,IA,JA)
 
    real(RP) :: fact_W0_mean, fact_W0
 
    integer :: k, i, j
    !---------------------------------------------------------------------------
 
    if ( param_atmos_phy_cp_wadapt ) then
       fact_w0_mean = 2.0_rp * max(cp_dtsec,time_dtsec) - time_dtsec
       fact_w0      = time_dtsec
    else ! w_time is tuning parameter
       fact_w0_mean = real(param_atmos_phy_cp_w_time,rp)
       fact_w0      = 1.0_rp
    endif
 
    do j = js, je
    do i = is, ie
    do k = ks, ke
       w0_mean(k,i,j) = ( w0_mean(k,i,j) * fact_w0_mean &
                        + w(k,i,j)       * fact_w0      ) / ( fact_w0_mean + fact_w0 )
    enddo
    enddo
    enddo
 
    return

References 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::ke, scale_atmos_grid_cartesc_index::ks, and scale_precision::rp.

Referenced by mod_atmos_phy_cp_driver::atmos_phy_cp_driver_calc_tendency().

Here is the caller graph for this function: