scale_atmos_phy_mp.F90

Go to the documentation of this file.

!-------------------------------------------------------------------------------
! Warning: This file was generated from atmos-physics/microphysics/scale_atmos_phy_mp.F90.erb.
!          Do not edit this file.
!-------------------------------------------------------------------------------
module scale_atmos_phy_mp
  !-----------------------------------------------------------------------------
  !
  !++ used modules
  !
  use scale_precision
  use scale_stdio
  use scale_prof
  use scale_grid_index
  use scale_tracer
  !-----------------------------------------------------------------------------
  implicit none
  private
  !-----------------------------------------------------------------------------
  !
  !++ Public procedure
  !
  public :: atmos_phy_mp_setup

  !-----------------------------------------------------------------------------
  !
  !++ Public parameters & variables
  !
  !-----------------------------------------------------------------------------
  !
  !++ Private procedure
  !
  !-----------------------------------------------------------------------------
  !
  !++ Private parameters & variables
  !
  abstract interface
     subroutine mp( &
          DENS,      &
          MOMZ,      &
          MOMX,      &
          MOMY,      &
          RHOT,      &
          QTRC,      &
          CCN,       &
          EVAPORATE, &
          SFLX_rain, &
          SFLX_snow  )
       use scale_precision
       use scale_grid_index
       use scale_tracer
       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,qa)
       real(RP), intent(in)    :: ccn(ka,ia,ja)
       real(RP), intent(out)   :: evaporate(ka,ia,ja)
       real(RP), intent(out)   :: sflx_rain(ia,ja)
       real(RP), intent(out)   :: sflx_snow(ia,ja)
     end subroutine mp
     subroutine cf( &
          cldfrac,       &
          QTRC,          &
          mask_criterion )
       use scale_precision
       use scale_grid_index
       use scale_tracer
       real(RP), intent(out) :: cldfrac(ka,ia,ja)
       real(RP), intent(in)  :: qtrc   (ka,ia,ja,qa)
       real(RP), intent(in)  :: mask_criterion
     end subroutine cf
     subroutine er( &
          Re,    &
          QTRC0, &
          DENS0, &
          TEMP0  )
       use scale_precision
       use scale_grid_index
       use scale_tracer
       real(RP), intent(out) :: re   (ka,ia,ja,mp_qa) ! effective radius          [cm]
       real(RP), intent(in)  :: qtrc0(ka,ia,ja,qa)    ! 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]
     end subroutine er
     subroutine mr( &
          Qe,    &
          QTRC0  )
       use scale_precision
       use scale_grid_index
       use scale_tracer
       real(RP), intent(out) :: qe   (ka,ia,ja,mp_qa) ! mixing ratio of each cateory [kg/kg]
       real(RP), intent(in)  :: qtrc0(ka,ia,ja,qa)    ! tracer mass concentration [kg/kg]
     end subroutine mr
  end interface
  procedure(mp), pointer :: atmos_phy_mp => null()
  procedure(cf), pointer :: atmos_phy_mp_cloudfraction => null()
  procedure(er), pointer :: atmos_phy_mp_effectiveradius => null()
  procedure(mr), pointer :: atmos_phy_mp_mixingratio => null()
  public :: atmos_phy_mp
  public :: atmos_phy_mp_cloudfraction
  public :: atmos_phy_mp_effectiveradius
  public :: atmos_phy_mp_mixingratio

  real(RP), pointer, public :: atmos_phy_mp_dens(:)

  !-----------------------------------------------------------------------------
contains
  !-----------------------------------------------------------------------------
  !-----------------------------------------------------------------------------
  subroutine atmos_phy_mp_setup( MP_TYPE )
    use scale_process, only: &
       prc_mpistop
#define EXTM(pre, name, post) pre ## name ## post
#define NAME(pre, name, post) EXTM(pre, name, post)
#ifdef MP
    use name(scale_atmos_phy_mp_, mp,), only: &
       name(atmos_phy_mp_, mp, _setup), &
       name(atmos_phy_mp_, mp,), &
       name(atmos_phy_mp_, mp, _cloudfraction), &
       name(atmos_phy_mp_, mp, _effectiveradius), &
       name(atmos_phy_mp_, mp, _mixingratio), &
       name(atmos_phy_mp_, mp, _dens) => atmos_phy_mp_dens
