module ATMOSPHERE / Physics Cumulus More...

Functions/Subroutines
subroutine, public	atmos_phy_cp_driver_setup
	Setup. More...

subroutine, public	atmos_phy_cp_driver_calc_tendency (update_flag)
	Driver. More...

Detailed Description

module ATMOSPHERE / Physics Cumulus

Description: Cumulus parameterization driver

Author: Team SCALE

NAMELIST

No namelist group

History Output

name	description	unit	variable
CBMFX	cloud base mass flux	kg/m2/s	MFLX_cloudbase
CUBASE	CP cloud base height	m	cloudbase
CUMFRC_DP	CP cloud fraction (deep)	1	cldfrac_dp
CUMFRC_SH	CP cloud fraction (shallow)	1	cldfrac_sh
CUMHGT	CP cloud top height	m	cloudtop
DENS_t_CP	tendency DENS in CP	kg/m3/s	DENS_t_CP
PREC_CP	surface precipitation rate by CP	kg/m2/s	SFLX_prec
QV_t_CP	tendency rho*QV in CP	kg/m3/s	RHOQV_t_CP
RAIN_CP	surface rain rate by CP	kg/m2/s	SFLX_rain
RHOT_t_CP	tendency RHOT in CP	K*kg/m3/s	RHOT_t_CP
SNOW_CP	surface snow rate by CP	kg/m2/s	SFLX_snow
kf_nca	advection or cumulus convection timescale for KF	s	kf_nca
{HYD_NAME}_t_CP	tendency rho*{HYD_NAME} in CP; {HYD_NAME} is QC, QR, QI, QS, QG, QH.	kg/m3/s	RHOHYD_t_CP
w0mean	running mean vertical wind velocity	kg/m2/s	w0mean

Function/Subroutine Documentation

◆ atmos_phy_cp_driver_setup()

subroutine, public mod_atmos_phy_cp_driver::atmos_phy_cp_driver_setup

Setup.

Definition at line 49 of file mod_atmos_phy_cp_driver.F90.

     use scale_prc, only: &
        prc_abort
     use scale_atmos_phy_cp_common, only: &
        atmos_phy_cp_common_setup
     use scale_atmos_phy_cp_kf, only: &
        atmos_phy_cp_kf_setup
     use mod_atmos_admin, only: &
        atmos_phy_cp_type, &
        atmos_sw_phy_cp
     use scale_time , only :&
        time_dtsec,             &
        time_dtsec_atmos_phy_cp
     use scale_atmos_grid_cartesc_real, only: &
        atmos_grid_cartesc_real_cz, &
        atmos_grid_cartesc_real_area
     use scale_atmos_hydrometeor, only: &
        atmos_hydrometeor_ice_phase
     implicit none
  
     logical :: warmrain
     !---------------------------------------------------------------------------
  
     log_newline
     log_info("ATMOS_PHY_CP_driver_setup",*) 'Setup'
  
     if ( atmos_sw_phy_cp ) then
  
        ! setup library component
        call atmos_phy_cp_common_setup
        select case ( atmos_phy_cp_type )
        case ( 'KF' )
           warmrain = ( .not. atmos_hydrometeor_ice_phase )
           call atmos_phy_cp_kf_setup( ka, ks, ke, ia, is, ie, ja, js, je, &
                                       atmos_grid_cartesc_real_cz, atmos_grid_cartesc_real_area, &
                                       warmrain                                                  )
        case default
           log_error("ATMOS_PHY_CP_driver_setup",*) 'ATMOS_PHY_CP_TYPE (', trim(atmos_phy_cp_type), ') is invalid. Check!'
           call prc_abort
        end select
  
     else
        log_info("ATMOS_PHY_CP_driver_setup",*) 'this component is never called.'
     endif
  
     return

Referenced by mod_atmos_driver::atmos_driver_setup().

Here is the call graph for this function:

Here is the caller graph for this function:

◆ atmos_phy_cp_driver_calc_tendency()

subroutine, public mod_atmos_phy_cp_driver::atmos_phy_cp_driver_calc_tendency ( logical, intent(in) update_flag )

Driver.

