scale_atmos_phy_ae_kajino13 Module Reference

module atmosphere / physics / aerosol / Kajino13 More...

Functions/Subroutines
subroutine, public	atmos_phy_ae_kajino13_tracer_setup (QA_AE)
	Tracer setup. More...
subroutine, public	atmos_phy_ae_kajino13_setup
	Setup. More...
subroutine, public	atmos_phy_ae_kajino13_tendency (KA, KS, KE, IA, IS, IE, JA, JS, JE, QA_AE, TEMP, PRES, QDRY, NREG, DENS, QV, QTRC, EMIT, dt, RHOQ_t_AE, CN, CCN)
	Aerosol Microphysics. More...
subroutine, public	atmos_phy_ae_kajino13_effective_radius (KA, IA, JA, QA_AE, QTRC, RH, Re)
	Calculate Effective Radius. More...
subroutine, public	atmos_phy_ae_kajino13_mkinit (KA, KS, KE, IA, IS, IE, JA, JS, JE, QA_AE, DENS, TEMP, PRES, QDRY, QV, m0_init, dg_init, sg_init, d_min_inp, d_max_inp, k_min_inp, k_max_inp, n_kap_inp, QTRC, CCN)
subroutine, public	atmos_phy_ae_kajino13_negative_fixer (KA, KS, KE, IA, IS, IE, JA, JS, JE, QA_AE, QTRC)

Variables
character(len=h_short), dimension(:), allocatable, public	atmos_phy_ae_kajino13_name
character(len=h_mid), dimension(:), allocatable, public	atmos_phy_ae_kajino13_desc
character(len=h_short), dimension(:), allocatable, public	atmos_phy_ae_kajino13_unit
real(rp), public	atmos_phy_ae_kajino13_h2so4dt = 5.E-6_RP
real(rp), public	atmos_phy_ae_kajino13_ocgasdt = 8.E-5_RP
real(rp), public	atmos_phy_ae_kajino13_c_kappa = 0.3_RP
logical, public	atmos_phy_ae_kajino13_flag_npf = .false.
logical, public	atmos_phy_ae_kajino13_flag_cond = .true.
logical, public	atmos_phy_ae_kajino13_flag_coag = .true.
logical, public	atmos_phy_ae_kajino13_flag_ccn_interactive = .true.
logical, public	atmos_phy_ae_kajino13_flag_regeneration = .true.
real(rp), public	atmos_phy_ae_kajino13_dg_reg = 5.E-7_RP
real(rp), public	atmos_phy_ae_kajino13_sg_reg = 1.6_RP
real(rp), public	atmos_phy_ae_kajino13_logk_aenucl = -12.4_RP

Detailed Description

module atmosphere / physics / aerosol / Kajino13

Description

kajino13 aerosol microphysics scheme

Author

Team SCALE

NAMELIST

• PARAM_ATMOS_PHY_AE_KAJINO13_TRACER

  name	type	default value	comment
  AE_CTG	integer	1
  NASIZ	integer, dimension(3)
  NAKAP	integer, dimension(3)

  
  
• PARAM_ATMOS_PHY_AE_KAJINO13

  name	type	default value	comment
  ATMOS_PHY_AE_KAJINO13_H2SO4DT	real(RP)	5.E-6_RP	h2so4 production rate (Temporal) [ug/m3/s]
  ATMOS_PHY_AE_KAJINO13_OCGASDT	real(RP)	8.E-5_RP	other condensational bas production rate 16*h2so4dt (see Kajino et al. 2013)

  
  

History Output

name	description	unit	variable
CGAS_emit	Emission ratio of Condensabule gas	ug/m3/s	EMIT
H2SO4_emit	Emission ratio of H2SO4 gas	ug/m3/s	EMIT
{ctg_name}_mass	Total mass mixing ratio of aerosol	kg/kg	total_aerosol_mass
{ctg_name}_mass_emit	Total mass mixing ratio of emitted aerosol	kg/kg	total_emit_aerosol_mass
{ctg_name}_number	Total number mixing ratio of aerosol	num/kg	total_aerosol_number
{ctg_name}_number_emit	Total number mixing ratio of emitted aerosol	num/kg	total_emit_aerosol_number

Function/Subroutine Documentation

atmos_phy_ae_kajino13_tracer_setup()

subroutine, public scale_atmos_phy_ae_kajino13::atmos_phy_ae_kajino13_tracer_setup

(

integer, intent(out)

QA_AE

)

Tracer setup.

Definition at line 183 of file scale_atmos_phy_ae_kajino13.F90.

References atmos_phy_ae_kajino13_desc, atmos_phy_ae_kajino13_name, atmos_phy_ae_kajino13_unit, scale_io::io_fid_conf, and scale_prc::prc_abort().

