scale_atmos_phy_mp_dry.F90

Go to the documentation of this file.

!-------------------------------------------------------------------------------
!-------------------------------------------------------------------------------
module scale_atmos_phy_mp_dry
  !-----------------------------------------------------------------------------
  !
  !++ used modules
  !
  use scale_precision
  use scale_stdio
  use scale_prof
  use scale_grid_index
  use scale_atmos_hydrometeor, only: &
     n_hyd
  !-----------------------------------------------------------------------------
  implicit none
  private
  !-----------------------------------------------------------------------------
  !
  !++ Public procedure
  !
  public :: atmos_phy_mp_dry_config
  public :: atmos_phy_mp_dry_setup
  public :: atmos_phy_mp_dry

  public :: atmos_phy_mp_dry_cloudfraction
  public :: atmos_phy_mp_dry_effectiveradius
  public :: atmos_phy_mp_dry_mixingratio

  !-----------------------------------------------------------------------------
  !
  !++ Public parameters & variables
  !
  integer, public, parameter :: qa_mp  = 0

  character(len=H_SHORT), public, target :: atmos_phy_mp_dry_name(qa_mp)
  character(len=H_MID)  , public, target :: atmos_phy_mp_dry_desc(qa_mp)
  character(len=H_SHORT), public, target :: atmos_phy_mp_dry_unit(qa_mp)

  real(RP), public, target :: atmos_phy_mp_dry_dens(n_hyd)

  !-----------------------------------------------------------------------------
  !
  !++ Private procedure
  !
  !-----------------------------------------------------------------------------
  !
  !++ Private parameters & variables
  !
  !-----------------------------------------------------------------------------
contains
  !-----------------------------------------------------------------------------
  subroutine atmos_phy_mp_dry_config( &
       MP_TYPE, &
       QA, QS   )
    use scale_process, only: &
       prc_mpistop
    implicit none

    character(len=*), intent(in)  :: mp_type
    integer,          intent(out) :: qa
    integer,          intent(out) :: qs
    !---------------------------------------------------------------------------

    if( io_l ) write(io_fid_log,*)
    if( io_l ) write(io_fid_log,*) '++++++ Module[Cloud Microphysics Tracer] / Categ[ATMOS PHYSICS] / Origin[SCALElib]'
    if( io_l ) write(io_fid_log,*) '*** No tracers for dummy process (dry Atmosphere)'

    if ( mp_type /= 'DRY' ) then
       write(*,*) 'xxx ATMOS_PHY_MP_TYPE is not DRY. Check!'
       call prc_mpistop
    endif

    qs = -1
    qa = qa_mp

    return
  end subroutine atmos_phy_mp_dry_config

  !-----------------------------------------------------------------------------
  subroutine atmos_phy_mp_dry_setup
    implicit none

    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,*) '*** dummy process (dry Atmosphere)'

    return
  end subroutine atmos_phy_mp_dry_setup

  !-----------------------------------------------------------------------------
  subroutine atmos_phy_mp_dry( &
       DENS,      &
       MOMZ,      &
       MOMX,      &
       MOMY,      &
       RHOT,      &
       QTRC,      &
       CCN ,      &
       EVAPORATE, &
       SFLX_rain, &
       SFLX_snow  )
    use scale_grid_index
    use scale_tracer, only: &
       qa
    use scale_const, only: &
       undef => const_undef
    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,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)
    !---------------------------------------------------------------------------

    if( io_l ) write(io_fid_log,*) '*** Atmos physics  step: Cloud microphysics(dummy)'

    atmos_phy_mp_dry_dens(:) = undef

    evaporate(:,:,:) = 0.0_rp
    sflx_rain(:,:) = 0.0_rp
    sflx_snow(:,:) = 0.0_rp

    return
  end subroutine atmos_phy_mp_dry

  !-----------------------------------------------------------------------------
  subroutine atmos_phy_mp_dry_cloudfraction( &
       cldfrac,       &
       QTRC,          &
       mask_criterion )
    use scale_grid_index
    use scale_tracer, only: &
       qa
    implicit none

    real(RP), intent(out) :: cldfrac(ka,ia,ja)
    real(RP), intent(in)  :: qtrc   (ka,ia,ja,qa)
    real(RP), intent(in)  :: mask_criterion ! not used
    !---------------------------------------------------------------------------

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

    return
  end subroutine atmos_phy_mp_dry_cloudfraction

  !-----------------------------------------------------------------------------
  subroutine atmos_phy_mp_dry_effectiveradius( &
       Re,    &
       QTRC0, &
       DENS0, &
       TEMP0  )
    use scale_grid_index
    use scale_tracer, only: &
       qa
    use scale_atmos_hydrometeor, only: &
       n_hyd
    implicit none

    real(RP), intent(out) :: re   (ka,ia,ja,n_hyd) ! 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]
    !---------------------------------------------------------------------------

    re = 0.0_rp

    return
  end subroutine atmos_phy_mp_dry_effectiveradius
  !-----------------------------------------------------------------------------
  subroutine atmos_phy_mp_dry_mixingratio( &
       Qe,   &
       QTRC0 )
    use scale_grid_index
    use scale_tracer, only: &
       qa
    use scale_atmos_hydrometeor, only: &
       n_hyd
    implicit none

    real(RP), intent(out) :: qe   (ka,ia,ja,n_hyd) ! mixing ratio of each cateory [kg/kg]
    real(RP), intent(in)  :: qtrc0(ka,ia,ja,qa)    ! tracer mass concentration [kg/kg]
    !---------------------------------------------------------------------------

    qe = 0.0_rp

    return
  end subroutine atmos_phy_mp_dry_mixingratio

end module scale_atmos_phy_mp_dry