scale_atmos_phy_ae_dummy.F90

Go to the documentation of this file.

!-------------------------------------------------------------------------------
!-------------------------------------------------------------------------------
module scale_atmos_phy_ae_dummy
  !-----------------------------------------------------------------------------
  !
  !++ 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_ae_dummy_config
  public :: atmos_phy_ae_dummy_setup
  public :: atmos_phy_ae_dummy

  public :: atmos_phy_ae_dummy_effectiveradius

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

  character(len=H_SHORT), public, target :: atmos_phy_ae_dummy_name(qa_ae)
  character(len=H_MID)  , public, target :: atmos_phy_ae_dummy_desc(qa_ae)
  character(len=H_SHORT), public, target :: atmos_phy_ae_dummy_unit(qa_ae)

  real(RP), public, target :: atmos_phy_ae_dummy_dens(qa_ae) ! hydrometeor density [kg/m3]=[g/L]

  !-----------------------------------------------------------------------------
  !
  !++ Private procedure
  !
  !-----------------------------------------------------------------------------
  !
  !++ Private parameters & variables
  !
  !-----------------------------------------------------------------------------
contains
  !-----------------------------------------------------------------------------
  subroutine atmos_phy_ae_dummy_config( &
       AE_TYPE, &
       QA_AE, QS_AE )
    use scale_process, only: &
       prc_mpistop
    implicit none
    character(len=*), intent(in)  :: ae_type
    integer,          intent(out) :: qa_ae
    integer,          intent(out) :: qs_ae
    !---------------------------------------------------------------------------

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

    if ( ae_type /= 'DUMMY' .AND. ae_type /= 'NONE' ) then
       write(*,*) 'xxx ATMOS_PHY_AE_TYPE is not DUMMY. Check!'
       call prc_mpistop
    endif

    qa_ae = 0
    qs_ae = -1

    return
  end subroutine atmos_phy_ae_dummy_config

  !-----------------------------------------------------------------------------
  subroutine atmos_phy_ae_dummy_setup
    implicit none

    if( io_l ) write(io_fid_log,*)
    if( io_l ) write(io_fid_log,*) '++++++ Module[Aerosol] / Categ[ATMOS PHYSICS] / Origin[SCALElib]'
    if( io_l ) write(io_fid_log,*) '*** dummy process'

    return
  end subroutine atmos_phy_ae_dummy_setup

  !-----------------------------------------------------------------------------
  subroutine atmos_phy_ae_dummy( &
       QQA,  &
       DENS, &
       MOMZ, &
       MOMX, &
       MOMY, &
       RHOT, &
       EMIT, &
       NREG, &
       QTRC, &
       CN,   &
       CCN,  &
       RHOQ_t_AE )
    use scale_grid_index
    use scale_tracer
    implicit none
    integer,  intent(in)    :: qqa
    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) :: emit(ka,ia,ja,qqa)
    real(RP), intent(in)    :: nreg(ka,ia,ja)
    real(RP), intent(inout) :: qtrc(ka,ia,ja,qa)
    real(RP), intent(out)   :: cn(ka,ia,ja)
    real(RP), intent(out)   :: ccn(ka,ia,ja)
    real(RP), intent(inout) :: rhoq_t_ae(ka,ia,ja,qa)

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

    cn(:,:,:) = 0.0_rp
    ccn(:,:,:) = 0.0_rp

    return
  end subroutine atmos_phy_ae_dummy

  !-----------------------------------------------------------------------------
  subroutine atmos_phy_ae_dummy_effectiveradius( &
       Re,   &
       QTRC, &
       RH    )
    use scale_grid_index
    use scale_tracer
    use scale_const, only: &
       undef => const_undef
    use scale_atmos_aerosol, only: &
       n_ae
    implicit none
    real(RP), intent(out) :: re  (ka,ia,ja,n_ae) ! effective radius
    real(RP), intent(in)  :: qtrc(ka,ia,ja,qa)   ! tracer mass concentration [kg/kg]
    real(RP), intent(in)  :: rh  (ka,ia,ja)      ! relative humidity         (0-1)
    integer :: iq
    !---------------------------------------------------------------------------

    do iq = 1, n_ae
       re(:,:,:,iq) = 0.0_rp
    end do

!    Re(:,:,:,I_ae_seasalt) = 2.E-4_RP
!    Re(:,:,:,I_ae_dust   ) = 4.E-6_RP
!    Re(:,:,:,I_ae_bc     ) = 4.E-8_RP
!    Re(:,:,:,I_ae_oc     ) = RH(:,:,:)
!    Re(:,:,:,I_ae_sulfate) = RH(:,:,:)

    return
  end subroutine atmos_phy_ae_dummy_effectiveradius

end module scale_atmos_phy_ae_dummy