chem/module_mosaic_init_aerpar.F

   1 module module_mosaic_init_aerpar
   2
   3   use module_data_mosaic_kind, only: r8
   4
   5   implicit none
   6
   7   private
   8
   9   public:: mosaic_init_aer_params
  10
  11   contains
  12
  13
  14 !-------------------------------------------------------------------------------
  15   subroutine mosaic_init_aer_params
  16     !BSINGH - All initialzations for Mosiac model
  17
  18     call load_mosaic_parameters
  19
  20   end subroutine mosaic_init_aer_params
  21
  22
  23 !-------------------------------------------------------------------------------
  24   subroutine load_mosaic_parameters
  25   !---------------------------------------------------------------------------------------!
  26   !BSINGH: load_mosaic_parameters subroutine is directly copied form the mosaic_box.25.f90
  27   !        code
  28   !---------------------------------------------------------------------------------------!
  29
  30   !---------------------------------------------------------------------------------------!
  31   ! Called only once per entire simulation to load gas and aerosol
  32   ! indices, parameters, physico-chemical constants, polynomial coeffs, etc.
  33   !
  34   ! author: Rahul A. Zaveri
  35   ! update: jan 2005
  36   !---------------------------------------------------------------------------------------!
  37
  38     use module_mosaic_soa_vbs, only: soa_vbs_load_params
  39
  40     use module_data_mosaic_aero, only: ipmcmos_aero, no_aerosol, all_solid, all_liquid,   &
  41          mixed, nelectrolyte, naercomp, naer, Ncation, Nanion,                            &
  42          ngas_aerchtot, ngas_volatile, nsalt,                                             &
  43          jsulf_poor_NUM, jsulf_rich_NUM, MDRH_T_NUM, d_mdrh_DIM2, phasestate, aer_name,   &
  44          gas_name, ename, jnh4so4, jlvcite, jnh4hso4, jnh4msa, jnh4no3, jnh4cl, jna2so4,  &
  45          jna3hso4, jnahso4, jnamsa, jnano3, jnacl, jcano3, jcacl2, jcamsa2, jh2so4, jmsa, &
  46          jhno3, jhcl, jhhso4, jcaso4, jcaco3, joc, jbc, join, jaro1, jaro2, jalk1, jole1, &
  47          japi1, japi2, jlim1, jlim2, jh2o, jc_h, jc_nh4, jc_na, jc_ca, ja_hso4, ja_so4,   &
  48          ja_no3, ja_cl, ja_msa, ih2so4_g, ihno3_g, ihcl_g, inh3_g, imsa_g, iaro1_g,       &
  49          iaro2_g, ialk1_g, iole1_g, iapi1_g, iapi2_g, ilim1_g, ilim2_g, iso4_a, ino3_a,   &
  50          icl_a, inh4_a, imsa_a, iaro1_a, iaro2_a, ialk1_a, iole1_a, iapi1_a, iapi2_a,     &
  51          ilim1_a, ilim2_a, ico3_a, ina_a, ica_a, ioin_a, ioc_a, ibc_a,                    &
  52          isoa_first, jsoa_first, nmax_ASTEM, b_mtem,                                      &
  53          zc, za, b_zsr, a_zsr, aw_min, mw_electrolyte, dens_electrolyte,                  &
  54          partial_molar_vol, MW_c, MW_a, mw_aer_mac,dens_aer_mac, kappa_aer_mac,           &
  55          dens_comp_a, mw_comp_a, ref_index_a, mw_gas, v_molar_gas,                        &
  56          rtol_mesa, jsalt_index, jsulf_poor, jsulf_rich, Nmax_mesa, d_mdrh, msoa_flag1,   &
  57          use_cam5mam_soa_params
  58
  59     ! local variables
  60     integer :: ia, iaer, ig, igas, ibin, ja, je, j_index
  61     logical :: use_mos31e_rz1_densities
  62     logical :: use_uniform_densities !BSINGH - 05/28/2013(RCE updates)
  63     logical :: use_sorgam_soa_species
  64     real(r8), dimension(nelectrolyte) :: G_MX,K_MX
  65
  66     !BSINGH - 05/28/2013(RCE updates)
  67     use_mos31e_rz1_densities = .true.
  68     if ( use_mos31e_rz1_densities ) then
  69        use_uniform_densities = .false.
  70     else
  71        use_uniform_densities = .true.
  72        if (ipmcmos_aero > 0) use_uniform_densities = .false.
  73     end if
  74     !BSINGH - 05/28/2013(RCE updates ENDS)
  75
  76     if (msoa_flag1 == 1) then
  77        use_sorgam_soa_species = .true.
  78     else
  79        use_sorgam_soa_species = .false.
  80     end if
  81
  82     ! rce 2013-07-31 -
  83     !   using a local saved variable like "first" no longer works
  84     !   the calling routine needs to determine if/when this routine is needed
  85     !   if(first)then
  86     !      first=.false.
  87
  88     !----------------------------------------------------------------
  89     ! control settings
  90     ! *** do not change mSIZE_FRAMEWORK here ***
  91     !      mSIZE_FRAMEWORK = mSECTIONAL        ! mMODAL or mSECTIONAL
  92     !      mDYNAMIC_SOLVER = mASTEM            ! mASTEM, mLSODES
  93     !      mGAS_AER_XFER   = mON               ! mON, mOFF
  94
  95     ! ASTEM parameters
  96     nmax_ASTEM      = 301              ! max number of time steps in ASTEM
  97     !      alpha_ASTEM     = 1.0               ! choose a value between 0.01 and 1.0
  98     !      rtol_eqb_ASTEM = 0.01               ! equilibrium tolerance in ASTEM
  99     !      ptol_mol_ASTEM = 0.01               ! mol percent tolerance in ASTEM
 100
 101     ! MESA parameters
 102     Nmax_MESA       = 80               ! max number of iterations in MESA_PTC
 103     rtol_mesa       = 0.01             ! MESA equilibrium tolerance
 104     !----------------------------------------------------------------
 105
 106     !
 107     ! set species indices
 108     !    ixxx_a are for aerosol     species - they apply to aer and total arrays
 109     !    ixxx_g are for gas         species - they apply to gas array
 110     !    jyyy   are for electrolyte species - they apply to electrolyte and comp_a arrays
 111     !
 112     ! *** CRITICAL "RULES" ABOUT THESE INDICES ***
 113     !
 114     ! *** THE CODE WILL FAIL IF THEY ARE NOT FOLLOWED  ***
 115     !
 116     ! for the volatile inorganic and organic species (1 thru ngas_volatile)
 117     !    the ixxx_a and ixxx_g must match (be indentical) for each aerosol-gas pair
 118     !
 119     ! for the soa species, the ordering of the ixxx_a, ixxx_g, and jyyy
 120     !    must be identical
 121     !
 122
 123     ! electrolyte indices (used for water content calculations)
 124     ! these indices are order sensitive
 125     ! inorganic species first
 126     jnh4so4    =  1    ! soluble
 127     jlvcite    =  2    ! soluble
 128     jnh4hso4   =  3    ! soluble
 129     jnh4msa    =  4    ! soluble: new
 130     jnh4no3    =  5    ! soluble
 131     jnh4cl     =  6    ! soluble
 132     jna2so4    =  7    ! soluble
 133     jna3hso4   =  8    ! soluble
 134     jnahso4    =  9    ! soluble
 135     jnamsa     = 10    ! soluble: new
 136     jnano3     = 11    ! soluble
 137     jnacl      = 12    ! soluble
 138     jcano3     = 13    ! soluble
 139     jcacl2     = 14    ! soluble
 140     jcamsa2    = 15    ! soluble     nsalt
 141     jh2so4     = 16    ! soluble
 142     jmsa       = 17    ! soluble
 143     jhno3      = 18    ! soluble
 144     jhcl       = 19    ! soluble
 145     jhhso4     = 20    ! soluble
 146     jcaso4     = 21    ! insoluble
 147     jcaco3     = 22    ! insoluble
 148     joc        = 23    ! insoluble - part of naercomp
 149     jbc        = 24    ! insoluble - part of naercomp
 150     join       = 25    ! insoluble - part of naercomp
 151
 152     ! aerosol and gas indices for inorganic species
 153     iso4_a     =  1 ;  ih2so4_g   =  1
 154     ino3_a     =  2 ;  ihno3_g    =  2
 155     icl_a      =  3 ;  ihcl_g     =  3
 156     inh4_a     =  4 ;  inh3_g     =  4
 157     imsa_a     =  5 ;  imsa_g     =  5
 158
 159     ia = imsa_a ; ig = imsa_g ; je = join
 160
 161     ! aerosol, gas, and electrolyte indices for sorgam secondary organic species
 162     if ( use_sorgam_soa_species ) then
 163        isoa_first = ia+1 ; jsoa_first = je+1
 164        iaro1_a    = ia+1 ;  iaro1_g    = ig+1 ;  jaro1      = je+1
 165        iaro2_a    = ia+2 ;  iaro2_g    = ig+2 ;  jaro2      = je+2
 166        ialk1_a    = ia+3 ;  ialk1_g    = ig+3 ;  jalk1      = je+3
 167        iole1_a    = ia+4 ;  iole1_g    = ig+4 ;  jole1      = je+4
 168        iapi1_a    = ia+5 ;  iapi1_g    = ig+5 ;  japi1      = je+5
 169        iapi2_a    = ia+6 ;  iapi2_g    = ig+6 ;  japi2      = je+6
 170        ilim1_a    = ia+7 ;  ilim1_g    = ig+7 ;  jlim1      = je+7
 171        ilim2_a    = ia+8 ;  ilim2_g    = ig+8 ;  jlim2      = je+8
 172        ia = ia+8 ; ig = ig+8 ; je = je+8
 173     else
 174        isoa_first = ia+1 ; jsoa_first = je+1
 175        iaro1_a    = -999888777 ;  iaro1_g    = -999888777 ;  jaro1      = -999888777
 176        iaro2_a    = -999888777 ;  iaro2_g    = -999888777 ;  jaro2      = -999888777
 177        ialk1_a    = -999888777 ;  ialk1_g    = -999888777 ;  jalk1      = -999888777
 178        iole1_a    = -999888777 ;  iole1_g    = -999888777 ;  jole1      = -999888777
 179        iapi1_a    = -999888777 ;  iapi1_g    = -999888777 ;  japi1      = -999888777
 180        iapi2_a    = -999888777 ;  iapi2_g    = -999888777 ;  japi2      = -999888777
 181        ilim1_a    = -999888777 ;  ilim1_g    = -999888777 ;  jlim1      = -999888777
 182        ilim2_a    = -999888777 ;  ilim2_g    = -999888777 ;  jlim2      = -999888777
 183        ia = ia+1 ; ig = ig+1 ; je = je+1
 184     end if
 185
 186     ! electrolyte indices for other species
 187     jh2o       = je+1    ! water - part of naercomp
 188
 189     ! aerosol indices for other species
 190     ico3_a     = ia+1
 191     ina_a      = ia+2
 192     ica_a      = ia+3
 193     ioin_a     = ia+4
 194     ioc_a      = ia+5
 195     ibc_a      = ia+6
 196
 197     ! gas indices for other species
 198     ! ico2_g   = 14  ! *** currently not used
 199                      ! *** if co3_a was treated as a volatile inorganic,
 200                      !     then ico3_a and ico2_g would have to be 6
 201
 202     ! local aerosol ions
 203     ! cations
 204     jc_h       =  1
 205     jc_nh4     =  2
 206     jc_na      =  3
 207     jc_ca      =  4
 208     !
 209     ! anions
 210     ja_hso4    =  1
 211     ja_so4     =  2
 212     ja_no3     =  3
 213     ja_cl      =  4
 214     ja_msa     =  5
 215     ! ja_co3   =  6   ! *** currently not used
 216
 217
 218     if (msoa_flag1 >= 1000) call soa_vbs_load_params( 1 )
 219
 220
 221     !--------------------------------------------------------------------
 222     ! phase state names
 223     phasestate(no_aerosol) = "NOAERO"
 224     phasestate(all_solid)  = "SOLID "
 225     phasestate(all_liquid) = "LIQUID"
 226     phasestate(mixed)      = "MIXED "
 227
 228     ! names of aer species
 229     do iaer = 1, naer
 230        write( aer_name(iaer), '(a,i4.4)' ) 'aer', iaer  ! default
 231     end do
 232
 233     aer_name(iso4_a) = "SO4"
 234     aer_name(ino3_a) = "NO3"
 235     aer_name(icl_a)  = "Cl "
 236     aer_name(inh4_a) = "NH4"
 237     aer_name(ioc_a)  = "OC "
 238     aer_name(imsa_a) = "MSA"
 239     aer_name(ico3_a) = "CO3"
 240     aer_name(ina_a)  = "Na "
 241     aer_name(ica_a)  = "Ca "
 242     aer_name(ibc_a)  = "BC "
 243     aer_name(ioin_a) = "OIN"
 244     if ( use_sorgam_soa_species ) then
 245     aer_name(iaro1_a)= "ARO1"
 246     aer_name(iaro2_a)= "ARO2"
 247     aer_name(ialk1_a)= "ALK1"
 248     aer_name(iole1_a)= "OLE1"
 249     aer_name(iapi1_a)= "API1"
 250     aer_name(iapi2_a)= "API2"
 251     aer_name(ilim1_a)= "LIM1"
 252     aer_name(ilim2_a)= "LIM2"
 253     end if
 254
 255     ! names of gas species
 256     do igas = 1, ngas_aerchtot
 257        write( gas_name(igas), '(a,i4.4)' ) 'gas', igas  ! default
 258     end do
 259
 260     gas_name(ih2so4_g) = "H2SO4"
 261     gas_name(ihno3_g)  = "HNO3 "
 262     gas_name(ihcl_g)   = "HCl  "
 263     gas_name(inh3_g)   = "NH3  "
 264     gas_name(imsa_g)   = "MSA  "
 265     if ( use_sorgam_soa_species ) then
 266     gas_name(iaro1_g)   = "ARO1 "
 267     gas_name(iaro2_g)   = "ARO2 "
 268     gas_name(ialk1_g)   = "ALK1 "
 269     gas_name(iole1_g)   = "OLE1 "
 270     gas_name(iapi1_g)   = "API1 "
 271     gas_name(iapi2_g)   = "API2 "
 272     gas_name(ilim1_g)   = "LIM1 "
 273     gas_name(ilim2_g)   = "LIM2 "
 274     end if
 275
 276     ! names of electrolytes
 277     ename(jnh4so4) = "AmSO4"
 278     ename(jlvcite) = "(NH4)3H(SO4)2"
 279     ename(jnh4hso4)= "NH4HSO4"
 280     ename(jnh4msa) = "CH3SO3NH4"
 281     ename(jnh4no3) = "NH4NO3"
 282     ename(jnh4cl)  = "NH4Cl"
 283     ename(jnacl)   = "NaCl"
 284     ename(jnano3)  = "NaNO3"
 285     ename(jna2so4) = "Na2SO4"
 286     ename(jna3hso4)= "Na3H(SO4)2"
 287     ename(jnamsa)  = "CH3SO3Na"
 288     ename(jnahso4) = "NaHSO4"
 289     ename(jcaso4)  = "CaSO4"
 290     ename(jcamsa2) = "(CH3SO3)2Ca"
 291     ename(jcano3)  = "Ca(NO3)2"
 292     ename(jcacl2)  = "CaCl2"
 293     ename(jcaco3)  = "CaCO3"
 294     ename(jh2so4)  = "H2SO4"
 295     ename(jhhso4)  = "HHSO4"
 296     ename(jhno3)   = "HNO3"
 297     ename(jhcl)    = "HCl"
 298     ename(jmsa)    = "CH3SO3H"
 299
 300     ! molecular weights of electrolytes
 301     mw_electrolyte(jnh4so4) = 132.0
 302     mw_electrolyte(jlvcite) = 247.0
 303     mw_electrolyte(jnh4hso4)= 115.0
 304     mw_electrolyte(jnh4msa) = 113.0
 305     mw_electrolyte(jnh4no3) = 80.0
 306     mw_electrolyte(jnh4cl)  = 53.5
 307     mw_electrolyte(jnacl)   = 58.5
 308     mw_electrolyte(jnano3)  = 85.0
 309     mw_electrolyte(jna2so4) = 142.0
 310     mw_electrolyte(jna3hso4)= 262.0
 311     mw_electrolyte(jnahso4) = 120.0
 312     mw_electrolyte(jnamsa)  = 118.0
 313     mw_electrolyte(jcaso4)  = 136.0
 314     mw_electrolyte(jcamsa2) = 230.0
 315     mw_electrolyte(jcano3)  = 164.0
 316     mw_electrolyte(jcacl2)  = 111.0
 317     mw_electrolyte(jcaco3)  = 100.0
 318     mw_electrolyte(jh2so4)  = 98.0
 319     mw_electrolyte(jhno3)   = 63.0
 320     mw_electrolyte(jhcl)    = 36.5
 321     mw_electrolyte(jmsa)    = 96.0
 322
 323
 324     ! molecular weights of ions [g/mol]
 325     MW_c(jc_h)  =  1.0
 326     MW_c(jc_nh4)= 18.0
 327     MW_c(jc_na) = 23.0
 328     MW_c(jc_ca) = 40.0
 329
 330     MW_a(ja_so4) = 96.0
 331     MW_a(ja_hso4)= 97.0
 332     MW_a(ja_no3) = 62.0
 333     MW_a(ja_cl)  = 35.5
 334     MW_a(ja_msa) = 95.0
 335
 336
 337     ! magnitude of the charges on ions
 338     zc(jc_h)   = 1
 339     zc(jc_nh4) = 1
 340     zc(jc_na)  = 1
 341     zc(jc_ca)  = 2
 342
 343     za(ja_hso4)= 1
 344     za(ja_so4) = 2
 345     za(ja_no3) = 1
 346     za(ja_cl)  = 1
 347     za(ja_msa) = 1
 348
 349
 350     ! densities of pure electrolytes in g/cc
 351     dens_electrolyte(jnh4so4)  = 1.8
 352     dens_electrolyte(jlvcite)  = 1.8
 353     dens_electrolyte(jnh4hso4) = 1.8
 354     dens_electrolyte(jnh4msa)  = 1.8 ! assumed same as nh4hso4
 355     dens_electrolyte(jnh4no3)  = 1.8
 356     dens_electrolyte(jnh4cl)   = 1.8
 357     dens_electrolyte(jnacl)    = 2.2
 358     dens_electrolyte(jnano3)   = 2.2
 359     dens_electrolyte(jna2so4)  = 2.2
 360     dens_electrolyte(jna3hso4) = 2.2
 361     dens_electrolyte(jnahso4)  = 2.2
 362     dens_electrolyte(jnamsa)   = 2.2 ! assumed same as nahso4
 363     dens_electrolyte(jcaso4)   = 2.6
 364     dens_electrolyte(jcamsa2)  = 2.6   ! assumed same as caso4
 365     dens_electrolyte(jcano3)   = 2.6
 366     dens_electrolyte(jcacl2)   = 2.6
 367     dens_electrolyte(jcaco3)   = 2.6
 368     dens_electrolyte(jh2so4)   = 1.8
 369     dens_electrolyte(jhhso4)   = 1.8
 370     dens_electrolyte(jhno3)    = 1.8
 371     dens_electrolyte(jhcl)     = 1.8
 372     dens_electrolyte(jmsa)     = 1.8 ! assumed same as h2so4
 373     if ( use_uniform_densities ) then!BSINGH - 05/28/2013(RCE updates)
 374        do je = 1, nelectrolyte
 375           dens_electrolyte(je) = 1.6
 376        enddo
 377     endif!BSINGH - 05/28/2013(RCE updates)
 378
 379     ! densities of compounds in g/cc
 380     dens_comp_a(jnh4so4)  = 1.8
 381     dens_comp_a(jlvcite)  = 1.8
 382     dens_comp_a(jnh4hso4) = 1.8
 383     dens_comp_a(jnh4msa)  = 1.8        ! assumed same as nh4hso4
 384     dens_comp_a(jnh4no3)  = 1.7
 385     dens_comp_a(jnh4cl)   = 1.5
 386     dens_comp_a(jnacl)    = 2.2
 387     dens_comp_a(jnano3)   = 2.2
 388     dens_comp_a(jna2so4)  = 2.2
 389     dens_comp_a(jna3hso4) = 2.2
 390     dens_comp_a(jnahso4)  = 2.2
 391     dens_comp_a(jnamsa)   = 2.2        ! assumed same as nahso4
 392     dens_comp_a(jcaso4)   = 2.6
 393     dens_comp_a(jcamsa2)  = 2.6        ! assumed same as caso4
 394     dens_comp_a(jcano3)   = 2.6
 395     dens_comp_a(jcacl2)   = 2.6
 396     dens_comp_a(jcaco3)   = 2.6
 397     dens_comp_a(jh2so4)   = 1.8
 398     dens_comp_a(jhhso4)   = 1.8
 399     dens_comp_a(jhno3)    = 1.8
 400     dens_comp_a(jhcl)     = 1.8
 401     dens_comp_a(jmsa)     = 1.8        ! assumed same as h2so4
 402     dens_comp_a(joc)      = 1.0
 403     dens_comp_a(jbc)      = 1.8
 404     dens_comp_a(join)     = 2.6
 405     if ( use_sorgam_soa_species ) then
 406     dens_comp_a(jaro1)    = 1.0
 407     dens_comp_a(jaro2)    = 1.0
 408     dens_comp_a(jalk1)    = 1.0
 409     dens_comp_a(jole1)    = 1.0
 410     dens_comp_a(japi1)    = 1.0
 411     dens_comp_a(japi2)    = 1.0
 412     dens_comp_a(jlim1)    = 1.0
 413     dens_comp_a(jlim2)    = 1.0
 414     end if
 415     dens_comp_a(jh2o)     = 1.0
 416     !BSINGH - 05/28/2013(RCE updates)
 417     ! following for comparison with mos31d_bs2 and. mos31e_rz1
 418     if ( use_mos31e_rz1_densities ) then
 419        dens_comp_a(joc)      = 1.4
 420        if ( use_sorgam_soa_species ) then
 421        dens_comp_a(jaro1)    = 1.4
 422        dens_comp_a(jaro2)    = 1.4
 423        dens_comp_a(jalk1)    = 1.4
 424        dens_comp_a(jole1)    = 1.4
 425        dens_comp_a(japi1)    = 1.4
 426        dens_comp_a(japi2)    = 1.4
 427        dens_comp_a(jlim1)    = 1.4
 428        dens_comp_a(jlim2)    = 1.4
 429        end if
 430     end if
 431
 432     if ( use_uniform_densities ) then
 433        !BSINGH - 05/28/2013(RCE updates ENDS)
 434        do je = 1, naercomp
 435           dens_comp_a(je) = 1.6
 436        enddo
 437        !BSINGH - 05/28/2013(RCE updates)
 438     endif
 439
 440     if (ipmcmos_aero > 0) then
 441        dens_comp_a(jnh4no3)  = 1.8
 442        dens_comp_a(jnh4cl)   = 1.8
 443        if ( use_sorgam_soa_species ) then
 444        dens_comp_a(jaro1)    = 1.4
 445        dens_comp_a(jaro2)    = 1.4
 446        dens_comp_a(jalk1)    = 1.4
 447        dens_comp_a(jole1)    = 1.4
 448        dens_comp_a(japi1)    = 1.4
 449        dens_comp_a(japi2)    = 1.4
 450        dens_comp_a(jlim1)    = 1.4
 451        dens_comp_a(jlim2)    = 1.4
 452        end if
 453     endif
 454     !BSINGH - 05/28/2013(RCE updates ENDS)
 455
 456     ! molecular weights of generic aerosol species
 457     mw_aer_mac(1:naer) = 200.0  ! default
 458
 459     mw_aer_mac(iso4_a) = 96.0
 460     mw_aer_mac(ino3_a) = 62.0
 461     mw_aer_mac(icl_a)  = 35.5
 462     mw_aer_mac(imsa_a) = 95.0  ! CH3SO3
 463     mw_aer_mac(ico3_a) = 60.0
 464     mw_aer_mac(inh4_a) = 18.0
 465     mw_aer_mac(ina_a)  = 23.0
 466     mw_aer_mac(ica_a)  = 40.0
 467     mw_aer_mac(ioin_a) = 1.0           ! not used
 468     mw_aer_mac(ibc_a)  = 1.0           ! not used
 469     mw_aer_mac(ioc_a)  = 1.0   ! 200 assumed for primary organics
 470     if ( use_sorgam_soa_species ) then
 471     mw_aer_mac(iaro1_a)= 150.0
 472     mw_aer_mac(iaro2_a)= 150.0
 473     mw_aer_mac(ialk1_a)= 140.0
 474     mw_aer_mac(iole1_a)= 140.0
 475     mw_aer_mac(iapi1_a)= 184.0
 476     mw_aer_mac(iapi2_a)= 184.0
 477     mw_aer_mac(ilim1_a)= 200.0
 478     mw_aer_mac(ilim2_a)= 200.0
 479     end if
 480
 481     ! molecular weights of compounds
 482     mw_comp_a(jnh4so4) = 132.0
 483     mw_comp_a(jlvcite) = 247.0
 484     mw_comp_a(jnh4hso4)= 115.0
 485     mw_comp_a(jnh4msa) = 113.0
 486     mw_comp_a(jnh4no3) =  80.0
 487     mw_comp_a(jnh4cl)  =  53.5
 488     mw_comp_a(jnacl)   =  58.5
 489     mw_comp_a(jnano3)  =  85.0
 490     mw_comp_a(jna2so4) = 142.0
 491     mw_comp_a(jna3hso4)= 262.0
 492     mw_comp_a(jnahso4) = 120.0
 493     mw_comp_a(jnamsa)  = 118.0
 494     mw_comp_a(jcaso4)  = 136.0
 495     mw_comp_a(jcamsa2) = 230.0
 496     mw_comp_a(jcano3)  = 164.0
 497     mw_comp_a(jcacl2)  = 111.0
 498     mw_comp_a(jcaco3)  = 100.0
 499     mw_comp_a(jh2so4)  =  98.0
 500     mw_comp_a(jhhso4)  =  98.0
 501     mw_comp_a(jhno3)   =  63.0
 502     mw_comp_a(jhcl)    =  36.5
 503     mw_comp_a(jmsa)    =  96.0
 504     mw_comp_a(joc)     =   1.0
 505     mw_comp_a(jbc)     =   1.0
 506     mw_comp_a(join)    =   1.0
 507     if ( use_sorgam_soa_species ) then
 508     mw_comp_a(jaro1)   = 150.0
 509     mw_comp_a(jaro2)   = 150.0
 510     mw_comp_a(jalk1)   = 140.0
 511     mw_comp_a(jole1)   = 140.0
 512     mw_comp_a(japi1)   = 184.0
 513     mw_comp_a(japi2)   = 184.0
 514     mw_comp_a(jlim1)   = 200.0
 515     mw_comp_a(jlim2)   = 200.0
 516     end if
 517     mw_comp_a(jh2o)    = 18.0
 518
 519     !BSINGH - 05/28/2013(RCE updates)
 520     ! partmc-2.2.1 jun-2012
 521     !#     dens (kg/m^3)   ions in soln (1)    molec wght (kg/mole)   kappa (1)
 522     ! SO4            1800                0                   96d-3      0.65
 523     ! NO3            1800                0                   62d-3      0.65
 524     ! Cl             2200                0                   35.5d-3    0.53
 525     ! NH4            1800                0                   18d-3      0.65
 526     ! MSA            1800                0                   95d-3      0.53
 527     ! ARO1           1400                0                   150d-3     0.1
 528     ! ARO2           1400                0                   150d-3     0.1
 529     ! ALK1           1400                0                   140d-3     0.1
 530     ! OLE1           1400                0                   140d-3     0.1
 531     ! API1           1400                0                   184d-3     0.1
 532     ! API2           1400                0                   184d-3     0.1
 533     ! LIM1           1400                0                   200d-3     0.1
 534     ! LIM2           1400                0                   200d-3     0.1
 535     ! CO3            2600                0                   60d-3      0.53
 536     ! Na             2200                0                   23d-3      0.53
 537     ! Ca             2600                0                   40d-3      0.53
 538     ! OIN            2600                0                   1d-3       0.1
 539     ! OC             1000                0                   1d-3       0.001
 540     ! BC             1800                0                   1d-3       0
 541     ! H2O            1000                0                   18d-3      0
 542     !BSINGH - 05/28/2013(RCE updates ENDS)
 543
 544     ! densities of generic aerosol species
 545     dens_aer_mac(1:naer) = 1.0   ! default
 546
 547     dens_aer_mac(iso4_a) = 1.8 ! used
 548     dens_aer_mac(ino3_a) = 1.8 ! used
 549     dens_aer_mac(icl_a)  = 2.2 ! used
 550     dens_aer_mac(imsa_a) = 1.8 ! used
 551     dens_aer_mac(ico3_a) = 2.6 ! used
 552     dens_aer_mac(inh4_a) = 1.8 ! used
 553     dens_aer_mac(ina_a)  = 2.2 ! used
 554     dens_aer_mac(ica_a)  = 2.6 ! used
 555     dens_aer_mac(ioin_a) = 2.6 ! used
 556     dens_aer_mac(ioc_a)  = 1.0 ! used
 557     dens_aer_mac(ibc_a)  = 1.8 ! used
 558     if ( use_sorgam_soa_species ) then
 559     dens_aer_mac(iaro1_a)= 1.0
 560     dens_aer_mac(iaro2_a)= 1.0
 561     dens_aer_mac(ialk1_a)= 1.0
 562     dens_aer_mac(iole1_a)= 1.0
 563     dens_aer_mac(iapi1_a)= 1.0
 564     dens_aer_mac(iapi2_a)= 1.0
 565     dens_aer_mac(ilim1_a)= 1.0
 566     dens_aer_mac(ilim2_a)= 1.0
 567     end if
 568     !BSINGH - 05/28/2013(RCE updates)
 569     ! following for comparison with mos31d_bs2 and. mos31e_rz1
 570     if ( use_mos31e_rz1_densities ) then
 571        dens_aer_mac(ioc_a)  = 1.4
 572        if ( use_sorgam_soa_species ) then
 573        dens_aer_mac(iaro1_a)= 1.4
 574        dens_aer_mac(iaro2_a)= 1.4
 575        dens_aer_mac(ialk1_a)= 1.4
 576        dens_aer_mac(iole1_a)= 1.4
 577        dens_aer_mac(iapi1_a)= 1.4
 578        dens_aer_mac(iapi2_a)= 1.4
 579        dens_aer_mac(ilim1_a)= 1.4
 580        dens_aer_mac(ilim2_a)= 1.4
 581        end if
 582     end if
 583
 584     if ( use_uniform_densities ) then
 585        !BSINGH - 05/28/2013(RCE updates ENDS)
 586
 587        do iaer = 1, naer
 588           dens_aer_mac(iaer) = 1.6
 589        enddo
 590     endif!BSINGH - 05/28/2013(RCE updates)
 591
 592     if (ipmcmos_aero > 0) then
 593        ! use partmc-mosaic densities
 594        dens_aer_mac(1:19) = (/ &
 595             1.80, 1.80, 2.20, 1.80, 1.80, 1.40, 1.40, 1.40, 1.40, 1.40, &
 596             1.40, 1.40, 1.40, 2.60, 2.20, 2.60, 2.60, 1.00, 1.80 /)!BSINGH - 05/28/2013(RCE updates)
 597        !           so4   no3   cl    nh4   msa   aro1  aro2  alk1  ole1  api1
 598        !           api2  lim1  lim2  co3   na    ca    oin   oc    bc
 599     end if
 600
 601     if ( use_cam5mam_soa_params > 0 ) then
 602        dens_aer_mac(ioc_a)   = 1.0
 603        dens_comp_a(joc)      = 1.0
 604        ! for oc, leave mw=1 because some of the mosaic code requires this
 605        if ( 1 <= ilim2_a .and. ilim2_a <= naer ) then
 606        dens_aer_mac(ilim2_a) = 1.0
 607        dens_comp_a(jlim2)    = 1.0
 608        mw_aer_mac(ilim2_a)   = 150.0
 609        mw_comp_a(jlim2)      = 150.0
 610        end if
 611     end if
 612
 613     ! kappa values (hygroscopicities) of generic aerosol species
 614     !
 615     ! for calculation of ccn properties, kappa of electrolytes
 616     !    should be used
 617     ! the multi-dimensional sectional code needs a "fixed" kappa
 618     !    for each generic aerosol species, just as the older
 619     !    1d sectional code needs a "fixed" dry density
 620     kappa_aer_mac(1:naer)  = 0.1  ! default
 621
 622     kappa_aer_mac(iso4_a)  = 0.65
 623     kappa_aer_mac(ino3_a)  = 0.65
 624     kappa_aer_mac(imsa_a)  = 0.65
 625     kappa_aer_mac(inh4_a)  = 0.65
 626     kappa_aer_mac(icl_a)   = 0.65
 627     kappa_aer_mac(ina_a)   = 0.65
 628     kappa_aer_mac(ico3_a)  = 0.001  ! ??
 629     kappa_aer_mac(ica_a)   = 0.001  ! ??
 630     kappa_aer_mac(ioin_a)  = 0.001
 631     kappa_aer_mac(ioc_a)   = 0.001
 632     kappa_aer_mac(ibc_a)   = 0.001
 633     if ( use_sorgam_soa_species ) then
 634     kappa_aer_mac(iaro1_a) = 0.1
 635     kappa_aer_mac(iaro2_a) = 0.1
 636     kappa_aer_mac(ialk1_a) = 0.1
 637     kappa_aer_mac(iole1_a) = 0.1
 638     kappa_aer_mac(iapi1_a) = 0.1
 639     kappa_aer_mac(iapi2_a) = 0.1
 640     kappa_aer_mac(ilim1_a) = 0.1
 641     kappa_aer_mac(ilim2_a) = 0.1
 642     end if
 643     !BSINGH - 05/28/2013(RCE updates)
 644     if (ipmcmos_aero > 0) then
 645        ! use partmc-mosaic kappas
 646        kappa_aer_mac(1:19) = (/ &
 647             0.65, 0.65, 0.53, 0.65, 0.53, 0.10, 0.10, 0.10, 0.10, 0.10, &
 648             0.10, 0.10, 0.10, 0.53, 0.53, 0.53, 0.10, 0.001, 0.0 /)
 649        !           so4   no3   cl    nh4   msa   aro1  aro2  alk1  ole1  api1
 650        !           api2  lim1  lim2  co3   na    ca    oin   oc    bc
 651     end if
 652     !BSINGH - 05/28/2013(RCE updates ENDS)
 653
 654     ! partial molar volumes of condensing gases
 655     partial_molar_vol(1:ngas_aerchtot)   = 200.0  ! default
 656
 657     partial_molar_vol(ih2so4_g) = 51.83
 658     partial_molar_vol(ihno3_g)  = 31.45
 659     partial_molar_vol(ihcl_g)   = 20.96
 660     partial_molar_vol(inh3_g)   = 24.03
 661     partial_molar_vol(imsa_g)   = 53.33
 662     if ( use_sorgam_soa_species ) then
 663     partial_molar_vol(iaro1_g)  = 150.0
 664     partial_molar_vol(iaro2_g)  = 150.0
 665     partial_molar_vol(ialk1_g)  = 140.0
 666     partial_molar_vol(iole1_g)  = 140.0
 667     partial_molar_vol(iapi1_g)  = 184.0
 668     partial_molar_vol(iapi2_g)  = 184.0
 669     partial_molar_vol(ilim1_g)  = 200.0
 670     partial_molar_vol(ilim2_g)  = 200.0
 671     end if
 672
 673     ! molecular weights of condensing gases
 674     mw_gas(1:ngas_aerchtot) = 200.0  ! default
 675     mw_gas(ih2so4_g) = 98.0
 676     mw_gas(ihno3_g)  = 63.0
 677     mw_gas(ihcl_g)   = 36.5
 678     mw_gas(inh3_g)   = 17.0
 679     mw_gas(imsa_g)   = 96.0
 680     if ( use_sorgam_soa_species ) then
 681     mw_gas(iaro1_g)  = 150.0
 682     mw_gas(iaro2_g)  = 150.0
 683     mw_gas(ialk1_g)  = 140.0
 684     mw_gas(iole1_g)  = 140.0
 685     mw_gas(iapi1_g)  = 184.0
 686     mw_gas(iapi2_g)  = 184.0
 687     mw_gas(ilim1_g)  = 200.0
 688     mw_gas(ilim2_g)  = 200.0
 689     end if
 690
 691     ! used to calculate diffusivities of condensing gases
 692     v_molar_gas(1:ngas_aerchtot) = 60.0  ! default
 693     v_molar_gas(ih2so4_g)= 42.88
 694     v_molar_gas(ihno3_g) = 24.11
 695     v_molar_gas(ihcl_g)  = 21.48
 696     v_molar_gas(inh3_g)  = 14.90
 697     v_molar_gas(imsa_g)  = 58.00
 698
 699     ! refractive index
 700     ref_index_a(jnh4so4) = cmplx(1.52,0.)
 701     ref_index_a(jlvcite) = cmplx(1.50,0.)
 702     ref_index_a(jnh4hso4)= cmplx(1.47,0.)
 703     ref_index_a(jnh4msa) = cmplx(1.50,0.)      ! assumed
 704     ref_index_a(jnh4no3) = cmplx(1.50,0.)
 705     ref_index_a(jnh4cl)  = cmplx(1.50,0.)
 706     ref_index_a(jnacl)   = cmplx(1.45,0.)
 707     ref_index_a(jnano3)  = cmplx(1.50,0.)
 708     ref_index_a(jna2so4) = cmplx(1.50,0.)
 709     ref_index_a(jna3hso4)= cmplx(1.50,0.)
 710     ref_index_a(jnahso4) = cmplx(1.50,0.)
 711     ref_index_a(jnamsa)  = cmplx(1.50,0.)      ! assumed
 712     ref_index_a(jcaso4)  = cmplx(1.56,0.006)
 713     ref_index_a(jcamsa2) = cmplx(1.56,0.006)   ! assumed
 714     ref_index_a(jcano3)  = cmplx(1.56,0.006)
 715     ref_index_a(jcacl2)  = cmplx(1.52,0.006)
 716     ref_index_a(jcaco3)  = cmplx(1.68,0.006)
 717     ref_index_a(jh2so4)  = cmplx(1.43,0.)
 718     ref_index_a(jhhso4)  = cmplx(1.43,0.)
 719     ref_index_a(jhno3)   = cmplx(1.50,0.)
 720     ref_index_a(jhcl)    = cmplx(1.50,0.)
 721     ref_index_a(jmsa)    = cmplx(1.43,0.)      ! assumed
 722     ref_index_a(joc)      = cmplx(1.45,0.)
 723     ref_index_a(jbc)      = cmplx(1.82,0.74)
 724     ref_index_a(join)    = cmplx(1.55,0.006)
 725     if ( use_sorgam_soa_species ) then
 726     ref_index_a(jaro1)   = cmplx(1.45,0.)
 727     ref_index_a(jaro2)   = cmplx(1.45,0.)
 728     ref_index_a(jalk1)   = cmplx(1.45,0.)
 729     ref_index_a(jole1)   = cmplx(1.45,0.)
 730     ref_index_a(japi1)   = cmplx(1.45,0.)
 731     ref_index_a(japi2)   = cmplx(1.45,0.)
 732     ref_index_a(jlim1)   = cmplx(1.45,0.)
 733     ref_index_a(jlim2)   = cmplx(1.45,0.)
 734     end if
 735     ref_index_a(jh2o)    = cmplx(1.33,0.)
 736
 737     ! jsalt_index
 738     jsalt_index(jnh4so4) = 5           ! AS
 739     jsalt_index(jlvcite) = 2           ! LV
 740     jsalt_index(jnh4hso4)= 1           ! AB
 741     jsalt_index(jnh4no3) = 2           ! AN
 742     jsalt_index(jnh4cl)  = 1           ! AC
 743     jsalt_index(jna2so4) = 60          ! SS
 744     jsalt_index(jnahso4) = 10          ! SB
 745     jsalt_index(jnano3)  = 40          ! SN
 746     jsalt_index(jnacl)   = 10          ! SC
 747     jsalt_index(jcano3)  = 120 ! CN
 748     jsalt_index(jcacl2)  = 80          ! CC
 749     jsalt_index(jnh4msa) = 0           ! AM    zero for now
 750     jsalt_index(jnamsa)  = 0           ! SM    zero for now
 751     jsalt_index(jcamsa2) = 0           ! CM    zero for now
 752
 753     ! Aerosol Indices
 754     !  AC = 1, AN = 2, AS = 5, SC = 10, SN = 40, SS = 60, CC = 80, CN = 120,
 755     !  AB = 1, LV = 2, SB = 10
 756     !
 757     ! SULFATE-POOR DOMAIN
 758     jsulf_poor(1)   =  1       !       AC
 759     jsulf_poor(2)   =  2       !       AN
 760     jsulf_poor(5)   =  3       !       AS
 761     jsulf_poor(10)  =  4       !       SC
 762     jsulf_poor(40)  =  5       !       SN
 763     jsulf_poor(60)  =  6       !       SS
 764     jsulf_poor(80)  =  7       !       CC
 765     jsulf_poor(120) =  8       !       CN
 766     jsulf_poor(3)   =  9       !       AN + AC
 767     jsulf_poor(6)   =  10      !       AS + AC
 768     jsulf_poor(7)   =  11      !       AS + AN
 769     jsulf_poor(8)   =          12      !       AS + AN + AC
 770     jsulf_poor(11)  =  13      !       SC + AC
 771     jsulf_poor(41)  =  14      !       SN + AC
 772     jsulf_poor(42)  =  15      !       SN + AN
 773     jsulf_poor(43)  =  16      !       SN + AN + AC
 774     jsulf_poor(50)  =  17      !       SN + SC
 775     jsulf_poor(51)  =  18      !       SN + SC + AC
 776     jsulf_poor(61)  =  19      !       SS + AC
 777     jsulf_poor(62)  =  20      !       SS + AN
 778     jsulf_poor(63)  =  21      !       SS + AN + AC
 779     jsulf_poor(65)  =  22      !       SS + AS
 780     jsulf_poor(66)  =  23      !       SS + AS + AC
 781     jsulf_poor(67)  =  24      !       SS + AS + AN
 782     jsulf_poor(68)  =  25      !       SS + AS + AN + AC
 783     jsulf_poor(70)  =  26      !       SS + SC
 784     jsulf_poor(71)  =  27      !       SS + SC + AC
 785     jsulf_poor(100) =  28      !       SS + SN
 786     jsulf_poor(101) =  29      !       SS + SN + AC
 787     jsulf_poor(102) =  30      !       SS + SN + AN
 788     jsulf_poor(103) =  31      !       SS + SN + AN + AC
 789     jsulf_poor(110) =  32      !       SS + SN + SC
 790     jsulf_poor(111) =  33      !       SS + SN + SC + AC
 791     jsulf_poor(81)  =  34      !       CC + AC
 792     jsulf_poor(90)  =  35      !       CC + SC
 793     jsulf_poor(91)  =  36      !       CC + SC + AC
 794     jsulf_poor(121) =  37      !       CN + AC
 795     jsulf_poor(122) =  38      !       CN + AN
 796     jsulf_poor(123) =  39      !       CN + AN + AC
 797     jsulf_poor(130) =  40      !       CN + SC
 798     jsulf_poor(131) =  41      !       CN + SC + AC
 799     jsulf_poor(160) =  42      !       CN + SN
 800     jsulf_poor(161) =  43      !       CN + SN + AC
 801     jsulf_poor(162) =  44      !       CN + SN + AN
 802     jsulf_poor(163) =  45      !       CN + SN + AN + AC
 803     jsulf_poor(170) =  46      !       CN + SN + SC
 804     jsulf_poor(171) =  47      !       CN + SN + SC + AC
 805     jsulf_poor(200) =  48      !       CN + CC
 806     jsulf_poor(201) =  49      !       CN + CC + AC
 807     jsulf_poor(210) =  50      !       CN + CC + SC
 808     jsulf_poor(211) =  51      !       CN + CC + SC + AC
 809     !
 810     ! SULFATE-RICH DOMAIN
 811     jsulf_rich(1)   =  52      !       AB
 812     jsulf_rich(2)   =  53      !       LV
 813     jsulf_rich(10)  =  54      !       SB
 814     jsulf_rich(3)   =  55      !       AB + LV
 815     jsulf_rich(7)   =  56      !       AS + LV
 816     jsulf_rich(70)  =  57      !       SS + SB
 817     jsulf_rich(62)  =  58      !       SS + LV
 818     jsulf_rich(67)  =  59      !       SS + AS + LV
 819     jsulf_rich(61)  =  60      !       SS + AB
 820     jsulf_rich(63)  =  61      !       SS + LV + AB
 821     jsulf_rich(11)  =  62      !       SB + AB
 822     jsulf_rich(71)  =  63      !       SS + SB + AB
 823     jsulf_rich(5)   =  3       !       AS
 824     jsulf_rich(60)  =  6       !       SS
 825     jsulf_rich(65)  =  22      !       SS + AS
 826
 827
 828
 829     !
 830     ! polynomial coefficients for binary molality (used in ZSR equation)
 831     !
 832     !
 833     ! a_zsr for aw < 0.97
 834     !
 835     ! (NH4)2SO4
 836     je = jnh4so4
 837     a_zsr(1,je)  =  1.30894
 838     a_zsr(2,je)  = -7.09922
 839     a_zsr(3,je)  =  20.62831
 840     a_zsr(4,je)  = -32.19965
 841     a_zsr(5,je)  =  25.17026
 842     a_zsr(6,je)  = -7.81632
 843     aw_min(je)   = 0.1
 844     !
 845     ! (NH4)3H(SO4)2
 846     je = jlvcite
 847     a_zsr(1,je)  =  1.10725
 848     a_zsr(2,je)  = -5.17978
 849     a_zsr(3,je)  =  12.29534
 850     a_zsr(4,je)  = -16.32545
 851     a_zsr(5,je)  =  11.29274
 852     a_zsr(6,je)  = -3.19164
 853     aw_min(je)   = 0.1
 854     !
 855     ! NH4HSO4
 856     je = jnh4hso4
 857     a_zsr(1,je)  =  1.15510
 858     a_zsr(2,je)  = -3.20815
 859     a_zsr(3,je)  =  2.71141
 860     a_zsr(4,je)  =  2.01155
 861     a_zsr(5,je)  = -4.71014
 862     a_zsr(6,je)  =  2.04616
 863     aw_min(je)   = 0.1
 864     !
 865     ! NH4MSA (assumed same as NH4HSO4)
 866     je = jnh4msa
 867     a_zsr(1,je)  =  1.15510
 868     a_zsr(2,je)  = -3.20815
 869     a_zsr(3,je)  =  2.71141
 870     a_zsr(4,je)  =  2.01155
 871     a_zsr(5,je)  = -4.71014
 872     a_zsr(6,je)  =  2.04616
 873     aw_min(je)   = 0.1
 874     !
 875     ! NH4NO3
 876     je = jnh4no3
 877     a_zsr(1,je)  =  0.43507
 878     a_zsr(2,je)  =  6.38220
 879     a_zsr(3,je)  = -30.19797
 880     a_zsr(4,je)  =  53.36470
 881     a_zsr(5,je)  = -43.44203
 882     a_zsr(6,je)  =  13.46158
 883     aw_min(je)   = 0.1
 884     !
 885     ! NH4Cl: revised on Nov 13, 2003. based on Chan and Ha (1999) JGR.
 886     je = jnh4cl
 887     a_zsr(1,je)  =  0.45309
 888     a_zsr(2,je)  =  2.65606
 889     a_zsr(3,je)  = -14.7730
 890     a_zsr(4,je)  =  26.2936
 891     a_zsr(5,je)  = -20.5735
 892     a_zsr(6,je)  =  5.94255
 893     aw_min(je)   = 0.1
 894     !
 895     ! NaCl
 896     je = jnacl
 897     a_zsr(1,je)  =  0.42922
 898     a_zsr(2,je)  = -1.17718
 899     a_zsr(3,je)  =  2.80208
 900     a_zsr(4,je)  = -4.51097
 901     a_zsr(5,je)  =  3.76963
 902     a_zsr(6,je)  = -1.31359
 903     aw_min(je)   = 0.1
 904     !
 905     ! NaNO3
 906     je = jnano3
 907     a_zsr(1,je)  =  1.34966
 908     a_zsr(2,je)  = -5.20116
 909     a_zsr(3,je)  =  11.49011
 910     a_zsr(4,je)  = -14.41380
 911     a_zsr(5,je)  =  9.07037
 912     a_zsr(6,je)  = -2.29769
 913     aw_min(je)   = 0.1
 914     !
 915     ! Na2SO4
 916     je = jna2so4
 917     a_zsr(1,je)  =  0.39888
 918     a_zsr(2,je)  = -1.27150
 919     a_zsr(3,je)  =  3.42792
 920     a_zsr(4,je)  = -5.92632
 921     a_zsr(5,je)  =  5.33351
 922     a_zsr(6,je)  = -1.96541
 923     aw_min(je)   = 0.1
 924     !
 925     ! Na3H(SO4)2  added on 1/14/2004
 926     je = jna3hso4
 927     a_zsr(1,je)  =  0.31480
 928     a_zsr(2,je)  = -1.01087
 929     a_zsr(3,je)  =  2.44029
 930     a_zsr(4,je)  = -3.66095
 931     a_zsr(5,je)  =  2.77632
 932     a_zsr(6,je)  = -0.86058
 933     aw_min(je)   = 0.1
 934     !
 935     ! NaHSO4
 936     je = jnahso4
 937     a_zsr(1,je)  =  0.62764
 938     a_zsr(2,je)  = -1.63520
 939     a_zsr(3,je)  =  4.62531
 940     a_zsr(4,je)  = -10.06925
 941     a_zsr(5,je)  =  10.33547
 942     a_zsr(6,je)  = -3.88729
 943     aw_min(je)   = 0.1
 944     !
 945     ! NaMSA (assumed same as NaHSO4)
 946     je = jnamsa
 947     a_zsr(1,je)  =  0.62764
 948     a_zsr(2,je)  = -1.63520
 949     a_zsr(3,je)  =  4.62531
 950     a_zsr(4,je)  = -10.06925
 951     a_zsr(5,je)  =  10.33547
 952     a_zsr(6,je)  = -3.88729
 953     aw_min(je)   = 0.1
 954     !
 955     ! Ca(NO3)2
 956     je = jcano3
 957     a_zsr(1,je)  =  0.38895
 958     a_zsr(2,je)  = -1.16013
 959     a_zsr(3,je)  =  2.16819
 960     a_zsr(4,je)  = -2.23079
 961     a_zsr(5,je)  =  1.00268
 962     a_zsr(6,je)  = -0.16923
 963     aw_min(je)   = 0.1
 964     !
 965     ! CaCl2: Kim and Seinfeld
 966     je = jcacl2
 967     a_zsr(1,je)  =  0.29891
 968     a_zsr(2,je)  = -1.31104
 969     a_zsr(3,je)  =  3.68759
 970     a_zsr(4,je)  = -5.81708
 971     a_zsr(5,je)  =  4.67520
 972     a_zsr(6,je)  = -1.53223
 973     aw_min(je)   = 0.1
 974     !
 975     ! H2SO4
 976     je = jh2so4
 977     a_zsr(1,je) =  0.32751
 978     a_zsr(2,je) = -1.00692
 979     a_zsr(3,je) =  2.59750
 980     a_zsr(4,je) = -4.40014
 981     a_zsr(5,je) =  3.88212
 982     a_zsr(6,je) = -1.39916
 983     aw_min(je)  = 0.1
 984     !
 985     ! MSA (assumed same as H2SO4)
 986     je = jmsa
 987     a_zsr(1,je) =  0.32751
 988     a_zsr(2,je) = -1.00692
 989     a_zsr(3,je) =  2.59750
 990     a_zsr(4,je) = -4.40014
 991     a_zsr(5,je) =  3.88212
 992     a_zsr(6,je) = -1.39916
 993     aw_min(je)  = 0.1
 994     !
 995     ! HHSO4
 996     je = jhhso4
 997     a_zsr(1,je) =  0.32751
 998     a_zsr(2,je) = -1.00692
 999     a_zsr(3,je) =  2.59750
