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_tracer_dry
  !-----------------------------------------------------------------------------
  implicit none
  private
  !-----------------------------------------------------------------------------
  !
  !++ Public procedure
  !
  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
  !
  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_dry_setup( MP_TYPE )
    use scale_process, only: &
       prc_mpistop
    implicit none

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

    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)'

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

    atmos_phy_mp_dens(:) = 0.0_rp

    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: &
       qad => qa
    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) 
    real(RP), intent(out)   :: SFLX_rain(ia,ja)
    real(RP), intent(out)   :: SFLX_snow(ia,ja)
    !---------------------------------------------------------------------------

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

    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     )
    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_dry_cloudfraction

  !-----------------------------------------------------------------------------
  subroutine atmos_phy_mp_dry_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          [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(:,:,:,:) = 8.e-6_rp * um2cm ! dummy

    return
  end subroutine atmos_phy_mp_dry_effectiveradius
  !-----------------------------------------------------------------------------
  subroutine atmos_phy_mp_dry_mixingratio( &
       Qe,    &
       QTRC0  )
    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_dry_mixingratio
  !-----------------------------------------------------------------------------

end module scale_atmos_phy_mp_dry