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_setup
  public :: atmos_phy_ae_dummy

  public :: atmos_phy_ae_dummy_effectiveradius

  !-----------------------------------------------------------------------------
  !
  !++ Public parameters & variables
  !
  real(RP), public, allocatable :: ae_dens(:) ! aerosol density [kg/m3]=[g/L]

  !-----------------------------------------------------------------------------
  !
  !++ Private procedure
  !
  !-----------------------------------------------------------------------------
  !
  !++ Private parameters & variables
  !
  !-----------------------------------------------------------------------------
contains
  !-----------------------------------------------------------------------------
  subroutine atmos_phy_ae_dummy_setup( AE_TYPE )
    use scale_process, only: &
       prc_mpistop
    implicit none

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

    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 aerosol process'
    if( io_l ) write(io_fid_log,*) '*** No namelists.'

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

    allocate( ae_dens(ae_qa) )
    ae_dens(:) = 0.0_rp

    return
  end subroutine atmos_phy_ae_dummy_setup

  !-----------------------------------------------------------------------------
  subroutine atmos_phy_ae_dummy( &
          DENS, &
          MOMZ, &
          MOMX, &
          MOMY, &
          RHOT, &
          EMIT, &
          NREG, &
          QTRC, &
          CN,   &
          CCN,  &
          RHOQ_t_AE )
    use scale_grid_index
    use scale_tracer
    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) :: EMIT(ka,ia,ja,qa_ae)
    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,*) '*** 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
    implicit none

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

    re(:,:,:,:) = undef

!    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