1000     a_zsr(4,je) = -4.40014
1001     a_zsr(5,je) =  3.88212
1002     a_zsr(6,je) = -1.39916
1003     aw_min(je)  = 1.0
1004     !
1005     ! HNO3
1006     je = jhno3
1007     a_zsr(1,je) =  0.75876
1008     a_zsr(2,je) = -3.31529
1009     a_zsr(3,je) =  9.26392
1010     a_zsr(4,je) = -14.89799
1011     a_zsr(5,je) =  12.08781
1012     a_zsr(6,je) = -3.89958
1013     aw_min(je)  = 0.1
1014     !
1015     ! HCl
1016     je = jhcl
1017     a_zsr(1,je) =  0.31133
1018     a_zsr(2,je) = -0.79688
1019     a_zsr(3,je) =  1.93995
1020     a_zsr(4,je) = -3.31582
1021     a_zsr(5,je) =  2.93513
1022     a_zsr(6,je) = -1.07268
1023     aw_min(je)  = 0.1
1024     !
1025     ! CaSO4
1026     je = jcaso4
1027     a_zsr(1,je)  =  0.0
1028     a_zsr(2,je)  =  0.0
1029     a_zsr(3,je)  =  0.0
1030     a_zsr(4,je)  =  0.0
1031     a_zsr(5,je)  =  0.0
1032     a_zsr(6,je)  =  0.0
1033     aw_min(je)   = 1.0
1034     !
1035     ! Ca(MSA)2 (assumed same as Ca(NO3)2)
1036     je = jcamsa2
1037     a_zsr(1,je)  =  0.38895
1038     a_zsr(2,je)  = -1.16013
1039     a_zsr(3,je)  =  2.16819
1040     a_zsr(4,je)  = -2.23079
1041     a_zsr(5,je)  =  1.00268
1042     a_zsr(6,je)  = -0.16923
1043     aw_min(je)   = 0.1
1044     !
1045     ! CaCO3
1046     je = jcaco3
1047     a_zsr(1,je)  =  0.0
1048     a_zsr(2,je)  =  0.0
1049     a_zsr(3,je)  =  0.0
1050     a_zsr(4,je)  =  0.0
1051     a_zsr(5,je)  =  0.0
1052     a_zsr(6,je)  =  0.0
1053     aw_min(je)   = 1.0
1054
1055
1056
1057     !-------------------------------------------
1058     ! b_zsr for aw => 0.97 to 0.99999
1059     !
1060     ! (NH4)2SO4
1061     b_zsr(jnh4so4)  = 28.0811
1062     !
1063     ! (NH4)3H(SO4)2
1064     b_zsr(jlvcite)  = 14.7178
1065     !
1066     ! NH4HSO4
1067     b_zsr(jnh4hso4) = 29.4779
1068     !
1069     ! NH4MSA
1070     b_zsr(jnh4msa)  = 29.4779 ! assumed same as NH4HSO4
1071     !
1072     ! NH4NO3
1073     b_zsr(jnh4no3)  = 33.4049
1074     !
1075     ! NH4Cl
1076     b_zsr(jnh4cl)   = 30.8888
1077     !
1078     ! NaCl
1079     b_zsr(jnacl)    = 29.8375
1080     !
1081     ! NaNO3
1082     b_zsr(jnano3)   = 32.2756
1083     !
1084     ! Na2SO4
1085     b_zsr(jna2so4)  = 27.6889
1086     !
1087     ! Na3H(SO4)2
1088     b_zsr(jna3hso4) = 14.2184
1089     !
1090     ! NaHSO4
1091     b_zsr(jnahso4)  = 28.3367
1092     !
1093     ! NaMSA
1094     b_zsr(jnamsa)   = 28.3367 ! assumed same as NaHSO4
1095     !
1096     ! Ca(NO3)2
1097     b_zsr(jcano3)   = 18.3661
1098     !
1099     ! CaCl2
1100     b_zsr(jcacl2)   = 20.8792
1101     !
1102     ! H2SO4
1103     b_zsr(jh2so4)   = 26.7347
1104     !
1105     ! HHSO4
1106     b_zsr(jhhso4)   = 26.7347
1107     !
1108     ! HNO3
1109     b_zsr(jhno3)    = 28.8257
1110     !
1111     ! HCl
1112     b_zsr(jhcl)     = 27.7108
1113     !
1114     ! MSA
1115     b_zsr(jmsa)     = 26.7347 ! assumed same as H2SO4
1116     !
1117     ! CaSO4
1118     b_zsr(jcaso4)   = 0.0
1119     !
1120     ! Ca(MSA)2
1121     b_zsr(jcamsa2)  = 18.3661 ! assumed same as Ca(NO3)2
1122     !
1123     ! CaCO3
1124     b_zsr(jcaco3)   = 0.0
1125
1126
1127
1128
1129
1130
1131
1132
1133
1134     !-------------------------------------------
1135     ! Li and Lu (2001) Surface tension model
1136     ! G_MX [mol/cm^2]; K_MX [-]
1137     !
1138     ! (NH4)2SO4
1139     G_MX(jnh4so4)  = -8.79e-7*1.e-4
1140     K_MX(jnh4so4)  =  3.84e+1
1141     !
1142     ! (NH4)3H(SO4)2
1143     G_MX(jlvcite)  = -8.79e-7*1.e-4    ! assumed same as (NH4)2SO4
1144     K_MX(jlvcite)  =  3.84e+1          ! assumed same as (NH4)2SO4
1145     !
1146     ! NH4HSO4
1147     G_MX(jnh4hso4) = -8.79e-7*1.e-4    ! assumed same as (NH4)2SO4
1148     K_MX(jnh4hso4) =  3.84e+1          ! assumed same as (NH4)2SO4
1149     !
1150     ! NH4MSA
1151     G_MX(jnh4msa)  = -8.79e-7*1.e-4    ! assumed same as (NH4)2SO4
1152     K_MX(jnh4msa)  =  3.84e+1          ! assumed same as (NH4)2SO4
1153     !
1154     ! NH4NO3
1155     G_MX(jnh4no3)  = -3.08e-6*1.e-4
1156     K_MX(jnh4no3)  =  4.89e-1
1157     !
1158     ! NH4Cl
1159     G_MX(jnh4cl)   = -1.01e-6*1.e-4
1160     K_MX(jnh4cl)   =  1.3
1161     !
1162     ! NaCl
1163     G_MX(jnacl)    = -1.05e-6*1.e-4
1164     K_MX(jnacl)    =  1.2
1165     !
1166     ! NaNO3
1167     G_MX(jnano3)   = -1.66e-6*1.e-4
1168     K_MX(jnano3)   =  1.25
1169     !
1170     ! Na2SO4
1171     G_MX(jna2so4)  = -8.37e-7*1.e-4
1172     K_MX(jna2so4)  =  7.57e+1
1173     !
1174     ! Na3H(SO4)2
1175     G_MX(jna3hso4) = -8.37e-7*1.e-4    ! assumed same as Na2SO4
1176     K_MX(jna3hso4) =  7.57e+1          ! assumed same as Na2SO4
1177     !
1178     ! NaHSO4
1179     G_MX(jnahso4)  = -8.37e-7*1.e-4    ! assumed same as Na2SO4
1180     K_MX(jnahso4)  =  7.57e+1          ! assumed same as Na2SO4
1181     !
1182     ! NaMSA
1183     G_MX(jnamsa)   = -8.37e-7*1.e-4
1184     K_MX(jnamsa)   =  7.57e+1
1185     !
1186     ! Ca(NO3)2
1187     G_MX(jcano3)   = -4.88e-7*1.e-4    ! assumed same as CaCl2
1188     K_MX(jcano3)   =  1.50e+1          ! assumed same as CaCl2
1189     !
1190     ! CaCl2
1191     G_MX(jcacl2)   = -4.88e-7*1.e-4
1192     K_MX(jcacl2)   =  1.50e+1
1193     !
1194     ! H2SO4
1195     G_MX(jh2so4)   = -6.75e-8*1.e-4
1196     K_MX(jh2so4)   =  1.65e+3
1197     !
1198     ! HHSO4
1199     G_MX(jh2so4)   = -6.75e-8*1.e-4    ! assumed same as H2SO4
1200     K_MX(jh2so4)   =  1.65e+3          ! assumed same as H2SO4
1201     !
1202     ! HNO3
1203     G_MX(jhno3)    =  8.05e-7*1.e-4
1204     K_MX(jhno3)    =  1.06e-1
1205     !
1206     ! HCl
1207     G_MX(jhcl)     =  4.12e-7*1.e-4
1208     K_MX(jhcl)     =  4.68e-3
1209     !&
1210
1211     ! MSA
1212     G_MX(jmsa)     =  8.05e-7*1.e-4    ! assumed same as HNO3
1213     K_MX(jmsa)     =  1.06e-1          ! assumed same as HNO3
1214     !
1215     ! CaSO4
1216     G_MX(jmsa)     =  0.0*1.e-4        ! assumed
1217     K_MX(jmsa)     =  0.0              ! assumed
1218     !
1219     ! Ca(MSA)2
1220     G_MX(jcamsa2)  =  0.0*1.e-4        ! assumed
1221     K_MX(jcamsa2)  =  0.0              ! assumed
1222     !
1223     ! CaCO3
1224     G_MX(jcaco3)   =  0.0*1.e-4        ! assumed
1225     K_MX(jcaco3)   =  0.0              ! assumed
1226
1227
1228
1229
1230
1231
1232
1233     !----------------------------------------------------------------
1234     ! parameters for MTEM mixing rule (Zaveri, Easter, and Wexler, 2005)
1235     ! log_gamZ(jA,jE)   A in E
1236     !----------------------------------------------------------------
1237     !
1238     b_mtem(:,:,:) = 0.0_r8 !BSINGH - Temporarily initialized, please modify if required *Ask dick about it* the code blows up if i initialize it with nan
1239     ! (NH4)2SO4 in E
1240     jA = jnh4so4
1241
1242     ! in (NH4)2SO4
1243     jE = jnh4so4
1244     b_mtem(1,jA,jE) = -2.94685
1245     b_mtem(2,jA,jE) = 17.3328
1246     b_mtem(3,jA,jE) = -64.8441
1247     b_mtem(4,jA,jE) = 122.7070
1248     b_mtem(5,jA,jE) = -114.4373
1249     b_mtem(6,jA,jE) = 41.6811
1250
1251     ! in NH4NO3
1252     jE = jnh4no3
1253     b_mtem(1,jA,jE) = -2.7503
1254     b_mtem(2,jA,jE) = 4.3806
1255     b_mtem(3,jA,jE) = -1.1110
1256     b_mtem(4,jA,jE) = -1.7005
1257     b_mtem(5,jA,jE) = -4.4207
1258     b_mtem(6,jA,jE) = 5.1990
1259
1260     ! in NH4Cl (revised on 11/15/2003)
1261     jE = jnh4cl
1262     b_mtem(1,jA,jE) = -2.06952
1263     b_mtem(2,jA,jE) = 7.1240
1264     b_mtem(3,jA,jE) = -24.4274
1265     b_mtem(4,jA,jE) = 51.1458
1266     b_mtem(5,jA,jE) = -54.2056
1267     b_mtem(6,jA,jE) = 22.0606
1268
1269     ! in Na2SO4
1270     jE = jna2so4
1271     b_mtem(1,jA,jE) = -2.17361
1272     b_mtem(2,jA,jE) = 15.9919
1273     b_mtem(3,jA,jE) = -69.0952
1274     b_mtem(4,jA,jE) = 139.8860
1275     b_mtem(5,jA,jE) = -134.9890
1276     b_mtem(6,jA,jE) = 49.8877
1277
1278     ! in NaNO3
1279     jE = jnano3
1280     b_mtem(1,jA,jE) = -4.4370
1281     b_mtem(2,jA,jE) = 24.0243
1282     b_mtem(3,jA,jE) = -76.2437
1283     b_mtem(4,jA,jE) = 128.6660
1284     b_mtem(5,jA,jE) = -110.0900
1285     b_mtem(6,jA,jE) = 37.7414
1286
1287     ! in NaCl
1288     jE = jnacl
1289     b_mtem(1,jA,jE) = -1.5394
1290     b_mtem(2,jA,jE) = 5.8671
1291     b_mtem(3,jA,jE) = -22.7726
1292     b_mtem(4,jA,jE) = 47.0547
1293     b_mtem(5,jA,jE) = -47.8266
1294     b_mtem(6,jA,jE) = 18.8489
1295
1296     ! in HNO3
1297     jE = jhno3
1298     b_mtem(1,jA,jE) = -0.35750
1299     b_mtem(2,jA,jE) = -3.82466
1300     b_mtem(3,jA,jE) = 4.55462
1301     b_mtem(4,jA,jE) = 5.05402
1302     b_mtem(5,jA,jE) = -14.7476
1303     b_mtem(6,jA,jE) = 8.8009
1304
1305     ! in HCl
1306     jE = jhcl
1307     b_mtem(1,jA,jE) = -2.15146
1308     b_mtem(2,jA,jE) = 5.50205
1309     b_mtem(3,jA,jE) = -19.1476
1310     b_mtem(4,jA,jE) = 39.1880
1311     b_mtem(5,jA,jE) = -39.9460
1312     b_mtem(6,jA,jE) = 16.0700
1313
1314     ! in H2SO4
1315     jE = jh2so4
1316     b_mtem(1,jA,jE) = -2.52604
1317     b_mtem(2,jA,jE) = 9.76022
1318     b_mtem(3,jA,jE) = -35.2540
1319     b_mtem(4,jA,jE) = 71.2981
1320     b_mtem(5,jA,jE) = -71.8207
1321     b_mtem(6,jA,jE) = 28.0758
1322
1323     !
1324     ! in NH4HSO4
1325     jE = jnh4hso4
1326     b_mtem(1,jA,jE) = -4.13219
1327     b_mtem(2,jA,jE) = 13.8863
1328     b_mtem(3,jA,jE) = -34.5387
1329     b_mtem(4,jA,jE) = 56.5012
1330     b_mtem(5,jA,jE) = -51.8702
1331     b_mtem(6,jA,jE) = 19.6232
1332
1333     !
1334     ! in (NH4)3H(SO4)2
1335     jE = jlvcite
1336     b_mtem(1,jA,jE) = -2.53482
1337     b_mtem(2,jA,jE) = 12.3333
1338     b_mtem(3,jA,jE) = -46.1020
1339     b_mtem(4,jA,jE) = 90.4775
1340     b_mtem(5,jA,jE) = -88.1254
1341     b_mtem(6,jA,jE) = 33.4715
1342
1343     !
1344     ! in NaHSO4
1345     jE = jnahso4
1346     b_mtem(1,jA,jE) = -3.23425
1347     b_mtem(2,jA,jE) = 18.7842
1348     b_mtem(3,jA,jE) = -78.7807
1349     b_mtem(4,jA,jE) = 161.517
1350     b_mtem(5,jA,jE) = -154.940
1351     b_mtem(6,jA,jE) = 56.2252
1352
1353     !
1354     ! in Na3H(SO4)2
1355     jE = jna3hso4
1356     b_mtem(1,jA,jE) = -1.25316
1357     b_mtem(2,jA,jE) = 7.40960
1358     b_mtem(3,jA,jE) = -34.8929
1359     b_mtem(4,jA,jE) = 72.8853
1360     b_mtem(5,jA,jE) = -72.4503
1361     b_mtem(6,jA,jE) = 27.7706
1362
1363
1364     !-----------------
1365     ! NH4NO3 in E
1366     jA = jnh4no3
1367
1368     ! in (NH4)2SO4
1369     jE = jnh4so4
1370     b_mtem(1,jA,jE) = -3.5201
1371     b_mtem(2,jA,jE) = 21.6584
1372     b_mtem(3,jA,jE) = -72.1499
1373     b_mtem(4,jA,jE) = 126.7000
1374     b_mtem(5,jA,jE) = -111.4550
1375     b_mtem(6,jA,jE) = 38.5677
1376
1377     ! in NH4NO3
1378     jE = jnh4no3
1379     b_mtem(1,jA,jE) = -2.2630
1380     b_mtem(2,jA,jE) = -0.1518
1381     b_mtem(3,jA,jE) = 17.0898
1382     b_mtem(4,jA,jE) = -36.7832
1383     b_mtem(5,jA,jE) = 29.8407
1384     b_mtem(6,jA,jE) = -7.9314
1385
1386     ! in NH4Cl (revised on 11/15/2003)
1387     jE = jnh4cl
1388     b_mtem(1,jA,jE) = -1.3851
1389     b_mtem(2,jA,jE) = -0.4462
1390     b_mtem(3,jA,jE) = 8.4567
1391     b_mtem(4,jA,jE) = -11.5988
1392     b_mtem(5,jA,jE) = 2.9802
1393     b_mtem(6,jA,jE) = 1.8132
1394
1395     ! in Na2SO4
1396     jE = jna2so4
1397     b_mtem(1,jA,jE) = -1.7602
1398     b_mtem(2,jA,jE) = 10.4044
1399     b_mtem(3,jA,jE) = -35.5894
1400     b_mtem(4,jA,jE) = 64.3584
1401     b_mtem(5,jA,jE) = -57.8931
1402     b_mtem(6,jA,jE) = 20.2141
1403
1404     ! in NaNO3
1405     jE = jnano3
1406     b_mtem(1,jA,jE) = -3.24346
1407     b_mtem(2,jA,jE) = 16.2794
1408     b_mtem(3,jA,jE) = -48.7601
1409     b_mtem(4,jA,jE) = 79.2246
1410     b_mtem(5,jA,jE) = -65.8169
1411     b_mtem(6,jA,jE) = 22.1500
1412
1413     ! in NaCl
1414     jE = jnacl
1415     b_mtem(1,jA,jE) = -1.75658
1416     b_mtem(2,jA,jE) = 7.71384
1417     b_mtem(3,jA,jE) = -22.7984
1418     b_mtem(4,jA,jE) = 39.1532
1419     b_mtem(5,jA,jE) = -34.6165
1420     b_mtem(6,jA,jE) = 12.1283
1421
1422     ! in Ca(NO3)2
1423     jE = jcano3
1424     b_mtem(1,jA,jE) = -0.97178
1425     b_mtem(2,jA,jE) = 6.61964
1426     b_mtem(3,jA,jE) = -26.2353
1427     b_mtem(4,jA,jE) = 50.5259
1428     b_mtem(5,jA,jE) = -47.6586
1429     b_mtem(6,jA,jE) = 17.5074
1430
1431     ! in CaCl2 added on 12/22/2003
1432     jE = jcacl2
1433     b_mtem(1,jA,jE) = -0.41515
1434     b_mtem(2,jA,jE) = 6.44101
1435     b_mtem(3,jA,jE) = -26.4473
1436     b_mtem(4,jA,jE) = 49.0718
1437     b_mtem(5,jA,jE) = -44.2631
1438     b_mtem(6,jA,jE) = 15.3771
1439
1440     ! in HNO3
1441     jE = jhno3
1442     b_mtem(1,jA,jE) = -1.20644
1443     b_mtem(2,jA,jE) = 5.70117
1444     b_mtem(3,jA,jE) = -18.2783
1445     b_mtem(4,jA,jE) = 31.7199
1446     b_mtem(5,jA,jE) = -27.8703
1447     b_mtem(6,jA,jE) = 9.7299
1448
1449     ! in HCl
1450     jE = jhcl
1451     b_mtem(1,jA,jE) = -0.680862
1452     b_mtem(2,jA,jE) = 3.59456
1453     b_mtem(3,jA,jE) = -10.7969
1454     b_mtem(4,jA,jE) = 17.8434
1455     b_mtem(5,jA,jE) = -15.3165
1456     b_mtem(6,jA,jE) = 5.17123
1457
1458
1459     !----------
1460     ! NH4Cl in E
1461     jA = jnh4cl
1462
1463     ! in (NH4)2SO4
1464     jE = jnh4so4
1465     b_mtem(1,jA,jE) = -2.8850
1466     b_mtem(2,jA,jE) = 20.6970
1467     b_mtem(3,jA,jE) = -70.6810
1468     b_mtem(4,jA,jE) = 124.3690
1469     b_mtem(5,jA,jE) = -109.2880
1470     b_mtem(6,jA,jE) = 37.5831
1471
1472     ! in NH4NO3
1473     jE = jnh4no3
1474     b_mtem(1,jA,jE) = -1.9386
1475     b_mtem(2,jA,jE) = 1.3238
1476     b_mtem(3,jA,jE) = 11.8500
1477     b_mtem(4,jA,jE) = -28.1168
1478     b_mtem(5,jA,jE) = 21.8543
1479     b_mtem(6,jA,jE) = -5.1671
1480
1481     ! in NH4Cl (revised on 11/15/2003)
1482     jE = jnh4cl
1483     b_mtem(1,jA,jE) = -0.9559
1484     b_mtem(2,jA,jE) = 0.8121
1485     b_mtem(3,jA,jE) = 4.3644
1486     b_mtem(4,jA,jE) = -8.9258
1487     b_mtem(5,jA,jE) = 4.2362
1488     b_mtem(6,jA,jE) = 0.2891
1489
1490     ! in Na2SO4
1491     jE = jna2so4
1492     b_mtem(1,jA,jE) = 0.0377
1493     b_mtem(2,jA,jE) = 6.0752
1494     b_mtem(3,jA,jE) = -30.8641
1495     b_mtem(4,jA,jE) = 63.3095
1496     b_mtem(5,jA,jE) = -61.0070
1497     b_mtem(6,jA,jE) = 22.1734
1498
1499     ! in NaNO3
1500     jE = jnano3
1501     b_mtem(1,jA,jE) = -1.8336
1502     b_mtem(2,jA,jE) = 12.8160
1503     b_mtem(3,jA,jE) = -42.3388
1504     b_mtem(4,jA,jE) = 71.1816
1505     b_mtem(5,jA,jE) = -60.5708
1506     b_mtem(6,jA,jE) = 20.5853
1507
1508     ! in NaCl
1509     jE = jnacl
1510     b_mtem(1,jA,jE) = -0.1429
1511     b_mtem(2,jA,jE) = 2.3561
1512     b_mtem(3,jA,jE) = -10.4425
1513     b_mtem(4,jA,jE) = 20.8951
1514     b_mtem(5,jA,jE) = -20.7739
1515     b_mtem(6,jA,jE) = 7.9355
1516
1517     ! in Ca(NO3)2
1518     jE = jcano3
1519     b_mtem(1,jA,jE) = 0.76235
1520     b_mtem(2,jA,jE) = 3.08323
1521     b_mtem(3,jA,jE) = -23.6772
1522     b_mtem(4,jA,jE) = 53.7415
1523     b_mtem(5,jA,jE) = -55.4043
1524     b_mtem(6,jA,jE) = 21.2944
1525
1526     ! in CaCl2 (revised on 11/27/2003)
1527     jE = jcacl2
1528     b_mtem(1,jA,jE) = 1.13864
1529     b_mtem(2,jA,jE) = -0.340539
1530     b_mtem(3,jA,jE) = -8.67025
1531     b_mtem(4,jA,jE) = 22.8008
1532     b_mtem(5,jA,jE) = -24.5181
1533     b_mtem(6,jA,jE) = 9.3663
1534
1535     ! in HNO3
1536     jE = jhno3
1537     b_mtem(1,jA,jE) = 2.42532
1538     b_mtem(2,jA,jE) = -14.1755
1539     b_mtem(3,jA,jE) = 38.804
1540     b_mtem(4,jA,jE) = -58.2437
1541     b_mtem(5,jA,jE) = 43.5431
1542     b_mtem(6,jA,jE) = -12.5824
1543
1544     ! in HCl
1545     jE = jhcl
1546     b_mtem(1,jA,jE) = 0.330337
1547     b_mtem(2,jA,jE) = 0.0778934
1548     b_mtem(3,jA,jE) = -2.30492
1549     b_mtem(4,jA,jE) = 4.73003
1550     b_mtem(5,jA,jE) = -4.80849
1551     b_mtem(6,jA,jE) = 1.78866
1552
1553
1554
1555     !----------
1556     ! Na2SO4 in E
1557     jA = jna2so4
1558
1559     ! in (NH4)2SO4
1560     jE = jnh4so4
1561     b_mtem(1,jA,jE) = -2.6982
1562     b_mtem(2,jA,jE) = 22.9875
1563     b_mtem(3,jA,jE) = -98.9840
1564     b_mtem(4,jA,jE) = 198.0180
1565     b_mtem(5,jA,jE) = -188.7270
1566     b_mtem(6,jA,jE) = 69.0548
1567
1568     ! in NH4NO3
1569     jE = jnh4no3
1570     b_mtem(1,jA,jE) = -2.4844
1571     b_mtem(2,jA,jE) = 6.5420
1572     b_mtem(3,jA,jE) = -9.8998
1573     b_mtem(4,jA,jE) = 11.3884
1574     b_mtem(5,jA,jE) = -13.6842
1575     b_mtem(6,jA,jE) = 7.7411
1576
1577     ! in NH4Cl (revised on 11/15/2003)
1578     jE = jnh4cl
1579     b_mtem(1,jA,jE) = -1.3325
1580     b_mtem(2,jA,jE) = 13.0406
1581     b_mtem(3,jA,jE) = -56.1935
1582     b_mtem(4,jA,jE) = 107.1170
1583     b_mtem(5,jA,jE) = -97.3721
1584     b_mtem(6,jA,jE) = 34.3763
1585
1586     ! in Na2SO4
1587     jE = jna2so4
1588     b_mtem(1,jA,jE) = -1.2832
1589     b_mtem(2,jA,jE) = 12.8526
1590     b_mtem(3,jA,jE) = -62.2087
1591     b_mtem(4,jA,jE) = 130.3876
1592     b_mtem(5,jA,jE) = -128.2627
1593     b_mtem(6,jA,jE) = 48.0340
1594
1595     ! in NaNO3
1596     jE = jnano3
1597     b_mtem(1,jA,jE) = -3.5384
1598     b_mtem(2,jA,jE) = 21.3758
1599     b_mtem(3,jA,jE) = -70.7638
1600     b_mtem(4,jA,jE) = 121.1580
1601     b_mtem(5,jA,jE) = -104.6230
1602     b_mtem(6,jA,jE) = 36.0557
1603
1604
1605     ! in NaCl
1606     jE = jnacl
1607     b_mtem(1,jA,jE) = 0.2175
1608     b_mtem(2,jA,jE) = -0.5648
1609     b_mtem(3,jA,jE) = -8.0288
1610     b_mtem(4,jA,jE) = 25.9734
1611     b_mtem(5,jA,jE) = -32.3577
1612     b_mtem(6,jA,jE) = 14.3924
1613
1614     ! in HNO3
1615     jE = jhno3
1616     b_mtem(1,jA,jE) = -0.309617
1617     b_mtem(2,jA,jE) = -1.82899
1618     b_mtem(3,jA,jE) = -1.5505
1619     b_mtem(4,jA,jE) = 13.3847
1620     b_mtem(5,jA,jE) = -20.1284
1621     b_mtem(6,jA,jE) = 9.93163
1622
1623     ! in HCl
1624     jE = jhcl
1625     b_mtem(1,jA,jE) = -0.259455
1626     b_mtem(2,jA,jE) = -0.819366
1627     b_mtem(3,jA,jE) = -4.28964
1628     b_mtem(4,jA,jE) = 16.4305
1629     b_mtem(5,jA,jE) = -21.8546
1630     b_mtem(6,jA,jE) = 10.3044
1631
1632     ! in H2SO4
1633     jE = jh2so4
1634     b_mtem(1,jA,jE) = -1.84257
1635     b_mtem(2,jA,jE) = 7.85788
1636     b_mtem(3,jA,jE) = -29.9275
1637     b_mtem(4,jA,jE) = 61.7515
1638     b_mtem(5,jA,jE) = -63.2308
1639     b_mtem(6,jA,jE) = 24.9542
1640
1641     ! in NH4HSO4
1642     jE = jnh4hso4
1643     b_mtem(1,jA,jE) = -1.05891
1644     b_mtem(2,jA,jE) = 2.84831
1645     b_mtem(3,jA,jE) = -21.1827
1646     b_mtem(4,jA,jE) = 57.5175
1647     b_mtem(5,jA,jE) = -64.8120
1648     b_mtem(6,jA,jE) = 26.1986
1649
1650     ! in (NH4)3H(SO4)2
1651     jE = jlvcite
1652     b_mtem(1,jA,jE) = -1.16584
1653     b_mtem(2,jA,jE) = 8.50075
1654     b_mtem(3,jA,jE) = -44.3420
1655     b_mtem(4,jA,jE) = 97.3974
1656     b_mtem(5,jA,jE) = -98.4549
1657     b_mtem(6,jA,jE) = 37.6104
1658
1659     ! in NaHSO4
1660     jE = jnahso4
1661     b_mtem(1,jA,jE) = -1.95805
1662     b_mtem(2,jA,jE) = 6.62417
1663     b_mtem(3,jA,jE) = -31.8072
1664     b_mtem(4,jA,jE) = 77.8603
1665     b_mtem(5,jA,jE) = -84.6458
1666     b_mtem(6,jA,jE) = 33.4963
1667
1668     ! in Na3H(SO4)2
1669     jE = jna3hso4
1670     b_mtem(1,jA,jE) = -0.36045
1671     b_mtem(2,jA,jE) = 3.55223
1672     b_mtem(3,jA,jE) = -24.0327
1673     b_mtem(4,jA,jE) = 54.4879
1674     b_mtem(5,jA,jE) = -56.6531
1675     b_mtem(6,jA,jE) = 22.4956
1676
1677
1678     !----------
1679     ! NaNO3 in E
1680     jA = jnano3
1681
1682     ! in (NH4)2SO4
1683     jE = jnh4so4
1684     b_mtem(1,jA,jE) = -2.5888
1685     b_mtem(2,jA,jE) = 17.6192
1686     b_mtem(3,jA,jE) = -63.2183
1687     b_mtem(4,jA,jE) = 115.3520
1688     b_mtem(5,jA,jE) = -104.0860
1689     b_mtem(6,jA,jE) = 36.7390
1690
1691     ! in NH4NO3
1692     jE = jnh4no3
1693
1694     b_mtem(1,jA,jE) = -2.0669
1695     b_mtem(2,jA,jE) = 1.4792
1696     b_mtem(3,jA,jE) = 10.5261
1697     b_mtem(4,jA,jE) = -27.0987
1698     b_mtem(5,jA,jE) = 23.0591
1699     b_mtem(6,jA,jE) = -6.0938
1700
1701     ! in NH4Cl (revised on 11/15/2003)
1702     jE = jnh4cl
1703     b_mtem(1,jA,jE) = -0.8325
1704     b_mtem(2,jA,jE) = 3.9933
1705     b_mtem(3,jA,jE) = -15.3789
1706     b_mtem(4,jA,jE) = 30.4050
1707     b_mtem(5,jA,jE) = -29.4204
1708     b_mtem(6,jA,jE) = 11.0597
1709
1710     ! in Na2SO4
1711     jE = jna2so4
1712     b_mtem(1,jA,jE) = -1.1233
1713     b_mtem(2,jA,jE) = 8.3998
1714     b_mtem(3,jA,jE) = -31.9002
1715     b_mtem(4,jA,jE) = 60.1450
1716     b_mtem(5,jA,jE) = -55.5503
1717     b_mtem(6,jA,jE) = 19.7757
1718
1719     ! in NaNO3
1720     jE = jnano3
1721     b_mtem(1,jA,jE) = -2.5386
1722     b_mtem(2,jA,jE) = 13.9039
1723     b_mtem(3,jA,jE) = -42.8467
1724     b_mtem(4,jA,jE) = 69.7442
1725     b_mtem(5,jA,jE) = -57.8988
1726     b_mtem(6,jA,jE) = 19.4635
1727
1728     ! in NaCl
1729     jE = jnacl
1730     b_mtem(1,jA,jE) = -0.4351
1731     b_mtem(2,jA,jE) = 2.8311
1732     b_mtem(3,jA,jE) = -11.4485
1733     b_mtem(4,jA,jE) = 22.7201
1734     b_mtem(5,jA,jE) = -22.4228
1735     b_mtem(6,jA,jE) = 8.5792
1736
1737     ! in Ca(NO3)2
1738     jE = jcano3
1739     b_mtem(1,jA,jE) = -0.72060
1740     b_mtem(2,jA,jE) = 5.64915
1741     b_mtem(3,jA,jE) = -23.5020
1742     b_mtem(4,jA,jE) = 46.0078
1743     b_mtem(5,jA,jE) = -43.8075
1744     b_mtem(6,jA,jE) = 16.1652
1745
1746     ! in CaCl2
1747     jE = jcacl2
1748
1749     b_mtem(1,jA,jE) = 0.003928
1750     b_mtem(2,jA,jE) = 3.54724
1751     b_mtem(3,jA,jE) = -18.6057
1752     b_mtem(4,jA,jE) = 38.1445
1753     b_mtem(5,jA,jE) = -36.7745
1754     b_mtem(6,jA,jE) = 13.4529
1755
1756     ! in HNO3
1757     jE = jhno3
1758     b_mtem(1,jA,jE) = -1.1712
1759     b_mtem(2,jA,jE) = 7.20907
1760     b_mtem(3,jA,jE) = -22.9215
1761     b_mtem(4,jA,jE) = 38.1257
1762     b_mtem(5,jA,jE) = -32.0759
1763     b_mtem(6,jA,jE) = 10.6443
1764
1765     ! in HCl
1766     jE = jhcl
1767     b_mtem(1,jA,jE) = 0.738022
1768     b_mtem(2,jA,jE) = -1.14313
1769     b_mtem(3,jA,jE) = 0.32251
1770     b_mtem(4,jA,jE) = 0.838679
1771     b_mtem(5,jA,jE) = -1.81747
1772     b_mtem(6,jA,jE) = 0.873986
1773
1774
1775     !----------
1776     ! NaCl in E
1777     jA = jnacl
1778
1779     ! in (NH4)2SO4
1780     jE = jnh4so4
1781     b_mtem(1,jA,jE) = -1.9525
1782     b_mtem(2,jA,jE) = 16.6433
1783     b_mtem(3,jA,jE) = -61.7090
1784     b_mtem(4,jA,jE) = 112.9910
1785     b_mtem(5,jA,jE) = -101.9370
1786     b_mtem(6,jA,jE) = 35.7760
1787
1788     ! in NH4NO3
1789     jE = jnh4no3
1790     b_mtem(1,jA,jE) = -1.7525
1791     b_mtem(2,jA,jE) = 3.0713
1792     b_mtem(3,jA,jE) = 4.8063
1793     b_mtem(4,jA,jE) = -17.5334
1794     b_mtem(5,jA,jE) = 14.2872
1795     b_mtem(6,jA,jE) = -3.0690
1796
1797     ! in NH4Cl (revised on 11/15/2003)
1798     jE = jnh4cl
1799     b_mtem(1,jA,jE) = -0.4021
1800     b_mtem(2,jA,jE) = 5.2399
1801     b_mtem(3,jA,jE) = -19.4278
1802     b_mtem(4,jA,jE) = 33.0027
1803     b_mtem(5,jA,jE) = -28.1020
1804     b_mtem(6,jA,jE) = 9.5159
1805
1806     ! in Na2SO4
1807     jE = jna2so4
1808     b_mtem(1,jA,jE) = 0.6692
1809     b_mtem(2,jA,jE) = 4.1207
1810     b_mtem(3,jA,jE) = -27.3314
1811     b_mtem(4,jA,jE) = 59.3112
1812     b_mtem(5,jA,jE) = -58.7998
1813     b_mtem(6,jA,jE) = 21.7674
1814
1815     ! in NaNO3
1816     jE = jnano3
1817     b_mtem(1,jA,jE) = -1.17444
1818     b_mtem(2,jA,jE) = 10.9927
1819     b_mtem(3,jA,jE) = -38.9013
1820     b_mtem(4,jA,jE) = 66.8521
1821     b_mtem(5,jA,jE) = -57.6564
1822     b_mtem(6,jA,jE) = 19.7296
1823
1824     ! in NaCl
1825     jE = jnacl
1826     b_mtem(1,jA,jE) = 1.17679
1827     b_mtem(2,jA,jE) = -2.5061
1828     b_mtem(3,jA,jE) = 0.8508
1829     b_mtem(4,jA,jE) = 4.4802
1830     b_mtem(5,jA,jE) = -8.4945
1831     b_mtem(6,jA,jE) = 4.3182
1832
1833     ! in Ca(NO3)2
1834     jE = jcano3
1835     b_mtem(1,jA,jE) = 1.01450
1836     b_mtem(2,jA,jE) = 2.10260
1837     b_mtem(3,jA,jE) = -20.9036
1838     b_mtem(4,jA,jE) = 49.1481
1839     b_mtem(5,jA,jE) = -51.4867
1840     b_mtem(6,jA,jE) = 19.9301
1841
1842     ! in CaCl2 (PSC92: revised on 11/27/2003)
1843     jE = jcacl2
1844     b_mtem(1,jA,jE) = 1.55463
1845     b_mtem(2,jA,jE) = -3.20122
1846     b_mtem(3,jA,jE) = -0.957075
1847     b_mtem(4,jA,jE) = 12.103
1848     b_mtem(5,jA,jE) = -17.221
1849     b_mtem(6,jA,jE) = 7.50264
1850
1851     ! in HNO3
1852     jE = jhno3
1853     b_mtem(1,jA,jE) = 2.46187
1854     b_mtem(2,jA,jE) = -12.6845
1855     b_mtem(3,jA,jE) = 34.2383
1856     b_mtem(4,jA,jE) = -51.9992
1857     b_mtem(5,jA,jE) = 39.4934
1858     b_mtem(6,jA,jE) = -11.7247
1859
1860     ! in HCl
1861     jE = jhcl
1862     b_mtem(1,jA,jE) = 1.74915
1863     b_mtem(2,jA,jE) = -4.65768
1864     b_mtem(3,jA,jE) = 8.80287
1865     b_mtem(4,jA,jE) = -12.2503
1866     b_mtem(5,jA,jE) = 8.668751
1867     b_mtem(6,jA,jE) = -2.50158
1868
1869
1870     !----------
1871     ! Ca(NO3)2 in E
1872     jA = jcano3
1873
1874     ! in NH4NO3
1875     jE = jnh4no3
1876     b_mtem(1,jA,jE) = -1.86260
1877     b_mtem(2,jA,jE) = 11.6178
1878     b_mtem(3,jA,jE) = -30.9069
1879     b_mtem(4,jA,jE) = 41.7578
1880     b_mtem(5,jA,jE) = -33.7338
1881     b_mtem(6,jA,jE) = 12.7541
1882
1883     ! in NH4Cl (revised on 11/15/2003)
1884     jE = jnh4cl
1885     b_mtem(1,jA,jE) = -1.1798
1886     b_mtem(2,jA,jE) = 25.9608
1887     b_mtem(3,jA,jE) = -98.9373
1888     b_mtem(4,jA,jE) = 160.2300
1889     b_mtem(5,jA,jE) = -125.9540
1890     b_mtem(6,jA,jE) = 39.5130
1891
1892     ! in NaNO3
1893     jE = jnano3
1894     b_mtem(1,jA,jE) = -1.44384
1895     b_mtem(2,jA,jE) = 13.6044
1896     b_mtem(3,jA,jE) = -54.4300
1897     b_mtem(4,jA,jE) = 100.582
1898     b_mtem(5,jA,jE) = -91.2364
1899     b_mtem(6,jA,jE) = 32.5970
1900
1901     ! in NaCl
1902     jE = jnacl
1903     b_mtem(1,jA,jE) = -0.099114
1904     b_mtem(2,jA,jE) = 2.84091
1905     b_mtem(3,jA,jE) = -16.9229
1906     b_mtem(4,jA,jE) = 37.4839
1907     b_mtem(5,jA,jE) = -39.5132
1908     b_mtem(6,jA,jE) = 15.8564
1909
1910     ! in Ca(NO3)2
1911     jE = jcano3
1912     b_mtem(1,jA,jE) = 0.055116
1913     b_mtem(2,jA,jE) = 4.58610
1914     b_mtem(3,jA,jE) = -27.6629
1915     b_mtem(4,jA,jE) = 60.8288
1916     b_mtem(5,jA,jE) = -61.4988
1917     b_mtem(6,jA,jE) = 23.3136
1918
1919     ! in CaCl2 (PSC92: revised on 11/27/2003)
1920     jE = jcacl2
1921     b_mtem(1,jA,jE) = 1.57155
1922     b_mtem(2,jA,jE) = -3.18486
1923     b_mtem(3,jA,jE) = -3.35758
1924     b_mtem(4,jA,jE) = 18.7501
1925     b_mtem(5,jA,jE) = -24.5604
1926     b_mtem(6,jA,jE) = 10.3798
1927
1928     ! in HNO3
1929     jE = jhno3
1930     b_mtem(1,jA,jE) = 1.04446
1931     b_mtem(2,jA,jE) = -3.19066
1932     b_mtem(3,jA,jE) = 2.44714
1933     b_mtem(4,jA,jE) = 2.07218
1934     b_mtem(5,jA,jE) = -6.43949
1935     b_mtem(6,jA,jE) = 3.66471
1936
1937     ! in HCl
1938     jE = jhcl
1939     b_mtem(1,jA,jE) = 1.05723
1940     b_mtem(2,jA,jE) = -1.46826
1941     b_mtem(3,jA,jE) = -1.0713
1942     b_mtem(4,jA,jE) = 4.64439
1943     b_mtem(5,jA,jE) = -6.32402
1944     b_mtem(6,jA,jE) = 2.78202
1945
1946
1947     !----------
1948     ! CaCl2 in E
1949     jA = jcacl2
1950
1951     ! in NH4NO3 (PSC92: revised on 12/22/2003)
1952     jE = jnh4no3
1953     b_mtem(1,jA,jE) = -1.43626
1954     b_mtem(2,jA,jE) = 13.6598
1955     b_mtem(3,jA,jE) = -38.2068
1956     b_mtem(4,jA,jE) = 53.9057
1957     b_mtem(5,jA,jE) = -44.9018
1958     b_mtem(6,jA,jE) = 16.6120
1959
1960     ! in NH4Cl (PSC92: revised on 11/27/2003)
1961     jE = jnh4cl
1962     b_mtem(1,jA,jE) = -0.603965
1963     b_mtem(2,jA,jE) = 27.6027
1964     b_mtem(3,jA,jE) = -104.258
1965     b_mtem(4,jA,jE) = 163.553
1966     b_mtem(5,jA,jE) = -124.076
1967     b_mtem(6,jA,jE) = 37.4153
1968
1969     ! in NaNO3 (PSC92: revised on 12/22/2003)
1970     jE = jnano3
1971     b_mtem(1,jA,jE) = 0.44648
1972     b_mtem(2,jA,jE) = 8.8850
1973     b_mtem(3,jA,jE) = -45.5232
1974     b_mtem(4,jA,jE) = 89.3263
1975     b_mtem(5,jA,jE) = -83.8604
1976     b_mtem(6,jA,jE) = 30.4069
1977
1978     ! in NaCl (PSC92: revised on 11/27/2003)
1979     jE = jnacl
1980     b_mtem(1,jA,jE) = 1.61927
1981     b_mtem(2,jA,jE) = 0.247547
1982     b_mtem(3,jA,jE) = -18.1252
1983     b_mtem(4,jA,jE) = 45.2479
1984     b_mtem(5,jA,jE) = -48.6072
1985     b_mtem(6,jA,jE) = 19.2784
1986
1987     ! in Ca(NO3)2 (PSC92: revised on 11/27/2003)
1988     jE = jcano3
1989     b_mtem(1,jA,jE) = 2.36667
1990     b_mtem(2,jA,jE) = -0.123309
1991     b_mtem(3,jA,jE) = -24.2723
1992     b_mtem(4,jA,jE) = 65.1486
1993     b_mtem(5,jA,jE) = -71.8504
1994     b_mtem(6,jA,jE) = 28.3696
1995
1996     ! in CaCl2 (PSC92: revised on 11/27/2003)
1997     jE = jcacl2
1998     b_mtem(1,jA,jE) = 3.64023
1999     b_mtem(2,jA,jE) = -12.1926
2000     b_mtem(3,jA,jE) = 20.2028
2001     b_mtem(4,jA,jE) = -16.0056
2002     b_mtem(5,jA,jE) = 1.52355
2003     b_mtem(6,jA,jE) = 2.44709
2004
2005     ! in HNO3
2006     jE = jhno3
2007     b_mtem(1,jA,jE) = 5.88794
2008     b_mtem(2,jA,jE) = -29.7083
2009     b_mtem(3,jA,jE) = 78.6309
2010     b_mtem(4,jA,jE) = -118.037
2011     b_mtem(5,jA,jE) = 88.932
2012     b_mtem(6,jA,jE) = -26.1407
2013
2014     ! in HCl
2015     jE = jhcl
2016     b_mtem(1,jA,jE) = 2.40628
2017     b_mtem(2,jA,jE) = -6.16566
2018     b_mtem(3,jA,jE) = 10.2851
2019     b_mtem(4,jA,jE) = -12.9035
2020     b_mtem(5,jA,jE) = 7.7441
2021     b_mtem(6,jA,jE) = -1.74821
2022
2023
2024     !----------
2025     ! HNO3 in E
2026     jA = jhno3
2027
2028     ! in (NH4)2SO4
2029     jE = jnh4so4
2030     b_mtem(1,jA,jE) = -3.57598
2031     b_mtem(2,jA,jE) = 21.5469
2032     b_mtem(3,jA,jE) = -77.4111
2033     b_mtem(4,jA,jE) = 144.136
2034     b_mtem(5,jA,jE) = -132.849
2035     b_mtem(6,jA,jE) = 47.9412
2036
2037     ! in NH4NO3
2038     jE = jnh4no3
2039     b_mtem(1,jA,jE) = -2.00209
2040     b_mtem(2,jA,jE) = -3.48399
2041     b_mtem(3,jA,jE) = 34.9906
2042     b_mtem(4,jA,jE) = -68.6653
2043     b_mtem(5,jA,jE) = 54.0992
2044     b_mtem(6,jA,jE) = -15.1343
2045
2046     ! in NH4Cl revised on 12/22/2003
2047     jE = jnh4cl
2048     b_mtem(1,jA,jE) = -0.63790
2049     b_mtem(2,jA,jE) = -1.67730
2050     b_mtem(3,jA,jE) = 10.1727
2051     b_mtem(4,jA,jE) = -14.9097
2052     b_mtem(5,jA,jE) = 7.67410
2053     b_mtem(6,jA,jE) = -0.79586
2054
2055     ! in NaCl
2056     jE = jnacl
2057     b_mtem(1,jA,jE) = 1.3446
2058     b_mtem(2,jA,jE) = -2.5578
2059     b_mtem(3,jA,jE) = 1.3464
2060     b_mtem(4,jA,jE) = 2.90537
2061     b_mtem(5,jA,jE) = -6.53014
2062     b_mtem(6,jA,jE) = 3.31339
2063
2064     ! in NaNO3
2065     jE = jnano3
2066     b_mtem(1,jA,jE) = -0.546636
2067     b_mtem(2,jA,jE) = 10.3127
2068     b_mtem(3,jA,jE) = -39.9603
2069     b_mtem(4,jA,jE) = 71.4609
2070     b_mtem(5,jA,jE) = -63.4958
2071     b_mtem(6,jA,jE) = 22.0679
2072
2073     ! in Na2SO4
2074     jE = jna2so4
2075     b_mtem(1,jA,jE) = 1.35059
2076     b_mtem(2,jA,jE) = 4.34557
2077     b_mtem(3,jA,jE) = -35.8425
2078     b_mtem(4,jA,jE) = 80.9868
2079     b_mtem(5,jA,jE) = -81.6544
2080     b_mtem(6,jA,jE) = 30.4841
2081
2082     ! in Ca(NO3)2
2083     jE = jcano3
2084     b_mtem(1,jA,jE) = 0.869414
2085     b_mtem(2,jA,jE) = 2.98486
2086     b_mtem(3,jA,jE) = -22.255
2087     b_mtem(4,jA,jE) = 50.1863
2088     b_mtem(5,jA,jE) = -51.214
2089     b_mtem(6,jA,jE) = 19.2235
2090
2091     ! in CaCl2 (KM) revised on 12/22/2003
2092     jE = jcacl2
2093     b_mtem(1,jA,jE) = 1.42800
2094     b_mtem(2,jA,jE) = -1.78959
2095     b_mtem(3,jA,jE) = -2.49075
2096     b_mtem(4,jA,jE) = 10.1877
2097     b_mtem(5,jA,jE) = -12.1948
2098     b_mtem(6,jA,jE) = 4.64475
2099
2100     ! in HNO3 (added on 12/06/2004)
2101     jE = jhno3
2102     b_mtem(1,jA,jE) = 0.22035
2103     b_mtem(2,jA,jE) = 2.94973
2104     b_mtem(3,jA,jE) = -12.1469
2105     b_mtem(4,jA,jE) = 20.4905
2106     b_mtem(5,jA,jE) = -17.3966
2107     b_mtem(6,jA,jE) = 5.70779
2108
2109     ! in HCl (added on 12/06/2004)
2110     jE = jhcl
2111     b_mtem(1,jA,jE) = 1.55503
2112     b_mtem(2,jA,jE) = -3.61226
2113     b_mtem(3,jA,jE) = 6.28265
2114     b_mtem(4,jA,jE) = -8.69575
2115     b_mtem(5,jA,jE) = 6.09372
2116     b_mtem(6,jA,jE) = -1.80898
2117
2118     ! in H2SO4
2119     jE = jh2so4
2120     b_mtem(1,jA,jE) = 1.10783
2121     b_mtem(2,jA,jE) = -1.3363
2122     b_mtem(3,jA,jE) = -1.83525
2123     b_mtem(4,jA,jE) = 7.47373
2124     b_mtem(5,jA,jE) = -9.72954
2125     b_mtem(6,jA,jE) = 4.12248
2126
2127     ! in NH4HSO4
2128     jE = jnh4hso4
2129     b_mtem(1,jA,jE) = -0.851026
2130     b_mtem(2,jA,jE) = 12.2515
2131     b_mtem(3,jA,jE) = -49.788
2132     b_mtem(4,jA,jE) = 91.6215
2133     b_mtem(5,jA,jE) = -81.4877
2134     b_mtem(6,jA,jE) = 28.0002
2135
2136     ! in (NH4)3H(SO4)2
2137     jE = jlvcite
2138     b_mtem(1,jA,jE) = -3.09464
2139     b_mtem(2,jA,jE) = 14.9303
2140     b_mtem(3,jA,jE) = -43.0454
2141     b_mtem(4,jA,jE) = 72.6695
2142     b_mtem(5,jA,jE) = -65.2140
2143     b_mtem(6,jA,jE) = 23.4814
2144
2145     ! in NaHSO4
2146     jE = jnahso4
2147     b_mtem(1,jA,jE) = 1.22973
2148     b_mtem(2,jA,jE) = 2.82702
2149     b_mtem(3,jA,jE) = -17.5869
2150     b_mtem(4,jA,jE) = 28.9564
2151     b_mtem(5,jA,jE) = -23.5814
2152     b_mtem(6,jA,jE) = 7.91153
2153
2154     ! in Na3H(SO4)2
2155     jE = jna3hso4
2156     b_mtem(1,jA,jE) = 1.64773
2157     b_mtem(2,jA,jE) = 0.94188
2158     b_mtem(3,jA,jE) = -19.1242
2159     b_mtem(4,jA,jE) = 46.9887
2160     b_mtem(5,jA,jE) = -50.9494
2161     b_mtem(6,jA,jE) = 20.2169
2162
2163
2164     !----------
2165     ! HCl in E
2166     jA = jhcl
2167
2168     ! in (NH4)2SO4
2169     jE = jnh4so4
2170     b_mtem(1,jA,jE) = -2.93783
2171     b_mtem(2,jA,jE) = 20.5546
2172     b_mtem(3,jA,jE) = -75.8548
2173     b_mtem(4,jA,jE) = 141.729
2174     b_mtem(5,jA,jE) = -130.697
2175     b_mtem(6,jA,jE) = 46.9905
2176
2177     ! in NH4NO3
2178     jE = jnh4no3
2179     b_mtem(1,jA,jE) = -1.69063
2180     b_mtem(2,jA,jE) = -1.85303
2181     b_mtem(3,jA,jE) = 29.0927
2182     b_mtem(4,jA,jE) = -58.7401
2183     b_mtem(5,jA,jE) = 44.999
2184     b_mtem(6,jA,jE) = -11.9988
2185
2186     ! in NH4Cl (revised on 11/15/2003)
2187     jE = jnh4cl
2188     b_mtem(1,jA,jE) = -0.2073
2189     b_mtem(2,jA,jE) = -0.4322
2190     b_mtem(3,jA,jE) = 6.1271
2191     b_mtem(4,jA,jE) = -12.3146
2192     b_mtem(5,jA,jE) = 8.9919
2193     b_mtem(6,jA,jE) = -2.3388
2194
2195     ! in NaCl
2196     jE = jnacl
2197     b_mtem(1,jA,jE) = 2.95913
2198     b_mtem(2,jA,jE) = -7.92254
2199     b_mtem(3,jA,jE) = 13.736
2200     b_mtem(4,jA,jE) = -15.433
2201     b_mtem(5,jA,jE) = 7.40386
2202     b_mtem(6,jA,jE) = -0.918641
2203
2204     ! in NaNO3
2205     jE = jnano3
2206     b_mtem(1,jA,jE) = 0.893272
2207     b_mtem(2,jA,jE) = 6.53768
2208     b_mtem(3,jA,jE) = -32.3458
2209     b_mtem(4,jA,jE) = 61.2834
2210     b_mtem(5,jA,jE) = -56.4446
2211     b_mtem(6,jA,jE) = 19.9202
2212
2213     ! in Na2SO4
2214     jE = jna2so4
2215     b_mtem(1,jA,jE) = 3.14484
2216     b_mtem(2,jA,jE) = 0.077019
2217     b_mtem(3,jA,jE) = -31.4199
2218     b_mtem(4,jA,jE) = 80.5865
2219     b_mtem(5,jA,jE) = -85.392
2220     b_mtem(6,jA,jE) = 32.6644
2221
2222     ! in Ca(NO3)2
2223     jE = jcano3
2224     b_mtem(1,jA,jE) = 2.60432
2225     b_mtem(2,jA,jE) = -0.55909
2226     b_mtem(3,jA,jE) = -19.6671
2227     b_mtem(4,jA,jE) = 53.3446
2228     b_mtem(5,jA,jE) = -58.9076
2229     b_mtem(6,jA,jE) = 22.9927
2230
2231     ! in CaCl2 (KM) revised on 3/13/2003 and again on 11/27/2003
2232     jE = jcacl2
2233     b_mtem(1,jA,jE) = 2.98036
2234     b_mtem(2,jA,jE) = -8.55365
2235     b_mtem(3,jA,jE) = 15.2108
2236     b_mtem(4,jA,jE) = -15.9359
2237     b_mtem(5,jA,jE) = 7.41772
2238     b_mtem(6,jA,jE) = -1.32143
2239
2240     ! in HNO3 (added on 12/06/2004)
2241     jE = jhno3
2242     b_mtem(1,jA,jE) = 3.8533
2243     b_mtem(2,jA,jE) = -16.9427
2244     b_mtem(3,jA,jE) = 45.0056
2245     b_mtem(4,jA,jE) = -69.6145
2246     b_mtem(5,jA,jE) = 54.1491
2247     b_mtem(6,jA,jE) = -16.6513
2248
2249     ! in HCl (added on 12/06/2004)
2250     jE = jhcl
2251     b_mtem(1,jA,jE) = 2.56665
2252     b_mtem(2,jA,jE) = -7.13585
2253     b_mtem(3,jA,jE) = 14.8103
2254     b_mtem(4,jA,jE) = -21.8881
2255     b_mtem(5,jA,jE) = 16.6808
2256     b_mtem(6,jA,jE) = -5.22091
2257
2258     ! in H2SO4
2259     jE = jh2so4
2260     b_mtem(1,jA,jE) = 2.50179
2261     b_mtem(2,jA,jE) = -6.69364
2262     b_mtem(3,jA,jE) = 11.6551
2263     b_mtem(4,jA,jE) = -13.6897
2264     b_mtem(5,jA,jE) = 7.36796
2265     b_mtem(6,jA,jE) = -1.33245
2266
2267     ! in NH4HSO4
2268     jE = jnh4hso4
2269     b_mtem(1,jA,jE) = 0.149955
2270     b_mtem(2,jA,jE) = 11.8213
2271     b_mtem(3,jA,jE) = -53.9164
2272     b_mtem(4,jA,jE) = 101.574
2273     b_mtem(5,jA,jE) = -91.4123
2274     b_mtem(6,jA,jE) = 31.5487
2275
2276     ! in (NH4)3H(SO4)2
2277     jE = jlvcite
2278     b_mtem(1,jA,jE) = -2.36927
2279     b_mtem(2,jA,jE) = 14.8359
2280     b_mtem(3,jA,jE) = -44.3443
2281     b_mtem(4,jA,jE) = 73.6229
2282     b_mtem(5,jA,jE) = -65.3366
2283     b_mtem(6,jA,jE) = 23.3250
2284
2285     ! in NaHSO4
2286     jE = jnahso4
2287     b_mtem(1,jA,jE) = 2.72993
2288     b_mtem(2,jA,jE) = -0.23406
2289     b_mtem(3,jA,jE) = -10.4103
2290     b_mtem(4,jA,jE) = 13.1586
2291     b_mtem(5,jA,jE) = -7.79925
2292     b_mtem(6,jA,jE) = 2.30843
2293
2294     ! in Na3H(SO4)2
2295     jE = jna3hso4
2296     b_mtem(1,jA,jE) = 3.51258
2297     b_mtem(2,jA,jE) = -3.95107
2298     b_mtem(3,jA,jE) = -11.0175
2299     b_mtem(4,jA,jE) = 38.8617
2300     b_mtem(5,jA,jE) = -48.1575
2301     b_mtem(6,jA,jE) = 20.4717
2302
2303
2304     !----------
2305     ! 2H.SO4 in E
2306     jA = jh2so4
2307
2308     ! in H2SO4
2309     jE = jh2so4
2310     b_mtem(1,jA,jE) = 0.76734
2311     b_mtem(2,jA,jE) = -1.12263
2312     b_mtem(3,jA,jE) = -9.08728
2313     b_mtem(4,jA,jE) = 30.3836
2314     b_mtem(5,jA,jE) = -38.4133
2315     b_mtem(6,jA,jE) = 17.0106
2316
2317     ! in NH4HSO4
2318     jE = jnh4hso4
2319     b_mtem(1,jA,jE) = -2.03879
2320     b_mtem(2,jA,jE) = 15.7033
2321     b_mtem(3,jA,jE) = -58.7363
2322     b_mtem(4,jA,jE) = 109.242
2323     b_mtem(5,jA,jE) = -102.237
2324     b_mtem(6,jA,jE) = 37.5350
2325
2326     ! in (NH4)3H(SO4)2
2327     jE = jlvcite
2328     b_mtem(1,jA,jE) = -3.10228
2329     b_mtem(2,jA,jE) = 16.6920
2330     b_mtem(3,jA,jE) = -59.1522
2331     b_mtem(4,jA,jE) = 113.487
2332     b_mtem(5,jA,jE) = -110.890
2333     b_mtem(6,jA,jE) = 42.4578
2334
2335     ! in (NH4)2SO4
2336     jE = jnh4so4
2337     b_mtem(1,jA,jE) = -3.43885
2338     b_mtem(2,jA,jE) = 21.0372
2339     b_mtem(3,jA,jE) = -84.7026
2340     b_mtem(4,jA,jE) = 165.324
2341     b_mtem(5,jA,jE) = -156.101
2342     b_mtem(6,jA,jE) = 57.3101
2343
2344     ! in NaHSO4
2345     jE = jnahso4
2346     b_mtem(1,jA,jE) = 0.33164
2347     b_mtem(2,jA,jE) = 6.55864
2348     b_mtem(3,jA,jE) = -33.5876
2349     b_mtem(4,jA,jE) = 65.1798
2350     b_mtem(5,jA,jE) = -63.2046
2351     b_mtem(6,jA,jE) = 24.1783
2352
2353     ! in Na3H(SO4)2
2354     jE = jna3hso4
2355     b_mtem(1,jA,jE) = 3.06830
2356     b_mtem(2,jA,jE) = -3.18408
2357     b_mtem(3,jA,jE) = -19.6332
2358     b_mtem(4,jA,jE) = 61.3657
2359     b_mtem(5,jA,jE) = -73.4438
2360     b_mtem(6,jA,jE) = 31.2334
2361
2362     ! in Na2SO4
2363     jE = jna2so4
2364     b_mtem(1,jA,jE) = 2.58649
2365     b_mtem(2,jA,jE) = 0.87921
2366     b_mtem(3,jA,jE) = -39.3023
2367     b_mtem(4,jA,jE) = 101.603
2368     b_mtem(5,jA,jE) = -109.469
2369     b_mtem(6,jA,jE) = 43.0188
2370
2371     ! in HNO3
2372     jE = jhno3
2373     b_mtem(1,jA,jE) = 1.54587
2374     b_mtem(2,jA,jE) = -7.50976
2375     b_mtem(3,jA,jE) = 12.8237
2376     b_mtem(4,jA,jE) = -10.1452
2377     b_mtem(5,jA,jE) = -0.541956
2378     b_mtem(6,jA,jE) = 3.34536
2379
2380     ! in HCl
2381     jE = jhcl
2382     b_mtem(1,jA,jE) = 0.829757
2383     b_mtem(2,jA,jE) = -4.11316
2384     b_mtem(3,jA,jE) = 3.67111
2385     b_mtem(4,jA,jE) = 3.6833
2386     b_mtem(5,jA,jE) = -11.2711
2387     b_mtem(6,jA,jE) = 6.71421
2388
2389
2390     !----------
2391     ! H.HSO4 in E
2392     jA = jhhso4
2393
2394     ! in H2SO4
2395     jE = jh2so4
2396     b_mtem(1,jA,jE) = 2.63953
2397     b_mtem(2,jA,jE) = -6.01532
2398     b_mtem(3,jA,jE) = 10.0204
2399     b_mtem(4,jA,jE) = -12.4840
2400     b_mtem(5,jA,jE) = 7.78853
2401     b_mtem(6,jA,jE) = -2.12638
2402
2403     ! in NH4HSO4
2404     jE = jnh4hso4
2405     b_mtem(1,jA,jE) = -0.77412
2406     b_mtem(2,jA,jE) = 14.1656
2407     b_mtem(3,jA,jE) = -53.4087
2408     b_mtem(4,jA,jE) = 93.2013
2409     b_mtem(5,jA,jE) = -80.5723
2410     b_mtem(6,jA,jE) = 27.1577
2411
2412     ! in (NH4)3H(SO4)2
2413     jE = jlvcite
2414     b_mtem(1,jA,jE) = -2.98882
2415     b_mtem(2,jA,jE) = 14.4436
2416     b_mtem(3,jA,jE) = -40.1774
2417     b_mtem(4,jA,jE) = 67.5937
2418     b_mtem(5,jA,jE) = -61.5040
2419     b_mtem(6,jA,jE) = 22.3695
2420
2421     ! in (NH4)2SO4
2422     jE = jnh4so4
2423     b_mtem(1,jA,jE) = -1.15502
2424     b_mtem(2,jA,jE) = 8.12309
2425     b_mtem(3,jA,jE) = -38.4726
2426     b_mtem(4,jA,jE) = 80.8861
2427     b_mtem(5,jA,jE) = -80.1644
2428     b_mtem(6,jA,jE) = 30.4717
2429
2430     ! in NaHSO4
2431     jE = jnahso4
2432     b_mtem(1,jA,jE) = 1.99641
2433     b_mtem(2,jA,jE) = -2.96061
2434     b_mtem(3,jA,jE) = 5.54778
2435     b_mtem(4,jA,jE) = -14.5488
2436     b_mtem(5,jA,jE) = 14.8492
2437     b_mtem(6,jA,jE) = -5.1389
2438
2439     ! in Na3H(SO4)2
2440     jE = jna3hso4
2441     b_mtem(1,jA,jE) = 2.23816
2442     b_mtem(2,jA,jE) = -3.20847
2443     b_mtem(3,jA,jE) = -4.82853
2444     b_mtem(4,jA,jE) = 20.9192
2445     b_mtem(5,jA,jE) = -27.2819
2446     b_mtem(6,jA,jE) = 11.8655
2447
2448     ! in Na2SO4
2449     jE = jna2so4
2450     b_mtem(1,jA,jE) = 2.56907
2451     b_mtem(2,jA,jE) = 1.13444
2452     b_mtem(3,jA,jE) = -34.6853
2453     b_mtem(4,jA,jE) = 87.9775
2454     b_mtem(5,jA,jE) = -93.2330
2455     b_mtem(6,jA,jE) = 35.9260
2456
2457     ! in HNO3
2458     jE = jhno3
2459     b_mtem(1,jA,jE) = 2.00024
2460     b_mtem(2,jA,jE) = -4.80868
2461     b_mtem(3,jA,jE) = 8.29222
2462     b_mtem(4,jA,jE) = -11.0849
2463     b_mtem(5,jA,jE) = 7.51262
2464     b_mtem(6,jA,jE) = -2.07654
2465
2466     ! in HCl
2467     jE = jhcl
2468     b_mtem(1,jA,jE) = 2.8009
2469     b_mtem(2,jA,jE) = -6.98416
2470     b_mtem(3,jA,jE) = 14.3146
2471     b_mtem(4,jA,jE) = -22.0068
2472     b_mtem(5,jA,jE) = 17.5557
2473     b_mtem(6,jA,jE) = -5.84917
2474
2475
2476     !----------
2477     ! NH4HSO4 in E
2478     jA = jnh4hso4
2479
2480     ! in H2SO4
2481     jE = jh2so4
2482     b_mtem(1,jA,jE) = 0.169160
2483     b_mtem(2,jA,jE) = 2.15094
2484     b_mtem(3,jA,jE) = -9.62904
2485     b_mtem(4,jA,jE) = 18.2631
2486     b_mtem(5,jA,jE) = -17.3333
2487     b_mtem(6,jA,jE) = 6.19835
2488
2489     ! in NH4HSO4
2490     jE = jnh4hso4
2491     b_mtem(1,jA,jE) = -2.34457
2492     b_mtem(2,jA,jE) = 12.8035
2493     b_mtem(3,jA,jE) = -35.2513
2494     b_mtem(4,jA,jE) = 53.6153
2495     b_mtem(5,jA,jE) = -42.7655
2496     b_mtem(6,jA,jE) = 13.7129
2497
2498     ! in (NH4)3H(SO4)2
2499     jE = jlvcite
2500     b_mtem(1,jA,jE) = -2.56109
2501     b_mtem(2,jA,jE) = 11.1414
2502     b_mtem(3,jA,jE) = -30.2361
2503     b_mtem(4,jA,jE) = 50.0320
2504     b_mtem(5,jA,jE) = -44.1586
2505     b_mtem(6,jA,jE) = 15.5393
2506
2507     ! in (NH4)2SO4
2508     jE = jnh4so4
2509     b_mtem(1,jA,jE) = -0.97315
2510     b_mtem(2,jA,jE) = 7.06295
2511     b_mtem(3,jA,jE) = -29.3032
2512     b_mtem(4,jA,jE) = 57.6101
2513     b_mtem(5,jA,jE) = -54.9020
2514     b_mtem(6,jA,jE) = 20.2222
2515
2516     ! in NaHSO4
2517     jE = jnahso4
2518     b_mtem(1,jA,jE) = -0.44450
2519     b_mtem(2,jA,jE) = 3.33451
2520     b_mtem(3,jA,jE) = -15.2791
2521     b_mtem(4,jA,jE) = 30.1413
2522     b_mtem(5,jA,jE) = -26.7710
2523     b_mtem(6,jA,jE) = 8.78462
2524
2525     ! in Na3H(SO4)2
2526     jE = jna3hso4
2527     b_mtem(1,jA,jE) = -0.99780
2528     b_mtem(2,jA,jE) = 4.69200
2529     b_mtem(3,jA,jE) = -16.1219
2530     b_mtem(4,jA,jE) = 29.3100
2531     b_mtem(5,jA,jE) = -26.3383
2532     b_mtem(6,jA,jE) = 9.20695
2533
2534     ! in Na2SO4
2535     jE = jna2so4
2536     b_mtem(1,jA,jE) = -0.52694
2537     b_mtem(2,jA,jE) = 7.02684
2538     b_mtem(3,jA,jE) = -33.7508
2539     b_mtem(4,jA,jE) = 70.0565
2540     b_mtem(5,jA,jE) = -68.3226
2541     b_mtem(6,jA,jE) = 25.2692
2542
2543     ! in HNO3
2544     jE = jhno3
2545     b_mtem(1,jA,jE) = 0.572926
2546     b_mtem(2,jA,jE) = -2.04791
2547     b_mtem(3,jA,jE) = 2.1134
2548     b_mtem(4,jA,jE) = 0.246654
2549     b_mtem(5,jA,jE) = -3.06019
2550     b_mtem(6,jA,jE) = 1.98126
2551
2552     ! in HCl
2553     jE = jhcl
2554     b_mtem(1,jA,jE) = 0.56514
2555     b_mtem(2,jA,jE) = 0.22287
2556     b_mtem(3,jA,jE) = -2.76973
2557     b_mtem(4,jA,jE) = 4.54444
2558     b_mtem(5,jA,jE) = -3.86549
2559     b_mtem(6,jA,jE) = 1.13441
2560
2561
2562     !----------
2563     ! (NH4)3H(SO4)2 in E
2564     jA = jlvcite
2565
2566     ! in H2SO4
2567     jE = jh2so4
2568     b_mtem(1,jA,jE) = -1.44811
2569     b_mtem(2,jA,jE) = 6.71815
2570     b_mtem(3,jA,jE) = -25.0141
2571     b_mtem(4,jA,jE) = 50.1109
2572     b_mtem(5,jA,jE) = -50.0561
2573     b_mtem(6,jA,jE) = 19.3370
2574
2575     ! in NH4HSO4
2576     jE = jnh4hso4
2577     b_mtem(1,jA,jE) = -3.41707
2578     b_mtem(2,jA,jE) = 13.4496
2579     b_mtem(3,jA,jE) = -34.8018
2580     b_mtem(4,jA,jE) = 55.2987
2581     b_mtem(5,jA,jE) = -48.1839
2582     b_mtem(6,jA,jE) = 17.2444
2583
2584     ! in (NH4)3H(SO4)2
2585     jE = jlvcite
2586     b_mtem(1,jA,jE) = -2.54479
2587     b_mtem(2,jA,jE) = 11.8501
2588     b_mtem(3,jA,jE) = -39.7286
2589     b_mtem(4,jA,jE) = 74.2479
2590     b_mtem(5,jA,jE) = -70.4934
2591     b_mtem(6,jA,jE) = 26.2836
2592
2593     ! in (NH4)2SO4
2594     jE = jnh4so4
2595     b_mtem(1,jA,jE) = -2.30561
2596     b_mtem(2,jA,jE) = 14.5806
2597     b_mtem(3,jA,jE) = -55.1238
2598     b_mtem(4,jA,jE) = 103.451
2599     b_mtem(5,jA,jE) = -95.2571
2600     b_mtem(6,jA,jE) = 34.2218
2601
2602     ! in NaHSO4
2603     jE = jnahso4
2604     b_mtem(1,jA,jE) = -2.20809
2605     b_mtem(2,jA,jE) = 13.6391
2606     b_mtem(3,jA,jE) = -57.8246
2607     b_mtem(4,jA,jE) = 117.907
2608     b_mtem(5,jA,jE) = -112.154
2609     b_mtem(6,jA,jE) = 40.3058
2610
2611     ! in Na3H(SO4)2
2612     jE = jna3hso4
2613     b_mtem(1,jA,jE) = -1.15099
2614     b_mtem(2,jA,jE) = 6.32269
2615     b_mtem(3,jA,jE) = -27.3860
2616     b_mtem(4,jA,jE) = 55.4592
2617     b_mtem(5,jA,jE) = -54.0100
2618     b_mtem(6,jA,jE) = 20.3469
2619
2620     ! in Na2SO4
2621     jE = jna2so4
2622     b_mtem(1,jA,jE) = -1.15678
2623     b_mtem(2,jA,jE) = 8.28718
2624     b_mtem(3,jA,jE) = -37.3231
2625     b_mtem(4,jA,jE) = 76.6124
2626     b_mtem(5,jA,jE) = -74.9307
2627     b_mtem(6,jA,jE) = 28.0559
2628
2629     ! in HNO3
2630     jE = jhno3
2631     b_mtem(1,jA,jE) = 0.01502
2632     b_mtem(2,jA,jE) = -3.1197
2633     b_mtem(3,jA,jE) = 3.61104
2634     b_mtem(4,jA,jE) = 3.05196
2635     b_mtem(5,jA,jE) = -9.98957
2636     b_mtem(6,jA,jE) = 6.04155
2637
2638     ! in HCl
2639     jE = jhcl
2640     b_mtem(1,jA,jE) = -1.06477
2641     b_mtem(2,jA,jE) = 3.38801
2642     b_mtem(3,jA,jE) = -12.5784
2643     b_mtem(4,jA,jE) = 25.2823
2644     b_mtem(5,jA,jE) = -25.4611
2645     b_mtem(6,jA,jE) = 10.0754
2646
2647
2648     !----------
2649     ! NaHSO4 in E
2650     jA = jnahso4
2651
2652     ! in H2SO4
2653     jE = jh2so4
2654     b_mtem(1,jA,jE) = 0.68259
2655     b_mtem(2,jA,jE) = 0.71468
2656     b_mtem(3,jA,jE) = -5.59003
2657     b_mtem(4,jA,jE) = 11.0089
2658     b_mtem(5,jA,jE) = -10.7983
2659     b_mtem(6,jA,jE) = 3.82335
2660
2661     ! in NH4HSO4
2662     jE = jnh4hso4
2663     b_mtem(1,jA,jE) = -0.03956
2664     b_mtem(2,jA,jE) = 4.52828
2665     b_mtem(3,jA,jE) = -25.2557
2666     b_mtem(4,jA,jE) = 54.4225
2667     b_mtem(5,jA,jE) = -52.5105
2668     b_mtem(6,jA,jE) = 18.6562
2669
2670     ! in (NH4)3H(SO4)2
2671     jE = jlvcite
2672     b_mtem(1,jA,jE) = -1.53503
2673     b_mtem(2,jA,jE) = 8.27608
2674     b_mtem(3,jA,jE) = -28.9539
2675     b_mtem(4,jA,jE) = 55.2876
2676     b_mtem(5,jA,jE) = -51.9563
2677     b_mtem(6,jA,jE) = 18.6576
2678
2679     ! in (NH4)2SO4
2680     jE = jnh4so4
2681     b_mtem(1,jA,jE) = -0.38793
2682     b_mtem(2,jA,jE) = 7.14680
2683     b_mtem(3,jA,jE) = -38.7201
2684     b_mtem(4,jA,jE) = 84.3965
2685     b_mtem(5,jA,jE) = -84.7453
2686     b_mtem(6,jA,jE) = 32.1283
2687
2688     ! in NaHSO4
2689     jE = jnahso4
2690     b_mtem(1,jA,jE) = -0.41982
2691     b_mtem(2,jA,jE) = 4.26491
2692     b_mtem(3,jA,jE) = -20.2351
2693     b_mtem(4,jA,jE) = 42.6764
2694     b_mtem(5,jA,jE) = -40.7503
2695     b_mtem(6,jA,jE) = 14.2868
2696
2697     ! in Na3H(SO4)2
2698     jE = jna3hso4
2699     b_mtem(1,jA,jE) = -0.32912
2700     b_mtem(2,jA,jE) = 1.80808
2701     b_mtem(3,jA,jE) = -8.01286
2702     b_mtem(4,jA,jE) = 15.5791
2703     b_mtem(5,jA,jE) = -14.5494
2704     b_mtem(6,jA,jE) = 5.27052
2705
2706     ! in Na2SO4
2707     jE = jna2so4
2708     b_mtem(1,jA,jE) = 0.10271
2709     b_mtem(2,jA,jE) = 5.09559
2710     b_mtem(3,jA,jE) = -30.3295
2711     b_mtem(4,jA,jE) = 66.2975
2712     b_mtem(5,jA,jE) = -66.3458
2713     b_mtem(6,jA,jE) = 24.9443
2714
2715     ! in HNO3
2716     jE = jhno3
2717     b_mtem(1,jA,jE) = 0.608309
2718     b_mtem(2,jA,jE) = -0.541905
2719     b_mtem(3,jA,jE) = -2.52084
2720     b_mtem(4,jA,jE) = 6.63297
2721     b_mtem(5,jA,jE) = -7.24599
2722     b_mtem(6,jA,jE) = 2.88811
2723
2724     ! in HCl
2725     jE = jhcl
2726     b_mtem(1,jA,jE) = 1.98399
2727     b_mtem(2,jA,jE) = -4.51562
2728     b_mtem(3,jA,jE) = 8.36059
2729     b_mtem(4,jA,jE) = -12.4948
2730     b_mtem(5,jA,jE) = 9.67514
2731     b_mtem(6,jA,jE) = -3.18004
2732
2733
2734     !----------
2735     ! Na3H(SO4)2 in E
2736     jA = jna3hso4
2737
2738     ! in H2SO4
2739     jE = jh2so4
2740     b_mtem(1,jA,jE) = -0.83214
2741     b_mtem(2,jA,jE) = 4.99572
2742     b_mtem(3,jA,jE) = -20.1697
2743     b_mtem(4,jA,jE) = 41.4066
2744     b_mtem(5,jA,jE) = -42.2119
2745     b_mtem(6,jA,jE) = 16.4855
2746
2747     ! in NH4HSO4
2748     jE = jnh4hso4
2749     b_mtem(1,jA,jE) = -0.65139
2750     b_mtem(2,jA,jE) = 3.52300
2751     b_mtem(3,jA,jE) = -22.8220
2752     b_mtem(4,jA,jE) = 56.2956
2753     b_mtem(5,jA,jE) = -59.9028
2754     b_mtem(6,jA,jE) = 23.1844
2755
2756     ! in (NH4)3H(SO4)2
2757     jE = jlvcite
2758     b_mtem(1,jA,jE) = -1.31331
2759     b_mtem(2,jA,jE) = 8.40835
2760     b_mtem(3,jA,jE) = -38.1757
2761     b_mtem(4,jA,jE) = 80.5312
2762     b_mtem(5,jA,jE) = -79.8346
2763     b_mtem(6,jA,jE) = 30.0219
2764
2765     ! in (NH4)2SO4
2766     jE = jnh4so4
2767     b_mtem(1,jA,jE) = -1.03054
2768     b_mtem(2,jA,jE) = 8.08155
2769     b_mtem(3,jA,jE) = -38.1046
2770     b_mtem(4,jA,jE) = 78.7168
2771     b_mtem(5,jA,jE) = -77.2263
2772     b_mtem(6,jA,jE) = 29.1521
2773
2774     ! in NaHSO4
2775     jE = jnahso4
2776     b_mtem(1,jA,jE) = -1.90695
2777     b_mtem(2,jA,jE) = 11.6241
2778     b_mtem(3,jA,jE) = -50.3175
2779     b_mtem(4,jA,jE) = 105.884
2780     b_mtem(5,jA,jE) = -103.258
2781     b_mtem(6,jA,jE) = 37.6588
2782
2783     ! in Na3H(SO4)2
2784     jE = jna3hso4
2785     b_mtem(1,jA,jE) = -0.34780
2786     b_mtem(2,jA,jE) = 2.85363
2787     b_mtem(3,jA,jE) = -17.6224
2788     b_mtem(4,jA,jE) = 38.9220
2789     b_mtem(5,jA,jE) = -39.8106
2790     b_mtem(6,jA,jE) = 15.6055
2791
2792     ! in Na2SO4
2793     jE = jna2so4
2794     b_mtem(1,jA,jE) = -0.75230
2795     b_mtem(2,jA,jE) = 10.0140
2796     b_mtem(3,jA,jE) = -50.5677
2797     b_mtem(4,jA,jE) = 106.941
2798     b_mtem(5,jA,jE) = -105.534
2799     b_mtem(6,jA,jE) = 39.5196
2800
2801     ! in HNO3
2802     jE = jhno3
2803     b_mtem(1,jA,jE) = 0.057456
2804     b_mtem(2,jA,jE) = -1.31264
2805     b_mtem(3,jA,jE) = -1.94662
2806     b_mtem(4,jA,jE) = 10.7024
2807     b_mtem(5,jA,jE) = -14.9946
2808     b_mtem(6,jA,jE) = 7.12161
2809
2810     ! in HCl
2811     jE = jhcl
2812     b_mtem(1,jA,jE) = 0.637894
2813     b_mtem(2,jA,jE) = -2.29719
2814     b_mtem(3,jA,jE) = 0.765361
2815     b_mtem(4,jA,jE) = 4.8748
2816     b_mtem(5,jA,jE) = -9.25978
2817     b_mtem(6,jA,jE) = 4.91773
2818     !
2819     !
2820     !
2821     !----------------------------------------------------------
2822     ! Coefficients for %MDRH(T) = d1 + d2*T + d3*T^2 + d4*T^3    (T in Kelvin)
2823     ! valid Temperature Range: 240 - 320 K
2824     !----------------------------------------------------------
2825     !
2826     ! SULFATE-POOR SYSTEMS
2827     ! AC
2828     j_index = 1
2829     d_mdrh(j_index,1) = -58.00268351
2830     d_mdrh(j_index,2) = 2.031077573
2831     d_mdrh(j_index,3) = -0.008281218
2832     d_mdrh(j_index,4) = 1.00447E-05
2833
2834     ! AN
2835     j_index = 2
2836     d_mdrh(j_index,1) = 1039.137773
2837     d_mdrh(j_index,2) = -11.47847095
2838     d_mdrh(j_index,3) = 0.047702786
2839     d_mdrh(j_index,4) = -6.77675E-05
2840
2841     ! AS
2842     j_index = 3
2843     d_mdrh(j_index,1) = 115.8366357
2844     d_mdrh(j_index,2) = 0.491881663
2845     d_mdrh(j_index,3) = -0.00422807
2846     d_mdrh(j_index,4) = 7.29274E-06
2847
2848     ! SC
2849     j_index = 4
2850     d_mdrh(j_index,1) = 253.2424151
2851     d_mdrh(j_index,2) = -1.429957864
2852     d_mdrh(j_index,3) = 0.003727554
2853     d_mdrh(j_index,4) = -3.13037E-06
2854
2855     ! SN
2856     j_index = 5
2857     d_mdrh(j_index,1) = -372.4306506
2858     d_mdrh(j_index,2) = 5.3955633
2859     d_mdrh(j_index,3) = -0.019804438
2860     d_mdrh(j_index,4) = 2.25662E-05
2861
2862     ! SS
2863     j_index = 6
2864     d_mdrh(j_index,1) = 286.1271416
2865     d_mdrh(j_index,2) = -1.670787758
2866     d_mdrh(j_index,3) = 0.004431373
2867     d_mdrh(j_index,4) = -3.57757E-06
2868
2869     ! CC
2870     j_index = 7
2871     d_mdrh(j_index,1) = -1124.07059
2872     d_mdrh(j_index,2) = 14.26364209
2873     d_mdrh(j_index,3) = -0.054816822
2874     d_mdrh(j_index,4) = 6.70107E-05
2875
2876     ! CN
2877     j_index = 8
2878     d_mdrh(j_index,1) = 1855.413934
2879     d_mdrh(j_index,2) = -20.29219473
2880     d_mdrh(j_index,3) = 0.07807482
2881     d_mdrh(j_index,4) = -1.017887858e-4
2882
2883     ! AN + AC
2884     j_index = 9
2885     d_mdrh(j_index,1) = 1761.176886
2886     d_mdrh(j_index,2) = -19.29811062
2887     d_mdrh(j_index,3) = 0.075676987
2888     d_mdrh(j_index,4) = -1.0116959e-4
2889
2890     ! AS + AC
2891     j_index = 10
2892     d_mdrh(j_index,1) = 122.1074303
2893     d_mdrh(j_index,2) = 0.429692122
2894     d_mdrh(j_index,3) = -0.003928277
2895     d_mdrh(j_index,4) = 6.43275E-06
2896
2897     ! AS + AN
2898     j_index = 11
2899     d_mdrh(j_index,1) = 2424.634678
2900     d_mdrh(j_index,2) = -26.54031307
2901     d_mdrh(j_index,3) = 0.101625387
2902     d_mdrh(j_index,4) = -1.31544547798e-4
2903
2904     ! AS + AN + AC
2905     j_index = 12
2906     d_mdrh(j_index,1) = 2912.082599
2907     d_mdrh(j_index,2) = -31.8894185
2908     d_mdrh(j_index,3) = 0.121185849
2909     d_mdrh(j_index,4) = -1.556534623e-4
2910
2911     ! SC + AC
2912     j_index = 13
2913     d_mdrh(j_index,1) = 172.2596493
2914     d_mdrh(j_index,2) = -0.511006195
2915     d_mdrh(j_index,3) = 4.27244597e-4
2916     d_mdrh(j_index,4) = 4.12797E-07
2917
2918     ! SN + AC
2919     j_index = 14
2920     d_mdrh(j_index,1) = 1596.184935
2921     d_mdrh(j_index,2) = -16.37945565
2922     d_mdrh(j_index,3) = 0.060281218
2923     d_mdrh(j_index,4) = -7.6161E-05
2924
2925     ! SN + AN
2926     j_index = 15
2927     d_mdrh(j_index,1) = 1916.072988
2928     d_mdrh(j_index,2) = -20.85594868
2929     d_mdrh(j_index,3) = 0.081140141
2930     d_mdrh(j_index,4) = -1.07954274796e-4
2931
2932     ! SN + AN + AC
2933     j_index = 16
2934     d_mdrh(j_index,1) = 1467.165935
2935     d_mdrh(j_index,2) = -16.01166196
2936     d_mdrh(j_index,3) = 0.063505582
2937     d_mdrh(j_index,4) = -8.66722E-05
2938
2939     ! SN + SC
2940     j_index = 17
2941     d_mdrh(j_index,1) = 158.447059
2942     d_mdrh(j_index,2) = -0.628167358
2943     d_mdrh(j_index,3) = 0.002014448
2944     d_mdrh(j_index,4) = -3.13037E-06
2945
2946     ! SN + SC + AC
2947     j_index = 18
2948     d_mdrh(j_index,1) = 1115.892468
2949     d_mdrh(j_index,2) = -11.76936534
2950     d_mdrh(j_index,3) = 0.045577399
2951     d_mdrh(j_index,4) = -6.05779E-05
2952
2953     ! SS + AC
2954     j_index = 19
2955     d_mdrh(j_index,1) = 269.5432407
2956     d_mdrh(j_index,2) = -1.319963885
2957     d_mdrh(j_index,3) = 0.002592363
2958     d_mdrh(j_index,4) = -1.44479E-06
2959
2960     ! SS + AN
2961     j_index = 20
2962     d_mdrh(j_index,1) = 2841.334784
2963     d_mdrh(j_index,2) = -31.1889487
2964     d_mdrh(j_index,3) = 0.118809274
2965     d_mdrh(j_index,4) = -1.53007e-4
2966
2967     ! SS + AN + AC
2968     j_index = 21
2969     d_mdrh(j_index,1) = 2199.36914
2970     d_mdrh(j_index,2) = -24.11926569
2971     d_mdrh(j_index,3) = 0.092932361
2972     d_mdrh(j_index,4) = -1.21774e-4
2973
2974     ! SS + AS
2975     j_index = 22
2976     d_mdrh(j_index,1) = 395.0051604
2977     d_mdrh(j_index,2) = -2.521101657
2978     d_mdrh(j_index,3) = 0.006139319
2979     d_mdrh(j_index,4) = -4.43756E-06
2980
2981     ! SS + AS + AC
2982     j_index = 23
2983     d_mdrh(j_index,1) = 386.5150675
2984     d_mdrh(j_index,2) = -2.4632138
2985     d_mdrh(j_index,3) = 0.006139319
2986     d_mdrh(j_index,4) = -4.98796E-06
2987
2988     ! SS + AS + AN
2989     j_index = 24
2990     d_mdrh(j_index,1) = 3101.538491
2991     d_mdrh(j_index,2) = -34.19978105
2992     d_mdrh(j_index,3) = 0.130118605
2993     d_mdrh(j_index,4) = -1.66873e-4
2994
2995     ! SS + AS + AN + AC
2996     j_index = 25
2997     d_mdrh(j_index,1) = 2307.579403
2998     d_mdrh(j_index,2) = -25.43136774
2999     d_mdrh(j_index,3) = 0.098064728
3000     d_mdrh(j_index,4) = -1.28301e-4
3001
3002     ! SS + SC
3003     j_index = 26
3004     d_mdrh(j_index,1) = 291.8309602
3005     d_mdrh(j_index,2) = -1.828912974
3006     d_mdrh(j_index,3) = 0.005053148
3007     d_mdrh(j_index,4) = -4.57516E-06
3008
3009     ! SS + SC + AC
3010     j_index = 27
3011     d_mdrh(j_index,1) = 188.3914345
3012     d_mdrh(j_index,2) = -0.631345031
3013     d_mdrh(j_index,3) = 0.000622807
3014     d_mdrh(j_index,4) = 4.47196E-07
3015
3016     ! SS + SN
3017     j_index = 28
3018     d_mdrh(j_index,1) = -167.1252839
3019     d_mdrh(j_index,2) = 2.969828002
3020     d_mdrh(j_index,3) = -0.010637255
3021     d_mdrh(j_index,4) = 1.13175E-05
3022
3023     ! SS + SN + AC
3024     j_index = 29
3025     d_mdrh(j_index,1) = 1516.782768
3026     d_mdrh(j_index,2) = -15.7922661
3027     d_mdrh(j_index,3) = 0.058942209
3028     d_mdrh(j_index,4) = -7.5301E-05
3029
3030     ! SS + SN + AN
3031     j_index = 30
3032     d_mdrh(j_index,1) = 1739.963163
3033     d_mdrh(j_index,2) = -19.06576022
3034     d_mdrh(j_index,3) = 0.07454963
3035     d_mdrh(j_index,4) = -9.94302E-05
3036
3037     ! SS + SN + AN + AC
3038     j_index = 31
3039     d_mdrh(j_index,1) = 2152.104877
3040     d_mdrh(j_index,2) = -23.74998008
3041     d_mdrh(j_index,3) = 0.092256654
3042     d_mdrh(j_index,4) = -1.21953e-4
3043
3044     ! SS + SN + SC
3045     j_index = 32
3046     d_mdrh(j_index,1) = 221.9976265
3047     d_mdrh(j_index,2) = -1.311331272
3048     d_mdrh(j_index,3) = 0.004406089
3049     d_mdrh(j_index,4) = -5.88235E-06
3050
3051     ! SS + SN + SC + AC
3052     j_index = 33
3053     d_mdrh(j_index,1) = 1205.645615
3054     d_mdrh(j_index,2) = -12.71353459
3055     d_mdrh(j_index,3) = 0.048803922
3056     d_mdrh(j_index,4) = -6.41899E-05
3057
3058     ! CC + AC
3059     j_index = 34
3060     d_mdrh(j_index,1) = 506.6737879
3061     d_mdrh(j_index,2) = -3.723520818
3062     d_mdrh(j_index,3) = 0.010814242
3063     d_mdrh(j_index,4) = -1.21087E-05
3064
3065     ! CC + SC
3066     j_index = 35
3067     d_mdrh(j_index,1) = -1123.523841
3068     d_mdrh(j_index,2) = 14.08345977
3069     d_mdrh(j_index,3) = -0.053687823
3070     d_mdrh(j_index,4) = 6.52219E-05
3071
3072     ! CC + SC + AC
3073     j_index = 36
3074     d_mdrh(j_index,1) = -1159.98607
3075     d_mdrh(j_index,2) = 14.44309169
3076     d_mdrh(j_index,3) = -0.054841073
3077     d_mdrh(j_index,4) = 6.64259E-05
3078
3079     ! CN + AC
3080     j_index = 37
3081     d_mdrh(j_index,1) = 756.0747916
3082     d_mdrh(j_index,2) = -8.546826257
3083     d_mdrh(j_index,3) = 0.035798677
3084     d_mdrh(j_index,4) = -5.06629E-05
3085
3086     ! CN + AN
3087     j_index = 38
3088     d_mdrh(j_index,1) = 338.668191
3089     d_mdrh(j_index,2) = -2.971223403
3090     d_mdrh(j_index,3) = 0.012294866
3091     d_mdrh(j_index,4) = -1.87558E-05
3092
3093     ! CN + AN + AC
3094     j_index = 39
3095     d_mdrh(j_index,1) = -53.18033508
3096     d_mdrh(j_index,2) = 0.663911748
3097     d_mdrh(j_index,3) = 9.16326e-4
3098     d_mdrh(j_index,4) = -6.70354E-06
3099
3100     ! CN + SC
3101     j_index = 40
3102     d_mdrh(j_index,1) = 3623.831129
3103     d_mdrh(j_index,2) = -39.27226457
3104     d_mdrh(j_index,3) = 0.144559515
3105     d_mdrh(j_index,4) = -1.78159e-4
3106
3107     ! CN + SC + AC
3108     j_index = 41
3109     d_mdrh(j_index,1) = 3436.656743
3110     d_mdrh(j_index,2) = -37.16192684
3111     d_mdrh(j_index,3) = 0.136641377
3112     d_mdrh(j_index,4) = -1.68262e-4
3113
3114     ! CN + SN
3115     j_index = 42
3116     d_mdrh(j_index,1) = 768.608476
3117     d_mdrh(j_index,2) = -8.051517149
3118     d_mdrh(j_index,3) = 0.032342332
3119     d_mdrh(j_index,4) = -4.52224E-05
3120
3121     ! CN + SN + AC
3122     j_index = 43
3123     d_mdrh(j_index,1) = 33.58027951
3124     d_mdrh(j_index,2) = -0.308772182
3125     d_mdrh(j_index,3) = 0.004713639
3126     d_mdrh(j_index,4) = -1.19658E-05
3127
3128     ! CN + SN + AN
3129     j_index = 44
3130     d_mdrh(j_index,1) = 57.80183041
3131     d_mdrh(j_index,2) = 0.215264604
3132     d_mdrh(j_index,3) = 4.11406e-4
3133     d_mdrh(j_index,4) = -4.30702E-06
3134
3135     ! CN + SN + AN + AC
3136     j_index = 45
3137     d_mdrh(j_index,1) = -234.368984
3138     d_mdrh(j_index,2) = 2.721045204
3139     d_mdrh(j_index,3) = -0.006688341
3140     d_mdrh(j_index,4) = 2.31729E-06
3141
3142     ! CN + SN + SC
3143     j_index = 46
3144     d_mdrh(j_index,1) = 3879.080557
3145     d_mdrh(j_index,2) = -42.13562874
3146     d_mdrh(j_index,3) = 0.155235005
3147     d_mdrh(j_index,4) = -1.91387e-4
3148
3149     ! CN + SN + SC + AC
3150     j_index = 47
3151     d_mdrh(j_index,1) = 3600.576985
3152     d_mdrh(j_index,2) = -39.0283489
3153     d_mdrh(j_index,3) = 0.143710316
3154     d_mdrh(j_index,4) = -1.77167e-4
3155
3156     ! CN + CC
3157     j_index = 48
3158     d_mdrh(j_index,1) = -1009.729826
3159     d_mdrh(j_index,2) = 12.9145339
3160     d_mdrh(j_index,3) = -0.049811146
3161     d_mdrh(j_index,4) = 6.09563E-05
3162
3163     ! CN + CC + AC
3164     j_index = 49
3165     d_mdrh(j_index,1) = -577.0919514
3166     d_mdrh(j_index,2) = 8.020324227
3167     d_mdrh(j_index,3) = -0.031469556
3168     d_mdrh(j_index,4) = 3.82181E-05
3169
3170     ! CN + CC + SC
3171     j_index = 50
3172     d_mdrh(j_index,1) = -728.9983499
3173     d_mdrh(j_index,2) = 9.849458215
3174     d_mdrh(j_index,3) = -0.03879257
3175     d_mdrh(j_index,4) = 4.78844E-05
3176
3177     ! CN + CC + SC + AC
3178     j_index = 51
3179     d_mdrh(j_index,1) = -803.7026845
3180     d_mdrh(j_index,2) = 10.61881494
3181     d_mdrh(j_index,3) = -0.041402993
3182     d_mdrh(j_index,4) = 5.08084E-05
3183
3184     !
3185     ! SULFATE-RICH SYSTEMS
3186     ! AB
3187     j_index = 52
3188     d_mdrh(j_index,1) = -493.6190458
3189     d_mdrh(j_index,2) = 6.747053851
3190     d_mdrh(j_index,3) = -0.026955267
3191     d_mdrh(j_index,4) = 3.45118E-05
3192
3193     ! LV
3194     j_index = 53
3195     d_mdrh(j_index,1) = 53.37874093
3196     d_mdrh(j_index,2) = 1.01368249
3197     d_mdrh(j_index,3) = -0.005887513
3198     d_mdrh(j_index,4) = 8.94393E-06
3199
3200     ! SB
3201     j_index = 54
3202     d_mdrh(j_index,1) = 206.619047
3203     d_mdrh(j_index,2) = -1.342735684
3204     d_mdrh(j_index,3) = 0.003197691
3205     d_mdrh(j_index,4) = -1.93603E-06
3206
3207     ! AB + LV
3208     j_index = 55
3209     d_mdrh(j_index,1) = -493.6190458
3210     d_mdrh(j_index,2) = 6.747053851
3211     d_mdrh(j_index,3) = -0.026955267
3212     d_mdrh(j_index,4) = 3.45118E-05
3213
3214     ! AS + LV
3215     j_index = 56
3216     d_mdrh(j_index,1) = 53.37874093
3217     d_mdrh(j_index,2) = 1.01368249
3218     d_mdrh(j_index,3) = -0.005887513
3219     d_mdrh(j_index,4) = 8.94393E-06
3220
3221     ! SS + SB
3222     j_index = 57
3223     d_mdrh(j_index,1) = 206.619047
3224     d_mdrh(j_index,2) = -1.342735684
3225     d_mdrh(j_index,3) = 0.003197691
3226     d_mdrh(j_index,4) = -1.93603E-06
3227
3228     ! SS + LV
3229     j_index = 58
3230     d_mdrh(j_index,1) = 41.7619047
3231     d_mdrh(j_index,2) = 1.303872053
3232     d_mdrh(j_index,3) = -0.007647908
3233     d_mdrh(j_index,4) = 1.17845E-05
3234
3235     ! SS + AS + LV
3236     j_index = 59
3237     d_mdrh(j_index,1) = 41.7619047
3238     d_mdrh(j_index,2) = 1.303872053
3239     d_mdrh(j_index,3) = -0.007647908
3240     d_mdrh(j_index,4) = 1.17845E-05
3241
3242     ! SS + AB
3243     j_index = 60
3244     d_mdrh(j_index,1) = -369.7142842
3245     d_mdrh(j_index,2) = 5.512878771
3246     d_mdrh(j_index,3) = -0.02301948
3247     d_mdrh(j_index,4) = 3.0303E-05
3248
3249     ! SS + LV + AB
3250     j_index = 61
3251     d_mdrh(j_index,1) = -369.7142842
3252     d_mdrh(j_index,2) = 5.512878771
3253     d_mdrh(j_index,3) = -0.02301948
3254     d_mdrh(j_index,4) = 3.0303E-05
3255
3256     ! SB + AB
3257     j_index = 62
3258     d_mdrh(j_index,1) = -162.8095232
3259     d_mdrh(j_index,2) = 2.399326592
3260     d_mdrh(j_index,3) = -0.009336219
3261     d_mdrh(j_index,4) = 1.17845E-05
3262
3263     ! SS + SB + AB
3264     j_index = 63
3265     d_mdrh(j_index,1) = -735.4285689
3266     d_mdrh(j_index,2) = 8.885521857
3267     d_mdrh(j_index,3) = -0.033488456
3268     d_mdrh(j_index,4) = 4.12458E-05
3269
3270     !   endif ! first
3271
3272
3273     if (msoa_flag1 >= 1000) call soa_vbs_load_params( 2 )
3274
3275
3276     return
3277   end subroutine load_mosaic_parameters
3278
3279
3280
3281 end module module_mosaic_init_aerpar