module ATMOSPHERE / Physics Cloud Microphysics More...

Functions/Subroutines
subroutine, public	atmos_phy_mp_kessler_setup (MP_TYPE)
	Setup. More...

subroutine, public	atmos_phy_mp_kessler (DENS, MOMZ, MOMX, MOMY, RHOT, QTRC, CCN, EVAPORATE, SFLX_rain, SFLX_snow)
	Cloud Microphysics. More...

subroutine, public	atmos_phy_mp_kessler_cloudfraction (cldfrac, QTRC, mask_criterion)
	Calculate Cloud Fraction. More...

subroutine, public	atmos_phy_mp_kessler_effectiveradius (Re, QTRC0, DENS0, TEMP0)
	Calculate Effective Radius. More...

subroutine, public	atmos_phy_mp_kessler_mixingratio (Qe, QTRC0)
	Calculate mixing ratio of each category. More...

Variables
real(rp), dimension(mp_qa), target, public	atmos_phy_mp_dens

Detailed Description

module ATMOSPHERE / Physics Cloud Microphysics

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

Author: Team SCALE

History

2012-01-14 (Y.Miyamoto) [new]
2012-03-23 (H.Yashiro) [mod] Explicit index parameter inclusion
2012-12-23 (H.Yashiro) [mod] Reconstruction
2015-09-08 (Y.Sato) [add] Add evaporated cloud number concentration

NAMELIST

PARAM_ATMOS_PHY_MP

name	type	default value	comment
MP_DOPRECIPITATION	logical	.true.	apply sedimentation (precipitation)?
MP_DONEGATIVE_FIXER	logical	.true.	apply negative fixer?
MP_NTMAX_SEDIMENTATION	integer	1	number of time step for sedimentation
MP_COUPLE_AEROSOL	logical	.false.	Consider CCN effect ?

History Output

name	description	unit	variable
Vterm_QR	terminal velocity of QR	m/s	vterm

Function/Subroutine Documentation

◆ atmos_phy_mp_kessler_setup()

subroutine, public scale_atmos_phy_mp_kessler::atmos_phy_mp_kessler_setup ( character(len=*), intent(in) MP_TYPE )

Setup.

Definition at line 80 of file scale_atmos_phy_mp_kessler.F90.

References atmos_phy_mp_dens, scale_const::const_dwatr, scale_grid::grid_cdz, scale_stdio::io_fid_conf, scale_stdio::io_fid_log, scale_stdio::io_l, scale_stdio::io_lnml, scale_grid_index::ka, scale_process::prc_mpistop(), and scale_time::time_dtsec_atmos_phy_mp.

