mod_atmos_phy_mp_driver Module Reference

module ATMOSPHERE / Physics Cloud Microphysics More...

Functions/Subroutines
subroutine, public	atmos_phy_mp_driver_config
	Config. More...
subroutine, public	atmos_phy_mp_driver_setup
	Setup. More...
subroutine, public	atmos_phy_mp_driver_resume
	resume More...
subroutine, public	atmos_phy_mp_driver (update_flag)
	Driver. More...

Detailed Description

module ATMOSPHERE / Physics Cloud Microphysics

Description

Cloud Microphysics driver

Author

Team SCALE

History

• 2013-12-06 (S.Nishizawa) [new]

NAMELIST

• No namelist group

History Output

name	description	unit	variable
DENS_t_MP	tendency DENS in MP	kg/m3/s	DENS_t_MP
EVAPORATE	evaporated cloud number	num/m3/s	EVAPORATE
MOMX_t_MP	tendency MOMX in MP	kg/m2/s2	MOMX_t_MP
MOMY_t_MP	tendency MOMY in MP	kg/m2/s2	MOMY_t_MP
MOMZ_t_MP	tendency MOMZ in MP	kg/m2/s2	MOMZ_t_MP
PREC_MP	surface precipitation rate by MP	kg/m2/s	precip
RAIN_MP	surface rain rate by MP	kg/m2/s	SFLX_rain
RHOT_t_MP	tendency RHOT in MP	K*kg/m3/s	RHOT_t_MP
Re_QC	effective radius cloud water	cm	Re
Re_QG	effective radius graupel	cm	Re
Re_QH	effective radius hail	cm	Re
Re_QI	effective radius ice water	cm	Re
Re_QR	effective radius rain	cm	Re
Re_QS	effective radius snow	cm	Re
SNOW_MP	surface snow rate by MP	kg/m2/s	SFLX_snow
TRACER_NAME_t_MP	tendency rho*TRACER_NAMEin MP; TRACER_NAME depends on the physics schemes, e.g., QV, QC, QR.	kg/m3/s	RHOQ_t_MP

Function/Subroutine Documentation

atmos_phy_mp_driver_config()

subroutine, public mod_atmos_phy_mp_driver::atmos_phy_mp_driver_config

(

)

Config.

Definition at line 53 of file mod_atmos_phy_mp_driver.f90.

References scale_atmos_phy_mp::atmos_phy_mp_config(), mod_atmos_admin::atmos_phy_mp_type, mod_atmos_admin::atmos_sw_phy_mp, scale_stdio::io_fid_log, and scale_stdio::io_l.

Referenced by mod_atmos_driver::atmos_driver_config().

    use scale_atmos_phy_mp, only: &
       atmos_phy_mp_config
    use mod_atmos_admin, only: &
       atmos_phy_mp_type, &
       atmos_sw_phy_mp
    implicit none
    !---------------------------------------------------------------------------

    if( io_l ) write(io_fid_log,*)
    if( io_l ) write(io_fid_log,*) '++++++ Module[CONFIG] / Categ[ATMOS PHY_MP] / Origin[SCALE-RM]'

    if ( atmos_sw_phy_mp ) then
       call atmos_phy_mp_config( atmos_phy_mp_type )
    end if

    return

Here is the call graph for this function:

Here is the caller graph for this function:

atmos_phy_mp_driver_setup()

subroutine, public mod_atmos_phy_mp_driver::atmos_phy_mp_driver_setup

(

)

Setup.

Definition at line 74 of file mod_atmos_phy_mp_driver.f90.

References scale_atmos_phy_mp::atmos_phy_mp_setup, mod_atmos_phy_mp_vars::atmos_phy_mp_sflx_rain, mod_atmos_phy_mp_vars::atmos_phy_mp_sflx_snow, mod_atmos_admin::atmos_phy_mp_type, mod_atmos_admin::atmos_sw_phy_mp, scale_stdio::io_fid_log, and scale_stdio::io_l.