Referenced by mod_atmos_phy_ae_driver::atmos_phy_ae_driver_tracer_setup().

    use scale_prc, only: &
       prc_abort

    integer, intent(out) :: qa_ae

    integer, allocatable :: aero_idx(:,:,:,:)
    integer :: ncat_max
    integer :: nasiz(3), nakap(3)
    character(len=H_SHORT) :: attribute, catego, aunit

    namelist / param_atmos_phy_ae_kajino13_tracer / &
       ae_ctg, &
       nasiz,  &
       nakap

    integer :: m, ierr, ik, ic, ia0, is0

    !---------------------------------------------------------------------------

    log_newline
    log_info("ATMOS_PHY_AE_kajino13_tracer_setup",*) 'Setup'

    ! note: tentatively, aerosol module should be called at all time. we need dummy subprogram.
    !if ( ATMOS_sw_phy_ae ) then

    log_warn("ATMOS_PHY_AE_kajino13_tracer_setup",*) '### Kajino13 scheme is still experimental'

    ncat_max = max( ic_mix, ic_sea, ic_dus )
   
    nasiz(:) = 64
    nakap(:) = 1
   
    rewind(io_fid_conf)
    read(io_fid_conf,nml=param_atmos_phy_ae_kajino13_tracer,iostat=ierr)

    if( ierr < 0 ) then !--- missing
       log_info("ATMOS_PHY_AE_kajino13_tracer_setup",*)  'Not found namelist. Default used.'
    elseif( ierr > 0 ) then !--- fatal error
       log_error("ATMOS_PHY_AE_kajino13_tracer_setup",*) 'Not appropriate names in namelist PARAM_ATMOS_PHY_AE_KAJINO13_TRACER, Check!'
       call prc_abort
    end if
   
    log_nml(param_atmos_phy_ae_kajino13_tracer)
   
    if( ae_ctg > ncat_max ) then
       log_error("ATMOS_PHY_AE_kajino13_tracer_setup",*) 'AE_CTG should be smaller than', ncat_max+1, 'stop'
       call prc_abort
    endif
   
    allocate( nsiz(ae_ctg) )
    allocate( nkap(ae_ctg) )
   
    nkap(1:ae_ctg) = nakap(1:ae_ctg)
    nsiz(1:ae_ctg) = nasiz(1:ae_ctg)

    if( maxval( nkap ) /= 1 .OR. minval( nkap ) /= 1 ) then
       log_error("ATMOS_PHY_AE_kajino13_tracer_setup",*) 'NKAP(:) /= 1 is not supported now, Stop!'
       call prc_abort
    end if
   
    ! do ia0 = 1, N_ATR
    do ic = 1, ae_ctg
    do ik = 1, nkap(ic)
    do is0 = 1, nsiz(ic)
       qa_ae   = qa_ae + n_atr
    enddo
    enddo
    enddo
    qa_ae = qa_ae + gas_ctg
   
    allocate( atmos_phy_ae_kajino13_name(qa_ae) )
    allocate( atmos_phy_ae_kajino13_desc(qa_ae) )
    allocate( atmos_phy_ae_kajino13_unit(qa_ae) )
   
    n_siz_max = 0
    n_kap_max = 0
    do ic = 1, ae_ctg
       n_siz_max = max(n_siz_max, nsiz(ic))
       n_kap_max = max(n_kap_max, nkap(ic))
    enddo
   
    allocate( aero_idx(n_atr,ae_ctg,n_kap_max,n_siz_max) )
    m = 0
   !       do ia0 = 1, N_ATR
    do ic = 1, ae_ctg
    do ik = 1, nkap(ic)
    do is0 = 1, nsiz(ic)
    do ia0 = 1, n_atr
       m = m+1
       aero_idx(ia0,ic,ik,is0) = m
    enddo
    enddo
    enddo
    enddo
   
    !-----------------------------------------------------------------------------
    !
    !++ calculate each category and aerosol
    !
    !-----------------------------------------------------------------------------
    ic = qa_ae-gas_ctg+ig_h2so4
    write(atmos_phy_ae_kajino13_unit(ic),'(a)')  'kg/kg'
    ic = qa_ae-gas_ctg+ig_cgas
    write(atmos_phy_ae_kajino13_unit(ic),'(a)')  'kg/kg'
   
   !do ia0 = 1, N_ATR
   !  if( ia0 == 1 ) then
   !     write(attribute,'(a)') "Number"
   !     write(aunit,'(a)') "num/kg"
   !  elseif( ia0 == 2 ) then
   !     write(attribute,'(a)') "Section"
   !     write(aunit,'(a)') "m2/kg"
   !  elseif( ia0 == 3 ) then
   !     write(attribute,'(a)') "Volume"
   !     write(aunit,'(a)') "m3/kg"
   !  elseif( ia0 == 4 ) then
   !     write(attribute,'(a)') "Mass"
   !     write(aunit,'(a)') "kg/kg"
   !  elseif( ia0 == 5 ) then
   !     write(attribute,'(a)') "kpXmass"
   !     write(aunit,'(a)') "kg/kg"
   !  endif
    do ic = 1, ae_ctg       !aerosol category
    do ik = 1, nkap(ic)   !kappa bin
    do is0 = 1, nsiz(ic)
    do ia0 = 1, n_atr
       if( ia0 == 1 ) then
          write(attribute,'(a)') "Number"
          write(aunit,'(a)') "num/kg"
       elseif( ia0 == 2 ) then
          write(attribute,'(a)') "Section"
          write(aunit,'(a)') "m2/kg"
       elseif( ia0 == 3 ) then
          write(attribute,'(a)') "Volume"
          write(aunit,'(a)') "m3/kg"
       elseif( ia0 == 4 ) then
          write(attribute,'(a)') "Mass"
          write(aunit,'(a)') "kg/kg"
       elseif( ia0 == 5 ) then
          write(attribute,'(a)') "kXm"
          write(aunit,'(a)') "kg/kg"
       endif
       if( ic == ic_mix ) then
          write(catego,'(a)') "Sulf_"
       elseif( ic == ic_sea ) then
          write(catego,'(a)') "Salt_"
       elseif( ic == ic_dus ) then
          write(catego,'(a)') "Dust_"
       endif
       write(atmos_phy_ae_kajino13_unit(aero_idx(ia0,ic,ik,is0)),'(a)')  trim(aunit)
       write(atmos_phy_ae_kajino13_name(aero_idx(ia0,ic,ik,is0)),'(a,a,i0)') trim(catego), trim(attribute), is0
       write(atmos_phy_ae_kajino13_desc(aero_idx(ia0,ic,ik,is0)),'(a,a,a,i0)') trim(attribute), ' mixing radio of ', trim(catego), is0
    enddo
    enddo
    enddo
    enddo
    ic = qa_ae-gas_ctg+ig_h2so4
    write(atmos_phy_ae_kajino13_name(ic),'(a)') 'H2SO4_Gas'
    ic = qa_ae-gas_ctg+ig_cgas
    write(atmos_phy_ae_kajino13_name(ic),'(a)') 'Condensable_GAS'

    ic = qa_ae-gas_ctg+ig_h2so4
    write(atmos_phy_ae_kajino13_desc(ic),'(a)') 'Mixing ratio of H2SO4 Gas'
    ic = qa_ae-gas_ctg+ig_cgas
    write(atmos_phy_ae_kajino13_desc(ic),'(a)') 'Mixing ratio of Condensable GAS'
   
    deallocate(aero_idx)

    return

Here is the call graph for this function:

Here is the caller graph for this function:

atmos_phy_ae_kajino13_setup()

subroutine, public scale_atmos_phy_ae_kajino13::atmos_phy_ae_kajino13_setup

(

)

Setup.

Definition at line 357 of file scale_atmos_phy_ae_kajino13.F90.

References atmos_phy_ae_kajino13_c_kappa, atmos_phy_ae_kajino13_dg_reg, atmos_phy_ae_kajino13_flag_ccn_interactive, atmos_phy_ae_kajino13_flag_coag, atmos_phy_ae_kajino13_flag_cond, atmos_phy_ae_kajino13_flag_npf, atmos_phy_ae_kajino13_flag_regeneration, atmos_phy_ae_kajino13_h2so4dt, atmos_phy_ae_kajino13_logk_aenucl, atmos_phy_ae_kajino13_ocgasdt, atmos_phy_ae_kajino13_sg_reg, float(), scale_io::io_fid_conf, and scale_prc::prc_abort().

Referenced by mod_atmos_phy_ae_driver::atmos_phy_ae_driver_setup(), and atmos_phy_ae_kajino13_mkinit().

    use scale_prc, only: &
       prc_abort
    implicit none

    real(RP), allocatable :: d_min_inp(:)
    real(RP), allocatable :: d_max_inp(:)
    real(RP), allocatable :: k_min_inp(:)
    real(RP), allocatable :: k_max_inp(:)
    integer , allocatable :: n_kap_inp(:)

    real(RP), parameter :: d_min_def = 1.e-9_rp ! default lower bound of 1st size bin
    real(RP), parameter :: d_max_def = 1.e-5_rp ! upper bound of last size bin
    integer , parameter :: n_kap_def = 1        ! number of kappa bins
    real(RP), parameter :: k_min_def = 0.e0_rp  ! lower bound of 1st kappa bin
    real(RP), parameter :: k_max_def = 1.e0_rp  ! upper bound of last kappa bin

    namelist / param_atmos_phy_ae_kajino13 / &
       atmos_phy_ae_kajino13_h2so4dt, &
       atmos_phy_ae_kajino13_ocgasdt, &
!       ATMOS_PHY_AE_KAJINO13_c_ratio, &
       atmos_phy_ae_kajino13_c_kappa, &
!       d_min_inp, &
!       d_max_inp, &
!       k_min_inp, &
!       k_max_inp, &
!       n_kap_inp, &
       atmos_phy_ae_kajino13_flag_npf, &
       atmos_phy_ae_kajino13_flag_cond, &
       atmos_phy_ae_kajino13_flag_coag, &
       atmos_phy_ae_kajino13_flag_ccn_interactive, &
       atmos_phy_ae_kajino13_flag_regeneration,    &
       atmos_phy_ae_kajino13_dg_reg, &
       atmos_phy_ae_kajino13_sg_reg, &
       atmos_phy_ae_kajino13_logk_aenucl, &
       atmos_phy_ae_kajino13_nbins_out

    integer :: it, ierr
    !---------------------------------------------------------------------------

    log_newline
    log_info("ATMOS_PHY_AE_kajino13_setup",*) 'Setup'
    log_info("ATMOS_PHY_AE_kajino13_setup",*) 'Kajino(2013) scheme'

    !--- setup parameter
    pi6   = pi / 6._rp              ! pi/6
    sixpi = 6._rp / pi              ! 6/pi
    forpi = 4._rp / pi              ! 4/pi


    n_ctg = ae_ctg
    allocate( rnum_out(atmos_phy_ae_kajino13_nbins_out) )
    allocate( n_siz(n_ctg) )
    allocate( d_min(n_ctg) )
    allocate( d_max(n_ctg) )
    allocate( n_kap(n_ctg) )
    allocate( k_min(n_ctg) )
    allocate( k_max(n_ctg) )
    allocate( d_min_inp(n_ctg) )
    allocate( d_max_inp(n_ctg) )
    allocate( n_kap_inp(n_ctg) )
    allocate( k_min_inp(n_ctg) )
    allocate( k_max_inp(n_ctg) )
    allocate( ctg_name(n_ctg) )

    n_siz(1:n_ctg) = nsiz(1:n_ctg)       ! number of size bins
    d_min(1:n_ctg) = d_min_def ! lower bound of 1st size bin
    d_max(1:n_ctg) = d_max_def ! upper bound of last size bin
    n_kap(1:n_ctg) = n_kap_def ! number of kappa bins
    k_min(1:n_ctg) = k_min_def ! lower bound of 1st kappa bin
    k_max(1:n_ctg) = k_max_def ! upper bound of last kappa bin

    do it = 1, n_ctg
     if( n_ctg == 1 ) then
       write(ctg_name(it),'(a)') "Sulfate"
     elseif( n_ctg == 2 ) then
       write(ctg_name(it),'(a)') "Seasalt"
     elseif( n_ctg == 3 ) then
       write(ctg_name(it),'(a)') "Dust"
     endif
    enddo

    !--- read namelist
    rewind(io_fid_conf)
    read(io_fid_conf,nml=param_atmos_phy_ae_kajino13,iostat=ierr)
    if( ierr < 0 ) then !--- missing
       log_info("ATMOS_PHY_AE_kajino13_setup",*) 'Not found namelist. Default used.'
    elseif( ierr > 0 ) then !--- fatal error
       log_error("ATMOS_PHY_AE_kajino13_setup",*) 'Not appropriate names in namelist PARAM_ATMOS_PHY_AE_KAJINO13. Check!'
       call prc_abort
    endif
    log_nml(param_atmos_phy_ae_kajino13)

    !--- now only the default setting is supported
!    n_siz(1:n_ctg) = NSIZ(1:n_ctg)       ! number of size bins
!    d_min(1:n_ctg) = d_min_inp(1:n_ctg)  ! lower bound of 1st size bin
!    d_max(1:n_ctg) = d_max_inp(1:n_ctg)  ! upper bound of last size bin
!    n_kap(1:n_ctg) = n_kap_inp(1:n_ctg)  ! number of kappa bins
!    k_min(1:n_ctg) = k_min_inp(1:n_ctg)  ! lower bound of 1st kappa bin
!    k_max(1:n_ctg) = k_max_inp(1:n_ctg)  ! upper bound of last kappa bin

    !--- diagnose parameters (n_trans, n_siz_max, n_kap_max)
    n_trans   = 0
    n_siz_max = 0
    n_kap_max = 0
    do ic = 1, n_ctg
      n_trans   = n_trans + n_siz(ic) * n_kap(ic) * n_atr
      n_siz_max = max(n_siz_max, n_siz(ic))
      n_kap_max = max(n_kap_max, n_kap(ic))
    enddo

    !--- bin settings
    allocate(d_lw(n_siz_max,n_ctg))
    allocate(d_ct(n_siz_max,n_ctg))
    allocate(d_up(n_siz_max,n_ctg))
    allocate(k_lw(n_kap_max,n_ctg))
    allocate(k_ct(n_kap_max,n_ctg))
    allocate(k_up(n_kap_max,n_ctg))
    d_lw(:,:) = 0.0_rp
    d_ct(:,:) = 0.0_rp
    d_up(:,:) = 0.0_rp
    k_lw(:,:) = 0.0_rp
    k_ct(:,:) = 0.0_rp
    k_up(:,:) = 0.0_rp

    do ic = 1, n_ctg
      dlogd = (log(d_max(ic)) - log(d_min(ic)))/float(n_siz(ic))
      do is0 = 1, n_siz(ic)  !size bin
        d_lw(is0,ic) = exp(log(d_min(ic))+dlogd* float(is0-1)      )
        d_ct(is0,ic) = exp(log(d_min(ic))+dlogd*(float(is0)-0.5_rp))
        d_up(is0,ic) = exp(log(d_min(ic))+dlogd* float(is0)        )
      enddo !is (1:n_siz(ic))

      dk  = (k_max(ic) - k_min(ic))/float(n_kap(ic))
      do ik = 1, n_kap(ic)  !size bin
        k_lw(ik,ic) = k_min(ic) + dk  * float(ik-1)
        k_ct(ik,ic) = k_min(ic) + dk  *(float(ik)-0.5_rp)
        k_up(ik,ic) = k_min(ic) + dk  * float(ik)
      enddo !ik (1:n_kap(ic))

    enddo !ic (1:n_ctg)

!find size bin of regenerated aerosols
    do is0 = 1, n_siz(ic_mix)
      if ( atmos_phy_ae_kajino13_dg_reg >= d_lw(is0,ic_mix) .AND. &
           atmos_phy_ae_kajino13_dg_reg < d_up(is0,ic_mix) ) then
       is0_reg = is0
      endif !d_lw < ATMOS_PHY_AE_KAJINO13_dg_reg < d_up
    enddo

    !--- coagulation rule
    !   [ NOTE: current version has one category and single
    !     hygroscopicity bins and thus does not consider inter-category
    !     nor inter-hygroscopicity-section coagulation ]
    mcomb = 0
    do ic = 1, n_ctg     !=1
    do ik = 1, n_kap(ic) !=1
      do is2 = 1, n_siz(ic)
      do is1 = 1, n_siz(ic)
        if ( d_ct(is2,ic) >= d_ct(is1,ic) ) then
          mcomb = mcomb + 1
        endif !d_ct(is2) >= d_ct(is1)
      enddo !is1 (1:n_siz(ic)  )
      enddo !is2 (1:n_siz(ic)  )
    enddo !ik(1:n_kap(ic))
    enddo !ic(1:n_ctg)

    allocate(is_i(mcomb))
    allocate(is_j(mcomb))
    allocate(is_k(mcomb))
    allocate(ik_i(mcomb))
    allocate(ik_j(mcomb))
    allocate(ik_k(mcomb))
    allocate(ic_i(mcomb))
    allocate(ic_j(mcomb))
    allocate(ic_k(mcomb))

    mc = 0
    do ic = 1, n_ctg     !=1
    do ik = 1, n_kap(ic) !=1
      do is2 = 1, n_siz(ic)
      do is1 = 1, n_siz(ic)
        if ( d_ct(is2,ic) >= d_ct(is1,ic) ) then
          mc = mc + 1
          is_i(mc) = is1
          ik_i(mc) = ik
          ic_i(mc) = ic
          is_j(mc) = is2
          ik_j(mc) = ik
          ic_j(mc) = ic
          is_k(mc) = is2
          ik_k(mc) = ik
          ic_k(mc) = ic
        endif !d_ct(is2) >= d_ct(is1)
      enddo !is1 (1:n_siz(ic)  )
      enddo !is2 (1:n_siz(ic)  )
    enddo !ik(1:n_kap(ic))
    enddo !ic(1:n_ctg)

    !--- gas concentration
!    conc_h2so4 = 0.0_RP
!    conc_cgas = 0.0_RP

    allocate( it_procs2trans(n_atr,n_siz_max,n_kap_max,n_ctg)  )
    allocate( ia_trans2procs(n_trans) )
    allocate( is_trans2procs(n_trans) )
    allocate( ik_trans2procs(n_trans) )
    allocate( ic_trans2procs(n_trans) )

    it_procs2trans(:,:,:,:)= -999
    ia_trans2procs(:)      = 0
    is_trans2procs(:)      = 0
    ik_trans2procs(:)      = 0
    ic_trans2procs(:)      = 0

    !--- get pointer for trans2procs, procs2trans
    it = 0
    do ic = 1, n_ctg       !aerosol category
    do ik = 1, n_kap(ic)   !kappa bin
    do is0 = 1, n_siz(ic)  !size bin
    do ia0 = 1, n_atr       !attributes
      it = it + 1
      it_procs2trans(ia0,is0,ik,ic)= it
      ia_trans2procs(it)         = ia0
      is_trans2procs(it)         = is0
      ik_trans2procs(it)         = ik
      ic_trans2procs(it)         = ic
    enddo !ia (1:n_atr_prog )
    enddo !is (1:n_siz(ic)  )
    enddo !ik (1:n_kap(ic)  )
    enddo !ic (1:n_ctg      )

    allocate( aerosol_procs(n_atr,n_siz_max,n_kap_max,n_ctg)  )
    allocate( aerosol_activ(n_atr,n_siz_max,n_kap_max,n_ctg)  )
    allocate( emis_procs(n_atr,n_siz_max,n_kap_max,n_ctg)  )
    allocate( emis_gas(gas_ctg)  )
    aerosol_procs(:,:,:,:) = 0.0_rp
    aerosol_activ(:,:,:,:) = 0.0_rp
    emis_procs(:,:,:,:) = 0.0_rp
    emis_gas(:)       = 0.0_rp

    return

Here is the call graph for this function:

Here is the caller graph for this function:

atmos_phy_ae_kajino13_tendency()

subroutine, public scale_atmos_phy_ae_kajino13::atmos_phy_ae_kajino13_tendency	(	integer, intent(in)	KA,
		integer, intent(in)	KS,
		integer, intent(in)	KE,
		integer, intent(in)	IA,
		integer, intent(in)	IS,
		integer, intent(in)	IE,
		integer, intent(in)	JA,
		integer, intent(in)	JS,
		integer, intent(in)	JE,
		integer, intent(in)	QA_AE,
		real(rp), dimension(ka,ia,ja), intent(in)	TEMP,
		real(rp), dimension(ka,ia,ja), intent(in)	PRES,
		real(rp), dimension(ka,ia,ja), intent(in)	QDRY,
		real(rp), dimension(ka,ia,ja), intent(in)	NREG,
		real(rp), dimension(ka,ia,ja), intent(in)	DENS,
		real(rp), dimension (ka,ia,ja), intent(in)	QV,
		real(rp), dimension(ka,ia,ja,qa_ae), intent(in)	QTRC,
		real(rp), dimension(ka,ia,ja,qa_ae), intent(in)	EMIT,
		real(dp), intent(in)	dt,
		real(rp), dimension(ka,ia,ja,qa_ae), intent(out)	RHOQ_t_AE,
		real(rp), dimension(ka,ia,ja), intent(out)	CN,
		real(rp), dimension(ka,ia,ja), intent(out)	CCN
	)

Aerosol Microphysics.

Definition at line 619 of file scale_atmos_phy_ae_kajino13.F90.

References atmos_phy_ae_kajino13_dg_reg, atmos_phy_ae_kajino13_flag_ccn_interactive, atmos_phy_ae_kajino13_flag_coag, atmos_phy_ae_kajino13_flag_cond, atmos_phy_ae_kajino13_flag_npf, atmos_phy_ae_kajino13_flag_regeneration, and atmos_phy_ae_kajino13_sg_reg.

Referenced by mod_atmos_phy_ae_driver::atmos_phy_ae_driver_calc_tendency().

    use scale_atmos_saturation, only: &
       pres2qsat_liq => atmos_saturation_pres2qsat_liq
    use scale_file_history, only: &
       file_history_in
    implicit none
    integer,  intent(in) :: ka, ks, ke
    integer,  intent(in) :: ia, is, ie
    integer,  intent(in) :: ja, js, je
    integer,  intent(in) :: qa_ae
    real(RP), intent(in) :: temp(ka,ia,ja)
    real(RP), intent(in) :: pres(ka,ia,ja)
    real(RP), intent(in) :: qdry(ka,ia,ja)
    real(RP), intent(in) :: nreg(ka,ia,ja)
    real(RP), intent(in) :: dens(ka,ia,ja)
    real(RP), intent(in) :: qv  (ka,ia,ja)
    real(RP), intent(in) :: qtrc(ka,ia,ja,qa_ae)
    real(RP), intent(in) :: emit(ka,ia,ja,qa_ae)
    real(DP), intent(in) :: dt

    real(RP), intent(out) :: rhoq_t_ae(ka,ia,ja,qa_ae)
    real(RP), intent(out) :: cn(ka,ia,ja)
    real(RP), intent(out) :: ccn(ka,ia,ja)

    real(RP) :: qtrc0(ka,ia,ja,qa_ae)
    real(RP) :: qtrc1(ka,ia,ja,qa_ae)

    !--- local
    real(RP) :: ssliq_ae(ka,ia,ja)
    real(RP) :: rrhog(ka,ia,ja)
!    real(RP) :: t_ccn, t_cn
    real(RP) :: qsat_tmp
    real(RP) :: m0_reg, m2_reg, m3_reg      !regenerated aerosols [m^k/m3]
    real(RP) :: ms_reg                      !regenerated aerosol mass [ug/m3]
    real(RP) :: reg_factor_m2,reg_factor_m3 !to save cpu time for moment conversion
    !--- aerosol variables
    real(RP) :: total_aerosol_mass(ka,ia,ja,n_ctg)
    real(RP) :: total_aerosol_number(ka,ia,ja,n_ctg)
    real(RP) :: total_emit_aerosol_mass(ka,ia,ja,n_ctg)
    real(RP) :: total_emit_aerosol_number(ka,ia,ja,n_ctg)
    !--- gas
!    real(RP),allocatable :: conc_h2so4(:,:,:,:)    !concentration [ug/m3]
!    real(RP) :: conc_gas(KA,IA,JA,GAS_CTG)    !concentration [ug/m3]
    real(RP) :: conc_gas(gas_ctg)    !concentration [ug/m3]
    integer :: i, j, k, iq, it

    log_progress(*) 'atmosphere / physics / aerosol / kajino13'

    aerosol_procs(:,:,:,:) = 0.0_rp
    aerosol_activ(:,:,:,:) = 0.0_rp
    emis_procs(:,:,:,:) = 0.0_rp
    emis_gas(:)       = 0.0_rp

    do iq = 1, qa_ae
    do j  = js, je
    do i  = is, ie
    do k  = ks, ke
       qtrc0(k,i,j,iq) = qtrc(k,i,j,iq) ! save
    enddo
    enddo
    enddo
    enddo

    reg_factor_m2 = atmos_phy_ae_kajino13_dg_reg**2._rp * exp( 2.0_rp *(log(atmos_phy_ae_kajino13_sg_reg)**2._rp)) !m0_reg to m2_reg
    reg_factor_m3 = atmos_phy_ae_kajino13_dg_reg**3._rp * exp( 4.5_rp *(log(atmos_phy_ae_kajino13_sg_reg)**2._rp)) !m0_reg to m3_reg

    !--- convert SCALE variable to zerochem variable

    do j = js, je
    do i = is, ie

       do k = ks, ke
          rrhog(k,i,j) = 1.0_rp / dens(k,i,j)
       enddo
       do k = ks, ke
          !--- calculate super saturation of water
          call pres2qsat_liq( temp(k,i,j), pres(k,i,j), qdry(k,i,j), & ! [IN]
                              qsat_tmp                               ) ! [OUT]
          ssliq_ae(k,i,j) = qv(k,i,j)/qsat_tmp - 1.0_rp
       enddo

    enddo
    enddo

  ! tiny number, tiny mass
    do j  = js, je
    do i  = is, ie
    do k  = ks, ke
      do ic = 1, n_ctg       !aerosol category
      do ik = 1, n_kap(ic)   !kappa bin
      do is0 = 1, n_siz(ic)   !size bin
        if (qtrc0(k,i,j,it_procs2trans(ia_m0,is0,ik,ic))*dens(k,i,j) < cleannumber) then
          do ia0 = 1, n_atr
            qtrc0(k,i,j,it_procs2trans(ia0,is0,ik,ic)) = 0._rp !to save cpu time and avoid underflow
          enddo !ia0 (1:n_atr      )
        endif
      enddo !is (1:n_siz(ic)  )
      enddo !ik (1:n_kap(ic)  )
      enddo !ic (1:n_ctg      )
    enddo
    enddo
    enddo

    do iq = 1, qa_ae
    do j  = js, je
    do i  = is, ie
    do k  = ks, ke
      qtrc1(k,i,j,iq) = qtrc0(k,i,j,iq) ! save
    enddo
    enddo
    enddo
    enddo

    !---- Calculate aerosol processs
    cn(:,:,:) = 0.0_rp
    ccn(:,:,:) = 0.0_rp
    do k = ks, ke
    do j = js, je
    do i = is, ie

       !--- aerosol_trans at initial time
       ! [xx/kg] -> [xx/m3]
       do it = 1, n_trans
            aerosol_procs(ia_trans2procs(it), &
                          is_trans2procs(it), &
                          ik_trans2procs(it), &
                          ic_trans2procs(it)) = qtrc1(k,i,j,it)*dens(k,i,j)
            emis_procs(ia_trans2procs(it), &
                          is_trans2procs(it), &
                          ik_trans2procs(it), &
                          ic_trans2procs(it)) = emit(k,i,j,it)*dens(k,i,j)
       enddo !it(1:n_trans)
       ! mixing ratio [kg/kg] -> concentration [ug/m3]
       conc_gas(1:gas_ctg) &
               = qtrc1(k,i,j,qa_ae-gas_ctg+1:qa_ae-gas_ctg+gas_ctg)*dens(k,i,j)*1.e+9_rp

       emis_gas(1:gas_ctg) = emit(k,i,j,qa_ae-gas_ctg+ig_h2so4:qa_ae-gas_ctg+ig_cgas)

       call aerosol_zerochem( &
         dt,                              & !--- in
         temp(k,i,j),                     & !--- in
         pres(k,i,j),                     & !--- in
         ssliq_ae(k,i,j),                 & !--- in
         atmos_phy_ae_kajino13_flag_npf,  & !--- in
         atmos_phy_ae_kajino13_flag_cond, & !--- in
         atmos_phy_ae_kajino13_flag_coag, & !--- in
         aerosol_procs,                   & !--- inout
         conc_gas,                        & !--- inout
         emis_procs,                      & !--- out
         emis_gas,                        & !--- out
         aerosol_activ                    ) !--- out

! aerosol loss due to activation to cloud droplets
       if (atmos_phy_ae_kajino13_flag_ccn_interactive) then
         do is0 = 1, n_siz(ic_mix)
         do ia0 = 1, n_atr       !attributes
           aerosol_activ(ia0,is0,ik_out,ic_mix) = min(max(0._rp, aerosol_activ(ia0,is0,ik_out,ic_mix)), &
                                                                 aerosol_procs(ia0,is0,ik_out,ic_mix) )

           aerosol_procs(ia0,is0,ik_out,ic_mix) = &
           aerosol_procs(ia0,is0,ik_out,ic_mix) - aerosol_activ(ia0,is0,ik_out,ic_mix)
         enddo
         enddo
       endif !flag_ccn_interactive

! aerosol regeneration due to evaporation of cloud droplets
! using prescribed size parameters and to internal mixture category (ic_mix)
       if (atmos_phy_ae_kajino13_flag_regeneration) then
         m0_reg = nreg(k,i,j)  !#/m3
!        m2_reg = m0_reg * ATMOS_PHY_AE_KAJINO13_dg_reg**2._RP * exp( 2.0_RP *(log(ATMOS_PHY_AE_KAJINO13_sg_reg)**2._RP)) !m2/m3
!        m3_reg = m0_reg * ATMOS_PHY_AE_KAJINO13_dg_reg**3._RP * exp( 4.5_RP *(log(ATMOS_PHY_AE_KAJINO13_sg_reg)**2._RP)) !m3/m3
         m2_reg = m0_reg * reg_factor_m2     !m2/m3
         m3_reg = m0_reg * reg_factor_m3     !m3/m3
         ms_reg = m3_reg * pi6 * conv_vl_ms  !ug/m3
         aerosol_procs(ia_m0,is0_reg,ik_out,ic_mix) = &
         aerosol_procs(ia_m0,is0_reg,ik_out,ic_mix) + m0_reg !#/m3
         aerosol_procs(ia_m2,is0_reg,ik_out,ic_mix) = &
         aerosol_procs(ia_m2,is0_reg,ik_out,ic_mix) + m2_reg !m2/m3
         aerosol_procs(ia_m3,is0_reg,ik_out,ic_mix) = &
         aerosol_procs(ia_m3,is0_reg,ik_out,ic_mix) + m3_reg !m3/m3
         aerosol_procs(ia_ms,is0_reg,ik_out,ic_mix) = &
         aerosol_procs(ia_ms,is0_reg,ik_out,ic_mix) + ms_reg !ug/m3
! additional attirbute to be added (ia_kp)
       endif !flag_regeneration

! diagnosed variables
       do is0 = 1, n_siz(ic_mix)
         ccn(k,i,j) = ccn(k,i,j) + aerosol_activ(ia_m0,is0,ik_out,ic_mix)
         cn(k,i,j) = cn(k,i,j) + aerosol_procs(ia_m0,is0,ik_out,ic_mix)
       enddo

!       call trans_ccn(aerosol_procs, aerosol_activ, t_ccn, t_cn,  &
!            n_ctg, n_kap_max, n_siz_max, N_ATR,         &
!            ic_mix, ia_m0, ia_m2, ia_m3, ik_out, n_siz, &
!            rnum_out, ATMOS_PHY_AE_KAJINO13_nbins_out)

!       CN(k,i,j) = t_cn
!       CCN(k,i,j) = t_ccn

       ! [xx/m3] -> [xx/kg]
       do ic = 1, n_ctg       !category
       do ik = 1, n_kap(ic)   !kappa bin
       do is0 = 1, n_siz(ic)   !size bin
       do ia0 = 1, n_atr       !attributes
          qtrc1(k,i,j,it_procs2trans(ia0,is0,ik,ic)) = aerosol_procs(ia0,is0,ik,ic) / dens(k,i,j)
       enddo !ia (1:N_ATR )
       enddo !is (1:n_siz(ic)  )
       enddo !ik (1:n_kap(ic)  )
       enddo !ic (1:n_ctg      )
       !  [ug/m3] -> mixing ratio [kg/kg]
       qtrc1(k,i,j,qa_ae-gas_ctg+1:qa_ae-gas_ctg+gas_ctg) = conc_gas(1:gas_ctg) / dens(k,i,j)*1.e-9_rp

       ! tiny number, tiny mass
       do ic = 1, n_ctg       !aerosol category
       do ik = 1, n_kap(ic)   !kappa bin
       do is0 = 1, n_siz(ic)   !size bin
         if (qtrc1(k,i,j,it_procs2trans(ia_m0,is0,ik,ic))*dens(k,i,j) < cleannumber) then
           do ia0 = 1, n_atr
             qtrc1(k,i,j,it_procs2trans(ia0,is0,ik,ic)) = 0._rp !to save cpu time and avoid underflow
           enddo !ia0 (1:n_atr      )
         endif
       enddo !is (1:n_siz(ic)  )
       enddo !ik (1:n_kap(ic)  )
       enddo !ic (1:n_ctg      )

       ! for history
       total_aerosol_mass(k,i,j,:) = 0.0_rp
       total_aerosol_number(k,i,j,:) = 0.0_rp
       total_emit_aerosol_mass(k,i,j,:) = 0.0_rp
       total_emit_aerosol_number(k,i,j,:) = 0.0_rp
       do ic = 1, n_ctg
       do ik = 1, n_kap(ic)
       do is0 = 1, n_siz(ic)
           total_aerosol_mass(k,i,j,ic) = total_aerosol_mass(k,i,j,ic) &
                                               + qtrc1(k,i,j,it_procs2trans(ia_ms,is0,ik,ic))
           total_aerosol_number(k,i,j,ic) = total_aerosol_number(k,i,j,ic) &
                                               + qtrc1(k,i,j,it_procs2trans(ia_m0,is0,ik,ic))
           total_emit_aerosol_mass(k,i,j,ic) = total_emit_aerosol_mass(k,i,j,ic) &
                                               + emit(k,i,j,it_procs2trans(ia_ms,is0,ik,ic))
           total_emit_aerosol_number(k,i,j,ic) = total_emit_aerosol_number(k,i,j,ic) &
                                               + emit(k,i,j,it_procs2trans(ia_m0,is0,ik,ic))
       enddo
       enddo
       enddo

    enddo
    enddo
    enddo

    do ic = 1, n_ctg
      call file_history_in( total_aerosol_mass(:,:,:,ic), trim(ctg_name(ic))//'_mass', 'Total mass mixing ratio of aerosol', 'kg/kg' )
      call file_history_in( total_aerosol_number(:,:,:,ic), trim(ctg_name(ic))//'_number', 'Total number mixing ratio of aerosol', 'num/kg' )
      call file_history_in( total_emit_aerosol_mass(:,:,:,ic), trim(ctg_name(ic))//'_mass_emit', 'Total mass mixing ratio of emitted aerosol', 'kg/kg' )
      call file_history_in( total_emit_aerosol_number(:,:,:,ic), trim(ctg_name(ic))//'_number_emit', 'Total number mixing ratio of emitted aerosol', 'num/kg' )
    enddo

    call file_history_in( emit(:,:,:,qa_ae-gas_ctg+ig_h2so4), 'H2SO4_emit', 'Emission ratio of H2SO4 gas', 'ug/m3/s' )
    call file_history_in( emit(:,:,:,qa_ae-gas_ctg+ig_cgas),  'CGAS_emit',  'Emission ratio of Condensabule gas', 'ug/m3/s' )


    do iq = 1, qa_ae
    do j  = js, je
    do i  = is, ie
    do k  = ks, ke
      rhoq_t_ae(k,i,j,iq) = ( qtrc1(k,i,j,iq) - qtrc0(k,i,j,iq) ) * dens(k,i,j) / dt
    enddo
    enddo
    enddo
    enddo

    return

Here is the caller graph for this function:

atmos_phy_ae_kajino13_effective_radius()

subroutine, public scale_atmos_phy_ae_kajino13::atmos_phy_ae_kajino13_effective_radius	(	integer, intent(in)	KA,
		integer, intent(in)	IA,
		integer, intent(in)	JA,
		integer, intent(in)	QA_AE,
		real(rp), dimension(ka,ia,ja,qa_ae), intent(in)	QTRC,
		real(rp), dimension (ka,ia,ja), intent(in)	RH,
		real(rp), dimension (ka,ia,ja,n_ae), intent(out)	Re
	)

Calculate Effective Radius.

Definition at line 897 of file scale_atmos_phy_ae_kajino13.F90.

References scale_atmos_aerosol::n_ae.

Referenced by mod_atmos_phy_ae_vars::atmos_phy_ae_vars_get_diagnostic().

    use scale_atmos_aerosol, only: &
       n_ae
    implicit none
    integer,  intent(in)  :: ka, ia, ja, qa_ae

    real(RP), intent(in)  :: qtrc(ka,ia,ja,qa_ae) ! tracer mass concentration [kg/kg]
    real(RP), intent(in)  :: rh  (ka,ia,ja)       ! relative humidity         (0-1)

    real(RP), intent(out) :: re  (ka,ia,ja,n_ae)  ! effective radius
    !---------------------------------------------------------------------------

    re(:,:,:,:) = 0.0_rp

!    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

Here is the caller graph for this function:

atmos_phy_ae_kajino13_mkinit()

subroutine, public scale_atmos_phy_ae_kajino13::atmos_phy_ae_kajino13_mkinit	(	integer, intent(in)	KA,
		integer, intent(in)	KS,
		integer, intent(in)	KE,
		integer, intent(in)	IA,
		integer, intent(in)	IS,
		integer, intent(in)	IE,
		integer, intent(in)	JA,
		integer, intent(in)	JS,
		integer, intent(in)	JE,
		integer, intent(in)	QA_AE,
		real(rp), dimension(ka,ia,ja), intent(in)	DENS,
		real(rp), dimension(ka,ia,ja), intent(in)	TEMP,
		real(rp), dimension(ka,ia,ja), intent(in)	PRES,
		real(rp), dimension(ka,ia,ja), intent(in)	QDRY,
		real(rp), dimension (ka,ia,ja), intent(in)	QV,
		real(rp), intent(in)	m0_init,
		real(rp), intent(in)	dg_init,
		real(rp), intent(in)	sg_init,
		real(rp), dimension(3), intent(in)	d_min_inp,
		real(rp), dimension(3), intent(in)	d_max_inp,
		real(rp), dimension(3), intent(in)	k_min_inp,
		real(rp), dimension(3), intent(in)	k_max_inp,
		integer, dimension(3), intent(in)	n_kap_inp,
		real(rp), dimension(ka,ia,ja,qa_ae), intent(out)	QTRC,
		real(rp), dimension (ka,ia,ja), intent(out)	CCN
	)

Definition at line 940 of file scale_atmos_phy_ae_kajino13.F90.

References atmos_phy_ae_kajino13_setup(), and scale_const::const_pi.

Referenced by mod_mkinit::rect_setup().

    use scale_const, only: &
       pi => const_pi
    use scale_atmos_saturation, only: &
       saturation_pres2qsat_liq => atmos_saturation_pres2qsat_liq
    implicit none

    integer,  intent(in)  :: ka, ks, ke
    integer,  intent(in)  :: ia, is, ie
    integer,  intent(in)  :: ja, js, je
    integer,  intent(in)  :: qa_ae

    real(RP), intent(in)  :: dens(ka,ia,ja)
    real(RP), intent(in)  :: temp(ka,ia,ja)
    real(RP), intent(in)  :: pres(ka,ia,ja)
    real(RP), intent(in)  :: qdry(ka,ia,ja)
    real(RP), intent(in)  :: qv  (ka,ia,ja)
    real(RP), intent(in)  :: m0_init             ! initial total num. conc. of modes (Atk,Acm,Cor) [#/m3]
    real(RP), intent(in)  :: dg_init             ! initial number equivalen diameters of modes     [m]
    real(RP), intent(in)  :: sg_init             ! initial standard deviation                      [-]
    real(RP), intent(in)  :: d_min_inp(3)
    real(RP), intent(in)  :: d_max_inp(3)
    real(RP), intent(in)  :: k_min_inp(3)
    real(RP), intent(in)  :: k_max_inp(3)
    integer,  intent(in)  :: n_kap_inp(3)

    real(RP), intent(out) :: qtrc(ka,ia,ja,qa_ae)
    real(RP), intent(out) :: ccn (ka,ia,ja)

    integer,  parameter :: ia_m0  = 1             ! 1. number conc        [#/m3]
    integer,  parameter :: ia_m2  = 2             ! 2. 2nd mom conc       [m2/m3]
    integer,  parameter :: ia_m3  = 3             ! 3. 3rd mom conc       [m3/m3]
    integer,  parameter :: ia_ms  = 4             ! 4. mass conc          [ug/m3]
    integer,  parameter :: ia_kp  = 5
    integer,  parameter :: ik_out = 1

    real(RP)            :: c_kappa     = 0.3_rp     ! hygroscopicity of condensable mass [-]
    real(RP), parameter :: cleannumber = 1.e-3_rp
    ! local variables
    real(RP), parameter :: rhod_ae    = 1.83_rp              ! particle density [g/cm3] sulfate assumed
    real(RP), parameter :: conv_vl_ms = rhod_ae/1.e-12_rp     ! M3(volume)[m3/m3] to mass[m3/m3]

    real(RP) :: m0t, dgt, sgt, m2t, m3t, mst
    !--- gas
    real(RP) :: conc_gas(gas_ctg)    !concentration [ug/m3]
    !--- bin  settings (lower bound, center, upper bound)
    real(RP), allocatable :: d_lw(:,:), d_ct(:,:), d_up(:,:)  !diameter [m]
    real(RP), allocatable :: k_lw(:,:), k_ct(:,:), k_up(:,:)  !kappa    [-]
    real(RP)  :: dlogd, dk                                    !delta log(D), delta K
    real(RP), allocatable :: d_min(:)              !lower bound of 1st size bin   (n_ctg)
    real(RP), allocatable :: d_max(:)              !upper bound of last size bin  (n_ctg)
    integer,  allocatable :: n_kap(:)              !number of kappa bins          (n_ctg)
    real(RP), allocatable :: k_min(:)              !lower bound of 1st kappa bin  (n_ctg)
    real(RP), allocatable :: k_max(:)

    real(RP) :: qsat_tmp, ssliq
    real(RP) :: pi6

    integer  :: ia0, ik, is0, ic, k, i, j, it
    !---------------------------------------------------------------------------


    call atmos_phy_ae_kajino13_setup


    pi6   = pi / 6._rp
    n_ctg = ae_ctg

    allocate( d_min(n_ctg) )
    allocate( d_max(n_ctg) )
    allocate( n_kap(n_ctg) )
    allocate( k_min(n_ctg) )
    allocate( k_max(n_ctg) )


    d_min(1:n_ctg) = d_min_inp(1:n_ctg)  ! lower bound of 1st size bin
    d_max(1:n_ctg) = d_max_inp(1:n_ctg)  ! upper bound of last size bin
    n_kap(1:n_ctg) = n_kap_inp(1:n_ctg)  ! number of kappa bins
    k_min(1:n_ctg) = k_min_inp(1:n_ctg)  ! lower bound of 1st kappa bin
    k_max(1:n_ctg) = k_max_inp(1:n_ctg)  ! upper bound of last kappa bin


    !bin setting
    allocate(d_lw(n_siz_max,n_ctg))
    allocate(d_ct(n_siz_max,n_ctg))
    allocate(d_up(n_siz_max,n_ctg))
    allocate(k_lw(n_kap_max,n_ctg))
    allocate(k_ct(n_kap_max,n_ctg))
    allocate(k_up(n_kap_max,n_ctg))
    d_lw(:,:) = 0._rp
    d_ct(:,:) = 0._rp
    d_up(:,:) = 0._rp
    k_lw(:,:) = 0._rp
    k_ct(:,:) = 0._rp
    k_up(:,:) = 0._rp

    do ic = 1, n_ctg

      dlogd = (log(d_max(ic)) - log(d_min(ic)))/real(NSIZ(ic),kind=rp)

      do is0 = 1, nsiz(ic)  !size bin
        d_lw(is0,ic) = exp(log(d_min(ic))+dlogd* real(is0-1,kind=RP)        )
        d_ct(is0,ic) = exp(log(d_min(ic))+dlogd*(real(is0  ,kind=rp)-0.5_rp))
        d_up(is0,ic) = exp(log(d_min(ic))+dlogd* real(is0  ,kind=RP)        )
      enddo !is (1:n_siz(ic))
      dk    = (k_max(ic) - k_min(ic))/real(n_kap(ic),kind=rp)
      do ik = 1, n_kap(ic)  !size bin
        k_lw(ik,ic) = k_min(ic) + dk  * real(ik-1,kind=rp)
        k_ct(ik,ic) = k_min(ic) + dk  *(real(ik  ,kind=rp)-0.5_rp)
        k_up(ik,ic) = k_min(ic) + dk  * real(ik  ,kind=rp)
      enddo !ik (1:n_kap(ic))

    enddo !ic (1:n_ctg)
!    ik  = 1       !only one kappa bin

    m0t = m0_init !total M0 [#/m3]
    dgt = dg_init ![m]
    sgt = sg_init ![-]

    if ( m0t <= cleannumber ) then
       m0t = cleannumber
       dgt = 0.1e-6_rp
       sgt = 1.3_rp
       log_warn("ATMOS_PHY_AE_kajino13_mkinit",*) 'Initial aerosol number is set as ', cleannumber, '[#/m3]'
    endif

    m2t = m0t*dgt**(2.d0) *dexp(2.0d0 *(dlog(real(sgt,kind=dp))**2.d0)) !total M2 [m2/m3]
    m3t = m0t*dgt**(3.d0) *dexp(4.5d0 *(dlog(real(sgt,kind=dp))**2.d0)) !total M3 [m3/m3]
    mst = m3t*pi6*conv_vl_ms                              !total Ms [ug/m3]

    do ic = 1, n_ctg
    !aerosol_procs initial condition
    do ik = 1, n_kap(ic)   !kappa bin
    do is0 = 1, nsiz(ic)
       if (dgt >= d_lw(is0,ic) .and. dgt < d_up(is0,ic)) then
         aerosol_procs(ia_m0,is0,ik,ic) = aerosol_procs(ia_m0,is0,ik,ic) + m0t          ![#/m3]
         aerosol_procs(ia_m2,is0,ik,ic) = aerosol_procs(ia_m2,is0,ik,ic) + m2t          ![m2/m3]
         aerosol_procs(ia_m3,is0,ik,ic) = aerosol_procs(ia_m3,is0,ik,ic) + m3t          ![m3/m3]
         aerosol_procs(ia_ms,is0,ik,ic) = aerosol_procs(ia_ms,is0,ik,ic) + mst*1.e-9_rp ![kg/m3]
       elseif (dgt < d_lw(1,ic)) then
         aerosol_procs(ia_m0,1 ,ik,ic) = aerosol_procs(ia_m0,1 ,ik,ic) + m0t          ![#/m3]
         aerosol_procs(ia_m2,1 ,ik,ic) = aerosol_procs(ia_m2,1 ,ik,ic) + m2t          ![m2/m3]
         aerosol_procs(ia_m3,1 ,ik,ic) = aerosol_procs(ia_m3,1 ,ik,ic) + m3t          ![m3/m3]
         aerosol_procs(ia_ms,1 ,ik,ic) = aerosol_procs(ia_ms,1 ,ik,ic) + mst*1.e-9_rp ![kg/m3]
       elseif (dgt >= d_up(nsiz(ic),ic)) then
         aerosol_procs(ia_m0,nsiz(ic),ik,ic) = aerosol_procs(ia_m0,nsiz(ic),ik,ic) + m0t          ![#/m3]
         aerosol_procs(ia_m2,nsiz(ic),ik,ic) = aerosol_procs(ia_m2,nsiz(ic),ik,ic) + m2t          ![m2/m3]
         aerosol_procs(ia_m3,nsiz(ic),ik,ic) = aerosol_procs(ia_m3,nsiz(ic),ik,ic) + m3t          ![m3/m3]
         aerosol_procs(ia_ms,nsiz(ic),ik,ic) = aerosol_procs(ia_ms,nsiz(ic),ik,ic) + mst*1.e-9_rp ![kg/m3]
       endif
    enddo
    enddo
    enddo

    conc_gas(:) = 0.0_rp
    do j = 1, ja
    do i = 1, ia
    do k = 1, ka
       it = 1
       do ic  = 1, n_ctg    ! category
       do ik  = 1, nkap(ic) ! kappa bin
       do is0 = 1, nsiz(ic) ! size bin
       do ia0 = 1, n_atr    ! attributes
          qtrc(k,i,j,it) = aerosol_procs(ia0,is0,ik,ic) / dens(k,i,j) !#,m2,m3,kg/m3 -> #,m2,m3,kg/kg
          it = it + 1
       enddo !ia0 (1:N_ATR )
       enddo !is (1:n_siz(ic)  )
       enddo !ik (1:n_kap(ic)  )
       enddo !ic (1:n_ctg      )

       do ic  = 1, gas_ctg ! GAS category
          qtrc(k,i,j,it) = conc_gas(ic) / dens(k,i,j) * 1.e-9_rp !mixing ratio [kg/kg]
          it = it + 1
       enddo
    enddo
    enddo
    enddo

    ccn(:,:,:) = 0.0_rp
    do j = js, je
    do i = is, ie
    do k = ks, ke
       !--- calculate super saturation of water
       call saturation_pres2qsat_liq( temp(k,i,j), pres(k,i,j), qdry(k,i,j), & ! [IN]
                                      qsat_tmp                               ) ! [OUT]
       ssliq = qv(k,i,j) / qsat_tmp - 1.0_rp

       call aerosol_activation( c_kappa, ssliq, temp(k,i,j), ia_m0, ia_m2, ia_m3,       &
                                n_atr, n_siz_max, n_kap_max, n_ctg, nsiz, n_kap, &
                                d_ct, aerosol_procs, aerosol_activ               )

       do is0 = 1, nsiz(ic_mix)
          ccn(k,i,j) = ccn(k,i,j) + aerosol_activ(ia_m0,is0,ik_out,ic_mix)
       enddo
    enddo
    enddo
    enddo

    return

Here is the call graph for this function:

Here is the caller graph for this function:

atmos_phy_ae_kajino13_negative_fixer()

subroutine, public scale_atmos_phy_ae_kajino13::atmos_phy_ae_kajino13_negative_fixer	(	integer, intent(in)	KA,
		integer, intent(in)	KS,
		integer, intent(in)	KE,
		integer, intent(in)	IA,
		integer, intent(in)	IS,
		integer, intent(in)	IE,
		integer, intent(in)	JA,
		integer, intent(in)	JS,
		integer, intent(in)	JE,
		integer, intent(in)	QA_AE,
		real(rp), dimension(ka,ia,ja,qa_ae), intent(inout)	QTRC
	)

Definition at line 1143 of file scale_atmos_phy_ae_kajino13.F90.

References atmos_phy_ae_kajino13_c_kappa, atmos_phy_ae_kajino13_h2so4dt, atmos_phy_ae_kajino13_logk_aenucl, atmos_phy_ae_kajino13_ocgasdt, float(), scale_atmos_grid_cartesc_index::ke, scale_atmos_grid_cartesc_index::ks, scale_prof::prof_rapend(), and scale_prof::prof_rapstart().

Referenced by mod_atmos_phy_ae_driver::atmos_phy_ae_driver_adjustment().

    integer, intent(in) :: ka, ks, ke
    integer, intent(in) :: ia, is, ie
    integer, intent(in) :: ja, js, je
    integer, intent(in) :: qa_ae

    real(RP), intent(inout) :: qtrc(ka,ia,ja,qa_ae)

    integer :: k ,i, j
    integer :: ic, ik, is0

    !--- Negative fixer
    do j  = js, je
    do i  = is, ie
    do k  = ks, ke
      do ic = 1, n_ctg       !aerosol category
      do ik = 1, n_kap(ic)   !kappa bin
      do is0 = 1, n_siz(ic)   !size bin
        if (qtrc(k,i,j,it_procs2trans(ia_m0,is0,ik,ic)) < 0.0_rp .or. &
            qtrc(k,i,j,it_procs2trans(ia_m2,is0,ik,ic)) < 0.0_rp .or. &
            qtrc(k,i,j,it_procs2trans(ia_m3,is0,ik,ic)) < 0.0_rp .or. &
            qtrc(k,i,j,it_procs2trans(ia_ms,is0,ik,ic)) < 0.0_rp .or. &
            qtrc(k,i,j,it_procs2trans(ia_kp,is0,ik,ic)) < 0.0_rp ) then
          qtrc(k,i,j,it_procs2trans(1:n_atr,is0,ik,ic)) = 0.0_rp
        endif
      enddo
      enddo
      enddo
    enddo
    enddo
    enddo

    return

Here is the call graph for this function:

Here is the caller graph for this function:

Variable Documentation

atmos_phy_ae_kajino13_name

character(len=h_short), dimension(:), allocatable, public scale_atmos_phy_ae_kajino13::atmos_phy_ae_kajino13_name

Definition at line 51 of file scale_atmos_phy_ae_kajino13.F90.

Referenced by mod_atmos_phy_ae_driver::atmos_phy_ae_driver_tracer_setup(), and atmos_phy_ae_kajino13_tracer_setup().

  character(len=H_SHORT), public, allocatable :: atmos_phy_ae_kajino13_name(:)

atmos_phy_ae_kajino13_desc

character(len=h_mid), dimension(:), allocatable, public scale_atmos_phy_ae_kajino13::atmos_phy_ae_kajino13_desc

Definition at line 52 of file scale_atmos_phy_ae_kajino13.F90.

Referenced by mod_atmos_phy_ae_driver::atmos_phy_ae_driver_tracer_setup(), and atmos_phy_ae_kajino13_tracer_setup().

  character(len=H_MID)  , public, allocatable :: atmos_phy_ae_kajino13_desc(:)

atmos_phy_ae_kajino13_unit

character(len=h_short), dimension(:), allocatable, public scale_atmos_phy_ae_kajino13::atmos_phy_ae_kajino13_unit

Definition at line 53 of file scale_atmos_phy_ae_kajino13.F90.

Referenced by mod_atmos_phy_ae_driver::atmos_phy_ae_driver_tracer_setup(), and atmos_phy_ae_kajino13_tracer_setup().

  character(len=H_SHORT), public, allocatable :: atmos_phy_ae_kajino13_unit(:)

atmos_phy_ae_kajino13_h2so4dt

real(rp), public scale_atmos_phy_ae_kajino13::atmos_phy_ae_kajino13_h2so4dt = 5.E-6_RP

Definition at line 55 of file scale_atmos_phy_ae_kajino13.F90.

Referenced by atmos_phy_ae_kajino13_negative_fixer(), and atmos_phy_ae_kajino13_setup().

  real(RP), public :: atmos_phy_ae_kajino13_h2so4dt = 5.e-6_rp   ! h2so4 production rate (Temporal)                [ug/m3/s]

atmos_phy_ae_kajino13_ocgasdt

real(rp), public scale_atmos_phy_ae_kajino13::atmos_phy_ae_kajino13_ocgasdt = 8.E-5_RP

Definition at line 56 of file scale_atmos_phy_ae_kajino13.F90.

Referenced by atmos_phy_ae_kajino13_negative_fixer(), and atmos_phy_ae_kajino13_setup().

  real(RP), public :: atmos_phy_ae_kajino13_ocgasdt = 8.e-5_rp   ! other condensational bas production rate 16*h2so4dt (see Kajino et al. 2013)

atmos_phy_ae_kajino13_c_kappa

real(rp), public scale_atmos_phy_ae_kajino13::atmos_phy_ae_kajino13_c_kappa = 0.3_RP

Definition at line 58 of file scale_atmos_phy_ae_kajino13.F90.

Referenced by atmos_phy_ae_kajino13_negative_fixer(), and atmos_phy_ae_kajino13_setup().

  real(RP), public :: atmos_phy_ae_kajino13_c_kappa = 0.3_rp     ! hygroscopicity of condensable mass              [-]

atmos_phy_ae_kajino13_flag_npf

logical, public scale_atmos_phy_ae_kajino13::atmos_phy_ae_kajino13_flag_npf = .false.

Definition at line 59 of file scale_atmos_phy_ae_kajino13.F90.

Referenced by atmos_phy_ae_kajino13_setup(), and atmos_phy_ae_kajino13_tendency().

  logical,  public :: atmos_phy_ae_kajino13_flag_npf = .false.

atmos_phy_ae_kajino13_flag_cond

logical, public scale_atmos_phy_ae_kajino13::atmos_phy_ae_kajino13_flag_cond = .true.

Definition at line 60 of file scale_atmos_phy_ae_kajino13.F90.

Referenced by atmos_phy_ae_kajino13_setup(), and atmos_phy_ae_kajino13_tendency().

  logical,  public :: atmos_phy_ae_kajino13_flag_cond = .true.

atmos_phy_ae_kajino13_flag_coag

logical, public scale_atmos_phy_ae_kajino13::atmos_phy_ae_kajino13_flag_coag = .true.

Definition at line 61 of file scale_atmos_phy_ae_kajino13.F90.

Referenced by atmos_phy_ae_kajino13_setup(), and atmos_phy_ae_kajino13_tendency().

  logical,  public :: atmos_phy_ae_kajino13_flag_coag = .true.

atmos_phy_ae_kajino13_flag_ccn_interactive

logical, public scale_atmos_phy_ae_kajino13::atmos_phy_ae_kajino13_flag_ccn_interactive = .true.

Definition at line 62 of file scale_atmos_phy_ae_kajino13.F90.

Referenced by atmos_phy_ae_kajino13_setup(), and atmos_phy_ae_kajino13_tendency().

  logical,  public :: atmos_phy_ae_kajino13_flag_ccn_interactive = .true.

atmos_phy_ae_kajino13_flag_regeneration

logical, public scale_atmos_phy_ae_kajino13::atmos_phy_ae_kajino13_flag_regeneration = .true.

Definition at line 63 of file scale_atmos_phy_ae_kajino13.F90.

Referenced by atmos_phy_ae_kajino13_setup(), and atmos_phy_ae_kajino13_tendency().

  logical,  public :: atmos_phy_ae_kajino13_flag_regeneration    = .true.

atmos_phy_ae_kajino13_dg_reg

real(rp), public scale_atmos_phy_ae_kajino13::atmos_phy_ae_kajino13_dg_reg = 5.E-7_RP

Definition at line 64 of file scale_atmos_phy_ae_kajino13.F90.

Referenced by atmos_phy_ae_kajino13_setup(), and atmos_phy_ae_kajino13_tendency().

  real(RP), public :: atmos_phy_ae_kajino13_dg_reg = 5.e-7_rp             ! dg of regenerated aerosol (droplet mode 500 nm) [m]

atmos_phy_ae_kajino13_sg_reg

real(rp), public scale_atmos_phy_ae_kajino13::atmos_phy_ae_kajino13_sg_reg = 1.6_RP

Definition at line 65 of file scale_atmos_phy_ae_kajino13.F90.

Referenced by atmos_phy_ae_kajino13_setup(), and atmos_phy_ae_kajino13_tendency().

  real(RP), public :: atmos_phy_ae_kajino13_sg_reg = 1.6_rp               ! sg of regenerated aerosol [-]

atmos_phy_ae_kajino13_logk_aenucl

real(rp), public scale_atmos_phy_ae_kajino13::atmos_phy_ae_kajino13_logk_aenucl = -12.4_RP

Definition at line 66 of file scale_atmos_phy_ae_kajino13.F90.

Referenced by atmos_phy_ae_kajino13_negative_fixer(), and atmos_phy_ae_kajino13_setup().

  real(RP), public :: atmos_phy_ae_kajino13_logk_aenucl = -12.4_rp          !constant coefficient for kinetic nucleation [-]