Referenced by scale_atmos_phy_mp::atmos_phy_mp_setup().

     use scale_process, only: &
        prc_mpistop
     use scale_const, only: &
        const_dwatr
     use scale_time, only: &
        time_dtsec_atmos_phy_mp
     use scale_grid, only: &
        cdz => grid_cdz
     implicit none
 
     character(len=*), intent(in) :: mp_type
 
     namelist / param_atmos_phy_mp / &
        mp_doprecipitation, &
        mp_donegative_fixer, &
        mp_ntmax_sedimentation, &
        mp_couple_aerosol
 
     real(RP), parameter :: max_term_vel = 10.0_rp  !-- terminal velocity for calculate dt of sedimentation
     integer :: nstep_max
     integer :: ierr
     !---------------------------------------------------------------------------
 
     if( io_l ) write(io_fid_log,*)
     if( io_l ) write(io_fid_log,*) '++++++ Module[Cloud Microphysics] / Categ[ATMOS PHYSICS] / Origin[SCALElib]'
     if( io_l ) write(io_fid_log,*) '*** KESSLER-type 1-moment bulk 3 category'
 
     if ( mp_type /= 'KESSLER' ) then
        write(*,*) 'xxx ATMOS_PHY_MP_TYPE is not KESSLER. Check!'
        call prc_mpistop
     endif
 
     if (      i_qv <= 0 &
          .OR. i_qc <= 0 &
          .OR. i_qr <= 0 ) then
        write(*,*) 'xxx KESSLER needs QV, QC, QR tracer. Check!'
        call prc_mpistop
     endif
 
     allocate( factor_vterm(ka) )
 
     !--- read namelist
     rewind(io_fid_conf)
     read(io_fid_conf,nml=param_atmos_phy_mp,iostat=ierr)
     if( ierr < 0 ) then !--- missing
        if( io_l ) write(io_fid_log,*) '*** Not found namelist. Default used.'
     elseif( ierr > 0 ) then !--- fatal error
        write(*,*) 'xxx Not appropriate names in namelist PARAM_ATMOS_PHY_MP. Check!'
        call prc_mpistop
     endif
     if( io_lnml ) write(io_fid_log,nml=param_atmos_phy_mp)
 
     if( mp_couple_aerosol ) then
        write(*,*) 'xxx MP_aerosol_couple should be .false. for KESSLER type MP!'
        call prc_mpistop
     endif
 
     if ( io_l ) write(io_fid_log,*)
     if ( io_l ) write(io_fid_log,*) '*** Enable negative fixer?                : ', mp_donegative_fixer
     if ( io_l ) write(io_fid_log,*) '*** Enable sedimentation (precipitation)? : ', mp_doprecipitation
 
     atmos_phy_mp_dens(i_mp_qc) = const_dwatr
     atmos_phy_mp_dens(i_mp_qr) = const_dwatr
 
     nstep_max = int( ( time_dtsec_atmos_phy_mp * max_term_vel ) / minval( cdz ) )
     mp_ntmax_sedimentation = max( mp_ntmax_sedimentation, nstep_max )
 
     mp_nstep_sedimentation  = mp_ntmax_sedimentation
     mp_rnstep_sedimentation = 1.0_rp / real(mp_ntmax_sedimentation,kind=rp)
     mp_dtsec_sedimentation  = time_dtsec_atmos_phy_mp * mp_rnstep_sedimentation
 
     if ( io_l ) write(io_fid_log,*)
     if ( io_l ) write(io_fid_log,*) '*** Timestep of sedimentation is divided into : ', mp_ntmax_sedimentation, ' step'
     if ( io_l ) write(io_fid_log,*) '*** DT of sedimentation is : ', mp_dtsec_sedimentation, '[s]'
 
     return

Here is the call graph for this function:

Here is the caller graph for this function:

◆ atmos_phy_mp_kessler()

subroutine, public scale_atmos_phy_mp_kessler::atmos_phy_mp_kessler	(	real(rp), dimension(ka,ia,ja), intent(inout)	DENS,
		real(rp), dimension(ka,ia,ja), intent(inout)	MOMZ,
		real(rp), dimension(ka,ia,ja), intent(inout)	MOMX,
		real(rp), dimension(ka,ia,ja), intent(inout)	MOMY,
		real(rp), dimension(ka,ia,ja), intent(inout)	RHOT,
		real(rp), dimension(ka,ia,ja,qad), intent(inout)	QTRC,
		real(rp), dimension(ka,ia,ja), intent(in)	CCN,
		real(rp), dimension(ka,ia,ja), intent(out)	EVAPORATE,
		real(rp), dimension(ia,ja), intent(out)	SFLX_rain,
		real(rp), dimension(ia,ja), intent(out)	SFLX_snow
	)

Cloud Microphysics.

Definition at line 171 of file scale_atmos_phy_mp_kessler.F90.

References scale_atmos_phy_mp_common::atmos_phy_mp_negative_fixer(), scale_atmos_phy_mp_common::atmos_phy_mp_precipitation(), scale_atmos_phy_mp_common::atmos_phy_mp_saturation_adjustment(), scale_atmos_refstate::atmos_refstate_dens, scale_comm::comm_horizontal_mean(), scale_const::const_lhv, scale_grid_index::ie, scale_stdio::io_fid_log, scale_stdio::io_l, scale_grid_index::is, scale_grid_index::je, scale_grid_index::js, scale_grid_index::ke, scale_grid_index::ks, scale_prof::prof_rapend(), scale_prof::prof_rapstart(), scale_tracer::qa, and scale_time::time_dtsec_atmos_phy_mp.