Referenced by mod_atmos_driver::atmos_driver_setup().

    use scale_atmos_phy_mp, only: &
       atmos_phy_mp_setup
    use mod_atmos_admin, only: &
       atmos_phy_mp_type, &
       atmos_sw_phy_mp
    use mod_atmos_phy_mp_vars, only: &
       sflx_rain => atmos_phy_mp_sflx_rain, &
       sflx_snow => atmos_phy_mp_sflx_snow
    implicit none
    !---------------------------------------------------------------------------

    if( io_l ) write(io_fid_log,*)
    if( io_l ) write(io_fid_log,*) '++++++ Module[DRIVER] / Categ[ATMOS PHY_MP] / Origin[SCALE-RM]'

    if ( atmos_sw_phy_mp ) then

       ! setup library component
       call atmos_phy_mp_setup

    else

       if( io_l ) write(io_fid_log,*) '*** this component is never called.'
       if( io_l ) write(io_fid_log,*) '*** SFLX_rain and SFLX_snow is set to zero.'
       sflx_rain(:,:) = 0.0_rp
       sflx_snow(:,:) = 0.0_rp

    endif

    return

Here is the caller graph for this function:

atmos_phy_mp_driver_resume()

subroutine, public mod_atmos_phy_mp_driver::atmos_phy_mp_driver_resume

(

)

resume

Definition at line 108 of file mod_atmos_phy_mp_driver.f90.

References atmos_phy_mp_driver(), mod_atmos_admin::atmos_sw_phy_mp, scale_prof::prof_rapend(), and scale_prof::prof_rapstart().

Referenced by mod_atmos_driver::atmos_driver_resume1().

    use mod_atmos_admin, only: &
       atmos_sw_phy_mp
    implicit none

    if ( atmos_sw_phy_mp ) then

       ! run once (only for the diagnostic value)
       call prof_rapstart('ATM_Microphysics', 1)
       call atmos_phy_mp_driver( update_flag = .true. )
       call prof_rapend  ('ATM_Microphysics', 1)

    end if

    return

Here is the call graph for this function:

Here is the caller graph for this function:

atmos_phy_mp_driver()

subroutine, public mod_atmos_phy_mp_driver::atmos_phy_mp_driver

(

logical, intent(in)

update_flag

)

Driver.

Definition at line 127 of file mod_atmos_phy_mp_driver.f90.

References mod_atmos_phy_ae_vars::atmos_phy_ae_ccn_t, scale_atmos_phy_mp::atmos_phy_mp, mod_atmos_phy_mp_vars::atmos_phy_mp_dens_t, scale_atmos_phy_mp::atmos_phy_mp_effectiveradius, mod_atmos_phy_mp_vars::atmos_phy_mp_evaporate, mod_atmos_phy_mp_vars::atmos_phy_mp_momx_t, mod_atmos_phy_mp_vars::atmos_phy_mp_momy_t, mod_atmos_phy_mp_vars::atmos_phy_mp_momz_t, mod_atmos_phy_mp_vars::atmos_phy_mp_rhoq_t, mod_atmos_phy_mp_vars::atmos_phy_mp_rhot_t, mod_atmos_phy_mp_vars::atmos_phy_mp_sflx_rain, mod_atmos_phy_mp_vars::atmos_phy_mp_sflx_snow, mod_atmos_vars::dens, mod_atmos_vars::dens_av, mod_atmos_vars::dens_tp, scale_atmos_hydrometeor::i_hc, scale_atmos_hydrometeor::i_hg, scale_atmos_hydrometeor::i_hh, scale_atmos_hydrometeor::i_hi, scale_atmos_hydrometeor::i_hr, scale_atmos_hydrometeor::i_hs, scale_grid_index::ia, scale_grid_index::ie, scale_grid_index::is, scale_grid_index::ja, scale_grid_index::je, scale_grid_index::js, scale_grid_index::ka, scale_grid_index::ke, scale_grid_index::ks, mod_atmos_vars::momx, mod_atmos_vars::momx_av, mod_atmos_vars::momx_tp, mod_atmos_vars::momy, mod_atmos_vars::momy_av, mod_atmos_vars::momy_tp, mod_atmos_vars::momz, mod_atmos_vars::momz_av, mod_atmos_vars::momz_tp, scale_atmos_hydrometeor::n_hyd, scale_tracer::qa, scale_atmos_phy_mp::qe_mp, scale_atmos_phy_mp::qs_mp, mod_atmos_vars::qtrc, mod_atmos_vars::qtrc_av, mod_atmos_vars::rhoq_tp, mod_atmos_vars::rhot, mod_atmos_vars::rhot_av, mod_atmos_vars::rhot_tp, scale_rm_statistics::statistics_checktotal, mod_atmos_vars::temp, scale_time::time_dtsec_atmos_phy_mp, and scale_tracer::tracer_name.