#else
    use scale_atmos_phy_mp_dry, only: &
       atmos_phy_mp_dry_setup, &
       atmos_phy_mp_dry, &
       atmos_phy_mp_dry_cloudfraction, &
       atmos_phy_mp_dry_effectiveradius, &
       atmos_phy_mp_dry_mixingratio, &
       atmos_phy_mp_dry_dens => atmos_phy_mp_dens
    use scale_atmos_phy_mp_kessler, only: &
       atmos_phy_mp_kessler_setup, &
       atmos_phy_mp_kessler, &
       atmos_phy_mp_kessler_cloudfraction, &
       atmos_phy_mp_kessler_effectiveradius, &
       atmos_phy_mp_kessler_mixingratio, &
       atmos_phy_mp_kessler_dens => atmos_phy_mp_dens
    use scale_atmos_phy_mp_tomita08, only: &
       atmos_phy_mp_tomita08_setup, &
       atmos_phy_mp_tomita08, &
       atmos_phy_mp_tomita08_cloudfraction, &
       atmos_phy_mp_tomita08_effectiveradius, &
       atmos_phy_mp_tomita08_mixingratio, &
       atmos_phy_mp_tomita08_dens => atmos_phy_mp_dens
    use scale_atmos_phy_mp_sn14, only: &
       atmos_phy_mp_sn14_setup, &
       atmos_phy_mp_sn14, &
       atmos_phy_mp_sn14_cloudfraction, &
       atmos_phy_mp_sn14_effectiveradius, &
       atmos_phy_mp_sn14_mixingratio, &
       atmos_phy_mp_sn14_dens => atmos_phy_mp_dens
    use scale_atmos_phy_mp_suzuki10, only: &
       atmos_phy_mp_suzuki10_setup, &
       atmos_phy_mp_suzuki10, &
       atmos_phy_mp_suzuki10_cloudfraction, &
       atmos_phy_mp_suzuki10_effectiveradius, &
       atmos_phy_mp_suzuki10_mixingratio, &
       atmos_phy_mp_suzuki10_dens => atmos_phy_mp_dens
    use scale_atmos_phy_mp_sdm, only: &
       atmos_phy_mp_sdm_setup, &
       atmos_phy_mp_sdm, &
       atmos_phy_mp_sdm_cloudfraction, &
       atmos_phy_mp_sdm_effectiveradius, &
       atmos_phy_mp_sdm_mixingratio, &
       atmos_phy_mp_sdm_dens => atmos_phy_mp_dens
#endif
    implicit none

    character(len=*), intent(in) :: MP_TYPE
    !---------------------------------------------------------------------------

#ifdef MP
    call name(atmos_phy_mp_, mp, _setup)( mp_type )
    atmos_phy_mp                 => name(atmos_phy_mp_, mp,)
    atmos_phy_mp_cloudfraction   => name(atmos_phy_mp_, mp, _cloudfraction)
    atmos_phy_mp_effectiveradius => name(atmos_phy_mp_, mp, _effectiveradius)
    atmos_phy_mp_mixingratio     => name(atmos_phy_mp_, mp, _mixingratio)
    atmos_phy_mp_dens            => name(atmos_phy_mp_, mp, _dens)
#else
    select case ( mp_type )
    case ( 'DRY' )
       call atmos_phy_mp_dry_setup( mp_type )
       atmos_phy_mp                 => atmos_phy_mp_dry
       atmos_phy_mp_cloudfraction   => atmos_phy_mp_dry_cloudfraction
       atmos_phy_mp_effectiveradius => atmos_phy_mp_dry_effectiveradius
       atmos_phy_mp_mixingratio     => atmos_phy_mp_dry_mixingratio
       atmos_phy_mp_dens            => atmos_phy_mp_dry_dens
    case ( 'KESSLER' )
       call atmos_phy_mp_kessler_setup( mp_type )
       atmos_phy_mp                 => atmos_phy_mp_kessler
       atmos_phy_mp_cloudfraction   => atmos_phy_mp_kessler_cloudfraction
       atmos_phy_mp_effectiveradius => atmos_phy_mp_kessler_effectiveradius
       atmos_phy_mp_mixingratio     => atmos_phy_mp_kessler_mixingratio
       atmos_phy_mp_dens            => atmos_phy_mp_kessler_dens
    case ( 'TOMITA08' )
       call atmos_phy_mp_tomita08_setup( mp_type )
       atmos_phy_mp                 => atmos_phy_mp_tomita08
       atmos_phy_mp_cloudfraction   => atmos_phy_mp_tomita08_cloudfraction
       atmos_phy_mp_effectiveradius => atmos_phy_mp_tomita08_effectiveradius
       atmos_phy_mp_mixingratio     => atmos_phy_mp_tomita08_mixingratio
       atmos_phy_mp_dens            => atmos_phy_mp_tomita08_dens
    case ( 'SN14' )
       call atmos_phy_mp_sn14_setup( mp_type )
       atmos_phy_mp                 => atmos_phy_mp_sn14
       atmos_phy_mp_cloudfraction   => atmos_phy_mp_sn14_cloudfraction
       atmos_phy_mp_effectiveradius => atmos_phy_mp_sn14_effectiveradius
       atmos_phy_mp_mixingratio     => atmos_phy_mp_sn14_mixingratio
       atmos_phy_mp_dens            => atmos_phy_mp_sn14_dens
    case ( 'SUZUKI10' )
       call atmos_phy_mp_suzuki10_setup( mp_type )
       atmos_phy_mp                 => atmos_phy_mp_suzuki10
       atmos_phy_mp_cloudfraction   => atmos_phy_mp_suzuki10_cloudfraction
       atmos_phy_mp_effectiveradius => atmos_phy_mp_suzuki10_effectiveradius
       atmos_phy_mp_mixingratio     => atmos_phy_mp_suzuki10_mixingratio
       atmos_phy_mp_dens            => atmos_phy_mp_suzuki10_dens
    case ( 'SDM' )
       call atmos_phy_mp_sdm_setup( mp_type )
       atmos_phy_mp                 => atmos_phy_mp_sdm
       atmos_phy_mp_cloudfraction   => atmos_phy_mp_sdm_cloudfraction
       atmos_phy_mp_effectiveradius => atmos_phy_mp_sdm_effectiveradius
       atmos_phy_mp_mixingratio     => atmos_phy_mp_sdm_mixingratio
       atmos_phy_mp_dens            => atmos_phy_mp_sdm_dens
    end select
#endif

    return
  end subroutine atmos_phy_mp_setup

end module scale_atmos_phy_mp