Referenced by scale_atmos_phy_mp::atmos_phy_mp_setup().

     use scale_grid_index
     use scale_comm, only: &
        comm_horizontal_mean
     use scale_history, only: &
        hist_in
     use scale_tracer, only: &
        qad => qa
     use scale_atmos_thermodyn, only: &
        thermodyn_rhoe        => atmos_thermodyn_rhoe,       &
        thermodyn_rhot        => atmos_thermodyn_rhot,       &
        thermodyn_temp_pres_e => atmos_thermodyn_temp_pres_e
     use scale_atmos_phy_mp_common, only: &
        mp_negative_fixer        => atmos_phy_mp_negative_fixer,       &
        mp_precipitation         => atmos_phy_mp_precipitation,        &
        mp_saturation_adjustment => atmos_phy_mp_saturation_adjustment
     implicit none
 
     real(RP), intent(inout) :: dens(ka,ia,ja)
     real(RP), intent(inout) :: momz(ka,ia,ja)
     real(RP), intent(inout) :: momx(ka,ia,ja)
     real(RP), intent(inout) :: momy(ka,ia,ja)
     real(RP), intent(inout) :: rhot(ka,ia,ja)
     real(RP), intent(inout) :: qtrc(ka,ia,ja,qad)
     real(RP), intent(in)    :: ccn(ka,ia,ja)
     real(RP), intent(out)   :: evaporate(ka,ia,ja)   !--- number of evaporated cloud [/m3]
     real(RP), intent(out)   :: sflx_rain(ia,ja)
     real(RP), intent(out)   :: sflx_snow(ia,ja)
 
     real(RP) :: rhoe_t(ka,ia,ja)
     real(RP) :: qtrc_t(ka,ia,ja,qad)
     real(RP) :: rhoe  (ka,ia,ja)
     real(RP) :: temp  (ka,ia,ja)
     real(RP) :: pres  (ka,ia,ja)
 
     real(RP) :: vterm   (ka,ia,ja,qad) ! terminal velocity of each tracer [m/s]
     real(RP) :: flx_rain(ka,ia,ja)
     real(RP) :: flx_snow(ka,ia,ja)
     real(RP) :: wflux_rain(ka,ia,ja)
     real(RP) :: wflux_snow(ka,ia,ja)
     integer  :: step
 
     real(RP) :: rho_prof(ka) ! averaged profile of rho
     integer  :: k, i, j
     !---------------------------------------------------------------------------
 
     if( io_l ) write(io_fid_log,*) '*** Physics step: Cloud microphysics(kessler)'
 
     if ( first ) then
        ! Calculate collection factor for terminal velocity of QR
        call comm_horizontal_mean( rho_prof(:), dens(:,:,:) )
        rho_prof(:) = rho_prof(:) * 1.e-3_rp ! [kg/m3]->[g/cc]
 
        do k = ks, ke
           factor_vterm(k) = sqrt( rho_prof(ks)/rho_prof(k) )
        enddo
        first = .false.
     endif
 
     evaporate(:,:,:)   = 0.0_rp
 
     if ( mp_donegative_fixer ) then
        call mp_negative_fixer( dens(:,:,:),  & ! [INOUT]
                                rhot(:,:,:),  & ! [INOUT]
                                qtrc(:,:,:,:) ) ! [INOUT]
     endif
 
     call thermodyn_rhoe( rhoe(:,:,:),  & ! [OUT]
                          rhot(:,:,:),  & ! [IN]
                          qtrc(:,:,:,:) ) ! [IN]
 
     !##### MP Main #####
     rhoe_t(:,:,:)   = 0.0_rp
     qtrc_t(:,:,:,:) = 0.0_rp
 
     call mp_kessler( rhoe_t(:,:,:),   & ! [INOUT]
                      qtrc_t(:,:,:,:), & ! [INOUT]
                      rhoe(:,:,:),   & ! [INOUT]
                      qtrc(:,:,:,:), & ! [INOUT]
                      dens(:,:,:)    ) ! [IN]
 
     if ( mp_doprecipitation ) then
 
        flx_rain(:,:,:) = 0.0_rp
        flx_snow(:,:,:) = 0.0_rp
 
        vterm(:,:,:,:) = 0.0_rp
 
        do step = 1, mp_nstep_sedimentation
 
          call mp_kessler_vterm( vterm(:,:,:,:), & ! [OUT]
                                 dens(:,:,:),   & ! [IN]
                                 qtrc(:,:,:,:)  ) ! [IN]
 
          call thermodyn_temp_pres_e( temp(:,:,:),  & ! [OUT]
                                      pres(:,:,:),  & ! [OUT]
                                      dens(:,:,:),  & ! [IN]
                                      rhoe(:,:,:),  & ! [IN]
                                      qtrc(:,:,:,:) ) ! [IN]
 
          call mp_precipitation( wflux_rain(:,:,:),     & ! [OUT]
                                 wflux_snow(:,:,:),     & ! [OUT]
                                 dens(:,:,:),       & ! [INOUT]
                                 momz(:,:,:),       & ! [INOUT]
                                 momx(:,:,:),       & ! [INOUT]
                                 momy(:,:,:),       & ! [INOUT]
                                 rhoe(:,:,:),       & ! [INOUT]
                                 qtrc(:,:,:,:),     & ! [INOUT]
                                 vterm(:,:,:,:),     & ! [IN]
                                 temp(:,:,:),       & ! [IN]
                                 mp_dtsec_sedimentation ) ! [IN]
 
          do j = js, je
          do i = is, ie
          do k = ks-1, ke
             flx_rain(k,i,j) = flx_rain(k,i,j) + wflux_rain(k,i,j) * mp_rnstep_sedimentation
             flx_snow(k,i,j) = flx_snow(k,i,j) + wflux_snow(k,i,j) * mp_rnstep_sedimentation
          enddo
          enddo
          enddo
 
        enddo
 
     else
        vterm(:,:,:,:) = 0.0_rp
 
        flx_rain(:,:,:) = 0.0_rp
        flx_snow(:,:,:) = 0.0_rp
     endif
 
     call mp_saturation_adjustment( rhoe_t(:,:,:),       & ! [INOUT]
                                    qtrc_t(:,:,:,:),     & ! [INOUT]
                                    rhoe(:,:,:),       & ! [INOUT]
                                    qtrc(:,:,:,:),     & ! [INOUT]
                                    dens(:,:,:),       & ! [IN]
                                    flag_liquid = .true. ) ! [IN]
 
     call hist_in( vterm(:,:,:,i_qr), 'Vterm_QR', 'terminal velocity of QR', 'm/s' )
 
     !##### END MP Main #####
 
     call thermodyn_rhot( rhot(:,:,:),  & ! [OUT]
                          rhoe(:,:,:),  & ! [IN]
                          qtrc(:,:,:,:) ) ! [IN]
 
     if ( mp_donegative_fixer ) then
        call mp_negative_fixer( dens(:,:,:),  & ! [INOUT]
                                rhot(:,:,:),  & ! [INOUT]
                                qtrc(:,:,:,:) ) ! [INOUT]
     endif
 
     sflx_rain(:,:) = flx_rain(ks-1,:,:)
     sflx_snow(:,:) = flx_snow(ks-1,:,:)
 
     return