Referenced by mod_atmos_driver::atmos_driver(), and atmos_phy_mp_driver_resume().

    use scale_time, only: &
       dt_mp => time_dtsec_atmos_phy_mp
    use scale_rm_statistics, only: &
       statistics_checktotal, &
       stat_total
    use scale_history, only: &
       hist_in
    use scale_atmos_hydrometeor, only: &
       n_hyd, &
       i_hc,  &
       i_hr,  &
       i_hi,  &
       i_hs,  &
       i_hg,  &
       i_hh
    use scale_atmos_phy_mp, only: &
       atmos_phy_mp,                 &
       atmos_phy_mp_effectiveradius, &
       qs_mp,                        &
       qe_mp
    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, &
       temp
    use mod_atmos_phy_mp_vars, only: &
       dens_t_mp => atmos_phy_mp_dens_t,    &
       momz_t_mp => atmos_phy_mp_momz_t,    &
       momx_t_mp => atmos_phy_mp_momx_t,    &
       momy_t_mp => atmos_phy_mp_momy_t,    &
       rhot_t_mp => atmos_phy_mp_rhot_t,    &
       rhoq_t_mp => atmos_phy_mp_rhoq_t,    &
       evaporate => atmos_phy_mp_evaporate, &
       sflx_rain => atmos_phy_mp_sflx_rain, &
       sflx_snow => atmos_phy_mp_sflx_snow
    use mod_atmos_phy_ae_vars, only: &
       ccn_t => atmos_phy_ae_ccn_t
    implicit none

    logical, intent(in) :: update_flag

    real(RP) :: DENS0(KA,IA,JA)
    real(RP) :: MOMZ0(KA,IA,JA)
    real(RP) :: MOMX0(KA,IA,JA)
    real(RP) :: MOMY0(KA,IA,JA)
    real(RP) :: RHOT0(KA,IA,JA)
    real(RP) :: QTRC0(KA,IA,JA,QA)
    real(RP) :: CCN(KA,IA,JA)

    real(RP) :: Re    (KA,IA,JA,N_HYD) ! effective radius
    real(RP) :: precip(IA,JA)
    real(RP) :: total ! dummy

    integer :: k, i, j, iq
    !---------------------------------------------------------------------------

    if ( update_flag ) then

!OCL XFILL
       !$omp parallel do default(none) &
       !$omp shared(JA,IA,KA,DENS0,MOMZ0,MOMX0,MOMY0,RHOT0,DENS,MOMZ,MOMX,MOMY,RHOT) &
       !$omp private(i,j,k) OMP_SCHEDULE_ collapse(2)
       do j  = 1, ja
       do i  = 1, ia
       do k  = 1, ka
          dens0(k,i,j) = dens(k,i,j) ! save
          momz0(k,i,j) = momz(k,i,j) ! save
          momx0(k,i,j) = momx(k,i,j) ! save
          momy0(k,i,j) = momy(k,i,j) ! save
          rhot0(k,i,j) = rhot(k,i,j) ! save
       enddo
       enddo
       enddo

