scale_atmos_phy_mp_sdm.F90

Go to the documentation of this file.

!-------------------------------------------------------------------------------
!-------------------------------------------------------------------------------
module scale_atmos_phy_mp_sdm
  !-----------------------------------------------------------------------------
  !
  !++ used modules
  !
  use scale_precision
  use scale_stdio
  use scale_prof
  use scale_grid_index

  use scale_tracer_sdm
  !-----------------------------------------------------------------------------
  implicit none
  private
  !-----------------------------------------------------------------------------
  !
  !++ Public procedure
  !
  public :: atmos_phy_mp_sdm_setup
  public :: atmos_phy_mp_sdm

  public :: atmos_phy_mp_sdm_cloudfraction
  public :: atmos_phy_mp_sdm_effectiveradius
  public :: atmos_phy_mp_sdm_mixingratio

  !-----------------------------------------------------------------------------
  !
  !++ Public parameters & variables
  !
  real(RP), public, target :: atmos_phy_mp_dens(mp_qa) ! hydrometeor density [kg/m3]=[g/L]

  !-----------------------------------------------------------------------------
  !
  !++ Private procedure
  !
  !-----------------------------------------------------------------------------
  !
  !++ Private parameters & variables
  !
  !-----------------------------------------------------------------------------
contains
  !-----------------------------------------------------------------------------
  !-----------------------------------------------------------------------------
  subroutine atmos_phy_mp_sdm_setup( MP_TYPE )
    use scale_process, only: &
       prc_mpistop
    implicit none
    character(len=H_SHORT), intent(in) :: MP_TYPE
    !---------------------------------------------------------------------------

    if( io_l ) write(io_fid_log,*)
    if( io_l ) write(io_fid_log,*) '*** SDM not supported.'
    if( io_l ) write(io_fid_log,*) '*** Please contact SCALE developers'
    call prc_mpistop

    return
  end subroutine atmos_phy_mp_sdm_setup

  !-----------------------------------------------------------------------------
  !-----------------------------------------------------------------------------
  subroutine atmos_phy_mp_sdm( &
       DENS,      &
       MOMZ,      &
       MOMX,      &
       MOMY,      &
       RHOT,      &
       QTRC,      &
       CCN,       &
       EVAPORATE, &
       SFLX_rain, &
       SFLX_snow  )
    use scale_grid_index
    use scale_tracer, only: &
       qad => qa
    use scale_process, only: &
       prc_mpistop
    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)   !---- evaporated cloud number concentration [/m3]
    real(RP), intent(out)   :: SFLX_rain(ia,ja)
    real(RP), intent(out)   :: SFLX_snow(ia,ja)
    !---------------------------------------------------------------------------

    if( io_l ) write(io_fid_log,*)
    if( io_l ) write(io_fid_log,*) '*** SDM not supported.'
    if( io_l ) write(io_fid_log,*) '*** Please contact SCALE developers'
    call prc_mpistop

    return
  end subroutine atmos_phy_mp_sdm

  !-----------------------------------------------------------------------------
  subroutine atmos_phy_mp_sdm_cloudfraction( &
       cldfrac, &
       QTRC     )
    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)
    !---------------------------------------------------------------------------

    cldfrac(:,:,:) = 0.0_rp ! dummy

    return
  end subroutine atmos_phy_mp_sdm_cloudfraction

  !-----------------------------------------------------------------------------
  subroutine atmos_phy_mp_sdm_effectiveradius( &
       Re,    &
       QTRC0, &
       DENS0, &
       TEMP0  )
    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
    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)        ! Temperatuer [K]
    !---------------------------------------------------------------------------

    re(:,:,:,:) = 8.e-6_rp ! dummy

    return
  end subroutine atmos_phy_mp_sdm_effectiveradius
  !-----------------------------------------------------------------------------
  subroutine atmos_phy_mp_sdm_mixingratio( &
       Qe,    &
       QTRC0  )
    use scale_const, only: &
       eps => const_eps
    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]

    integer  :: ihydro
    !---------------------------------------------------------------------------

    do ihydro = 1, mp_qa
       qe(:,:,:,ihydro) = 8.e-6_rp ! dummy
    enddo

    return

  end subroutine atmos_phy_mp_sdm_mixingratio
  !-----------------------------------------------------------------------------

end module scale_atmos_phy_mp_sdm