Here is the call graph for this function:

Here is the caller graph for this function:

◆ atmos_phy_mp_kessler_cloudfraction()

subroutine, public scale_atmos_phy_mp_kessler::atmos_phy_mp_kessler_cloudfraction	(	real(rp), dimension(ka,ia,ja), intent(out)	cldfrac,
		real(rp), dimension (ka,ia,ja,qad), intent(in)	QTRC,
		real(rp), intent(in)	mask_criterion
	)

Calculate Cloud Fraction.

Definition at line 517 of file scale_atmos_phy_mp_kessler.F90.

References scale_grid_index::ie, scale_grid_index::is, scale_grid_index::je, scale_grid_index::js, scale_grid_index::ke, scale_grid_index::ks, and scale_tracer::qa.

Referenced by scale_atmos_phy_mp::atmos_phy_mp_setup().

     use scale_grid_index
     use scale_tracer, only: &
        qad => qa
     implicit none
 
     real(RP), intent(out) :: cldfrac(ka,ia,ja)
     real(RP), intent(in)  :: qtrc   (ka,ia,ja,qad)
     real(RP), intent(in)  :: mask_criterion
 
     real(RP) :: qhydro
     integer  :: k, i, j, iq
     !---------------------------------------------------------------------------
 
     do j  = js, je
     do i  = is, ie
     do k  = ks, ke
        qhydro = 0.0_rp
        do iq = 1, mp_qa
           qhydro = qhydro + qtrc(k,i,j,i_mp2all(iq))
        enddo
        cldfrac(k,i,j) = 0.5_rp + sign(0.5_rp,qhydro-mask_criterion)
     enddo
     enddo
     enddo
 
     return