!OCL XFILL
       do iq = 1, qa
       do j  = 1, ja
       do i  = 1, ia
       do k  = 1, ka
          qtrc0(k,i,j,iq) = qtrc(k,i,j,iq) ! save
       enddo
       enddo
       enddo
       enddo

       ccn(:,:,:) = ccn_t(:,:,:) * dt_mp

       call atmos_phy_mp( dens0(:,:,:),   & ! [INOUT]
                          momz0(:,:,:),   & ! [INOUT]
                          momx0(:,:,:),   & ! [INOUT]
                          momy0(:,:,:),   & ! [INOUT]
                          rhot0(:,:,:),   & ! [INOUT]
                          qtrc0(:,:,:,:), & ! [INOUT]
                          ccn(:,:,:),   & ! [IN]
                          evaporate(:,:,:),   & ! [OUT]
                          sflx_rain(:,:),     & ! [OUT]
                          sflx_snow(:,:)      ) ! [OUT]

!OCL XFILL
       !$omp parallel do default(none) private(i,j,k) OMP_SCHEDULE_ collapse(2) &
       !$omp shared(JS,JE,IS,IE,KS,KE,DENS_t_MP,DENS0,DENS,MOMZ_t_MP,MOMZ0,MOMZ,MOMX_t_MP,MOMX0) &
       !$omp shared(MOMX,MOMY_t_MP,MOMY0,MOMY,RHOT_t_MP,RHOT0,RHOT,dt_MP)
       do j  = js, je
       do i  = is, ie
       do k  = ks, ke
          dens_t_mp(k,i,j) = ( dens0(k,i,j) - dens(k,i,j) ) / dt_mp
          momz_t_mp(k,i,j) = ( momz0(k,i,j) - momz(k,i,j) ) / dt_mp
          momx_t_mp(k,i,j) = ( momx0(k,i,j) - momx(k,i,j) ) / dt_mp
          momy_t_mp(k,i,j) = ( momy0(k,i,j) - momy(k,i,j) ) / dt_mp
          rhot_t_mp(k,i,j) = ( rhot0(k,i,j) - rhot(k,i,j) ) / dt_mp
       enddo
       enddo
       enddo

!OCL XFILL
       do iq = qs_mp, qe_mp
       !$omp parallel do default(none) &
       !$omp shared(JS,JE,IS,IE,KS,KE,RHOQ_t_MP,iq,QTRC0,QTRC,DENS0,DENS,dt_MP) &
       !$omp private(i,j,k) OMP_SCHEDULE_ collapse(2)
       do j  = js, je
       do i  = is, ie
       do k  = ks, ke
          rhoq_t_mp(k,i,j,iq) = ( qtrc0(k,i,j,iq) * dens0(k,i,j) &
                                - qtrc(k,i,j,iq) * dens(k,i,j) ) / dt_mp
       enddo
       enddo
       enddo
       enddo

