scale_atmos_phy_ae.F90

Go to the documentation of this file.

!-------------------------------------------------------------------------------
!-------------------------------------------------------------------------------
module scale_atmos_phy_ae
  !-----------------------------------------------------------------------------
  !
  !++ 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_setup

  !-----------------------------------------------------------------------------
  !
  !++ Public parameters & variables
  !
  abstract interface
     subroutine ae( &
          DENS, &
          MOMZ, &
          MOMX, &
          MOMY, &
          RHOT, &
          EMIT, &
          NREG, &
          QTRC, &
          CN,   &
          CCN,  &
          RHOQ_t_AE )
       use scale_precision
       use scale_grid_index
       use scale_tracer
       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)
     end subroutine ae

     subroutine er( RE, QTRC, RH )
       use scale_precision
       use scale_grid_index
       use scale_tracer
       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]
     end subroutine er
  end interface

  procedure(ae), pointer :: atmos_phy_ae => null()
  procedure(er), pointer :: atmos_phy_ae_effectiveradius => null()
  public :: atmos_phy_ae
  public :: atmos_phy_ae_effectiveradius

  real(RP), public, pointer :: ae_dens(:) ! aerosol density [kg/m3]=[g/L]
  !-----------------------------------------------------------------------------
  !
  !++ Private procedure
  !
  !-----------------------------------------------------------------------------
  !
  !++ Private parameters & variables
  !
  !-----------------------------------------------------------------------------
contains
  !-----------------------------------------------------------------------------
  subroutine atmos_phy_ae_setup( AE_TYPE )
    use scale_process, only: &
       prc_mpistop
#define EXTM(pre, name, post) pre ## name ## post
#define NAME(pre, name, post) EXTM(pre, name, post)
#ifdef AE
    use name(scale_atmos_phy_ae_, ae,), only: &
       name(atmos_phy_ae_, ae, _setup), &
       name(atmos_phy_ae_, ae,)
       name(atmos_phy_ae_, ae, _effectiveradius)
#else
    use scale_atmos_phy_ae_dummy, only: &
       atmos_phy_ae_dummy_setup, &
       atmos_phy_ae_dummy, &
       atmos_phy_ae_dummy_effectiveradius
    use scale_atmos_phy_ae_kajino13, only: &
       atmos_phy_ae_kajino13_setup, &
       atmos_phy_ae_kajino13, &
       atmos_phy_ae_kajino13_effectiveradius
#endif
    implicit none

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

    select case( ae_type )
    case ( 'DUMMY', 'NONE' )
       call atmos_phy_ae_dummy_setup( ae_type )
       atmos_phy_ae                 => atmos_phy_ae_dummy
       atmos_phy_ae_effectiveradius => atmos_phy_ae_dummy_effectiveradius
    case ( 'KAJINO13' )
       call atmos_phy_ae_kajino13_setup( ae_type )
       atmos_phy_ae                 => atmos_phy_ae_kajino13
       atmos_phy_ae_effectiveradius => atmos_phy_ae_kajino13_effectiveradius
       write(*,*) '### aerosol type(', ae_type, '). is not recommended in current version!'
    case default
       write(*,*) 'xxx invalid aerosol type(', ae_type, '). CHECK!'
       call prc_mpistop
    end select

    return
  end subroutine atmos_phy_ae_setup

end module scale_atmos_phy_ae