Here is the caller graph for this function:

◆ atmos_phy_mp_kessler_effectiveradius()

subroutine, public scale_atmos_phy_mp_kessler::atmos_phy_mp_kessler_effectiveradius	(	real(rp), dimension (ka,ia,ja,mp_qad), intent(out)	Re,
		real(rp), dimension(ka,ia,ja,qad), intent(in)	QTRC0,
		real(rp), dimension(ka,ia,ja), intent(in)	DENS0,
		real(rp), dimension(ka,ia,ja), intent(in)	TEMP0
	)

Calculate Effective Radius.

Definition at line 552 of file scale_atmos_phy_mp_kessler.F90.

References scale_tracer::mp_qa, and scale_tracer::qa.

Referenced by scale_atmos_phy_mp::atmos_phy_mp_setup().

     use scale_grid_index
     use scale_tracer, only: &
        qad => qa, &
        mp_qad => mp_qa
     implicit none
 
     real(RP), intent(out) :: re   (ka,ia,ja,mp_qad) ! effective radius          [cm]
     real(RP), intent(in)  :: qtrc0(ka,ia,ja,qad)    ! tracer mass concentration [kg/kg]
     real(RP), intent(in)  :: dens0(ka,ia,ja)        ! density                   [kg/m3]
     real(RP), intent(in)  :: temp0(ka,ia,ja)        ! temperature               [K]
 
     real(RP), parameter :: um2cm = 100.0_rp
     !---------------------------------------------------------------------------
 
     re(:,:,:,i_mp_qc) =   8.e-6_rp * um2cm
     re(:,:,:,i_mp_qr) = 100.e-6_rp * um2cm
 
     return

Here is the caller graph for this function:

◆ atmos_phy_mp_kessler_mixingratio()

subroutine, public scale_atmos_phy_mp_kessler::atmos_phy_mp_kessler_mixingratio	(	real(rp), dimension (ka,ia,ja,mp_qad), intent(out)	Qe,
		real(rp), dimension(ka,ia,ja,qad), intent(in)	QTRC0
	)

Calculate mixing ratio of each category.

Definition at line 577 of file scale_atmos_phy_mp_kessler.F90.

References scale_tracer::mp_qa, and scale_tracer::qa.

Referenced by scale_atmos_phy_mp::atmos_phy_mp_setup().

     use scale_grid_index
     use scale_tracer, only: &
        qad => qa, &
        mp_qad => mp_qa
     implicit none
 
     real(RP), intent(out) :: qe   (ka,ia,ja,mp_qad) ! mixing ratio of each cateory [kg/kg]
     real(RP), intent(in)  :: qtrc0(ka,ia,ja,qad)    ! tracer mass concentration [kg/kg]
     !---------------------------------------------------------------------------
 
     qe(:,:,:,i_mp_qc) = qtrc0(:,:,:,i_qc)
     qe(:,:,:,i_mp_qr) = qtrc0(:,:,:,i_qr)
 
     return