!OCL XFILL
       do j  = js, je
       do i  = is, ie
          precip(i,j) = sflx_rain(i,j) + sflx_snow(i,j)
       end do
       end do

       call hist_in( sflx_rain(:,:),   'RAIN_MP',   'surface rain rate by MP',          'kg/m2/s',  nohalo=.true. )
       call hist_in( sflx_snow(:,:),   'SNOW_MP',   'surface snow rate by MP',          'kg/m2/s',  nohalo=.true. )
       call hist_in( precip(:,:),   'PREC_MP',   'surface precipitation rate by MP', 'kg/m2/s',  nohalo=.true. )
       call hist_in( evaporate(:,:,:), 'EVAPORATE', 'evaporated cloud number',          'num/m3/s', nohalo=.true. )

       call hist_in( dens_t_mp(:,:,:), 'DENS_t_MP', 'tendency DENS in MP', 'kg/m3/s'  , nohalo=.true. )
       call hist_in( momz_t_mp(:,:,:), 'MOMZ_t_MP', 'tendency MOMZ in MP', 'kg/m2/s2' , nohalo=.true. )
       call hist_in( momx_t_mp(:,:,:), 'MOMX_t_MP', 'tendency MOMX in MP', 'kg/m2/s2' , nohalo=.true. )
       call hist_in( momy_t_mp(:,:,:), 'MOMY_t_MP', 'tendency MOMY in MP', 'kg/m2/s2' , nohalo=.true. )
       call hist_in( rhot_t_mp(:,:,:), 'RHOT_t_MP', 'tendency RHOT in MP', 'K*kg/m3/s', nohalo=.true. )

       do iq = qs_mp, qe_mp
          call hist_in( rhoq_t_mp(:,:,:,iq), trim(tracer_name(iq))//'_t_MP', &
                        'tendency rho*'//trim(tracer_name(iq))//'in MP', 'kg/m3/s', nohalo=.true. )
       enddo

       call atmos_phy_mp_effectiveradius( re(:,:,:,:), & ! [OUT]
                                          qtrc(:,:,:,:), & ! [IN]
                                          dens(:,:,:),   & ! [IN]
                                          temp(:,:,:)    ) ! [IN]

       call hist_in( re(:,:,:,i_hc), 'Re_QC', 'effective radius cloud water', 'cm', nohalo=.true. )
       call hist_in( re(:,:,:,i_hr), 'Re_QR', 'effective radius rain',        'cm', nohalo=.true. )
       call hist_in( re(:,:,:,i_hi), 'Re_QI', 'effective radius ice water',   'cm', nohalo=.true. )
       call hist_in( re(:,:,:,i_hs), 'Re_QS', 'effective radius snow',        'cm', nohalo=.true. )
       call hist_in( re(:,:,:,i_hg), 'Re_QG', 'effective radius graupel',     'cm', nohalo=.true. )
       call hist_in( re(:,:,:,i_hh), 'Re_QH', 'effective radius hail',        'cm', nohalo=.true. )
    endif

    !$omp parallel do default(none) private(i,j,k) OMP_SCHEDULE_ collapse(2) &
    !$omp shared(JS,JE,IS,IE,KS,KE,DENS_t,DENS_t_MP,MOMZ_t,MOMZ_t_MP,MOMX_t,MOMX_t_MP,MOMY_t) &
    !$omp shared(MOMY_t_MP,RHOT_t,RHOT_t_MP)
    do j = js, je
    do i = is, ie
    do k = ks, ke
       dens_t(k,i,j) = dens_t(k,i,j) + dens_t_mp(k,i,j)
       momz_t(k,i,j) = momz_t(k,i,j) + momz_t_mp(k,i,j)
       momx_t(k,i,j) = momx_t(k,i,j) + momx_t_mp(k,i,j)
       momy_t(k,i,j) = momy_t(k,i,j) + momy_t_mp(k,i,j)
       rhot_t(k,i,j) = rhot_t(k,i,j) + rhot_t_mp(k,i,j)
    enddo
    enddo
    enddo

    do iq = qs_mp, qe_mp
    !$omp parallel do default(none) private(i,j,k) OMP_SCHEDULE_ &
    !$omp shared(JS,JE,IS,IE,KS,KE,RHOQ_t,iq,RHOQ_t_MP)
    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_mp(k,i,j,iq)
    enddo
    enddo
    enddo
    enddo

    if ( statistics_checktotal ) then
       call stat_total( total, dens_t_mp(:,:,:), 'DENS_t_MP' )
       call stat_total( total, momz_t_mp(:,:,:), 'MOMZ_t_MP' )
       call stat_total( total, momx_t_mp(:,:,:), 'MOMX_t_MP' )
       call stat_total( total, momy_t_mp(:,:,:), 'MOMY_t_MP' )
       call stat_total( total, rhot_t_mp(:,:,:), 'RHOT_t_MP' )

       do iq = qs_mp, qe_mp
          call stat_total( total, rhoq_t_mp(:,:,:,iq), trim(tracer_name(iq))//'_t_MP' )
       enddo
    endif

    return

Here is the caller graph for this function: