module_data_mosaic_asect.F

!**********************************************************************************
! This computer software was prepared by Battelle Memorial Institute, hereinafter
! the Contractor, under Contract No. DE-AC05-76RL0 1830 with the Department of
! Energy (DOE). NEITHER THE GOVERNMENT NOR THE CONTRACTOR MAKES ANY WARRANTY,
! EXPRESS OR IMPLIED, OR ASSUMES ANY LIABILITY FOR THE USE OF THIS SOFTWARE.
!
! MOSAIC module: see module_mosaic_driver.F for references and terms of use
!**********************************************************************************
module module_data_mosaic_asect

   implicit none

!-----------------------------------------------------------------------
!
!   The variables in this module provide a means of organizing and accessing
!   aerosol species in the "chem" array by their chemical component,
!   size bin (or mode), "type", and "phase"
!
!   Their purpose is to allow flexible coding of process modules,
!   compared to "hard-coding" using the chem array p_xxx indices
!   (e.g., p_so4_a01, p_so4_a02, ...; p_num_a01, ...)
!
!-----------------------------------------------------------------------
!
!   rce & sg 2004-dec-03 - added phase and type capability,
!        which changed this module almost completely
!
!-----------------------------------------------------------------------
!
!   maxd_atype = maximum allowable number of aerosol types
!   maxd_asize = maximum allowable number of aerosol size bins
!   maxd_acomp = maximum allowable number of chemical components
!        in each aerosol size bin
!   maxd_aphase = maximum allowable number of aerosol phases
!        (gas, cloud, ice, rain, ...)
!
!   ntype_aer = number of aerosol types
!        The aerosol type will allow treatment of an externally mixed
!        aerosol.  The current MOSAIC code has only 1 type, with the implicit
!        assumption of internal mixing.  Eventually, multiple types
!        could treat fresh primary BC/OC, fresh SO4 from nucleation,
!        aged BC/OC/SO4/... mixture, soil dust, sea salt, ...
!
!   nphase_aer = number of aerosol phases
!
!   ai_phase = phase (p) index for interstitial (unactivated) aerosol particles
!   cw_phase = phase (p) index for aerosol particles in cloud water
!   ci_phase = phase (p) index for aerosol particles in cloud ice
!   rn_phase = phase (p) index for aerosol particles in rain
!   sn_phase = phase (p) index for aerosol particles in snow
!   gr_phase = phase (p) index for aerosol particles in graupel
!   [Note:  the value of "xx_phase" will be between 1 and nphase_aer
!        for phases that are active in a simulation.  The others
!        will have non-positive values.]
!
!   nsize_aer(t) = number of aerosol size bins for aerosol type t
!
!   ncomp_aer(t) = number of "regular" chemical components for aerosol type t
!   ncomp_plustracer_aer(t) = number of "regular" plus "tracer"
!        chemical components for aerosol type t
!   [Note:  only "regular" components are used for calculating
!        aerosol physical (mass, volume) and chemical properties.
!        "Tracer" components are optional, and can be used to track source
!        regions, source mechanisms, etc.]
!   [Note:  for aerosol type t, all phases have the same number of size
!        bins, and all size bins have the same number of
!        both regular and tracer components.]
!
!   ntot_mastercomp_aer = number of aerosol chemical components defined
!        in the "master component list".
!   [Note:  each aerosol type will use some but not necessarily all
!        of the components in the "master component list".]
!
!   mastercompptr_aer(c,t) = the position/index/i.d. in the
!       "master component list" for chemical component c of aerosol type t.
!        (1=sulfate, others to be defined by user.)
!
!   massptr_aer(c,s,t,p) = the position/index in the chem array for mixing-
!        ratio for chemical component c, size bin s, type t, and phase p.
!
!   lptr_so4_aer(s,t,p) = the position/index in the chem array for mixing-
!        ratio for sulfate for aerosol size bin s, type t, and phase p
!   (similar lptr's are defined for no3, cl, msa, co3,
!        nh4, na, ca, oin, oc, bc, ...)
!   [Note:  the massptr_aer allow you to loop over all species of
!        an aerosol type.  The lptr_so4_aer, etc., allow you to access
!        a specific chemical component.]
!
!   waterptr_aer(s,t) = the position/index in the chem array for mixing-
!        ratio of aerosol water content for size bin s, type t.
!        [Note:  water content is only carried for the interstitial aerosol
!        phase, so there is no p dimension.]
!
!   hyswptr_aer(s,t) = the position/index in the chem array for mixing-
!        ratio of aerosol "hysteresis water" content for size bin s, type t.
!        This is used to determine if aerosol is in the dry or wet state, when
!        the ambient RH is between the crystallization and deliquescence RH.
!        [Note:  hysteresis water content is only carried for the
!        interstitial aerosol phase, so there is no p dimension.]
!
!   numptr_aer(s,t,p) = the position/index in the chem array for mixing-
!        ratio of particle number for size bin s, type t, and phase p.
!
!   mprognum_aer(s,t,p) - if positive, number mixing-ratio for size s, type t,
!       and phase p will be prognosed.  Otherwise, it is diagnosed using
!        mass mixing-ratio add assumed/prescribed size.
!
!        mixing ratio (moles-water/mole-air) for water
!       associated with aerosol size bin s and type t
!
!
!   mastercompindx_so4_aer = the position/index in the
!       "master component list" for sulfate.
!   (similar lptr's are defined for no3, cl, msa, co3,
!        nh4, na, ca, oin, oc, bc, ...)
!   [Note:  the mastercompindx_xxx_aer are used primarily in
!        initialization routines, and generally aren't needed elsewhere.]
!
!-----------------------------------------------------------------------
!
!   dens_mastercomp_aer(mc) = dry density (g/cm^3) of component mc
!        of the master component list.
!   dens_aer(c,t) = dry density (g/cm^3) of aerosol chemical component
!        c of type t
!   [Note:  dens_aer(c,t) == dens_mastercomp_aer(mastercompptr_aer(c,t))
!        The dens_mastercomp_aer is used in some initialization routines.
!        The dens_aer is used in most other places because of convenience.]
!
!   mw_mastercomp_aer(mc) = molecular weight (g/mole) of component mc
!        of the master component list.
!   mw_aer(c,t) = molecular weight (g/mole) of aerosol chemical component
!        c of type t
!   [Note:  mw_aer(c,t) == mw_mastercomp_aer(mastercompptr_aer(c,t)) ]
!
!   name_mastercomp_aer(mc) = name of component mc of the
!        master component list (e.g., "sulfate", "nitrate", ...).
!   name_aer(c,t) = molecular weight (g/mole) of aerosol chemical component
!        c of type t
!   [Note:  name_aer(c,t) == name_mastercomp_aer(mastercompptr_aer(c,t)) ]
!
!   hygro_mastercomp_aer(mc) = bulk hygroscopicity (--) at dilute conditions
!        (RH near 100%) of component mc of the master component list.
!   hygro_aer(c,t) = bulk hygroscopicity (--) at dilute conditions
!        (RH near 100%) of aerosol chemical component c of type t
!   [For definition of bulk hygroscopicity,
!        see Abdul-Razzak and Ghan, 2004, J Geophys Res, V105, p. 6837-6844.]
!   [Note:  hygro_aer(c,t) == hygro_mastercomp_aer(mastercompptr_aer(c,t)) ]
!
!-----------------------------------------------------------------------
!
!   volumlo_sect(s,t) = 1-particle volume (cm^3) at lower boundary of section m
!   volumhi_sect(s,t) = 1-particle volume (cm^3) at upper boundary of section m
!   volumcen_sect(s,t)= 1-particle volume (cm^3) at "center" of section m
!
!   dlo_sect(s,t) = 1-particle diameter (cm) at lower boundary of section m
!   dhi_sect(s,t) = 1-particle diameter (cm) at upper boundary of section m
!   dcen_sect(s,t) = 1-particle diameter (cm) at "center" section m
!
!   [Note:  the "center" values are defined as follows:
!       volumcen_sect == 0.5*(volumlo_sect + volumhi_sect)
!                     == (pi/6) * (dcen_sect**3) ]
!
!-----------------------------------------------------------------------
!
!   msectional - if positive, each aerosol size bin is a section.
!   if equals 10, use jacobson moving center
!   if equals 20, use tzivion mass-number advection
!        if zero/negative, each size bin is a mode (aitken, accumulation, ...)
!
!   maerosolincw - if positive, both unactivated/interstitial and activated
!       aerosol species are simulated.  if zero/negative, only the
!       unactivated are simulated.  [maerosolincw>0 only when cw_phase>0]
!
!   maerocoag - if positive, aerosol coagulation is done.
!        If zero/negative, it is skipped.
!        (This is not yet implemented in WRF-Chem.)
!
!   maerchem - if positive, aerosol gas-particle condensation/evaporation
!        of inorganic species is done.  If zero/negative, it is skipped.
!        (This is not yet implemented in WRF-Chem.)
!
!   maerchem_boxtest_output - if positive, "boxtest" output is done from
!        the aerchemistry routine.  If zero/negative, it is skipped.
!        (This is not yet implemented in WRF-Chem.)
!
!   maeroptical - if positive, aerosol optical properties are calculated.
!        If zero/negative, it is skipped.
!        (This is not yet implemented in WRF-Chem.)
!
!-----------------------------------------------------------------------

   integer, parameter :: maxd_atype = 1
   integer, parameter :: maxd_asize = 8
!        integer, parameter :: maxd_acomp = 120 ! for additional SOA species changed by Manish Shrivastava on 01/25/10
   integer, parameter :: maxd_acomp = 400 ! for additional SOA species changed by Manish Shrivastava on 01/25/10
   integer, parameter :: maxd_aphase = 2

   integer, save :: ai_phase = -999888777
   integer, save :: cw_phase = -999888777
   integer, save :: ci_phase = -999888777
   integer, save :: rn_phase = -999888777
   integer, save :: sn_phase = -999888777
   integer, save :: gr_phase = -999888777

   integer, save :: ntype_aer = 0 ! number of types
   integer, save :: ntot_mastercomp_aer = 0 ! number of master components
   integer, save :: nphase_aer = 0 ! number of phases

   integer, save :: &
      nsize_aer(maxd_atype), & ! number of size bins
      ncomp_aer(maxd_atype), & ! number of chemical components
      ncomp_plustracer_aer(maxd_atype), &
      mastercompptr_aer(maxd_acomp, maxd_atype), & !  mastercomp index
      massptr_aer(maxd_acomp, maxd_asize, maxd_atype, maxd_aphase), &
      ! index for mixing ratio
      waterptr_aer(maxd_asize, maxd_atype), & ! index for aerosol water
      hyswptr_aer(maxd_asize, maxd_atype), &
      numptr_aer(maxd_asize, maxd_atype, maxd_aphase), &
      ! index for the number mixing ratio
      mprognum_aer(maxd_asize, maxd_atype, maxd_aphase)

!   these indices give the location in the "mastercomp list" of
!   the different aerosol chemical (or tracer) components
   integer, save :: mastercompindx_so4_aer = -999888777
   integer, save :: mastercompindx_no3_aer = -999888777
   integer, save :: mastercompindx_cl_aer = -999888777
   integer, save :: mastercompindx_msa_aer = -999888777
   integer, save :: mastercompindx_co3_aer = -999888777
   integer, save :: mastercompindx_nh4_aer = -999888777
   integer, save :: mastercompindx_na_aer = -999888777
   integer, save :: mastercompindx_ca_aer = -999888777
   integer, save :: mastercompindx_oin_aer = -999888777
   integer, save :: mastercompindx_oc_aer = -999888777
   integer, save :: mastercompindx_bc_aer = -999888777
   integer, save :: mastercompindx_pcg1_b_c_aer = -999888777
   integer, save :: mastercompindx_pcg2_b_c_aer = -999888777
   integer, save :: mastercompindx_pcg3_b_c_aer = -999888777
   integer, save :: mastercompindx_pcg4_b_c_aer = -999888777
   integer, save :: mastercompindx_pcg5_b_c_aer = -999888777
   integer, save :: mastercompindx_pcg6_b_c_aer = -999888777
   integer, save :: mastercompindx_pcg7_b_c_aer = -999888777
   integer, save :: mastercompindx_pcg8_b_c_aer = -999888777
   integer, save :: mastercompindx_pcg9_b_c_aer = -999888777
   integer, save :: mastercompindx_pcg1_b_o_aer = -999888777
   integer, save :: mastercompindx_pcg2_b_o_aer = -999888777
   integer, save :: mastercompindx_pcg3_b_o_aer = -999888777
   integer, save :: mastercompindx_pcg4_b_o_aer = -999888777
   integer, save :: mastercompindx_pcg5_b_o_aer = -999888777
   integer, save :: mastercompindx_pcg6_b_o_aer = -999888777
   integer, save :: mastercompindx_pcg7_b_o_aer = -999888777
   integer, save :: mastercompindx_pcg8_b_o_aer = -999888777
   integer, save :: mastercompindx_pcg9_b_o_aer = -999888777
   integer, save :: mastercompindx_opcg1_b_c_aer = -999888777
   integer, save :: mastercompindx_opcg2_b_c_aer = -999888777
   integer, save :: mastercompindx_opcg3_b_c_aer = -999888777
   integer, save :: mastercompindx_opcg4_b_c_aer = -999888777
   integer, save :: mastercompindx_opcg5_b_c_aer = -999888777
   integer, save :: mastercompindx_opcg6_b_c_aer = -999888777
   integer, save :: mastercompindx_opcg7_b_c_aer = -999888777
   integer, save :: mastercompindx_opcg8_b_c_aer = -999888777
   integer, save :: mastercompindx_opcg1_b_o_aer = -999888777
   integer, save :: mastercompindx_opcg2_b_o_aer = -999888777
   integer, save :: mastercompindx_opcg3_b_o_aer = -999888777
   integer, save :: mastercompindx_opcg4_b_o_aer = -999888777
   integer, save :: mastercompindx_opcg5_b_o_aer = -999888777
   integer, save :: mastercompindx_opcg6_b_o_aer = -999888777
   integer, save :: mastercompindx_opcg7_b_o_aer = -999888777
   integer, save :: mastercompindx_opcg8_b_o_aer = -999888777
   integer, save :: mastercompindx_pcg1_f_c_aer = -999888777
   integer, save :: mastercompindx_pcg2_f_c_aer = -999888777
   integer, save :: mastercompindx_pcg3_f_c_aer = -999888777
   integer, save :: mastercompindx_pcg4_f_c_aer = -999888777
   integer, save :: mastercompindx_pcg5_f_c_aer = -999888777
   integer, save :: mastercompindx_pcg6_f_c_aer = -999888777
   integer, save :: mastercompindx_pcg7_f_c_aer = -999888777
   integer, save :: mastercompindx_pcg8_f_c_aer = -999888777
   integer, save :: mastercompindx_pcg9_f_c_aer = -999888777
   integer, save :: mastercompindx_pcg1_f_o_aer = -999888777
   integer, save :: mastercompindx_pcg2_f_o_aer = -999888777
   integer, save :: mastercompindx_pcg3_f_o_aer = -999888777
   integer, save :: mastercompindx_pcg4_f_o_aer = -999888777
   integer, save :: mastercompindx_pcg5_f_o_aer = -999888777
   integer, save :: mastercompindx_pcg6_f_o_aer = -999888777
   integer, save :: mastercompindx_pcg7_f_o_aer = -999888777
   integer, save :: mastercompindx_pcg8_f_o_aer = -999888777
   integer, save :: mastercompindx_pcg9_f_o_aer = -999888777
   integer, save :: mastercompindx_opcg1_f_c_aer = -999888777
   integer, save :: mastercompindx_opcg2_f_c_aer = -999888777
   integer, save :: mastercompindx_opcg3_f_c_aer = -999888777
   integer, save :: mastercompindx_opcg4_f_c_aer = -999888777
   integer, save :: mastercompindx_opcg5_f_c_aer = -999888777
   integer, save :: mastercompindx_opcg6_f_c_aer = -999888777
   integer, save :: mastercompindx_opcg7_f_c_aer = -999888777
   integer, save :: mastercompindx_opcg8_f_c_aer = -999888777
   integer, save :: mastercompindx_opcg1_f_o_aer = -999888777
   integer, save :: mastercompindx_opcg2_f_o_aer = -999888777
   integer, save :: mastercompindx_opcg3_f_o_aer = -999888777
   integer, save :: mastercompindx_opcg4_f_o_aer = -999888777
   integer, save :: mastercompindx_opcg5_f_o_aer = -999888777
   integer, save :: mastercompindx_opcg6_f_o_aer = -999888777
   integer, save :: mastercompindx_opcg7_f_o_aer = -999888777
   integer, save :: mastercompindx_opcg8_f_o_aer = -999888777
   integer, save :: mastercompindx_smpa_aer = -999888777
   integer, save :: mastercompindx_smpbb_aer = -999888777

   integer, save :: mastercompindx_glysoa_r1_aer = -999888777
   integer, save :: mastercompindx_glysoa_r2_aer = -999888777
   integer, save :: mastercompindx_glysoa_oh_aer = -999888777
   integer, save :: mastercompindx_glysoa_nh4_aer = -999888777
   integer, save :: mastercompindx_glysoa_sfc_aer = -999888777

   integer, save :: mastercompindx_ant1_c_aer = -999888777
   integer, save :: mastercompindx_ant2_c_aer = -999888777
   integer, save :: mastercompindx_ant3_c_aer = -999888777
   integer, save :: mastercompindx_ant4_c_aer = -999888777
   integer, save :: mastercompindx_ant1_o_aer = -999888777
   integer, save :: mastercompindx_ant2_o_aer = -999888777
   integer, save :: mastercompindx_ant3_o_aer = -999888777
   integer, save :: mastercompindx_ant4_o_aer = -999888777
   integer, save :: mastercompindx_biog1_c_aer = -999888777
   integer, save :: mastercompindx_biog2_c_aer = -999888777
   integer, save :: mastercompindx_biog3_c_aer = -999888777
   integer, save :: mastercompindx_biog4_c_aer = -999888777
   integer, save :: mastercompindx_biog1_o_aer = -999888777
   integer, save :: mastercompindx_biog2_o_aer = -999888777
   integer, save :: mastercompindx_biog3_o_aer = -999888777
   integer, save :: mastercompindx_biog4_o_aer = -999888777

   integer, save :: mastercompindx_asoaX_aer = -999888777
   integer, save :: mastercompindx_asoa1_aer = -999888777
   integer, save :: mastercompindx_asoa2_aer = -999888777
   integer, save :: mastercompindx_asoa3_aer = -999888777
   integer, save :: mastercompindx_asoa4_aer = -999888777
   integer, save :: mastercompindx_bsoaX_aer = -999888777
   integer, save :: mastercompindx_bsoa1_aer = -999888777
   integer, save :: mastercompindx_bsoa2_aer = -999888777
   integer, save :: mastercompindx_bsoa3_aer = -999888777
   integer, save :: mastercompindx_bsoa4_aer = -999888777

   real, save :: &
      dens_aer(maxd_acomp, maxd_atype), &
      dens_mastercomp_aer(maxd_acomp), &
      mw_mastercomp_aer(maxd_acomp), &
      mw_aer(maxd_acomp, maxd_atype), &
      hygro_mastercomp_aer(maxd_acomp), &
      hygro_aer(maxd_acomp, maxd_atype)

   real, save :: &
      volumcen_sect(maxd_asize, maxd_atype), &
      volumlo_sect(maxd_asize, maxd_atype), &
      volumhi_sect(maxd_asize, maxd_atype), &
      dcen_sect(maxd_asize, maxd_atype), &
      dlo_sect(maxd_asize, maxd_atype), &
      dhi_sect(maxd_asize, maxd_atype), &
      sigmag_aer(maxd_asize, maxd_atype)

   character*20, save :: &
      name_mastercomp_aer(maxd_acomp), &
      name_aer(maxd_acomp, maxd_atype)

   integer, save :: &
      lptr_so4_aer(maxd_asize, maxd_atype, maxd_aphase), &
      lptr_msa_aer(maxd_asize, maxd_atype, maxd_aphase), &
      lptr_no3_aer(maxd_asize, maxd_atype, maxd_aphase), &
      lptr_cl_aer(maxd_asize, maxd_atype, maxd_aphase), &
      lptr_co3_aer(maxd_asize, maxd_atype, maxd_aphase), &
      lptr_nh4_aer(maxd_asize, maxd_atype, maxd_aphase), &
      lptr_na_aer(maxd_asize, maxd_atype, maxd_aphase), &
      lptr_ca_aer(maxd_asize, maxd_atype, maxd_aphase), &
      lptr_oin_aer(maxd_asize, maxd_atype, maxd_aphase), &
      lptr_oc_aer(maxd_asize, maxd_atype, maxd_aphase), &
      lptr_bc_aer(maxd_asize, maxd_atype, maxd_aphase), &
      lptr_pcg1_b_c_aer(maxd_asize, maxd_atype, maxd_aphase), &
      lptr_pcg2_b_c_aer(maxd_asize, maxd_atype, maxd_aphase), &
      lptr_pcg3_b_c_aer(maxd_asize, maxd_atype, maxd_aphase), &
      lptr_pcg4_b_c_aer(maxd_asize, maxd_atype, maxd_aphase), &
      lptr_pcg5_b_c_aer(maxd_asize, maxd_atype, maxd_aphase), &
      lptr_pcg6_b_c_aer(maxd_asize, maxd_atype, maxd_aphase), &
      lptr_pcg7_b_c_aer(maxd_asize, maxd_atype, maxd_aphase), &
      lptr_pcg8_b_c_aer(maxd_asize, maxd_atype, maxd_aphase), &
      lptr_pcg9_b_c_aer(maxd_asize, maxd_atype, maxd_aphase), &
      lptr_pcg1_b_o_aer(maxd_asize, maxd_atype, maxd_aphase), &
      lptr_pcg2_b_o_aer(maxd_asize, maxd_atype, maxd_aphase), &
      lptr_pcg3_b_o_aer(maxd_asize, maxd_atype, maxd_aphase), &
      lptr_pcg4_b_o_aer(maxd_asize, maxd_atype, maxd_aphase), &
      lptr_pcg5_b_o_aer(maxd_asize, maxd_atype, maxd_aphase), &
      lptr_pcg6_b_o_aer(maxd_asize, maxd_atype, maxd_aphase), &
      lptr_pcg7_b_o_aer(maxd_asize, maxd_atype, maxd_aphase), &
      lptr_pcg8_b_o_aer(maxd_asize, maxd_atype, maxd_aphase), &
      lptr_pcg9_b_o_aer(maxd_asize, maxd_atype, maxd_aphase), &
      lptr_opcg1_b_c_aer(maxd_asize, maxd_atype, maxd_aphase), &
      lptr_opcg2_b_c_aer(maxd_asize, maxd_atype, maxd_aphase), &
      lptr_opcg3_b_c_aer(maxd_asize, maxd_atype, maxd_aphase), &
      lptr_opcg4_b_c_aer(maxd_asize, maxd_atype, maxd_aphase), &
      lptr_opcg5_b_c_aer(maxd_asize, maxd_atype, maxd_aphase), &
      lptr_opcg6_b_c_aer(maxd_asize, maxd_atype, maxd_aphase), &
      lptr_opcg7_b_c_aer(maxd_asize, maxd_atype, maxd_aphase), &
      lptr_opcg8_b_c_aer(maxd_asize, maxd_atype, maxd_aphase), &
      lptr_opcg1_b_o_aer(maxd_asize, maxd_atype, maxd_aphase), &
      lptr_opcg2_b_o_aer(maxd_asize, maxd_atype, maxd_aphase), &
      lptr_opcg3_b_o_aer(maxd_asize, maxd_atype, maxd_aphase), &
      lptr_opcg4_b_o_aer(maxd_asize, maxd_atype, maxd_aphase), &
      lptr_opcg5_b_o_aer(maxd_asize, maxd_atype, maxd_aphase), &
      lptr_opcg6_b_o_aer(maxd_asize, maxd_atype, maxd_aphase), &
      lptr_opcg7_b_o_aer(maxd_asize, maxd_atype, maxd_aphase), &
      lptr_opcg8_b_o_aer(maxd_asize, maxd_atype, maxd_aphase), &
      lptr_pcg1_f_c_aer(maxd_asize, maxd_atype, maxd_aphase), &
      lptr_pcg2_f_c_aer(maxd_asize, maxd_atype, maxd_aphase), &
      lptr_pcg3_f_c_aer(maxd_asize, maxd_atype, maxd_aphase), &
      lptr_pcg4_f_c_aer(maxd_asize, maxd_atype, maxd_aphase), &
      lptr_pcg5_f_c_aer(maxd_asize, maxd_atype, maxd_aphase), &
      lptr_pcg6_f_c_aer(maxd_asize, maxd_atype, maxd_aphase), &
      lptr_pcg7_f_c_aer(maxd_asize, maxd_atype, maxd_aphase), &
      lptr_pcg8_f_c_aer(maxd_asize, maxd_atype, maxd_aphase), &
      lptr_pcg9_f_c_aer(maxd_asize, maxd_atype, maxd_aphase), &
      lptr_pcg1_f_o_aer(maxd_asize, maxd_atype, maxd_aphase), &
      lptr_pcg2_f_o_aer(maxd_asize, maxd_atype, maxd_aphase), &
      lptr_pcg3_f_o_aer(maxd_asize, maxd_atype, maxd_aphase), &
      lptr_pcg4_f_o_aer(maxd_asize, maxd_atype, maxd_aphase), &
      lptr_pcg5_f_o_aer(maxd_asize, maxd_atype, maxd_aphase), &
      lptr_pcg6_f_o_aer(maxd_asize, maxd_atype, maxd_aphase), &
      lptr_pcg7_f_o_aer(maxd_asize, maxd_atype, maxd_aphase), &
      lptr_pcg8_f_o_aer(maxd_asize, maxd_atype, maxd_aphase), &
      lptr_pcg9_f_o_aer(maxd_asize, maxd_atype, maxd_aphase), &
      lptr_opcg1_f_c_aer(maxd_asize, maxd_atype, maxd_aphase), &
      lptr_opcg2_f_c_aer(maxd_asize, maxd_atype, maxd_aphase), &
      lptr_opcg3_f_c_aer(maxd_asize, maxd_atype, maxd_aphase), &
      lptr_opcg4_f_c_aer(maxd_asize, maxd_atype, maxd_aphase), &
      lptr_opcg5_f_c_aer(maxd_asize, maxd_atype, maxd_aphase), &
      lptr_opcg6_f_c_aer(maxd_asize, maxd_atype, maxd_aphase), &
      lptr_opcg7_f_c_aer(maxd_asize, maxd_atype, maxd_aphase), &
      lptr_opcg8_f_c_aer(maxd_asize, maxd_atype, maxd_aphase), &
      lptr_opcg1_f_o_aer(maxd_asize, maxd_atype, maxd_aphase), &
      lptr_opcg2_f_o_aer(maxd_asize, maxd_atype, maxd_aphase), &
      lptr_opcg3_f_o_aer(maxd_asize, maxd_atype, maxd_aphase), &
      lptr_opcg4_f_o_aer(maxd_asize, maxd_atype, maxd_aphase), &
      lptr_opcg5_f_o_aer(maxd_asize, maxd_atype, maxd_aphase), &
      lptr_opcg6_f_o_aer(maxd_asize, maxd_atype, maxd_aphase), &
      lptr_opcg7_f_o_aer(maxd_asize, maxd_atype, maxd_aphase), &
      lptr_opcg8_f_o_aer(maxd_asize, maxd_atype, maxd_aphase), &
      lptr_smpa_aer(maxd_asize, maxd_atype, maxd_aphase), &
      lptr_smpbb_aer(maxd_asize, maxd_atype, maxd_aphase), &
      lptr_glysoa_r1_aer(maxd_asize, maxd_atype, maxd_aphase), &
      lptr_glysoa_r2_aer(maxd_asize, maxd_atype, maxd_aphase), &
      lptr_glysoa_oh_aer(maxd_asize, maxd_atype, maxd_aphase), &
      lptr_glysoa_nh4_aer(maxd_asize, maxd_atype, maxd_aphase), &
      lptr_glysoa_sfc_aer(maxd_asize, maxd_atype, maxd_aphase), &
      lptr_ant1_c_aer(maxd_asize, maxd_atype, maxd_aphase), &
      lptr_ant2_c_aer(maxd_asize, maxd_atype, maxd_aphase), &
      lptr_ant3_c_aer(maxd_asize, maxd_atype, maxd_aphase), &
      lptr_ant4_c_aer(maxd_asize, maxd_atype, maxd_aphase), &
      lptr_ant1_o_aer(maxd_asize, maxd_atype, maxd_aphase), &
      lptr_ant2_o_aer(maxd_asize, maxd_atype, maxd_aphase), &
      lptr_ant3_o_aer(maxd_asize, maxd_atype, maxd_aphase), &
      lptr_ant4_o_aer(maxd_asize, maxd_atype, maxd_aphase), &
      lptr_biog1_c_aer(maxd_asize, maxd_atype, maxd_aphase), &
      lptr_biog2_c_aer(maxd_asize, maxd_atype, maxd_aphase), &
      lptr_biog3_c_aer(maxd_asize, maxd_atype, maxd_aphase), &
      lptr_biog4_c_aer(maxd_asize, maxd_atype, maxd_aphase), &
      lptr_biog1_o_aer(maxd_asize, maxd_atype, maxd_aphase), &
      lptr_biog2_o_aer(maxd_asize, maxd_atype, maxd_aphase), &
      lptr_biog3_o_aer(maxd_asize, maxd_atype, maxd_aphase), &
      lptr_biog4_o_aer(maxd_asize, maxd_atype, maxd_aphase), &
      lptr_asoaX_aer(maxd_asize, maxd_atype, maxd_aphase), &
      lptr_asoa1_aer(maxd_asize, maxd_atype, maxd_aphase), &
      lptr_asoa2_aer(maxd_asize, maxd_atype, maxd_aphase), &
      lptr_asoa3_aer(maxd_asize, maxd_atype, maxd_aphase), &
      lptr_asoa4_aer(maxd_asize, maxd_atype, maxd_aphase), &
      lptr_bsoaX_aer(maxd_asize, maxd_atype, maxd_aphase), &
      lptr_bsoa1_aer(maxd_asize, maxd_atype, maxd_aphase), &
      lptr_bsoa2_aer(maxd_asize, maxd_atype, maxd_aphase), &
      lptr_bsoa3_aer(maxd_asize, maxd_atype, maxd_aphase), &
      lptr_bsoa4_aer(maxd_asize, maxd_atype, maxd_aphase)

! rce 11-sep-2004 - eliminated all of the "..._wrfch" pointers
!    so now there is only one set of pointers ("..._amode")
! sg/rce nov-2004 - totally new pointer system  - "..._aer"

!   molecular weights (g/mol)
   real, parameter :: mw_so4_aer = 96.066
   real, parameter :: mw_no3_aer = 62.007
   real, parameter :: mw_cl_aer = 35.450
   real, parameter :: mw_msa_aer = 96.109
   real, parameter :: mw_co3_aer = 60.007
   real, parameter :: mw_nh4_aer = 18.042
   real, parameter :: mw_na_aer = 22.990
   real, parameter :: mw_ca_aer = 40.080
   real, parameter :: mw_oin_aer = 1.0
   real, parameter :: mw_oc_aer = 250.0
   real, parameter :: mw_bc_aer = 1.0
   real, parameter :: mw_water_aer = 18.016
   real, parameter :: mw_pcg1_b_c_aer = 250.0
   real, parameter :: mw_pcg2_b_c_aer = 250.0
   real, parameter :: mw_pcg3_b_c_aer = 250.0
   real, parameter :: mw_pcg4_b_c_aer = 250.0
   real, parameter :: mw_pcg5_b_c_aer = 250.0
   real, parameter :: mw_pcg6_b_c_aer = 250.0
   real, parameter :: mw_pcg7_b_c_aer = 250.0
   real, parameter :: mw_pcg8_b_c_aer = 250.0
   real, parameter :: mw_pcg9_b_c_aer = 250.0
   real, parameter :: mw_pcg1_b_o_aer = 250.0
   real, parameter :: mw_pcg2_b_o_aer = 250.0
   real, parameter :: mw_pcg3_b_o_aer = 250.0
   real, parameter :: mw_pcg4_b_o_aer = 250.0
   real, parameter :: mw_pcg5_b_o_aer = 250.0
   real, parameter :: mw_pcg6_b_o_aer = 250.0
   real, parameter :: mw_pcg7_b_o_aer = 250.0
   real, parameter :: mw_pcg8_b_o_aer = 250.0
   real, parameter :: mw_pcg9_b_o_aer = 250.0
   real, parameter :: mw_opcg1_b_c_aer = 250.0
   real, parameter :: mw_opcg2_b_c_aer = 250.0
   real, parameter :: mw_opcg3_b_c_aer = 250.0
   real, parameter :: mw_opcg4_b_c_aer = 250.0
   real, parameter :: mw_opcg5_b_c_aer = 250.0
   real, parameter :: mw_opcg6_b_c_aer = 250.0
   real, parameter :: mw_opcg7_b_c_aer = 250.0
   real, parameter :: mw_opcg8_b_c_aer = 250.0
   real, parameter :: mw_opcg1_b_o_aer = 250.0
   real, parameter :: mw_opcg2_b_o_aer = 250.0
   real, parameter :: mw_opcg3_b_o_aer = 250.0
   real, parameter :: mw_opcg4_b_o_aer = 250.0
   real, parameter :: mw_opcg5_b_o_aer = 250.0
   real, parameter :: mw_opcg6_b_o_aer = 250.0
   real, parameter :: mw_opcg7_b_o_aer = 250.0
   real, parameter :: mw_opcg8_b_o_aer = 250.0
   real, parameter :: mw_pcg1_f_c_aer = 250.0
   real, parameter :: mw_pcg2_f_c_aer = 250.0
   real, parameter :: mw_pcg3_f_c_aer = 250.0
   real, parameter :: mw_pcg4_f_c_aer = 250.0
   real, parameter :: mw_pcg5_f_c_aer = 250.0
   real, parameter :: mw_pcg6_f_c_aer = 250.0
   real, parameter :: mw_pcg7_f_c_aer = 250.0
   real, parameter :: mw_pcg8_f_c_aer = 250.0
   real, parameter :: mw_pcg9_f_c_aer = 250.0
   real, parameter :: mw_pcg1_f_o_aer = 250.0
   real, parameter :: mw_pcg2_f_o_aer = 250.0
   real, parameter :: mw_pcg3_f_o_aer = 250.0
   real, parameter :: mw_pcg4_f_o_aer = 250.0
   real, parameter :: mw_pcg5_f_o_aer = 250.0
   real, parameter :: mw_pcg6_f_o_aer = 250.0
   real, parameter :: mw_pcg7_f_o_aer = 250.0
   real, parameter :: mw_pcg8_f_o_aer = 250.0
   real, parameter :: mw_pcg9_f_o_aer = 250.0
   real, parameter :: mw_opcg1_f_c_aer = 250.0
   real, parameter :: mw_opcg2_f_c_aer = 250.0
   real, parameter :: mw_opcg3_f_c_aer = 250.0
   real, parameter :: mw_opcg4_f_c_aer = 250.0
   real, parameter :: mw_opcg5_f_c_aer = 250.0
   real, parameter :: mw_opcg6_f_c_aer = 250.0
   real, parameter :: mw_opcg7_f_c_aer = 250.0
   real, parameter :: mw_opcg8_f_c_aer = 250.0
   real, parameter :: mw_opcg1_f_o_aer = 250.0
   real, parameter :: mw_opcg2_f_o_aer = 250.0
   real, parameter :: mw_opcg3_f_o_aer = 250.0
   real, parameter :: mw_opcg4_f_o_aer = 250.0
   real, parameter :: mw_opcg5_f_o_aer = 250.0
   real, parameter :: mw_opcg6_f_o_aer = 250.0
   real, parameter :: mw_opcg7_f_o_aer = 250.0
   real, parameter :: mw_opcg8_f_o_aer = 250.0
   real, parameter :: mw_smpa_aer = 250
   real, parameter :: mw_smpbb_aer = 250

   real, parameter :: mw_glysoa_r1_aer = 250.0
   real, parameter :: mw_glysoa_r2_aer = 250.0
   real, parameter :: mw_glysoa_oh_aer = 250.0
   real, parameter :: mw_glysoa_nh4_aer = 250.0
   real, parameter :: mw_glysoa_sfc_aer = 250.0

   real, parameter :: mw_ant1_c_aer = 250
   real, parameter :: mw_ant2_c_aer = 250
   real, parameter :: mw_ant3_c_aer = 250
   real, parameter :: mw_ant4_c_aer = 250
   real, parameter :: mw_ant1_o_aer = 250
   real, parameter :: mw_ant2_o_aer = 250
   real, parameter :: mw_ant3_o_aer = 250
   real, parameter :: mw_ant4_o_aer = 250
   real, parameter :: mw_biog1_c_aer = 250.0
   real, parameter :: mw_biog2_c_aer = 250.0
   real, parameter :: mw_biog3_c_aer = 250.0
   real, parameter :: mw_biog4_c_aer = 250.0
   real, parameter :: mw_biog1_o_aer = 250.0
   real, parameter :: mw_biog2_o_aer = 250.0
   real, parameter :: mw_biog3_o_aer = 250.0
   real, parameter :: mw_biog4_o_aer = 250.0

   real, parameter :: mw_asoaX_aer = 250.0
   real, parameter :: mw_asoa1_aer = 250.0
   real, parameter :: mw_asoa2_aer = 250.0
   real, parameter :: mw_asoa3_aer = 250.0
   real, parameter :: mw_asoa4_aer = 250.0
   real, parameter :: mw_bsoaX_aer = 250.0
   real, parameter :: mw_bsoa1_aer = 250.0
   real, parameter :: mw_bsoa2_aer = 250.0
   real, parameter :: mw_bsoa3_aer = 250.0
   real, parameter :: mw_bsoa4_aer = 250.0

!   dry densities (g/cm3)
   real, parameter :: dens_so4_aer = 1.80
   real, parameter :: dens_no3_aer = 1.80
   real, parameter :: dens_cl_aer = 2.20
   real, parameter :: dens_msa_aer = 1.80
   real, parameter :: dens_co3_aer = 2.60
   real, parameter :: dens_nh4_aer = 1.80
   real, parameter :: dens_na_aer = 2.20
   real, parameter :: dens_ca_aer = 2.60
   real, parameter :: dens_oin_aer = 2.60
   real, parameter :: dens_oc_aer = 1.00
   real, parameter :: dens_bc_aer = 1.70
   real, parameter :: dens_pcg1_b_c_aer = 1.0
   real, parameter :: dens_pcg2_b_c_aer = 1.0
   real, parameter :: dens_pcg3_b_c_aer = 1.0
   real, parameter :: dens_pcg4_b_c_aer = 1.0
   real, parameter :: dens_pcg5_b_c_aer = 1.0
   real, parameter :: dens_pcg6_b_c_aer = 1.0
   real, parameter :: dens_pcg7_b_c_aer = 1.0
   real, parameter :: dens_pcg8_b_c_aer = 1.0
   real, parameter :: dens_pcg9_b_c_aer = 1.0
   real, parameter :: dens_pcg1_b_o_aer = 1.0
   real, parameter :: dens_pcg2_b_o_aer = 1.0
   real, parameter :: dens_pcg3_b_o_aer = 1.0
   real, parameter :: dens_pcg4_b_o_aer = 1.0
   real, parameter :: dens_pcg5_b_o_aer = 1.0
   real, parameter :: dens_pcg6_b_o_aer = 1.0
   real, parameter :: dens_pcg7_b_o_aer = 1.0
   real, parameter :: dens_pcg8_b_o_aer = 1.0
   real, parameter :: dens_pcg9_b_o_aer = 1.0
   real, parameter :: dens_opcg1_b_c_aer = 1.0
   real, parameter :: dens_opcg2_b_c_aer = 1.0
   real, parameter :: dens_opcg3_b_c_aer = 1.0
   real, parameter :: dens_opcg4_b_c_aer = 1.0
   real, parameter :: dens_opcg5_b_c_aer = 1.0
   real, parameter :: dens_opcg6_b_c_aer = 1.0
   real, parameter :: dens_opcg7_b_c_aer = 1.0
   real, parameter :: dens_opcg8_b_c_aer = 1.0
   real, parameter :: dens_opcg1_b_o_aer = 1.0
   real, parameter :: dens_opcg2_b_o_aer = 1.0
   real, parameter :: dens_opcg3_b_o_aer = 1.0
   real, parameter :: dens_opcg4_b_o_aer = 1.0
   real, parameter :: dens_opcg5_b_o_aer = 1.0
   real, parameter :: dens_opcg6_b_o_aer = 1.0
   real, parameter :: dens_opcg7_b_o_aer = 1.0
   real, parameter :: dens_opcg8_b_o_aer = 1.0
   real, parameter :: dens_pcg1_f_c_aer = 1.0
   real, parameter :: dens_pcg2_f_c_aer = 1.0
   real, parameter :: dens_pcg3_f_c_aer = 1.0
   real, parameter :: dens_pcg4_f_c_aer = 1.0
   real, parameter :: dens_pcg5_f_c_aer = 1.0
   real, parameter :: dens_pcg6_f_c_aer = 1.0
   real, parameter :: dens_pcg7_f_c_aer = 1.0
   real, parameter :: dens_pcg8_f_c_aer = 1.0
   real, parameter :: dens_pcg9_f_c_aer = 1.0
   real, parameter :: dens_pcg1_f_o_aer = 1.0
   real, parameter :: dens_pcg2_f_o_aer = 1.0
   real, parameter :: dens_pcg3_f_o_aer = 1.0
   real, parameter :: dens_pcg4_f_o_aer = 1.0
   real, parameter :: dens_pcg5_f_o_aer = 1.0
   real, parameter :: dens_pcg6_f_o_aer = 1.0
   real, parameter :: dens_pcg7_f_o_aer = 1.0
   real, parameter :: dens_pcg8_f_o_aer = 1.0
   real, parameter :: dens_pcg9_f_o_aer = 1.0
   real, parameter :: dens_opcg1_f_c_aer = 1.0
   real, parameter :: dens_opcg2_f_c_aer = 1.0
   real, parameter :: dens_opcg3_f_c_aer = 1.0
   real, parameter :: dens_opcg4_f_c_aer = 1.0
   real, parameter :: dens_opcg5_f_c_aer = 1.0
   real, parameter :: dens_opcg6_f_c_aer = 1.0
   real, parameter :: dens_opcg7_f_c_aer = 1.0
   real, parameter :: dens_opcg8_f_c_aer = 1.0
   real, parameter :: dens_opcg1_f_o_aer = 1.0
   real, parameter :: dens_opcg2_f_o_aer = 1.0
   real, parameter :: dens_opcg3_f_o_aer = 1.0
   real, parameter :: dens_opcg4_f_o_aer = 1.0
   real, parameter :: dens_opcg5_f_o_aer = 1.0
   real, parameter :: dens_opcg6_f_o_aer = 1.0
   real, parameter :: dens_opcg7_f_o_aer = 1.0
   real, parameter :: dens_opcg8_f_o_aer = 1.0
   real, parameter :: dens_smpa_aer = 1.0
   real, parameter :: dens_smpbb_aer = 1.0

   real, parameter :: dens_glysoa_r1_aer = 1.0
   real, parameter :: dens_glysoa_r2_aer = 1.0
   real, parameter :: dens_glysoa_oh_aer = 1.0
   real, parameter :: dens_glysoa_nh4_aer = 1.0
   real, parameter :: dens_glysoa_sfc_aer = 1.0

   real, parameter :: dens_ant1_c_aer = 1.0
   real, parameter :: dens_ant2_c_aer = 1.0
   real, parameter :: dens_ant3_c_aer = 1.0
   real, parameter :: dens_ant4_c_aer = 1.0
   real, parameter :: dens_ant1_o_aer = 1.0
   real, parameter :: dens_ant2_o_aer = 1.0
   real, parameter :: dens_ant3_o_aer = 1.0
   real, parameter :: dens_ant4_o_aer = 1.0
   real, parameter :: dens_biog1_c_aer = 1.0
   real, parameter :: dens_biog2_c_aer = 1.0
   real, parameter :: dens_biog3_c_aer = 1.0
   real, parameter :: dens_biog4_c_aer = 1.0
   real, parameter :: dens_biog1_o_aer = 1.0
   real, parameter :: dens_biog2_o_aer = 1.0
   real, parameter :: dens_biog3_o_aer = 1.0
   real, parameter :: dens_biog4_o_aer = 1.0

! dens based on Murphy and Pandis, ES&T, 2009
   real, parameter :: dens_asoaX_aer = 1.5
   real, parameter :: dens_asoa1_aer = 1.5
   real, parameter :: dens_asoa2_aer = 1.5
   real, parameter :: dens_asoa3_aer = 1.5
   real, parameter :: dens_asoa4_aer = 1.5
   real, parameter :: dens_bsoaX_aer = 1.5
   real, parameter :: dens_bsoa1_aer = 1.5
   real, parameter :: dens_bsoa2_aer = 1.5
   real, parameter :: dens_bsoa3_aer = 1.5
   real, parameter :: dens_bsoa4_aer = 1.5

!   water density (g/cm3)
!        real, parameter :: dens_water_asize  = 1.0
   real, parameter :: dens_water_aer = 1.0

!   hygroscopicities (dimensionless)
   real, parameter :: hygro_so4_aer = 0.5
   real, parameter :: hygro_no3_aer = 0.5
   real, parameter :: hygro_ca_aer = 0.1
   real, parameter :: hygro_co3_aer = 0.1
   real, parameter :: hygro_nh4_aer = 0.5
   real, parameter :: hygro_msa_aer = 0.58
   real, parameter :: hygro_cl_aer = 1.16
   real, parameter :: hygro_na_aer = 1.16
   real, parameter :: hygro_oin_aer = 0.14
   real, parameter :: hygro_oc_aer = 0.14
   real, parameter :: hygro_bc_aer = 1.e-6
   real, parameter :: hygro_pcg1_b_c_aer = 0.14
   real, parameter :: hygro_pcg2_b_c_aer = 0.14
   real, parameter :: hygro_pcg3_b_c_aer = 0.14
   real, parameter :: hygro_pcg4_b_c_aer = 0.14
   real, parameter :: hygro_pcg5_b_c_aer = 0.14
   real, parameter :: hygro_pcg6_b_c_aer = 0.14
   real, parameter :: hygro_pcg7_b_c_aer = 0.14
   real, parameter :: hygro_pcg8_b_c_aer = 0.14
   real, parameter :: hygro_pcg9_b_c_aer = 0.14
   real, parameter :: hygro_pcg1_b_o_aer = 0.14
   real, parameter :: hygro_pcg2_b_o_aer = 0.14
   real, parameter :: hygro_pcg3_b_o_aer = 0.14
   real, parameter :: hygro_pcg4_b_o_aer = 0.14
   real, parameter :: hygro_pcg5_b_o_aer = 0.14
   real, parameter :: hygro_pcg6_b_o_aer = 0.14
   real, parameter :: hygro_pcg7_b_o_aer = 0.14
   real, parameter :: hygro_pcg8_b_o_aer = 0.14
   real, parameter :: hygro_pcg9_b_o_aer = 0.14
   real, parameter :: hygro_opcg1_b_c_aer = 0.14
   real, parameter :: hygro_opcg2_b_c_aer = 0.14
   real, parameter :: hygro_opcg3_b_c_aer = 0.14
   real, parameter :: hygro_opcg4_b_c_aer = 0.14
   real, parameter :: hygro_opcg5_b_c_aer = 0.14
   real, parameter :: hygro_opcg6_b_c_aer = 0.14
   real, parameter :: hygro_opcg7_b_c_aer = 0.14
   real, parameter :: hygro_opcg8_b_c_aer = 0.14
   real, parameter :: hygro_opcg1_b_o_aer = 0.14
   real, parameter :: hygro_opcg2_b_o_aer = 0.14
   real, parameter :: hygro_opcg3_b_o_aer = 0.14
   real, parameter :: hygro_opcg4_b_o_aer = 0.14
   real, parameter :: hygro_opcg5_b_o_aer = 0.14
   real, parameter :: hygro_opcg6_b_o_aer = 0.14
   real, parameter :: hygro_opcg7_b_o_aer = 0.14
   real, parameter :: hygro_opcg8_b_o_aer = 0.14
   real, parameter :: hygro_pcg1_f_c_aer = 0.14
   real, parameter :: hygro_pcg2_f_c_aer = 0.14
   real, parameter :: hygro_pcg3_f_c_aer = 0.14
   real, parameter :: hygro_pcg4_f_c_aer = 0.14
   real, parameter :: hygro_pcg5_f_c_aer = 0.14
   real, parameter :: hygro_pcg6_f_c_aer = 0.14
   real, parameter :: hygro_pcg7_f_c_aer = 0.14
   real, parameter :: hygro_pcg8_f_c_aer = 0.14
   real, parameter :: hygro_pcg9_f_c_aer = 0.14
   real, parameter :: hygro_pcg1_f_o_aer = 0.14
   real, parameter :: hygro_pcg2_f_o_aer = 0.14
   real, parameter :: hygro_pcg3_f_o_aer = 0.14
   real, parameter :: hygro_pcg4_f_o_aer = 0.14
   real, parameter :: hygro_pcg5_f_o_aer = 0.14
   real, parameter :: hygro_pcg6_f_o_aer = 0.14
   real, parameter :: hygro_pcg7_f_o_aer = 0.14
   real, parameter :: hygro_pcg8_f_o_aer = 0.14
   real, parameter :: hygro_pcg9_f_o_aer = 0.14
   real, parameter :: hygro_opcg1_f_c_aer = 0.14
   real, parameter :: hygro_opcg2_f_c_aer = 0.14
   real, parameter :: hygro_opcg3_f_c_aer = 0.14
   real, parameter :: hygro_opcg4_f_c_aer = 0.14
   real, parameter :: hygro_opcg5_f_c_aer = 0.14
   real, parameter :: hygro_opcg6_f_c_aer = 0.14
   real, parameter :: hygro_opcg7_f_c_aer = 0.14
   real, parameter :: hygro_opcg8_f_c_aer = 0.14
   real, parameter :: hygro_opcg1_f_o_aer = 0.14
   real, parameter :: hygro_opcg2_f_o_aer = 0.14
   real, parameter :: hygro_opcg3_f_o_aer = 0.14
   real, parameter :: hygro_opcg4_f_o_aer = 0.14
   real, parameter :: hygro_opcg5_f_o_aer = 0.14
   real, parameter :: hygro_opcg6_f_o_aer = 0.14
   real, parameter :: hygro_opcg7_f_o_aer = 0.14
   real, parameter :: hygro_opcg8_f_o_aer = 0.14
   real, parameter :: hygro_smpa_aer = 0.14
   real, parameter :: hygro_smpbb_aer = 0.14

   real, parameter :: hygro_glysoa_r1_aer = 0.14
   real, parameter :: hygro_glysoa_r2_aer = 0.14
   real, parameter :: hygro_glysoa_oh_aer = 0.14
   real, parameter :: hygro_glysoa_nh4_aer = 0.14
   real, parameter :: hygro_glysoa_sfc_aer = 0.14

   real, parameter :: hygro_ant1_c_aer = 0.14
   real, parameter :: hygro_ant2_c_aer = 0.14
   real, parameter :: hygro_ant3_c_aer = 0.14
   real, parameter :: hygro_ant4_c_aer = 0.14
   real, parameter :: hygro_ant1_o_aer = 0.14
   real, parameter :: hygro_ant2_o_aer = 0.14
   real, parameter :: hygro_ant3_o_aer = 0.14
   real, parameter :: hygro_ant4_o_aer = 0.14
   real, parameter :: hygro_biog1_c_aer = 0.14
   real, parameter :: hygro_biog2_c_aer = 0.14
   real, parameter :: hygro_biog3_c_aer = 0.14
   real, parameter :: hygro_biog4_c_aer = 0.14
   real, parameter :: hygro_biog1_o_aer = 0.14
   real, parameter :: hygro_biog2_o_aer = 0.14
   real, parameter :: hygro_biog3_o_aer = 0.14
   real, parameter :: hygro_biog4_o_aer = 0.14

   real, parameter :: hygro_asoaX_aer = 0.14
   real, parameter :: hygro_asoa1_aer = 0.14
   real, parameter :: hygro_asoa2_aer = 0.14
   real, parameter :: hygro_asoa3_aer = 0.14
   real, parameter :: hygro_asoa4_aer = 0.14
   real, parameter :: hygro_bsoaX_aer = 0.14
   real, parameter :: hygro_bsoa1_aer = 0.14
   real, parameter :: hygro_bsoa2_aer = 0.14
   real, parameter :: hygro_bsoa3_aer = 0.14
   real, parameter :: hygro_bsoa4_aer = 0.14

   integer, save :: &
      msectional, maerosolincw, &
      maerocoag, maerchem, maeroptical, maerchem_boxtest_output

!-----------------------------------------------------------------------
!   the following arrays are used during "column calculations" within
!        mosaic modules
!
!   aqvoldry_sub(n,k,m) = dry-volume (cm^3-aerosol/mole-air) for mode n,
!                        layer k, subarea m
!   aqmassdry_sub(n,k,m) = dry-mass (g-aerosol/mole-air) for mode n, ...
!   adrydens_sub(n,k,m) = dry-density (g-aerosol/cm^3-aerosol) for mode n, ...
!                        == amassdry_sub/avoldry_sub
!   awetdens_sub(n,k,m) = wet-density (g-aerosol/cm^3-aerosol) for mode n, ...
!
!   admeandry_sub(n,k,m) = current mean dry-diameter (cm)
!        for unactivated aerosol in mode n, layer k, subarea m.
!        (Used in sectional code)
!   admeanwet_sub(n,k,m) = current mean wet-diameter (cm)
!        for unactivated aerosol in mode n, layer k, subarea m.
!        (Used in sectional code)
!
!   awetdens_sfc(n,i,j) = wet-density for mode n and k=1.  This value is
!        saved after the second pass thru subr. aerosol_wetsize
!        for use in subr. bounds
!   admeanwet_sfc(n,i,j) = mean wet-diameter for mode n and k=1.
!        This value is saved after the second pass thru
!        subr. aerosol_wetsize for use in subr. bounds
!
!
!   following are used in aerosol growth routine and are locally
!   defined for the current layer (k) and subarea (m)
!        drymass_pregrow(n) = dry-mass (g/mole-air) for section n
!                        before the aerosol growth
!        drymass_aftgrow(n) = dry-mass (g/mole-air) for section n
!                        after the growth but before inter-section transfer
!        drydens_pregrow(n) = dry-density (g/cm3) for section n
!                        before the aerosol growth
!        drydens_aftgrow(n) = dry-density (g/cm3) for section n
!                        after the growth but before inter-section transfer
!
!-----------------------------------------------------------------------

!        integer, parameter :: imaxd_asize = 78
!        integer, parameter :: jmaxd_asize = 28
   integer, parameter :: kmaxd_asize = 100
   integer, parameter :: nsubareamaxd_asize = 5

! rce 22-jul-2006 - added maxd_atype dimension to following 6 arrays
!    and changed avoldry_sub,amassdry_sub to aqvoldry_sub,aqmassdry_sub
   real, save :: aqvoldry_sub(maxd_asize, maxd_atype, kmaxd_asize, nsubareamaxd_asize)
   real, save :: aqmassdry_sub(maxd_asize, maxd_atype, kmaxd_asize, nsubareamaxd_asize)
   real, save :: adrydens_sub(maxd_asize, maxd_atype, kmaxd_asize, nsubareamaxd_asize)
   real, save :: awetdens_sub(maxd_asize, maxd_atype, kmaxd_asize, nsubareamaxd_asize)
   real, save :: admeandry_sub(maxd_asize, maxd_atype, kmaxd_asize, nsubareamaxd_asize)
   real, save :: admeanwet_sub(maxd_asize, maxd_atype, kmaxd_asize, nsubareamaxd_asize)

!        real, save :: awetdens_sfc(maxd_asize,imaxd_asize,jmaxd_asize)
!        real, save :: admeanwet_sfc(maxd_asize,imaxd_asize,jmaxd_asize)

   real, save :: drymass_pregrow(maxd_asize, maxd_atype)
   real, save :: drydens_pregrow(maxd_asize, maxd_atype)
   real, save :: drymass_aftgrow(maxd_asize, maxd_atype)
   real, save :: drydens_aftgrow(maxd_asize, maxd_atype)

! table lookup of aerosol impaction/interception
! scavenging rates
   real dlndg_nimptblgrow
   integer nimptblgrow_mind, nimptblgrow_maxd
   parameter(nimptblgrow_mind=-7, nimptblgrow_maxd=12)
   real scavimptblnum(4, nimptblgrow_mind:nimptblgrow_maxd, maxd_asize, maxd_atype), &
      scavimptblvol(4, nimptblgrow_mind:nimptblgrow_maxd, maxd_asize, maxd_atype)

end module module_data_mosaic_asect