Definition at line 99 of file mod_atmos_phy_cp_driver.F90.

     use scale_const, only: &
        pre00 => const_pre00
     use scale_statistics, only: &
        statistics_checktotal, &
        statistics_total
     use scale_atmos_grid_cartesc_real, only: &
        atmos_grid_cartesc_real_vol,       &
        atmos_grid_cartesc_real_totvol,    &
        atmos_grid_cartesc_real_volwxy,    &
        atmos_grid_cartesc_real_totvolwxy, &
        atmos_grid_cartesc_real_volzuy,    &
        atmos_grid_cartesc_real_totvolzuy, &
        atmos_grid_cartesc_real_volzxv,    &
        atmos_grid_cartesc_real_totvolzxv
     use scale_file_history, only: &
        file_history_in
     use scale_time , only :&
        time_dtsec, &
        time_dtsec_atmos_phy_cp
     use scale_atmos_grid_cartesc_real, only: &
        fz => atmos_grid_cartesc_real_fz
     use scale_atmos_hydrometeor, only: &
        hyd_name, &
        cv_water
     use scale_atmos_phy_cp_kf, only: &
        atmos_phy_cp_kf_tendency
     use scale_atmos_phy_cp_common, only: &
        atmos_phy_cp_common_wmean
     use mod_atmos_phy_mp_vars, only: &
        qs_mp, &
        qe_mp
     use mod_atmos_phy_mp_driver, only: &
        atmos_phy_mp_driver_qhyd2qtrc
     use mod_atmos_admin, only: &
        atmos_phy_cp_type
     use mod_atmos_vars, only: &
        dens   => dens_av, &
        momz   => momz_av, &
        momx   => momx_av, &
        momy   => momy_av, &
        rhot   => rhot_av, &
        qtrc   => qtrc_av, &
        dens_t => dens_tp, &
        momz_t => momz_tp, &
        momx_t => momx_tp, &
        momy_t => momy_tp, &
        rhot_t => rhot_tp, &
        rhoq_t => rhoq_tp, &
        u,                 &
        v,                 &
        w,                 &
        temp,              &
        pres,              &
        qdry,              &
        qv
     use scale_atmos_hydrometeor, only: &
        n_hyd
     use mod_atmos_phy_cp_vars, only: &
        dens_t_cp      => atmos_phy_cp_dens_t,         &
        rhot_t_cp      => atmos_phy_cp_rhot_t,         &
        rhoqv_t_cp     => atmos_phy_cp_rhoqv_t,        &
        rhohyd_t_cp    => atmos_phy_cp_rhohyd_t,       &
        mflx_cloudbase => atmos_phy_cp_mflx_cloudbase, &
        sflx_rain      => atmos_phy_cp_sflx_rain,      &  ! convective rain [kg/m2/s]
        sflx_snow      => atmos_phy_cp_sflx_snow,      &  ! convective snow [kg/m2/s]
        sflx_engi      => atmos_phy_cp_sflx_engi,      &  ! internal energy flux [J/m2/s]
        cloudtop       => atmos_phy_cp_cloudtop,       &  ! cloud top height [m]
        cloudbase      => atmos_phy_cp_cloudbase,      &  ! cloud base height [m]
        cldfrac_dp     => atmos_phy_cp_cldfrac_dp,     &  ! cloud fraction (deep convection) (0-1)
        cldfrac_sh     => atmos_phy_cp_cldfrac_sh,     &  ! cloud fraction (shallow convection) (0-1)
        w0mean         => atmos_phy_cp_w0mean,         &  ! running mean vertical wind velocity [m/s]
        kf_nca         => atmos_phy_cp_kf_nca             ! advection/cumulus convection timescale/dt for KF [step]
     implicit none
  
     logical, intent(in) :: update_flag
  
     real(RP) :: RHOQ_t_CP(KA,IA,JA,QS_MP:QE_MP)
  
     real(RP) :: SFLX_prec(IA,JA)
  
     integer  :: k, i, j, iq
     !---------------------------------------------------------------------------
  
     ! temporal running mean of vertical velocity
     call atmos_phy_cp_common_wmean( ka, ks, ke, ia, is, ie, ja, js, je, &
                                     w(:,:,:),                            & ! [IN]
                                     time_dtsec, time_dtsec_atmos_phy_cp, & ! [IN]
                                     w0mean(:,:,:)                        ) ! [INOUT]
     call file_history_in( w0mean(:,:,:), 'w0mean', 'running mean vertical wind velocity', 'kg/m2/s', fill_halo=.true. )
  
     if ( update_flag ) then ! update
        select case ( atmos_phy_cp_type )
        case ( 'KF' )
           call atmos_phy_cp_kf_tendency( ka, ks, ke, ia, is, ie, ja, js, je, &
                                          dens(:,:,:),                              & ! [IN]
                                          u(:,:,:), v(:,:,:),                       & ! [IN]
                                          rhot(:,:,:), temp(:,:,:), pres(:,:,:),    & ! [IN]
                                          qdry(:,:,:), qv(:,:,:),                   & ! [IN]
                                          w0mean(:,:,:),                            & ! [IN]
                                          fz(:,:,:),                                & ! [IN]
                                          time_dtsec_atmos_phy_cp,                  & ! [IN]
                                          dens_t_cp(:,:,:),                         & ! [INOUT]
                                          rhot_t_cp(:,:,:),                         & ! [INOUT]
                                          rhoqv_t_cp(:,:,:), rhohyd_t_cp(:,:,:,:),  & ! [INOUT]
                                          kf_nca(:,:),                              & ! [INOUT]
                                          sflx_rain(:,:), sflx_snow(:,:),           & ! [INOUT]
                                          sflx_engi(:,:),                           & ! [INOUT]
                                          cloudtop(:,:), cloudbase(:,:),            & ! [INOUT]
                                          cldfrac_dp(:,:,:), cldfrac_sh(:,:,:)      ) ! [INOUT]
  
        end select
  
 !OCL XFILL
        !$omp parallel do
        do j  = js, je
        do i  = is, ie
           mflx_cloudbase(i,j) = 0.0_rp
        enddo
        enddo
  
        ! diagnose tendency of number concentration
  
        !$omp parallel do
        do j = js, je
        do i = is, ie
           sflx_prec(i,j) = sflx_rain(i,j) + sflx_snow(i,j)
        end do
        end do
  
        call file_history_in( mflx_cloudbase(:,:),   'CBMFX',     'cloud base mass flux',             'kg/m2/s', fill_halo=.true. )
        call file_history_in( sflx_rain(:,:),   'RAIN_CP',   'surface rain rate by CP',          'kg/m2/s', fill_halo=.true. )
        call file_history_in( sflx_snow(:,:),   'SNOW_CP',   'surface snow rate by CP',          'kg/m2/s', fill_halo=.true. )
        call file_history_in( sflx_prec(:,:),   'PREC_CP',   'surface precipitation rate by CP', 'kg/m2/s', fill_halo=.true. )
        call file_history_in( cloudtop(:,:),   'CUMHGT',    'CP cloud top height',              'm',       fill_halo=.true. )
        call file_history_in( cloudbase(:,:),   'CUBASE',    'CP cloud base height',             'm',       fill_halo=.true. )
        call file_history_in( cldfrac_dp(:,:,:), 'CUMFRC_DP', 'CP cloud fraction (deep)',         '1',       fill_halo=.true. )
        call file_history_in( cldfrac_sh(:,:,:), 'CUMFRC_SH', 'CP cloud fraction (shallow)',      '1',       fill_halo=.true. )
        call file_history_in( kf_nca(:,:),   'kf_nca',    'advection or cumulus convection timescale for KF', 's', fill_halo=.true. )
  
        call file_history_in( dens_t_cp(:,:,:), 'DENS_t_CP', 'tendency DENS in CP', 'kg/m3/s'  , fill_halo=.true. )
        call file_history_in( rhot_t_cp(:,:,:), 'RHOT_t_CP', 'tendency RHOT in CP', 'K*kg/m3/s', fill_halo=.true. )
  
        call file_history_in( rhoqv_t_cp(:,:,:), 'QV_t_CP',  'tendency rho*QV in CP', 'kg/m3/s', fill_halo=.true. )
        do iq = 1, n_hyd
           call file_history_in( rhohyd_t_cp(:,:,:,iq), trim(hyd_name(iq))//'_t_CP', &
                                 'tendency rho*'//trim(hyd_name(iq))//' in CP', 'kg/m3/s', fill_halo=.true. )
        enddo
  
     endif ! update
  
     !$omp parallel do private(i,j,k) OMP_SCHEDULE_ collapse(2)
     do j = jsb, jeb
     do i = isb, ieb
     do k = ks, ke
        dens_t(k,i,j) = dens_t(k,i,j) + dens_t_cp(k,i,j)
        rhot_t(k,i,j) = rhot_t(k,i,j) + rhot_t_cp(k,i,j)
     enddo
     enddo
     enddo
  
     call atmos_phy_mp_driver_qhyd2qtrc( ka, ks, ke, ia, is, ie, ja, js, je, &
                                         rhoqv_t_cp(:,:,:), rhohyd_t_cp(:,:,:,:), & ! [IN]
                                         rhoq_t_cp(:,:,:,qs_mp:qe_mp)             ) ! [OUT]
  
     do iq = qs_mp, qe_mp
     !$omp parallel do private(i,j,k) OMP_SCHEDULE_ collapse(3)
     do j  = js, je
     do i  = is, ie
     do k  = ks, ke
        rhoq_t(k,i,j,iq) = rhoq_t(k,i,j,iq) + rhoq_t_cp(k,i,j,iq)
     enddo
     enddo
     enddo
     enddo
  
     if ( statistics_checktotal ) then
        call statistics_total( ka, ks, ke, ia, is, ie, ja, js, je, &
                               dens_t_cp(:,:,:), 'DENS_t_CP',       &
                               atmos_grid_cartesc_real_vol(:,:,:),  &
                               atmos_grid_cartesc_real_totvol       )
        call statistics_total( ka, ks, ke, ia, is, ie, ja, js, je, &
                               rhot_t_cp(:,:,:), 'RHOT_t_CP',       &
                               atmos_grid_cartesc_real_vol(:,:,:),  &
                               atmos_grid_cartesc_real_totvol       )
  
        do iq = qs_mp, qe_mp
           call statistics_total( ka, ks, ke, ia, is, ie, ja, js, je, &
                                  rhoq_t_cp(:,:,:,iq), trim(tracer_name(iq))//'_t_CP', &
                                  atmos_grid_cartesc_real_vol(:,:,:),                  &
                                  atmos_grid_cartesc_real_totvol                       )
        enddo
     endif
  
     return

Referenced by mod_atmos_driver::atmos_driver_calc_tendency().