module module_mosaic_init_aerpar use module_data_mosaic_kind, only: r8 implicit none private public:: mosaic_init_aer_params contains !------------------------------------------------------------------------------- subroutine mosaic_init_aer_params !BSINGH - All initialzations for Mosiac model call load_mosaic_parameters end subroutine mosaic_init_aer_params !------------------------------------------------------------------------------- subroutine load_mosaic_parameters !---------------------------------------------------------------------------------------! !BSINGH: load_mosaic_parameters subroutine is directly copied form the mosaic_box.25.f90 ! code !---------------------------------------------------------------------------------------! !---------------------------------------------------------------------------------------! ! Called only once per entire simulation to load gas and aerosol ! indices, parameters, physico-chemical constants, polynomial coeffs, etc. ! ! author: Rahul A. Zaveri ! update: jan 2005 !---------------------------------------------------------------------------------------! use module_mosaic_soa_vbs, only: soa_vbs_load_params use module_data_mosaic_aero, only: ipmcmos_aero, no_aerosol, all_solid, all_liquid, & mixed, nelectrolyte, naercomp, naer, Ncation, Nanion, & ngas_aerchtot, ngas_volatile, nsalt, & jsulf_poor_NUM, jsulf_rich_NUM, MDRH_T_NUM, d_mdrh_DIM2, phasestate, aer_name, & gas_name, ename, jnh4so4, jlvcite, jnh4hso4, jnh4msa, jnh4no3, jnh4cl, jna2so4, & jna3hso4, jnahso4, jnamsa, jnano3, jnacl, jcano3, jcacl2, jcamsa2, jh2so4, jmsa, & jhno3, jhcl, jhhso4, jcaso4, jcaco3, joc, jbc, join, jaro1, jaro2, jalk1, jole1, & japi1, japi2, jlim1, jlim2, jh2o, jc_h, jc_nh4, jc_na, jc_ca, ja_hso4, ja_so4, & ja_no3, ja_cl, ja_msa, ih2so4_g, ihno3_g, ihcl_g, inh3_g, imsa_g, iaro1_g, & iaro2_g, ialk1_g, iole1_g, iapi1_g, iapi2_g, ilim1_g, ilim2_g, iso4_a, ino3_a, & icl_a, inh4_a, imsa_a, iaro1_a, iaro2_a, ialk1_a, iole1_a, iapi1_a, iapi2_a, & ilim1_a, ilim2_a, ico3_a, ina_a, ica_a, ioin_a, ioc_a, ibc_a, & isoa_first, jsoa_first, nmax_ASTEM, b_mtem, & zc, za, b_zsr, a_zsr, aw_min, mw_electrolyte, dens_electrolyte, & partial_molar_vol, MW_c, MW_a, mw_aer_mac,dens_aer_mac, kappa_aer_mac, & dens_comp_a, mw_comp_a, ref_index_a, mw_gas, v_molar_gas, & rtol_mesa, jsalt_index, jsulf_poor, jsulf_rich, Nmax_mesa, d_mdrh, msoa_flag1, & use_cam5mam_soa_params ! local variables integer :: ia, iaer, ig, igas, ibin, ja, je, j_index logical :: use_mos31e_rz1_densities logical :: use_uniform_densities !BSINGH - 05/28/2013(RCE updates) logical :: use_sorgam_soa_species real(r8), dimension(nelectrolyte) :: G_MX,K_MX !BSINGH - 05/28/2013(RCE updates) use_mos31e_rz1_densities = .true. if ( use_mos31e_rz1_densities ) then use_uniform_densities = .false. else use_uniform_densities = .true. if (ipmcmos_aero > 0) use_uniform_densities = .false. end if !BSINGH - 05/28/2013(RCE updates ENDS) if (msoa_flag1 == 1) then use_sorgam_soa_species = .true. else use_sorgam_soa_species = .false. end if ! rce 2013-07-31 - ! using a local saved variable like "first" no longer works ! the calling routine needs to determine if/when this routine is needed ! if(first)then ! first=.false. !---------------------------------------------------------------- ! control settings ! *** do not change mSIZE_FRAMEWORK here *** ! mSIZE_FRAMEWORK = mSECTIONAL ! mMODAL or mSECTIONAL ! mDYNAMIC_SOLVER = mASTEM ! mASTEM, mLSODES ! mGAS_AER_XFER = mON ! mON, mOFF ! ASTEM parameters nmax_ASTEM = 301 ! max number of time steps in ASTEM ! alpha_ASTEM = 1.0 ! choose a value between 0.01 and 1.0 ! rtol_eqb_ASTEM = 0.01 ! equilibrium tolerance in ASTEM ! ptol_mol_ASTEM = 0.01 ! mol percent tolerance in ASTEM ! MESA parameters Nmax_MESA = 80 ! max number of iterations in MESA_PTC rtol_mesa = 0.01 ! MESA equilibrium tolerance !---------------------------------------------------------------- ! ! set species indices ! ixxx_a are for aerosol species - they apply to aer and total arrays ! ixxx_g are for gas species - they apply to gas array ! jyyy are for electrolyte species - they apply to electrolyte and comp_a arrays ! ! *** CRITICAL "RULES" ABOUT THESE INDICES *** ! ! *** THE CODE WILL FAIL IF THEY ARE NOT FOLLOWED *** ! ! for the volatile inorganic and organic species (1 thru ngas_volatile) ! the ixxx_a and ixxx_g must match (be indentical) for each aerosol-gas pair ! ! for the soa species, the ordering of the ixxx_a, ixxx_g, and jyyy ! must be identical ! ! electrolyte indices (used for water content calculations) ! these indices are order sensitive ! inorganic species first jnh4so4 = 1 ! soluble jlvcite = 2 ! soluble jnh4hso4 = 3 ! soluble jnh4msa = 4 ! soluble: new jnh4no3 = 5 ! soluble jnh4cl = 6 ! soluble jna2so4 = 7 ! soluble jna3hso4 = 8 ! soluble jnahso4 = 9 ! soluble jnamsa = 10 ! soluble: new jnano3 = 11 ! soluble jnacl = 12 ! soluble jcano3 = 13 ! soluble jcacl2 = 14 ! soluble jcamsa2 = 15 ! soluble nsalt jh2so4 = 16 ! soluble jmsa = 17 ! soluble jhno3 = 18 ! soluble jhcl = 19 ! soluble jhhso4 = 20 ! soluble jcaso4 = 21 ! insoluble jcaco3 = 22 ! insoluble joc = 23 ! insoluble - part of naercomp jbc = 24 ! insoluble - part of naercomp join = 25 ! insoluble - part of naercomp ! aerosol and gas indices for inorganic species iso4_a = 1 ; ih2so4_g = 1 ino3_a = 2 ; ihno3_g = 2 icl_a = 3 ; ihcl_g = 3 inh4_a = 4 ; inh3_g = 4 imsa_a = 5 ; imsa_g = 5 ia = imsa_a ; ig = imsa_g ; je = join ! aerosol, gas, and electrolyte indices for sorgam secondary organic species if ( use_sorgam_soa_species ) then isoa_first = ia+1 ; jsoa_first = je+1 iaro1_a = ia+1 ; iaro1_g = ig+1 ; jaro1 = je+1 iaro2_a = ia+2 ; iaro2_g = ig+2 ; jaro2 = je+2 ialk1_a = ia+3 ; ialk1_g = ig+3 ; jalk1 = je+3 iole1_a = ia+4 ; iole1_g = ig+4 ; jole1 = je+4 iapi1_a = ia+5 ; iapi1_g = ig+5 ; japi1 = je+5 iapi2_a = ia+6 ; iapi2_g = ig+6 ; japi2 = je+6 ilim1_a = ia+7 ; ilim1_g = ig+7 ; jlim1 = je+7 ilim2_a = ia+8 ; ilim2_g = ig+8 ; jlim2 = je+8 ia = ia+8 ; ig = ig+8 ; je = je+8 else isoa_first = ia+1 ; jsoa_first = je+1 iaro1_a = -999888777 ; iaro1_g = -999888777 ; jaro1 = -999888777 iaro2_a = -999888777 ; iaro2_g = -999888777 ; jaro2 = -999888777 ialk1_a = -999888777 ; ialk1_g = -999888777 ; jalk1 = -999888777 iole1_a = -999888777 ; iole1_g = -999888777 ; jole1 = -999888777 iapi1_a = -999888777 ; iapi1_g = -999888777 ; japi1 = -999888777 iapi2_a = -999888777 ; iapi2_g = -999888777 ; japi2 = -999888777 ilim1_a = -999888777 ; ilim1_g = -999888777 ; jlim1 = -999888777 ilim2_a = -999888777 ; ilim2_g = -999888777 ; jlim2 = -999888777 ia = ia+1 ; ig = ig+1 ; je = je+1 end if ! electrolyte indices for other species jh2o = je+1 ! water - part of naercomp ! aerosol indices for other species ico3_a = ia+1 ina_a = ia+2 ica_a = ia+3 ioin_a = ia+4 ioc_a = ia+5 ibc_a = ia+6 ! gas indices for other species ! ico2_g = 14 ! *** currently not used ! *** if co3_a was treated as a volatile inorganic, ! then ico3_a and ico2_g would have to be 6 ! local aerosol ions ! cations jc_h = 1 jc_nh4 = 2 jc_na = 3 jc_ca = 4 ! ! anions ja_hso4 = 1 ja_so4 = 2 ja_no3 = 3 ja_cl = 4 ja_msa = 5 ! ja_co3 = 6 ! *** currently not used if (msoa_flag1 >= 1000) call soa_vbs_load_params( 1 ) !-------------------------------------------------------------------- ! phase state names phasestate(no_aerosol) = "NOAERO" phasestate(all_solid) = "SOLID " phasestate(all_liquid) = "LIQUID" phasestate(mixed) = "MIXED " ! names of aer species do iaer = 1, naer write( aer_name(iaer), '(a,i4.4)' ) 'aer', iaer ! default end do aer_name(iso4_a) = "SO4" aer_name(ino3_a) = "NO3" aer_name(icl_a) = "Cl " aer_name(inh4_a) = "NH4" aer_name(ioc_a) = "OC " aer_name(imsa_a) = "MSA" aer_name(ico3_a) = "CO3" aer_name(ina_a) = "Na " aer_name(ica_a) = "Ca " aer_name(ibc_a) = "BC " aer_name(ioin_a) = "OIN" if ( use_sorgam_soa_species ) then aer_name(iaro1_a)= "ARO1" aer_name(iaro2_a)= "ARO2" aer_name(ialk1_a)= "ALK1" aer_name(iole1_a)= "OLE1" aer_name(iapi1_a)= "API1" aer_name(iapi2_a)= "API2" aer_name(ilim1_a)= "LIM1" aer_name(ilim2_a)= "LIM2" end if ! names of gas species do igas = 1, ngas_aerchtot write( gas_name(igas), '(a,i4.4)' ) 'gas', igas ! default end do gas_name(ih2so4_g) = "H2SO4" gas_name(ihno3_g) = "HNO3 " gas_name(ihcl_g) = "HCl " gas_name(inh3_g) = "NH3 " gas_name(imsa_g) = "MSA " if ( use_sorgam_soa_species ) then gas_name(iaro1_g) = "ARO1 " gas_name(iaro2_g) = "ARO2 " gas_name(ialk1_g) = "ALK1 " gas_name(iole1_g) = "OLE1 " gas_name(iapi1_g) = "API1 " gas_name(iapi2_g) = "API2 " gas_name(ilim1_g) = "LIM1 " gas_name(ilim2_g) = "LIM2 " end if ! names of electrolytes ename(jnh4so4) = "AmSO4" ename(jlvcite) = "(NH4)3H(SO4)2" ename(jnh4hso4)= "NH4HSO4" ename(jnh4msa) = "CH3SO3NH4" ename(jnh4no3) = "NH4NO3" ename(jnh4cl) = "NH4Cl" ename(jnacl) = "NaCl" ename(jnano3) = "NaNO3" ename(jna2so4) = "Na2SO4" ename(jna3hso4)= "Na3H(SO4)2" ename(jnamsa) = "CH3SO3Na" ename(jnahso4) = "NaHSO4" ename(jcaso4) = "CaSO4" ename(jcamsa2) = "(CH3SO3)2Ca" ename(jcano3) = "Ca(NO3)2" ename(jcacl2) = "CaCl2" ename(jcaco3) = "CaCO3" ename(jh2so4) = "H2SO4" ename(jhhso4) = "HHSO4" ename(jhno3) = "HNO3" ename(jhcl) = "HCl" ename(jmsa) = "CH3SO3H" ! molecular weights of electrolytes mw_electrolyte(jnh4so4) = 132.0 mw_electrolyte(jlvcite) = 247.0 mw_electrolyte(jnh4hso4)= 115.0 mw_electrolyte(jnh4msa) = 113.0 mw_electrolyte(jnh4no3) = 80.0 mw_electrolyte(jnh4cl) = 53.5 mw_electrolyte(jnacl) = 58.5 mw_electrolyte(jnano3) = 85.0 mw_electrolyte(jna2so4) = 142.0 mw_electrolyte(jna3hso4)= 262.0 mw_electrolyte(jnahso4) = 120.0 mw_electrolyte(jnamsa) = 118.0 mw_electrolyte(jcaso4) = 136.0 mw_electrolyte(jcamsa2) = 230.0 mw_electrolyte(jcano3) = 164.0 mw_electrolyte(jcacl2) = 111.0 mw_electrolyte(jcaco3) = 100.0 mw_electrolyte(jh2so4) = 98.0 mw_electrolyte(jhno3) = 63.0 mw_electrolyte(jhcl) = 36.5 mw_electrolyte(jmsa) = 96.0 ! molecular weights of ions [g/mol] MW_c(jc_h) = 1.0 MW_c(jc_nh4)= 18.0 MW_c(jc_na) = 23.0 MW_c(jc_ca) = 40.0 MW_a(ja_so4) = 96.0 MW_a(ja_hso4)= 97.0 MW_a(ja_no3) = 62.0 MW_a(ja_cl) = 35.5 MW_a(ja_msa) = 95.0 ! magnitude of the charges on ions zc(jc_h) = 1 zc(jc_nh4) = 1 zc(jc_na) = 1 zc(jc_ca) = 2 za(ja_hso4)= 1 za(ja_so4) = 2 za(ja_no3) = 1 za(ja_cl) = 1 za(ja_msa) = 1 ! densities of pure electrolytes in g/cc dens_electrolyte(jnh4so4) = 1.8 dens_electrolyte(jlvcite) = 1.8 dens_electrolyte(jnh4hso4) = 1.8 dens_electrolyte(jnh4msa) = 1.8 ! assumed same as nh4hso4 dens_electrolyte(jnh4no3) = 1.8 dens_electrolyte(jnh4cl) = 1.8 dens_electrolyte(jnacl) = 2.2 dens_electrolyte(jnano3) = 2.2 dens_electrolyte(jna2so4) = 2.2 dens_electrolyte(jna3hso4) = 2.2 dens_electrolyte(jnahso4) = 2.2 dens_electrolyte(jnamsa) = 2.2 ! assumed same as nahso4 dens_electrolyte(jcaso4) = 2.6 dens_electrolyte(jcamsa2) = 2.6 ! assumed same as caso4 dens_electrolyte(jcano3) = 2.6 dens_electrolyte(jcacl2) = 2.6 dens_electrolyte(jcaco3) = 2.6 dens_electrolyte(jh2so4) = 1.8 dens_electrolyte(jhhso4) = 1.8 dens_electrolyte(jhno3) = 1.8 dens_electrolyte(jhcl) = 1.8 dens_electrolyte(jmsa) = 1.8 ! assumed same as h2so4 if ( use_uniform_densities ) then!BSINGH - 05/28/2013(RCE updates) do je = 1, nelectrolyte dens_electrolyte(je) = 1.6 enddo endif!BSINGH - 05/28/2013(RCE updates) ! densities of compounds in g/cc dens_comp_a(jnh4so4) = 1.8 dens_comp_a(jlvcite) = 1.8 dens_comp_a(jnh4hso4) = 1.8 dens_comp_a(jnh4msa) = 1.8 ! assumed same as nh4hso4 dens_comp_a(jnh4no3) = 1.7 dens_comp_a(jnh4cl) = 1.5 dens_comp_a(jnacl) = 2.2 dens_comp_a(jnano3) = 2.2 dens_comp_a(jna2so4) = 2.2 dens_comp_a(jna3hso4) = 2.2 dens_comp_a(jnahso4) = 2.2 dens_comp_a(jnamsa) = 2.2 ! assumed same as nahso4 dens_comp_a(jcaso4) = 2.6 dens_comp_a(jcamsa2) = 2.6 ! assumed same as caso4 dens_comp_a(jcano3) = 2.6 dens_comp_a(jcacl2) = 2.6 dens_comp_a(jcaco3) = 2.6 dens_comp_a(jh2so4) = 1.8 dens_comp_a(jhhso4) = 1.8 dens_comp_a(jhno3) = 1.8 dens_comp_a(jhcl) = 1.8 dens_comp_a(jmsa) = 1.8 ! assumed same as h2so4 dens_comp_a(joc) = 1.0 dens_comp_a(jbc) = 1.8 dens_comp_a(join) = 2.6 if ( use_sorgam_soa_species ) then dens_comp_a(jaro1) = 1.0 dens_comp_a(jaro2) = 1.0 dens_comp_a(jalk1) = 1.0 dens_comp_a(jole1) = 1.0 dens_comp_a(japi1) = 1.0 dens_comp_a(japi2) = 1.0 dens_comp_a(jlim1) = 1.0 dens_comp_a(jlim2) = 1.0 end if dens_comp_a(jh2o) = 1.0 !BSINGH - 05/28/2013(RCE updates) ! following for comparison with mos31d_bs2 and. mos31e_rz1 if ( use_mos31e_rz1_densities ) then dens_comp_a(joc) = 1.4 if ( use_sorgam_soa_species ) then dens_comp_a(jaro1) = 1.4 dens_comp_a(jaro2) = 1.4 dens_comp_a(jalk1) = 1.4 dens_comp_a(jole1) = 1.4 dens_comp_a(japi1) = 1.4 dens_comp_a(japi2) = 1.4 dens_comp_a(jlim1) = 1.4 dens_comp_a(jlim2) = 1.4 end if end if if ( use_uniform_densities ) then !BSINGH - 05/28/2013(RCE updates ENDS) do je = 1, naercomp dens_comp_a(je) = 1.6 enddo !BSINGH - 05/28/2013(RCE updates) endif if (ipmcmos_aero > 0) then dens_comp_a(jnh4no3) = 1.8 dens_comp_a(jnh4cl) = 1.8 if ( use_sorgam_soa_species ) then dens_comp_a(jaro1) = 1.4 dens_comp_a(jaro2) = 1.4 dens_comp_a(jalk1) = 1.4 dens_comp_a(jole1) = 1.4 dens_comp_a(japi1) = 1.4 dens_comp_a(japi2) = 1.4 dens_comp_a(jlim1) = 1.4 dens_comp_a(jlim2) = 1.4 end if endif !BSINGH - 05/28/2013(RCE updates ENDS) ! molecular weights of generic aerosol species mw_aer_mac(1:naer) = 200.0 ! default mw_aer_mac(iso4_a) = 96.0 mw_aer_mac(ino3_a) = 62.0 mw_aer_mac(icl_a) = 35.5 mw_aer_mac(imsa_a) = 95.0 ! CH3SO3 mw_aer_mac(ico3_a) = 60.0 mw_aer_mac(inh4_a) = 18.0 mw_aer_mac(ina_a) = 23.0 mw_aer_mac(ica_a) = 40.0 mw_aer_mac(ioin_a) = 1.0 ! not used mw_aer_mac(ibc_a) = 1.0 ! not used mw_aer_mac(ioc_a) = 1.0 ! 200 assumed for primary organics if ( use_sorgam_soa_species ) then mw_aer_mac(iaro1_a)= 150.0 mw_aer_mac(iaro2_a)= 150.0 mw_aer_mac(ialk1_a)= 140.0 mw_aer_mac(iole1_a)= 140.0 mw_aer_mac(iapi1_a)= 184.0 mw_aer_mac(iapi2_a)= 184.0 mw_aer_mac(ilim1_a)= 200.0 mw_aer_mac(ilim2_a)= 200.0 end if ! molecular weights of compounds mw_comp_a(jnh4so4) = 132.0 mw_comp_a(jlvcite) = 247.0 mw_comp_a(jnh4hso4)= 115.0 mw_comp_a(jnh4msa) = 113.0 mw_comp_a(jnh4no3) = 80.0 mw_comp_a(jnh4cl) = 53.5 mw_comp_a(jnacl) = 58.5 mw_comp_a(jnano3) = 85.0 mw_comp_a(jna2so4) = 142.0 mw_comp_a(jna3hso4)= 262.0 mw_comp_a(jnahso4) = 120.0 mw_comp_a(jnamsa) = 118.0 mw_comp_a(jcaso4) = 136.0 mw_comp_a(jcamsa2) = 230.0 mw_comp_a(jcano3) = 164.0 mw_comp_a(jcacl2) = 111.0 mw_comp_a(jcaco3) = 100.0 mw_comp_a(jh2so4) = 98.0 mw_comp_a(jhhso4) = 98.0 mw_comp_a(jhno3) = 63.0 mw_comp_a(jhcl) = 36.5 mw_comp_a(jmsa) = 96.0 mw_comp_a(joc) = 1.0 mw_comp_a(jbc) = 1.0 mw_comp_a(join) = 1.0 if ( use_sorgam_soa_species ) then mw_comp_a(jaro1) = 150.0 mw_comp_a(jaro2) = 150.0 mw_comp_a(jalk1) = 140.0 mw_comp_a(jole1) = 140.0 mw_comp_a(japi1) = 184.0 mw_comp_a(japi2) = 184.0 mw_comp_a(jlim1) = 200.0 mw_comp_a(jlim2) = 200.0 end if mw_comp_a(jh2o) = 18.0 !BSINGH - 05/28/2013(RCE updates) ! partmc-2.2.1 jun-2012 !# dens (kg/m^3) ions in soln (1) molec wght (kg/mole) kappa (1) ! SO4 1800 0 96d-3 0.65 ! NO3 1800 0 62d-3 0.65 ! Cl 2200 0 35.5d-3 0.53 ! NH4 1800 0 18d-3 0.65 ! MSA 1800 0 95d-3 0.53 ! ARO1 1400 0 150d-3 0.1 ! ARO2 1400 0 150d-3 0.1 ! ALK1 1400 0 140d-3 0.1 ! OLE1 1400 0 140d-3 0.1 ! API1 1400 0 184d-3 0.1 ! API2 1400 0 184d-3 0.1 ! LIM1 1400 0 200d-3 0.1 ! LIM2 1400 0 200d-3 0.1 ! CO3 2600 0 60d-3 0.53 ! Na 2200 0 23d-3 0.53 ! Ca 2600 0 40d-3 0.53 ! OIN 2600 0 1d-3 0.1 ! OC 1000 0 1d-3 0.001 ! BC 1800 0 1d-3 0 ! H2O 1000 0 18d-3 0 !BSINGH - 05/28/2013(RCE updates ENDS) ! densities of generic aerosol species dens_aer_mac(1:naer) = 1.0 ! default dens_aer_mac(iso4_a) = 1.8 ! used dens_aer_mac(ino3_a) = 1.8 ! used dens_aer_mac(icl_a) = 2.2 ! used dens_aer_mac(imsa_a) = 1.8 ! used dens_aer_mac(ico3_a) = 2.6 ! used dens_aer_mac(inh4_a) = 1.8 ! used dens_aer_mac(ina_a) = 2.2 ! used dens_aer_mac(ica_a) = 2.6 ! used dens_aer_mac(ioin_a) = 2.6 ! used dens_aer_mac(ioc_a) = 1.0 ! used dens_aer_mac(ibc_a) = 1.8 ! used if ( use_sorgam_soa_species ) then dens_aer_mac(iaro1_a)= 1.0 dens_aer_mac(iaro2_a)= 1.0 dens_aer_mac(ialk1_a)= 1.0 dens_aer_mac(iole1_a)= 1.0 dens_aer_mac(iapi1_a)= 1.0 dens_aer_mac(iapi2_a)= 1.0 dens_aer_mac(ilim1_a)= 1.0 dens_aer_mac(ilim2_a)= 1.0 end if !BSINGH - 05/28/2013(RCE updates) ! following for comparison with mos31d_bs2 and. mos31e_rz1 if ( use_mos31e_rz1_densities ) then dens_aer_mac(ioc_a) = 1.4 if ( use_sorgam_soa_species ) then dens_aer_mac(iaro1_a)= 1.4 dens_aer_mac(iaro2_a)= 1.4 dens_aer_mac(ialk1_a)= 1.4 dens_aer_mac(iole1_a)= 1.4 dens_aer_mac(iapi1_a)= 1.4 dens_aer_mac(iapi2_a)= 1.4 dens_aer_mac(ilim1_a)= 1.4 dens_aer_mac(ilim2_a)= 1.4 end if end if if ( use_uniform_densities ) then !BSINGH - 05/28/2013(RCE updates ENDS) do iaer = 1, naer dens_aer_mac(iaer) = 1.6 enddo endif!BSINGH - 05/28/2013(RCE updates) if (ipmcmos_aero > 0) then ! use partmc-mosaic densities dens_aer_mac(1:19) = (/ & 1.80, 1.80, 2.20, 1.80, 1.80, 1.40, 1.40, 1.40, 1.40, 1.40, & 1.40, 1.40, 1.40, 2.60, 2.20, 2.60, 2.60, 1.00, 1.80 /)!BSINGH - 05/28/2013(RCE updates) ! so4 no3 cl nh4 msa aro1 aro2 alk1 ole1 api1 ! api2 lim1 lim2 co3 na ca oin oc bc end if if ( use_cam5mam_soa_params > 0 ) then dens_aer_mac(ioc_a) = 1.0 dens_comp_a(joc) = 1.0 ! for oc, leave mw=1 because some of the mosaic code requires this if ( 1 <= ilim2_a .and. ilim2_a <= naer ) then dens_aer_mac(ilim2_a) = 1.0 dens_comp_a(jlim2) = 1.0 mw_aer_mac(ilim2_a) = 150.0 mw_comp_a(jlim2) = 150.0 end if end if ! kappa values (hygroscopicities) of generic aerosol species ! ! for calculation of ccn properties, kappa of electrolytes ! should be used ! the multi-dimensional sectional code needs a "fixed" kappa ! for each generic aerosol species, just as the older ! 1d sectional code needs a "fixed" dry density kappa_aer_mac(1:naer) = 0.1 ! default kappa_aer_mac(iso4_a) = 0.65 kappa_aer_mac(ino3_a) = 0.65 kappa_aer_mac(imsa_a) = 0.65 kappa_aer_mac(inh4_a) = 0.65 kappa_aer_mac(icl_a) = 0.65 kappa_aer_mac(ina_a) = 0.65 kappa_aer_mac(ico3_a) = 0.001 ! ?? kappa_aer_mac(ica_a) = 0.001 ! ?? kappa_aer_mac(ioin_a) = 0.001 kappa_aer_mac(ioc_a) = 0.001 kappa_aer_mac(ibc_a) = 0.001 if ( use_sorgam_soa_species ) then kappa_aer_mac(iaro1_a) = 0.1 kappa_aer_mac(iaro2_a) = 0.1 kappa_aer_mac(ialk1_a) = 0.1 kappa_aer_mac(iole1_a) = 0.1 kappa_aer_mac(iapi1_a) = 0.1 kappa_aer_mac(iapi2_a) = 0.1 kappa_aer_mac(ilim1_a) = 0.1 kappa_aer_mac(ilim2_a) = 0.1 end if !BSINGH - 05/28/2013(RCE updates) if (ipmcmos_aero > 0) then ! use partmc-mosaic kappas kappa_aer_mac(1:19) = (/ & 0.65, 0.65, 0.53, 0.65, 0.53, 0.10, 0.10, 0.10, 0.10, 0.10, & 0.10, 0.10, 0.10, 0.53, 0.53, 0.53, 0.10, 0.001, 0.0 /) ! so4 no3 cl nh4 msa aro1 aro2 alk1 ole1 api1 ! api2 lim1 lim2 co3 na ca oin oc bc end if !BSINGH - 05/28/2013(RCE updates ENDS) ! partial molar volumes of condensing gases partial_molar_vol(1:ngas_aerchtot) = 200.0 ! default partial_molar_vol(ih2so4_g) = 51.83 partial_molar_vol(ihno3_g) = 31.45 partial_molar_vol(ihcl_g) = 20.96 partial_molar_vol(inh3_g) = 24.03 partial_molar_vol(imsa_g) = 53.33 if ( use_sorgam_soa_species ) then partial_molar_vol(iaro1_g) = 150.0 partial_molar_vol(iaro2_g) = 150.0 partial_molar_vol(ialk1_g) = 140.0 partial_molar_vol(iole1_g) = 140.0 partial_molar_vol(iapi1_g) = 184.0 partial_molar_vol(iapi2_g) = 184.0 partial_molar_vol(ilim1_g) = 200.0 partial_molar_vol(ilim2_g) = 200.0 end if ! molecular weights of condensing gases mw_gas(1:ngas_aerchtot) = 200.0 ! default mw_gas(ih2so4_g) = 98.0 mw_gas(ihno3_g) = 63.0 mw_gas(ihcl_g) = 36.5 mw_gas(inh3_g) = 17.0 mw_gas(imsa_g) = 96.0 if ( use_sorgam_soa_species ) then mw_gas(iaro1_g) = 150.0 mw_gas(iaro2_g) = 150.0 mw_gas(ialk1_g) = 140.0 mw_gas(iole1_g) = 140.0 mw_gas(iapi1_g) = 184.0 mw_gas(iapi2_g) = 184.0 mw_gas(ilim1_g) = 200.0 mw_gas(ilim2_g) = 200.0 end if ! used to calculate diffusivities of condensing gases v_molar_gas(1:ngas_aerchtot) = 60.0 ! default v_molar_gas(ih2so4_g)= 42.88 v_molar_gas(ihno3_g) = 24.11 v_molar_gas(ihcl_g) = 21.48 v_molar_gas(inh3_g) = 14.90 v_molar_gas(imsa_g) = 58.00 ! refractive index ref_index_a(jnh4so4) = cmplx(1.52,0.) ref_index_a(jlvcite) = cmplx(1.50,0.) ref_index_a(jnh4hso4)= cmplx(1.47,0.) ref_index_a(jnh4msa) = cmplx(1.50,0.) ! assumed ref_index_a(jnh4no3) = cmplx(1.50,0.) ref_index_a(jnh4cl) = cmplx(1.50,0.) ref_index_a(jnacl) = cmplx(1.45,0.) ref_index_a(jnano3) = cmplx(1.50,0.) ref_index_a(jna2so4) = cmplx(1.50,0.) ref_index_a(jna3hso4)= cmplx(1.50,0.) ref_index_a(jnahso4) = cmplx(1.50,0.) ref_index_a(jnamsa) = cmplx(1.50,0.) ! assumed ref_index_a(jcaso4) = cmplx(1.56,0.006) ref_index_a(jcamsa2) = cmplx(1.56,0.006) ! assumed ref_index_a(jcano3) = cmplx(1.56,0.006) ref_index_a(jcacl2) = cmplx(1.52,0.006) ref_index_a(jcaco3) = cmplx(1.68,0.006) ref_index_a(jh2so4) = cmplx(1.43,0.) ref_index_a(jhhso4) = cmplx(1.43,0.) ref_index_a(jhno3) = cmplx(1.50,0.) ref_index_a(jhcl) = cmplx(1.50,0.) ref_index_a(jmsa) = cmplx(1.43,0.) ! assumed ref_index_a(joc) = cmplx(1.45,0.) ref_index_a(jbc) = cmplx(1.82,0.74) ref_index_a(join) = cmplx(1.55,0.006) if ( use_sorgam_soa_species ) then ref_index_a(jaro1) = cmplx(1.45,0.) ref_index_a(jaro2) = cmplx(1.45,0.) ref_index_a(jalk1) = cmplx(1.45,0.) ref_index_a(jole1) = cmplx(1.45,0.) ref_index_a(japi1) = cmplx(1.45,0.) ref_index_a(japi2) = cmplx(1.45,0.) ref_index_a(jlim1) = cmplx(1.45,0.) ref_index_a(jlim2) = cmplx(1.45,0.) end if ref_index_a(jh2o) = cmplx(1.33,0.) ! jsalt_index jsalt_index(jnh4so4) = 5 ! AS jsalt_index(jlvcite) = 2 ! LV jsalt_index(jnh4hso4)= 1 ! AB jsalt_index(jnh4no3) = 2 ! AN jsalt_index(jnh4cl) = 1 ! AC jsalt_index(jna2so4) = 60 ! SS jsalt_index(jnahso4) = 10 ! SB jsalt_index(jnano3) = 40 ! SN jsalt_index(jnacl) = 10 ! SC jsalt_index(jcano3) = 120 ! CN jsalt_index(jcacl2) = 80 ! CC jsalt_index(jnh4msa) = 0 ! AM zero for now jsalt_index(jnamsa) = 0 ! SM zero for now jsalt_index(jcamsa2) = 0 ! CM zero for now ! Aerosol Indices ! AC = 1, AN = 2, AS = 5, SC = 10, SN = 40, SS = 60, CC = 80, CN = 120, ! AB = 1, LV = 2, SB = 10 ! ! SULFATE-POOR DOMAIN jsulf_poor(1) = 1 ! AC jsulf_poor(2) = 2 ! AN jsulf_poor(5) = 3 ! AS jsulf_poor(10) = 4 ! SC jsulf_poor(40) = 5 ! SN jsulf_poor(60) = 6 ! SS jsulf_poor(80) = 7 ! CC jsulf_poor(120) = 8 ! CN jsulf_poor(3) = 9 ! AN + AC jsulf_poor(6) = 10 ! AS + AC jsulf_poor(7) = 11 ! AS + AN jsulf_poor(8) = 12 ! AS + AN + AC jsulf_poor(11) = 13 ! SC + AC jsulf_poor(41) = 14 ! SN + AC jsulf_poor(42) = 15 ! SN + AN jsulf_poor(43) = 16 ! SN + AN + AC jsulf_poor(50) = 17 ! SN + SC jsulf_poor(51) = 18 ! SN + SC + AC jsulf_poor(61) = 19 ! SS + AC jsulf_poor(62) = 20 ! SS + AN jsulf_poor(63) = 21 ! SS + AN + AC jsulf_poor(65) = 22 ! SS + AS jsulf_poor(66) = 23 ! SS + AS + AC jsulf_poor(67) = 24 ! SS + AS + AN jsulf_poor(68) = 25 ! SS + AS + AN + AC jsulf_poor(70) = 26 ! SS + SC jsulf_poor(71) = 27 ! SS + SC + AC jsulf_poor(100) = 28 ! SS + SN jsulf_poor(101) = 29 ! SS + SN + AC jsulf_poor(102) = 30 ! SS + SN + AN jsulf_poor(103) = 31 ! SS + SN + AN + AC jsulf_poor(110) = 32 ! SS + SN + SC jsulf_poor(111) = 33 ! SS + SN + SC + AC jsulf_poor(81) = 34 ! CC + AC jsulf_poor(90) = 35 ! CC + SC jsulf_poor(91) = 36 ! CC + SC + AC jsulf_poor(121) = 37 ! CN + AC jsulf_poor(122) = 38 ! CN + AN jsulf_poor(123) = 39 ! CN + AN + AC jsulf_poor(130) = 40 ! CN + SC jsulf_poor(131) = 41 ! CN + SC + AC jsulf_poor(160) = 42 ! CN + SN jsulf_poor(161) = 43 ! CN + SN + AC jsulf_poor(162) = 44 ! CN + SN + AN jsulf_poor(163) = 45 ! CN + SN + AN + AC jsulf_poor(170) = 46 ! CN + SN + SC jsulf_poor(171) = 47 ! CN + SN + SC + AC jsulf_poor(200) = 48 ! CN + CC jsulf_poor(201) = 49 ! CN + CC + AC jsulf_poor(210) = 50 ! CN + CC + SC jsulf_poor(211) = 51 ! CN + CC + SC + AC ! ! SULFATE-RICH DOMAIN jsulf_rich(1) = 52 ! AB jsulf_rich(2) = 53 ! LV jsulf_rich(10) = 54 ! SB jsulf_rich(3) = 55 ! AB + LV jsulf_rich(7) = 56 ! AS + LV jsulf_rich(70) = 57 ! SS + SB jsulf_rich(62) = 58 ! SS + LV jsulf_rich(67) = 59 ! SS + AS + LV jsulf_rich(61) = 60 ! SS + AB jsulf_rich(63) = 61 ! SS + LV + AB jsulf_rich(11) = 62 ! SB + AB jsulf_rich(71) = 63 ! SS + SB + AB jsulf_rich(5) = 3 ! AS jsulf_rich(60) = 6 ! SS jsulf_rich(65) = 22 ! SS + AS ! ! polynomial coefficients for binary molality (used in ZSR equation) ! ! ! a_zsr for aw < 0.97 ! ! (NH4)2SO4 je = jnh4so4 a_zsr(1,je) = 1.30894 a_zsr(2,je) = -7.09922 a_zsr(3,je) = 20.62831 a_zsr(4,je) = -32.19965 a_zsr(5,je) = 25.17026 a_zsr(6,je) = -7.81632 aw_min(je) = 0.1 ! ! (NH4)3H(SO4)2 je = jlvcite a_zsr(1,je) = 1.10725 a_zsr(2,je) = -5.17978 a_zsr(3,je) = 12.29534 a_zsr(4,je) = -16.32545 a_zsr(5,je) = 11.29274 a_zsr(6,je) = -3.19164 aw_min(je) = 0.1 ! ! NH4HSO4 je = jnh4hso4 a_zsr(1,je) = 1.15510 a_zsr(2,je) = -3.20815 a_zsr(3,je) = 2.71141 a_zsr(4,je) = 2.01155 a_zsr(5,je) = -4.71014 a_zsr(6,je) = 2.04616 aw_min(je) = 0.1 ! ! NH4MSA (assumed same as NH4HSO4) je = jnh4msa a_zsr(1,je) = 1.15510 a_zsr(2,je) = -3.20815 a_zsr(3,je) = 2.71141 a_zsr(4,je) = 2.01155 a_zsr(5,je) = -4.71014 a_zsr(6,je) = 2.04616 aw_min(je) = 0.1 ! ! NH4NO3 je = jnh4no3 a_zsr(1,je) = 0.43507 a_zsr(2,je) = 6.38220 a_zsr(3,je) = -30.19797 a_zsr(4,je) = 53.36470 a_zsr(5,je) = -43.44203 a_zsr(6,je) = 13.46158 aw_min(je) = 0.1 ! ! NH4Cl: revised on Nov 13, 2003. based on Chan and Ha (1999) JGR. je = jnh4cl a_zsr(1,je) = 0.45309 a_zsr(2,je) = 2.65606 a_zsr(3,je) = -14.7730 a_zsr(4,je) = 26.2936 a_zsr(5,je) = -20.5735 a_zsr(6,je) = 5.94255 aw_min(je) = 0.1 ! ! NaCl je = jnacl a_zsr(1,je) = 0.42922 a_zsr(2,je) = -1.17718 a_zsr(3,je) = 2.80208 a_zsr(4,je) = -4.51097 a_zsr(5,je) = 3.76963 a_zsr(6,je) = -1.31359 aw_min(je) = 0.1 ! ! NaNO3 je = jnano3 a_zsr(1,je) = 1.34966 a_zsr(2,je) = -5.20116 a_zsr(3,je) = 11.49011 a_zsr(4,je) = -14.41380 a_zsr(5,je) = 9.07037 a_zsr(6,je) = -2.29769 aw_min(je) = 0.1 ! ! Na2SO4 je = jna2so4 a_zsr(1,je) = 0.39888 a_zsr(2,je) = -1.27150 a_zsr(3,je) = 3.42792 a_zsr(4,je) = -5.92632 a_zsr(5,je) = 5.33351 a_zsr(6,je) = -1.96541 aw_min(je) = 0.1 ! ! Na3H(SO4)2 added on 1/14/2004 je = jna3hso4 a_zsr(1,je) = 0.31480 a_zsr(2,je) = -1.01087 a_zsr(3,je) = 2.44029 a_zsr(4,je) = -3.66095 a_zsr(5,je) = 2.77632 a_zsr(6,je) = -0.86058 aw_min(je) = 0.1 ! ! NaHSO4 je = jnahso4 a_zsr(1,je) = 0.62764 a_zsr(2,je) = -1.63520 a_zsr(3,je) = 4.62531 a_zsr(4,je) = -10.06925 a_zsr(5,je) = 10.33547 a_zsr(6,je) = -3.88729 aw_min(je) = 0.1 ! ! NaMSA (assumed same as NaHSO4) je = jnamsa a_zsr(1,je) = 0.62764 a_zsr(2,je) = -1.63520 a_zsr(3,je) = 4.62531 a_zsr(4,je) = -10.06925 a_zsr(5,je) = 10.33547 a_zsr(6,je) = -3.88729 aw_min(je) = 0.1 ! ! Ca(NO3)2 je = jcano3 a_zsr(1,je) = 0.38895 a_zsr(2,je) = -1.16013 a_zsr(3,je) = 2.16819 a_zsr(4,je) = -2.23079 a_zsr(5,je) = 1.00268 a_zsr(6,je) = -0.16923 aw_min(je) = 0.1 ! ! CaCl2: Kim and Seinfeld je = jcacl2 a_zsr(1,je) = 0.29891 a_zsr(2,je) = -1.31104 a_zsr(3,je) = 3.68759 a_zsr(4,je) = -5.81708 a_zsr(5,je) = 4.67520 a_zsr(6,je) = -1.53223 aw_min(je) = 0.1 ! ! H2SO4 je = jh2so4 a_zsr(1,je) = 0.32751 a_zsr(2,je) = -1.00692 a_zsr(3,je) = 2.59750 a_zsr(4,je) = -4.40014 a_zsr(5,je) = 3.88212 a_zsr(6,je) = -1.39916 aw_min(je) = 0.1 ! ! MSA (assumed same as H2SO4) je = jmsa a_zsr(1,je) = 0.32751 a_zsr(2,je) = -1.00692 a_zsr(3,je) = 2.59750 a_zsr(4,je) = -4.40014 a_zsr(5,je) = 3.88212 a_zsr(6,je) = -1.39916 aw_min(je) = 0.1 ! ! HHSO4 je = jhhso4 a_zsr(1,je) = 0.32751 a_zsr(2,je) = -1.00692 a_zsr(3,je) = 2.59750 a_zsr(4,je) = -4.40014 a_zsr(5,je) = 3.88212 a_zsr(6,je) = -1.39916 aw_min(je) = 1.0 ! ! HNO3 je = jhno3 a_zsr(1,je) = 0.75876 a_zsr(2,je) = -3.31529 a_zsr(3,je) = 9.26392 a_zsr(4,je) = -14.89799 a_zsr(5,je) = 12.08781 a_zsr(6,je) = -3.89958 aw_min(je) = 0.1 ! ! HCl je = jhcl a_zsr(1,je) = 0.31133 a_zsr(2,je) = -0.79688 a_zsr(3,je) = 1.93995 a_zsr(4,je) = -3.31582 a_zsr(5,je) = 2.93513 a_zsr(6,je) = -1.07268 aw_min(je) = 0.1 ! ! CaSO4 je = jcaso4 a_zsr(1,je) = 0.0 a_zsr(2,je) = 0.0 a_zsr(3,je) = 0.0 a_zsr(4,je) = 0.0 a_zsr(5,je) = 0.0 a_zsr(6,je) = 0.0 aw_min(je) = 1.0 ! ! Ca(MSA)2 (assumed same as Ca(NO3)2) je = jcamsa2 a_zsr(1,je) = 0.38895 a_zsr(2,je) = -1.16013 a_zsr(3,je) = 2.16819 a_zsr(4,je) = -2.23079 a_zsr(5,je) = 1.00268 a_zsr(6,je) = -0.16923 aw_min(je) = 0.1 ! ! CaCO3 je = jcaco3 a_zsr(1,je) = 0.0 a_zsr(2,je) = 0.0 a_zsr(3,je) = 0.0 a_zsr(4,je) = 0.0 a_zsr(5,je) = 0.0 a_zsr(6,je) = 0.0 aw_min(je) = 1.0 !------------------------------------------- ! b_zsr for aw => 0.97 to 0.99999 ! ! (NH4)2SO4 b_zsr(jnh4so4) = 28.0811 ! ! (NH4)3H(SO4)2 b_zsr(jlvcite) = 14.7178 ! ! NH4HSO4 b_zsr(jnh4hso4) = 29.4779 ! ! NH4MSA b_zsr(jnh4msa) = 29.4779 ! assumed same as NH4HSO4 ! ! NH4NO3 b_zsr(jnh4no3) = 33.4049 ! ! NH4Cl b_zsr(jnh4cl) = 30.8888 ! ! NaCl b_zsr(jnacl) = 29.8375 ! ! NaNO3 b_zsr(jnano3) = 32.2756 ! ! Na2SO4 b_zsr(jna2so4) = 27.6889 ! ! Na3H(SO4)2 b_zsr(jna3hso4) = 14.2184 ! ! NaHSO4 b_zsr(jnahso4) = 28.3367 ! ! NaMSA b_zsr(jnamsa) = 28.3367 ! assumed same as NaHSO4 ! ! Ca(NO3)2 b_zsr(jcano3) = 18.3661 ! ! CaCl2 b_zsr(jcacl2) = 20.8792 ! ! H2SO4 b_zsr(jh2so4) = 26.7347 ! ! HHSO4 b_zsr(jhhso4) = 26.7347 ! ! HNO3 b_zsr(jhno3) = 28.8257 ! ! HCl b_zsr(jhcl) = 27.7108 ! ! MSA b_zsr(jmsa) = 26.7347 ! assumed same as H2SO4 ! ! CaSO4 b_zsr(jcaso4) = 0.0 ! ! Ca(MSA)2 b_zsr(jcamsa2) = 18.3661 ! assumed same as Ca(NO3)2 ! ! CaCO3 b_zsr(jcaco3) = 0.0 !------------------------------------------- ! Li and Lu (2001) Surface tension model ! G_MX [mol/cm^2]; K_MX [-] ! ! (NH4)2SO4 G_MX(jnh4so4) = -8.79e-7*1.e-4 K_MX(jnh4so4) = 3.84e+1 ! ! (NH4)3H(SO4)2 G_MX(jlvcite) = -8.79e-7*1.e-4 ! assumed same as (NH4)2SO4 K_MX(jlvcite) = 3.84e+1 ! assumed same as (NH4)2SO4 ! ! NH4HSO4 G_MX(jnh4hso4) = -8.79e-7*1.e-4 ! assumed same as (NH4)2SO4 K_MX(jnh4hso4) = 3.84e+1 ! assumed same as (NH4)2SO4 ! ! NH4MSA G_MX(jnh4msa) = -8.79e-7*1.e-4 ! assumed same as (NH4)2SO4 K_MX(jnh4msa) = 3.84e+1 ! assumed same as (NH4)2SO4 ! ! NH4NO3 G_MX(jnh4no3) = -3.08e-6*1.e-4 K_MX(jnh4no3) = 4.89e-1 ! ! NH4Cl G_MX(jnh4cl) = -1.01e-6*1.e-4 K_MX(jnh4cl) = 1.3 ! ! NaCl G_MX(jnacl) = -1.05e-6*1.e-4 K_MX(jnacl) = 1.2 ! ! NaNO3 G_MX(jnano3) = -1.66e-6*1.e-4 K_MX(jnano3) = 1.25 ! ! Na2SO4 G_MX(jna2so4) = -8.37e-7*1.e-4 K_MX(jna2so4) = 7.57e+1 ! ! Na3H(SO4)2 G_MX(jna3hso4) = -8.37e-7*1.e-4 ! assumed same as Na2SO4 K_MX(jna3hso4) = 7.57e+1 ! assumed same as Na2SO4 ! ! NaHSO4 G_MX(jnahso4) = -8.37e-7*1.e-4 ! assumed same as Na2SO4 K_MX(jnahso4) = 7.57e+1 ! assumed same as Na2SO4 ! ! NaMSA G_MX(jnamsa) = -8.37e-7*1.e-4 K_MX(jnamsa) = 7.57e+1 ! ! Ca(NO3)2 G_MX(jcano3) = -4.88e-7*1.e-4 ! assumed same as CaCl2 K_MX(jcano3) = 1.50e+1 ! assumed same as CaCl2 ! ! CaCl2 G_MX(jcacl2) = -4.88e-7*1.e-4 K_MX(jcacl2) = 1.50e+1 ! ! H2SO4 G_MX(jh2so4) = -6.75e-8*1.e-4 K_MX(jh2so4) = 1.65e+3 ! ! HHSO4 G_MX(jh2so4) = -6.75e-8*1.e-4 ! assumed same as H2SO4 K_MX(jh2so4) = 1.65e+3 ! assumed same as H2SO4 ! ! HNO3 G_MX(jhno3) = 8.05e-7*1.e-4 K_MX(jhno3) = 1.06e-1 ! ! HCl G_MX(jhcl) = 4.12e-7*1.e-4 K_MX(jhcl) = 4.68e-3 !& ! MSA G_MX(jmsa) = 8.05e-7*1.e-4 ! assumed same as HNO3 K_MX(jmsa) = 1.06e-1 ! assumed same as HNO3 ! ! CaSO4 G_MX(jmsa) = 0.0*1.e-4 ! assumed K_MX(jmsa) = 0.0 ! assumed ! ! Ca(MSA)2 G_MX(jcamsa2) = 0.0*1.e-4 ! assumed K_MX(jcamsa2) = 0.0 ! assumed ! ! CaCO3 G_MX(jcaco3) = 0.0*1.e-4 ! assumed K_MX(jcaco3) = 0.0 ! assumed !---------------------------------------------------------------- ! parameters for MTEM mixing rule (Zaveri, Easter, and Wexler, 2005) ! log_gamZ(jA,jE) A in E !---------------------------------------------------------------- ! 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 ! (NH4)2SO4 in E jA = jnh4so4 ! in (NH4)2SO4 jE = jnh4so4 b_mtem(1,jA,jE) = -2.94685 b_mtem(2,jA,jE) = 17.3328 b_mtem(3,jA,jE) = -64.8441 b_mtem(4,jA,jE) = 122.7070 b_mtem(5,jA,jE) = -114.4373 b_mtem(6,jA,jE) = 41.6811 ! in NH4NO3 jE = jnh4no3 b_mtem(1,jA,jE) = -2.7503 b_mtem(2,jA,jE) = 4.3806 b_mtem(3,jA,jE) = -1.1110 b_mtem(4,jA,jE) = -1.7005 b_mtem(5,jA,jE) = -4.4207 b_mtem(6,jA,jE) = 5.1990 ! in NH4Cl (revised on 11/15/2003) jE = jnh4cl b_mtem(1,jA,jE) = -2.06952 b_mtem(2,jA,jE) = 7.1240 b_mtem(3,jA,jE) = -24.4274 b_mtem(4,jA,jE) = 51.1458 b_mtem(5,jA,jE) = -54.2056 b_mtem(6,jA,jE) = 22.0606 ! in Na2SO4 jE = jna2so4 b_mtem(1,jA,jE) = -2.17361 b_mtem(2,jA,jE) = 15.9919 b_mtem(3,jA,jE) = -69.0952 b_mtem(4,jA,jE) = 139.8860 b_mtem(5,jA,jE) = -134.9890 b_mtem(6,jA,jE) = 49.8877 ! in NaNO3 jE = jnano3 b_mtem(1,jA,jE) = -4.4370 b_mtem(2,jA,jE) = 24.0243 b_mtem(3,jA,jE) = -76.2437 b_mtem(4,jA,jE) = 128.6660 b_mtem(5,jA,jE) = -110.0900 b_mtem(6,jA,jE) = 37.7414 ! in NaCl jE = jnacl b_mtem(1,jA,jE) = -1.5394 b_mtem(2,jA,jE) = 5.8671 b_mtem(3,jA,jE) = -22.7726 b_mtem(4,jA,jE) = 47.0547 b_mtem(5,jA,jE) = -47.8266 b_mtem(6,jA,jE) = 18.8489 ! in HNO3 jE = jhno3 b_mtem(1,jA,jE) = -0.35750 b_mtem(2,jA,jE) = -3.82466 b_mtem(3,jA,jE) = 4.55462 b_mtem(4,jA,jE) = 5.05402 b_mtem(5,jA,jE) = -14.7476 b_mtem(6,jA,jE) = 8.8009 ! in HCl jE = jhcl b_mtem(1,jA,jE) = -2.15146 b_mtem(2,jA,jE) = 5.50205 b_mtem(3,jA,jE) = -19.1476 b_mtem(4,jA,jE) = 39.1880 b_mtem(5,jA,jE) = -39.9460 b_mtem(6,jA,jE) = 16.0700 ! in H2SO4 jE = jh2so4 b_mtem(1,jA,jE) = -2.52604 b_mtem(2,jA,jE) = 9.76022 b_mtem(3,jA,jE) = -35.2540 b_mtem(4,jA,jE) = 71.2981 b_mtem(5,jA,jE) = -71.8207 b_mtem(6,jA,jE) = 28.0758 ! ! in NH4HSO4 jE = jnh4hso4 b_mtem(1,jA,jE) = -4.13219 b_mtem(2,jA,jE) = 13.8863 b_mtem(3,jA,jE) = -34.5387 b_mtem(4,jA,jE) = 56.5012 b_mtem(5,jA,jE) = -51.8702 b_mtem(6,jA,jE) = 19.6232 ! ! in (NH4)3H(SO4)2 jE = jlvcite b_mtem(1,jA,jE) = -2.53482 b_mtem(2,jA,jE) = 12.3333 b_mtem(3,jA,jE) = -46.1020 b_mtem(4,jA,jE) = 90.4775 b_mtem(5,jA,jE) = -88.1254 b_mtem(6,jA,jE) = 33.4715 ! ! in NaHSO4 jE = jnahso4 b_mtem(1,jA,jE) = -3.23425 b_mtem(2,jA,jE) = 18.7842 b_mtem(3,jA,jE) = -78.7807 b_mtem(4,jA,jE) = 161.517 b_mtem(5,jA,jE) = -154.940 b_mtem(6,jA,jE) = 56.2252 ! ! in Na3H(SO4)2 jE = jna3hso4 b_mtem(1,jA,jE) = -1.25316 b_mtem(2,jA,jE) = 7.40960 b_mtem(3,jA,jE) = -34.8929 b_mtem(4,jA,jE) = 72.8853 b_mtem(5,jA,jE) = -72.4503 b_mtem(6,jA,jE) = 27.7706 !----------------- ! NH4NO3 in E jA = jnh4no3 ! in (NH4)2SO4 jE = jnh4so4 b_mtem(1,jA,jE) = -3.5201 b_mtem(2,jA,jE) = 21.6584 b_mtem(3,jA,jE) = -72.1499 b_mtem(4,jA,jE) = 126.7000 b_mtem(5,jA,jE) = -111.4550 b_mtem(6,jA,jE) = 38.5677 ! in NH4NO3 jE = jnh4no3 b_mtem(1,jA,jE) = -2.2630 b_mtem(2,jA,jE) = -0.1518 b_mtem(3,jA,jE) = 17.0898 b_mtem(4,jA,jE) = -36.7832 b_mtem(5,jA,jE) = 29.8407 b_mtem(6,jA,jE) = -7.9314 ! in NH4Cl (revised on 11/15/2003) jE = jnh4cl b_mtem(1,jA,jE) = -1.3851 b_mtem(2,jA,jE) = -0.4462 b_mtem(3,jA,jE) = 8.4567 b_mtem(4,jA,jE) = -11.5988 b_mtem(5,jA,jE) = 2.9802 b_mtem(6,jA,jE) = 1.8132 ! in Na2SO4 jE = jna2so4 b_mtem(1,jA,jE) = -1.7602 b_mtem(2,jA,jE) = 10.4044 b_mtem(3,jA,jE) = -35.5894 b_mtem(4,jA,jE) = 64.3584 b_mtem(5,jA,jE) = -57.8931 b_mtem(6,jA,jE) = 20.2141 ! in NaNO3 jE = jnano3 b_mtem(1,jA,jE) = -3.24346 b_mtem(2,jA,jE) = 16.2794 b_mtem(3,jA,jE) = -48.7601 b_mtem(4,jA,jE) = 79.2246 b_mtem(5,jA,jE) = -65.8169 b_mtem(6,jA,jE) = 22.1500 ! in NaCl jE = jnacl b_mtem(1,jA,jE) = -1.75658 b_mtem(2,jA,jE) = 7.71384 b_mtem(3,jA,jE) = -22.7984 b_mtem(4,jA,jE) = 39.1532 b_mtem(5,jA,jE) = -34.6165 b_mtem(6,jA,jE) = 12.1283 ! in Ca(NO3)2 jE = jcano3 b_mtem(1,jA,jE) = -0.97178 b_mtem(2,jA,jE) = 6.61964 b_mtem(3,jA,jE) = -26.2353 b_mtem(4,jA,jE) = 50.5259 b_mtem(5,jA,jE) = -47.6586 b_mtem(6,jA,jE) = 17.5074 ! in CaCl2 added on 12/22/2003 jE = jcacl2 b_mtem(1,jA,jE) = -0.41515 b_mtem(2,jA,jE) = 6.44101 b_mtem(3,jA,jE) = -26.4473 b_mtem(4,jA,jE) = 49.0718 b_mtem(5,jA,jE) = -44.2631 b_mtem(6,jA,jE) = 15.3771 ! in HNO3 jE = jhno3 b_mtem(1,jA,jE) = -1.20644 b_mtem(2,jA,jE) = 5.70117 b_mtem(3,jA,jE) = -18.2783 b_mtem(4,jA,jE) = 31.7199 b_mtem(5,jA,jE) = -27.8703 b_mtem(6,jA,jE) = 9.7299 ! in HCl jE = jhcl b_mtem(1,jA,jE) = -0.680862 b_mtem(2,jA,jE) = 3.59456 b_mtem(3,jA,jE) = -10.7969 b_mtem(4,jA,jE) = 17.8434 b_mtem(5,jA,jE) = -15.3165 b_mtem(6,jA,jE) = 5.17123 !---------- ! NH4Cl in E jA = jnh4cl ! in (NH4)2SO4 jE = jnh4so4 b_mtem(1,jA,jE) = -2.8850 b_mtem(2,jA,jE) = 20.6970 b_mtem(3,jA,jE) = -70.6810 b_mtem(4,jA,jE) = 124.3690 b_mtem(5,jA,jE) = -109.2880 b_mtem(6,jA,jE) = 37.5831 ! in NH4NO3 jE = jnh4no3 b_mtem(1,jA,jE) = -1.9386 b_mtem(2,jA,jE) = 1.3238 b_mtem(3,jA,jE) = 11.8500 b_mtem(4,jA,jE) = -28.1168 b_mtem(5,jA,jE) = 21.8543 b_mtem(6,jA,jE) = -5.1671 ! in NH4Cl (revised on 11/15/2003) jE = jnh4cl b_mtem(1,jA,jE) = -0.9559 b_mtem(2,jA,jE) = 0.8121 b_mtem(3,jA,jE) = 4.3644 b_mtem(4,jA,jE) = -8.9258 b_mtem(5,jA,jE) = 4.2362 b_mtem(6,jA,jE) = 0.2891 ! in Na2SO4 jE = jna2so4 b_mtem(1,jA,jE) = 0.0377 b_mtem(2,jA,jE) = 6.0752 b_mtem(3,jA,jE) = -30.8641 b_mtem(4,jA,jE) = 63.3095 b_mtem(5,jA,jE) = -61.0070 b_mtem(6,jA,jE) = 22.1734 ! in NaNO3 jE = jnano3 b_mtem(1,jA,jE) = -1.8336 b_mtem(2,jA,jE) = 12.8160 b_mtem(3,jA,jE) = -42.3388 b_mtem(4,jA,jE) = 71.1816 b_mtem(5,jA,jE) = -60.5708 b_mtem(6,jA,jE) = 20.5853 ! in NaCl jE = jnacl b_mtem(1,jA,jE) = -0.1429 b_mtem(2,jA,jE) = 2.3561 b_mtem(3,jA,jE) = -10.4425 b_mtem(4,jA,jE) = 20.8951 b_mtem(5,jA,jE) = -20.7739 b_mtem(6,jA,jE) = 7.9355 ! in Ca(NO3)2 jE = jcano3 b_mtem(1,jA,jE) = 0.76235 b_mtem(2,jA,jE) = 3.08323 b_mtem(3,jA,jE) = -23.6772 b_mtem(4,jA,jE) = 53.7415 b_mtem(5,jA,jE) = -55.4043 b_mtem(6,jA,jE) = 21.2944 ! in CaCl2 (revised on 11/27/2003) jE = jcacl2 b_mtem(1,jA,jE) = 1.13864 b_mtem(2,jA,jE) = -0.340539 b_mtem(3,jA,jE) = -8.67025 b_mtem(4,jA,jE) = 22.8008 b_mtem(5,jA,jE) = -24.5181 b_mtem(6,jA,jE) = 9.3663 ! in HNO3 jE = jhno3 b_mtem(1,jA,jE) = 2.42532 b_mtem(2,jA,jE) = -14.1755 b_mtem(3,jA,jE) = 38.804 b_mtem(4,jA,jE) = -58.2437 b_mtem(5,jA,jE) = 43.5431 b_mtem(6,jA,jE) = -12.5824 ! in HCl jE = jhcl b_mtem(1,jA,jE) = 0.330337 b_mtem(2,jA,jE) = 0.0778934 b_mtem(3,jA,jE) = -2.30492 b_mtem(4,jA,jE) = 4.73003 b_mtem(5,jA,jE) = -4.80849 b_mtem(6,jA,jE) = 1.78866 !---------- ! Na2SO4 in E jA = jna2so4 ! in (NH4)2SO4 jE = jnh4so4 b_mtem(1,jA,jE) = -2.6982 b_mtem(2,jA,jE) = 22.9875 b_mtem(3,jA,jE) = -98.9840 b_mtem(4,jA,jE) = 198.0180 b_mtem(5,jA,jE) = -188.7270 b_mtem(6,jA,jE) = 69.0548 ! in NH4NO3 jE = jnh4no3 b_mtem(1,jA,jE) = -2.4844 b_mtem(2,jA,jE) = 6.5420 b_mtem(3,jA,jE) = -9.8998 b_mtem(4,jA,jE) = 11.3884 b_mtem(5,jA,jE) = -13.6842 b_mtem(6,jA,jE) = 7.7411 ! in NH4Cl (revised on 11/15/2003) jE = jnh4cl b_mtem(1,jA,jE) = -1.3325 b_mtem(2,jA,jE) = 13.0406 b_mtem(3,jA,jE) = -56.1935 b_mtem(4,jA,jE) = 107.1170 b_mtem(5,jA,jE) = -97.3721 b_mtem(6,jA,jE) = 34.3763 ! in Na2SO4 jE = jna2so4 b_mtem(1,jA,jE) = -1.2832 b_mtem(2,jA,jE) = 12.8526 b_mtem(3,jA,jE) = -62.2087 b_mtem(4,jA,jE) = 130.3876 b_mtem(5,jA,jE) = -128.2627 b_mtem(6,jA,jE) = 48.0340 ! in NaNO3 jE = jnano3 b_mtem(1,jA,jE) = -3.5384 b_mtem(2,jA,jE) = 21.3758 b_mtem(3,jA,jE) = -70.7638 b_mtem(4,jA,jE) = 121.1580 b_mtem(5,jA,jE) = -104.6230 b_mtem(6,jA,jE) = 36.0557 ! in NaCl jE = jnacl b_mtem(1,jA,jE) = 0.2175 b_mtem(2,jA,jE) = -0.5648 b_mtem(3,jA,jE) = -8.0288 b_mtem(4,jA,jE) = 25.9734 b_mtem(5,jA,jE) = -32.3577 b_mtem(6,jA,jE) = 14.3924 ! in HNO3 jE = jhno3 b_mtem(1,jA,jE) = -0.309617 b_mtem(2,jA,jE) = -1.82899 b_mtem(3,jA,jE) = -1.5505 b_mtem(4,jA,jE) = 13.3847 b_mtem(5,jA,jE) = -20.1284 b_mtem(6,jA,jE) = 9.93163 ! in HCl jE = jhcl b_mtem(1,jA,jE) = -0.259455 b_mtem(2,jA,jE) = -0.819366 b_mtem(3,jA,jE) = -4.28964 b_mtem(4,jA,jE) = 16.4305 b_mtem(5,jA,jE) = -21.8546 b_mtem(6,jA,jE) = 10.3044 ! in H2SO4 jE = jh2so4 b_mtem(1,jA,jE) = -1.84257 b_mtem(2,jA,jE) = 7.85788 b_mtem(3,jA,jE) = -29.9275 b_mtem(4,jA,jE) = 61.7515 b_mtem(5,jA,jE) = -63.2308 b_mtem(6,jA,jE) = 24.9542 ! in NH4HSO4 jE = jnh4hso4 b_mtem(1,jA,jE) = -1.05891 b_mtem(2,jA,jE) = 2.84831 b_mtem(3,jA,jE) = -21.1827 b_mtem(4,jA,jE) = 57.5175 b_mtem(5,jA,jE) = -64.8120 b_mtem(6,jA,jE) = 26.1986 ! in (NH4)3H(SO4)2 jE = jlvcite b_mtem(1,jA,jE) = -1.16584 b_mtem(2,jA,jE) = 8.50075 b_mtem(3,jA,jE) = -44.3420 b_mtem(4,jA,jE) = 97.3974 b_mtem(5,jA,jE) = -98.4549 b_mtem(6,jA,jE) = 37.6104 ! in NaHSO4 jE = jnahso4 b_mtem(1,jA,jE) = -1.95805 b_mtem(2,jA,jE) = 6.62417 b_mtem(3,jA,jE) = -31.8072 b_mtem(4,jA,jE) = 77.8603 b_mtem(5,jA,jE) = -84.6458 b_mtem(6,jA,jE) = 33.4963 ! in Na3H(SO4)2 jE = jna3hso4 b_mtem(1,jA,jE) = -0.36045 b_mtem(2,jA,jE) = 3.55223 b_mtem(3,jA,jE) = -24.0327 b_mtem(4,jA,jE) = 54.4879 b_mtem(5,jA,jE) = -56.6531 b_mtem(6,jA,jE) = 22.4956 !---------- ! NaNO3 in E jA = jnano3 ! in (NH4)2SO4 jE = jnh4so4 b_mtem(1,jA,jE) = -2.5888 b_mtem(2,jA,jE) = 17.6192 b_mtem(3,jA,jE) = -63.2183 b_mtem(4,jA,jE) = 115.3520 b_mtem(5,jA,jE) = -104.0860 b_mtem(6,jA,jE) = 36.7390 ! in NH4NO3 jE = jnh4no3 b_mtem(1,jA,jE) = -2.0669 b_mtem(2,jA,jE) = 1.4792 b_mtem(3,jA,jE) = 10.5261 b_mtem(4,jA,jE) = -27.0987 b_mtem(5,jA,jE) = 23.0591 b_mtem(6,jA,jE) = -6.0938 ! in NH4Cl (revised on 11/15/2003) jE = jnh4cl b_mtem(1,jA,jE) = -0.8325 b_mtem(2,jA,jE) = 3.9933 b_mtem(3,jA,jE) = -15.3789 b_mtem(4,jA,jE) = 30.4050 b_mtem(5,jA,jE) = -29.4204 b_mtem(6,jA,jE) = 11.0597 ! in Na2SO4 jE = jna2so4 b_mtem(1,jA,jE) = -1.1233 b_mtem(2,jA,jE) = 8.3998 b_mtem(3,jA,jE) = -31.9002 b_mtem(4,jA,jE) = 60.1450 b_mtem(5,jA,jE) = -55.5503 b_mtem(6,jA,jE) = 19.7757 ! in NaNO3 jE = jnano3 b_mtem(1,jA,jE) = -2.5386 b_mtem(2,jA,jE) = 13.9039 b_mtem(3,jA,jE) = -42.8467 b_mtem(4,jA,jE) = 69.7442 b_mtem(5,jA,jE) = -57.8988 b_mtem(6,jA,jE) = 19.4635 ! in NaCl jE = jnacl b_mtem(1,jA,jE) = -0.4351 b_mtem(2,jA,jE) = 2.8311 b_mtem(3,jA,jE) = -11.4485 b_mtem(4,jA,jE) = 22.7201 b_mtem(5,jA,jE) = -22.4228 b_mtem(6,jA,jE) = 8.5792 ! in Ca(NO3)2 jE = jcano3 b_mtem(1,jA,jE) = -0.72060 b_mtem(2,jA,jE) = 5.64915 b_mtem(3,jA,jE) = -23.5020 b_mtem(4,jA,jE) = 46.0078 b_mtem(5,jA,jE) = -43.8075 b_mtem(6,jA,jE) = 16.1652 ! in CaCl2 jE = jcacl2 b_mtem(1,jA,jE) = 0.003928 b_mtem(2,jA,jE) = 3.54724 b_mtem(3,jA,jE) = -18.6057 b_mtem(4,jA,jE) = 38.1445 b_mtem(5,jA,jE) = -36.7745 b_mtem(6,jA,jE) = 13.4529 ! in HNO3 jE = jhno3 b_mtem(1,jA,jE) = -1.1712 b_mtem(2,jA,jE) = 7.20907 b_mtem(3,jA,jE) = -22.9215 b_mtem(4,jA,jE) = 38.1257 b_mtem(5,jA,jE) = -32.0759 b_mtem(6,jA,jE) = 10.6443 ! in HCl jE = jhcl b_mtem(1,jA,jE) = 0.738022 b_mtem(2,jA,jE) = -1.14313 b_mtem(3,jA,jE) = 0.32251 b_mtem(4,jA,jE) = 0.838679 b_mtem(5,jA,jE) = -1.81747 b_mtem(6,jA,jE) = 0.873986 !---------- ! NaCl in E jA = jnacl ! in (NH4)2SO4 jE = jnh4so4 b_mtem(1,jA,jE) = -1.9525 b_mtem(2,jA,jE) = 16.6433 b_mtem(3,jA,jE) = -61.7090 b_mtem(4,jA,jE) = 112.9910 b_mtem(5,jA,jE) = -101.9370 b_mtem(6,jA,jE) = 35.7760 ! in NH4NO3 jE = jnh4no3 b_mtem(1,jA,jE) = -1.7525 b_mtem(2,jA,jE) = 3.0713 b_mtem(3,jA,jE) = 4.8063 b_mtem(4,jA,jE) = -17.5334 b_mtem(5,jA,jE) = 14.2872 b_mtem(6,jA,jE) = -3.0690 ! in NH4Cl (revised on 11/15/2003) jE = jnh4cl b_mtem(1,jA,jE) = -0.4021 b_mtem(2,jA,jE) = 5.2399 b_mtem(3,jA,jE) = -19.4278 b_mtem(4,jA,jE) = 33.0027 b_mtem(5,jA,jE) = -28.1020 b_mtem(6,jA,jE) = 9.5159 ! in Na2SO4 jE = jna2so4 b_mtem(1,jA,jE) = 0.6692 b_mtem(2,jA,jE) = 4.1207 b_mtem(3,jA,jE) = -27.3314 b_mtem(4,jA,jE) = 59.3112 b_mtem(5,jA,jE) = -58.7998 b_mtem(6,jA,jE) = 21.7674 ! in NaNO3 jE = jnano3 b_mtem(1,jA,jE) = -1.17444 b_mtem(2,jA,jE) = 10.9927 b_mtem(3,jA,jE) = -38.9013 b_mtem(4,jA,jE) = 66.8521 b_mtem(5,jA,jE) = -57.6564 b_mtem(6,jA,jE) = 19.7296 ! in NaCl jE = jnacl b_mtem(1,jA,jE) = 1.17679 b_mtem(2,jA,jE) = -2.5061 b_mtem(3,jA,jE) = 0.8508 b_mtem(4,jA,jE) = 4.4802 b_mtem(5,jA,jE) = -8.4945 b_mtem(6,jA,jE) = 4.3182 ! in Ca(NO3)2 jE = jcano3 b_mtem(1,jA,jE) = 1.01450 b_mtem(2,jA,jE) = 2.10260 b_mtem(3,jA,jE) = -20.9036 b_mtem(4,jA,jE) = 49.1481 b_mtem(5,jA,jE) = -51.4867 b_mtem(6,jA,jE) = 19.9301 ! in CaCl2 (PSC92: revised on 11/27/2003) jE = jcacl2 b_mtem(1,jA,jE) = 1.55463 b_mtem(2,jA,jE) = -3.20122 b_mtem(3,jA,jE) = -0.957075 b_mtem(4,jA,jE) = 12.103 b_mtem(5,jA,jE) = -17.221 b_mtem(6,jA,jE) = 7.50264 ! in HNO3 jE = jhno3 b_mtem(1,jA,jE) = 2.46187 b_mtem(2,jA,jE) = -12.6845 b_mtem(3,jA,jE) = 34.2383 b_mtem(4,jA,jE) = -51.9992 b_mtem(5,jA,jE) = 39.4934 b_mtem(6,jA,jE) = -11.7247 ! in HCl jE = jhcl b_mtem(1,jA,jE) = 1.74915 b_mtem(2,jA,jE) = -4.65768 b_mtem(3,jA,jE) = 8.80287 b_mtem(4,jA,jE) = -12.2503 b_mtem(5,jA,jE) = 8.668751 b_mtem(6,jA,jE) = -2.50158 !---------- ! Ca(NO3)2 in E jA = jcano3 ! in NH4NO3 jE = jnh4no3 b_mtem(1,jA,jE) = -1.86260 b_mtem(2,jA,jE) = 11.6178 b_mtem(3,jA,jE) = -30.9069 b_mtem(4,jA,jE) = 41.7578 b_mtem(5,jA,jE) = -33.7338 b_mtem(6,jA,jE) = 12.7541 ! in NH4Cl (revised on 11/15/2003) jE = jnh4cl b_mtem(1,jA,jE) = -1.1798 b_mtem(2,jA,jE) = 25.9608 b_mtem(3,jA,jE) = -98.9373 b_mtem(4,jA,jE) = 160.2300 b_mtem(5,jA,jE) = -125.9540 b_mtem(6,jA,jE) = 39.5130 ! in NaNO3 jE = jnano3 b_mtem(1,jA,jE) = -1.44384 b_mtem(2,jA,jE) = 13.6044 b_mtem(3,jA,jE) = -54.4300 b_mtem(4,jA,jE) = 100.582 b_mtem(5,jA,jE) = -91.2364 b_mtem(6,jA,jE) = 32.5970 ! in NaCl jE = jnacl b_mtem(1,jA,jE) = -0.099114 b_mtem(2,jA,jE) = 2.84091 b_mtem(3,jA,jE) = -16.9229 b_mtem(4,jA,jE) = 37.4839 b_mtem(5,jA,jE) = -39.5132 b_mtem(6,jA,jE) = 15.8564 ! in Ca(NO3)2 jE = jcano3 b_mtem(1,jA,jE) = 0.055116 b_mtem(2,jA,jE) = 4.58610 b_mtem(3,jA,jE) = -27.6629 b_mtem(4,jA,jE) = 60.8288 b_mtem(5,jA,jE) = -61.4988 b_mtem(6,jA,jE) = 23.3136 ! in CaCl2 (PSC92: revised on 11/27/2003) jE = jcacl2 b_mtem(1,jA,jE) = 1.57155 b_mtem(2,jA,jE) = -3.18486 b_mtem(3,jA,jE) = -3.35758 b_mtem(4,jA,jE) = 18.7501 b_mtem(5,jA,jE) = -24.5604 b_mtem(6,jA,jE) = 10.3798 ! in HNO3 jE = jhno3 b_mtem(1,jA,jE) = 1.04446 b_mtem(2,jA,jE) = -3.19066 b_mtem(3,jA,jE) = 2.44714 b_mtem(4,jA,jE) = 2.07218 b_mtem(5,jA,jE) = -6.43949 b_mtem(6,jA,jE) = 3.66471 ! in HCl jE = jhcl b_mtem(1,jA,jE) = 1.05723 b_mtem(2,jA,jE) = -1.46826 b_mtem(3,jA,jE) = -1.0713 b_mtem(4,jA,jE) = 4.64439 b_mtem(5,jA,jE) = -6.32402 b_mtem(6,jA,jE) = 2.78202 !---------- ! CaCl2 in E jA = jcacl2 ! in NH4NO3 (PSC92: revised on 12/22/2003) jE = jnh4no3 b_mtem(1,jA,jE) = -1.43626 b_mtem(2,jA,jE) = 13.6598 b_mtem(3,jA,jE) = -38.2068 b_mtem(4,jA,jE) = 53.9057 b_mtem(5,jA,jE) = -44.9018 b_mtem(6,jA,jE) = 16.6120 ! in NH4Cl (PSC92: revised on 11/27/2003) jE = jnh4cl b_mtem(1,jA,jE) = -0.603965 b_mtem(2,jA,jE) = 27.6027 b_mtem(3,jA,jE) = -104.258 b_mtem(4,jA,jE) = 163.553 b_mtem(5,jA,jE) = -124.076 b_mtem(6,jA,jE) = 37.4153 ! in NaNO3 (PSC92: revised on 12/22/2003) jE = jnano3 b_mtem(1,jA,jE) = 0.44648 b_mtem(2,jA,jE) = 8.8850 b_mtem(3,jA,jE) = -45.5232 b_mtem(4,jA,jE) = 89.3263 b_mtem(5,jA,jE) = -83.8604 b_mtem(6,jA,jE) = 30.4069 ! in NaCl (PSC92: revised on 11/27/2003) jE = jnacl b_mtem(1,jA,jE) = 1.61927 b_mtem(2,jA,jE) = 0.247547 b_mtem(3,jA,jE) = -18.1252 b_mtem(4,jA,jE) = 45.2479 b_mtem(5,jA,jE) = -48.6072 b_mtem(6,jA,jE) = 19.2784 ! in Ca(NO3)2 (PSC92: revised on 11/27/2003) jE = jcano3 b_mtem(1,jA,jE) = 2.36667 b_mtem(2,jA,jE) = -0.123309 b_mtem(3,jA,jE) = -24.2723 b_mtem(4,jA,jE) = 65.1486 b_mtem(5,jA,jE) = -71.8504 b_mtem(6,jA,jE) = 28.3696 ! in CaCl2 (PSC92: revised on 11/27/2003) jE = jcacl2 b_mtem(1,jA,jE) = 3.64023 b_mtem(2,jA,jE) = -12.1926 b_mtem(3,jA,jE) = 20.2028 b_mtem(4,jA,jE) = -16.0056 b_mtem(5,jA,jE) = 1.52355 b_mtem(6,jA,jE) = 2.44709 ! in HNO3 jE = jhno3 b_mtem(1,jA,jE) = 5.88794 b_mtem(2,jA,jE) = -29.7083 b_mtem(3,jA,jE) = 78.6309 b_mtem(4,jA,jE) = -118.037 b_mtem(5,jA,jE) = 88.932 b_mtem(6,jA,jE) = -26.1407 ! in HCl jE = jhcl b_mtem(1,jA,jE) = 2.40628 b_mtem(2,jA,jE) = -6.16566 b_mtem(3,jA,jE) = 10.2851 b_mtem(4,jA,jE) = -12.9035 b_mtem(5,jA,jE) = 7.7441 b_mtem(6,jA,jE) = -1.74821 !---------- ! HNO3 in E jA = jhno3 ! in (NH4)2SO4 jE = jnh4so4 b_mtem(1,jA,jE) = -3.57598 b_mtem(2,jA,jE) = 21.5469 b_mtem(3,jA,jE) = -77.4111 b_mtem(4,jA,jE) = 144.136 b_mtem(5,jA,jE) = -132.849 b_mtem(6,jA,jE) = 47.9412 ! in NH4NO3 jE = jnh4no3 b_mtem(1,jA,jE) = -2.00209 b_mtem(2,jA,jE) = -3.48399 b_mtem(3,jA,jE) = 34.9906 b_mtem(4,jA,jE) = -68.6653 b_mtem(5,jA,jE) = 54.0992 b_mtem(6,jA,jE) = -15.1343 ! in NH4Cl revised on 12/22/2003 jE = jnh4cl b_mtem(1,jA,jE) = -0.63790 b_mtem(2,jA,jE) = -1.67730 b_mtem(3,jA,jE) = 10.1727 b_mtem(4,jA,jE) = -14.9097 b_mtem(5,jA,jE) = 7.67410 b_mtem(6,jA,jE) = -0.79586 ! in NaCl jE = jnacl b_mtem(1,jA,jE) = 1.3446 b_mtem(2,jA,jE) = -2.5578 b_mtem(3,jA,jE) = 1.3464 b_mtem(4,jA,jE) = 2.90537 b_mtem(5,jA,jE) = -6.53014 b_mtem(6,jA,jE) = 3.31339 ! in NaNO3 jE = jnano3 b_mtem(1,jA,jE) = -0.546636 b_mtem(2,jA,jE) = 10.3127 b_mtem(3,jA,jE) = -39.9603 b_mtem(4,jA,jE) = 71.4609 b_mtem(5,jA,jE) = -63.4958 b_mtem(6,jA,jE) = 22.0679 ! in Na2SO4 jE = jna2so4 b_mtem(1,jA,jE) = 1.35059 b_mtem(2,jA,jE) = 4.34557 b_mtem(3,jA,jE) = -35.8425 b_mtem(4,jA,jE) = 80.9868 b_mtem(5,jA,jE) = -81.6544 b_mtem(6,jA,jE) = 30.4841 ! in Ca(NO3)2 jE = jcano3 b_mtem(1,jA,jE) = 0.869414 b_mtem(2,jA,jE) = 2.98486 b_mtem(3,jA,jE) = -22.255 b_mtem(4,jA,jE) = 50.1863 b_mtem(5,jA,jE) = -51.214 b_mtem(6,jA,jE) = 19.2235 ! in CaCl2 (KM) revised on 12/22/2003 jE = jcacl2 b_mtem(1,jA,jE) = 1.42800 b_mtem(2,jA,jE) = -1.78959 b_mtem(3,jA,jE) = -2.49075 b_mtem(4,jA,jE) = 10.1877 b_mtem(5,jA,jE) = -12.1948 b_mtem(6,jA,jE) = 4.64475 ! in HNO3 (added on 12/06/2004) jE = jhno3 b_mtem(1,jA,jE) = 0.22035 b_mtem(2,jA,jE) = 2.94973 b_mtem(3,jA,jE) = -12.1469 b_mtem(4,jA,jE) = 20.4905 b_mtem(5,jA,jE) = -17.3966 b_mtem(6,jA,jE) = 5.70779 ! in HCl (added on 12/06/2004) jE = jhcl b_mtem(1,jA,jE) = 1.55503 b_mtem(2,jA,jE) = -3.61226 b_mtem(3,jA,jE) = 6.28265 b_mtem(4,jA,jE) = -8.69575 b_mtem(5,jA,jE) = 6.09372 b_mtem(6,jA,jE) = -1.80898 ! in H2SO4 jE = jh2so4 b_mtem(1,jA,jE) = 1.10783 b_mtem(2,jA,jE) = -1.3363 b_mtem(3,jA,jE) = -1.83525 b_mtem(4,jA,jE) = 7.47373 b_mtem(5,jA,jE) = -9.72954 b_mtem(6,jA,jE) = 4.12248 ! in NH4HSO4 jE = jnh4hso4 b_mtem(1,jA,jE) = -0.851026 b_mtem(2,jA,jE) = 12.2515 b_mtem(3,jA,jE) = -49.788 b_mtem(4,jA,jE) = 91.6215 b_mtem(5,jA,jE) = -81.4877 b_mtem(6,jA,jE) = 28.0002 ! in (NH4)3H(SO4)2 jE = jlvcite b_mtem(1,jA,jE) = -3.09464 b_mtem(2,jA,jE) = 14.9303 b_mtem(3,jA,jE) = -43.0454 b_mtem(4,jA,jE) = 72.6695 b_mtem(5,jA,jE) = -65.2140 b_mtem(6,jA,jE) = 23.4814 ! in NaHSO4 jE = jnahso4 b_mtem(1,jA,jE) = 1.22973 b_mtem(2,jA,jE) = 2.82702 b_mtem(3,jA,jE) = -17.5869 b_mtem(4,jA,jE) = 28.9564 b_mtem(5,jA,jE) = -23.5814 b_mtem(6,jA,jE) = 7.91153 ! in Na3H(SO4)2 jE = jna3hso4 b_mtem(1,jA,jE) = 1.64773 b_mtem(2,jA,jE) = 0.94188 b_mtem(3,jA,jE) = -19.1242 b_mtem(4,jA,jE) = 46.9887 b_mtem(5,jA,jE) = -50.9494 b_mtem(6,jA,jE) = 20.2169 !---------- ! HCl in E jA = jhcl ! in (NH4)2SO4 jE = jnh4so4 b_mtem(1,jA,jE) = -2.93783 b_mtem(2,jA,jE) = 20.5546 b_mtem(3,jA,jE) = -75.8548 b_mtem(4,jA,jE) = 141.729 b_mtem(5,jA,jE) = -130.697 b_mtem(6,jA,jE) = 46.9905 ! in NH4NO3 jE = jnh4no3 b_mtem(1,jA,jE) = -1.69063 b_mtem(2,jA,jE) = -1.85303 b_mtem(3,jA,jE) = 29.0927 b_mtem(4,jA,jE) = -58.7401 b_mtem(5,jA,jE) = 44.999 b_mtem(6,jA,jE) = -11.9988 ! in NH4Cl (revised on 11/15/2003) jE = jnh4cl b_mtem(1,jA,jE) = -0.2073 b_mtem(2,jA,jE) = -0.4322 b_mtem(3,jA,jE) = 6.1271 b_mtem(4,jA,jE) = -12.3146 b_mtem(5,jA,jE) = 8.9919 b_mtem(6,jA,jE) = -2.3388 ! in NaCl jE = jnacl b_mtem(1,jA,jE) = 2.95913 b_mtem(2,jA,jE) = -7.92254 b_mtem(3,jA,jE) = 13.736 b_mtem(4,jA,jE) = -15.433 b_mtem(5,jA,jE) = 7.40386 b_mtem(6,jA,jE) = -0.918641 ! in NaNO3 jE = jnano3 b_mtem(1,jA,jE) = 0.893272 b_mtem(2,jA,jE) = 6.53768 b_mtem(3,jA,jE) = -32.3458 b_mtem(4,jA,jE) = 61.2834 b_mtem(5,jA,jE) = -56.4446 b_mtem(6,jA,jE) = 19.9202 ! in Na2SO4 jE = jna2so4 b_mtem(1,jA,jE) = 3.14484 b_mtem(2,jA,jE) = 0.077019 b_mtem(3,jA,jE) = -31.4199 b_mtem(4,jA,jE) = 80.5865 b_mtem(5,jA,jE) = -85.392 b_mtem(6,jA,jE) = 32.6644 ! in Ca(NO3)2 jE = jcano3 b_mtem(1,jA,jE) = 2.60432 b_mtem(2,jA,jE) = -0.55909 b_mtem(3,jA,jE) = -19.6671 b_mtem(4,jA,jE) = 53.3446 b_mtem(5,jA,jE) = -58.9076 b_mtem(6,jA,jE) = 22.9927 ! in CaCl2 (KM) revised on 3/13/2003 and again on 11/27/2003 jE = jcacl2 b_mtem(1,jA,jE) = 2.98036 b_mtem(2,jA,jE) = -8.55365 b_mtem(3,jA,jE) = 15.2108 b_mtem(4,jA,jE) = -15.9359 b_mtem(5,jA,jE) = 7.41772 b_mtem(6,jA,jE) = -1.32143 ! in HNO3 (added on 12/06/2004) jE = jhno3 b_mtem(1,jA,jE) = 3.8533 b_mtem(2,jA,jE) = -16.9427 b_mtem(3,jA,jE) = 45.0056 b_mtem(4,jA,jE) = -69.6145 b_mtem(5,jA,jE) = 54.1491 b_mtem(6,jA,jE) = -16.6513 ! in HCl (added on 12/06/2004) jE = jhcl b_mtem(1,jA,jE) = 2.56665 b_mtem(2,jA,jE) = -7.13585 b_mtem(3,jA,jE) = 14.8103 b_mtem(4,jA,jE) = -21.8881 b_mtem(5,jA,jE) = 16.6808 b_mtem(6,jA,jE) = -5.22091 ! in H2SO4 jE = jh2so4 b_mtem(1,jA,jE) = 2.50179 b_mtem(2,jA,jE) = -6.69364 b_mtem(3,jA,jE) = 11.6551 b_mtem(4,jA,jE) = -13.6897 b_mtem(5,jA,jE) = 7.36796 b_mtem(6,jA,jE) = -1.33245 ! in NH4HSO4 jE = jnh4hso4 b_mtem(1,jA,jE) = 0.149955 b_mtem(2,jA,jE) = 11.8213 b_mtem(3,jA,jE) = -53.9164 b_mtem(4,jA,jE) = 101.574 b_mtem(5,jA,jE) = -91.4123 b_mtem(6,jA,jE) = 31.5487 ! in (NH4)3H(SO4)2 jE = jlvcite b_mtem(1,jA,jE) = -2.36927 b_mtem(2,jA,jE) = 14.8359 b_mtem(3,jA,jE) = -44.3443 b_mtem(4,jA,jE) = 73.6229 b_mtem(5,jA,jE) = -65.3366 b_mtem(6,jA,jE) = 23.3250 ! in NaHSO4 jE = jnahso4 b_mtem(1,jA,jE) = 2.72993 b_mtem(2,jA,jE) = -0.23406 b_mtem(3,jA,jE) = -10.4103 b_mtem(4,jA,jE) = 13.1586 b_mtem(5,jA,jE) = -7.79925 b_mtem(6,jA,jE) = 2.30843 ! in Na3H(SO4)2 jE = jna3hso4 b_mtem(1,jA,jE) = 3.51258 b_mtem(2,jA,jE) = -3.95107 b_mtem(3,jA,jE) = -11.0175 b_mtem(4,jA,jE) = 38.8617 b_mtem(5,jA,jE) = -48.1575 b_mtem(6,jA,jE) = 20.4717 !---------- ! 2H.SO4 in E jA = jh2so4 ! in H2SO4 jE = jh2so4 b_mtem(1,jA,jE) = 0.76734 b_mtem(2,jA,jE) = -1.12263 b_mtem(3,jA,jE) = -9.08728 b_mtem(4,jA,jE) = 30.3836 b_mtem(5,jA,jE) = -38.4133 b_mtem(6,jA,jE) = 17.0106 ! in NH4HSO4 jE = jnh4hso4 b_mtem(1,jA,jE) = -2.03879 b_mtem(2,jA,jE) = 15.7033 b_mtem(3,jA,jE) = -58.7363 b_mtem(4,jA,jE) = 109.242 b_mtem(5,jA,jE) = -102.237 b_mtem(6,jA,jE) = 37.5350 ! in (NH4)3H(SO4)2 jE = jlvcite b_mtem(1,jA,jE) = -3.10228 b_mtem(2,jA,jE) = 16.6920 b_mtem(3,jA,jE) = -59.1522 b_mtem(4,jA,jE) = 113.487 b_mtem(5,jA,jE) = -110.890 b_mtem(6,jA,jE) = 42.4578 ! in (NH4)2SO4 jE = jnh4so4 b_mtem(1,jA,jE) = -3.43885 b_mtem(2,jA,jE) = 21.0372 b_mtem(3,jA,jE) = -84.7026 b_mtem(4,jA,jE) = 165.324 b_mtem(5,jA,jE) = -156.101 b_mtem(6,jA,jE) = 57.3101 ! in NaHSO4 jE = jnahso4 b_mtem(1,jA,jE) = 0.33164 b_mtem(2,jA,jE) = 6.55864 b_mtem(3,jA,jE) = -33.5876 b_mtem(4,jA,jE) = 65.1798 b_mtem(5,jA,jE) = -63.2046 b_mtem(6,jA,jE) = 24.1783 ! in Na3H(SO4)2 jE = jna3hso4 b_mtem(1,jA,jE) = 3.06830 b_mtem(2,jA,jE) = -3.18408 b_mtem(3,jA,jE) = -19.6332 b_mtem(4,jA,jE) = 61.3657 b_mtem(5,jA,jE) = -73.4438 b_mtem(6,jA,jE) = 31.2334 ! in Na2SO4 jE = jna2so4 b_mtem(1,jA,jE) = 2.58649 b_mtem(2,jA,jE) = 0.87921 b_mtem(3,jA,jE) = -39.3023 b_mtem(4,jA,jE) = 101.603 b_mtem(5,jA,jE) = -109.469 b_mtem(6,jA,jE) = 43.0188 ! in HNO3 jE = jhno3 b_mtem(1,jA,jE) = 1.54587 b_mtem(2,jA,jE) = -7.50976 b_mtem(3,jA,jE) = 12.8237 b_mtem(4,jA,jE) = -10.1452 b_mtem(5,jA,jE) = -0.541956 b_mtem(6,jA,jE) = 3.34536 ! in HCl jE = jhcl b_mtem(1,jA,jE) = 0.829757 b_mtem(2,jA,jE) = -4.11316 b_mtem(3,jA,jE) = 3.67111 b_mtem(4,jA,jE) = 3.6833 b_mtem(5,jA,jE) = -11.2711 b_mtem(6,jA,jE) = 6.71421 !---------- ! H.HSO4 in E jA = jhhso4 ! in H2SO4 jE = jh2so4 b_mtem(1,jA,jE) = 2.63953 b_mtem(2,jA,jE) = -6.01532 b_mtem(3,jA,jE) = 10.0204 b_mtem(4,jA,jE) = -12.4840 b_mtem(5,jA,jE) = 7.78853 b_mtem(6,jA,jE) = -2.12638 ! in NH4HSO4 jE = jnh4hso4 b_mtem(1,jA,jE) = -0.77412 b_mtem(2,jA,jE) = 14.1656 b_mtem(3,jA,jE) = -53.4087 b_mtem(4,jA,jE) = 93.2013 b_mtem(5,jA,jE) = -80.5723 b_mtem(6,jA,jE) = 27.1577 ! in (NH4)3H(SO4)2 jE = jlvcite b_mtem(1,jA,jE) = -2.98882 b_mtem(2,jA,jE) = 14.4436 b_mtem(3,jA,jE) = -40.1774 b_mtem(4,jA,jE) = 67.5937 b_mtem(5,jA,jE) = -61.5040 b_mtem(6,jA,jE) = 22.3695 ! in (NH4)2SO4 jE = jnh4so4 b_mtem(1,jA,jE) = -1.15502 b_mtem(2,jA,jE) = 8.12309 b_mtem(3,jA,jE) = -38.4726 b_mtem(4,jA,jE) = 80.8861 b_mtem(5,jA,jE) = -80.1644 b_mtem(6,jA,jE) = 30.4717 ! in NaHSO4 jE = jnahso4 b_mtem(1,jA,jE) = 1.99641 b_mtem(2,jA,jE) = -2.96061 b_mtem(3,jA,jE) = 5.54778 b_mtem(4,jA,jE) = -14.5488 b_mtem(5,jA,jE) = 14.8492 b_mtem(6,jA,jE) = -5.1389 ! in Na3H(SO4)2 jE = jna3hso4 b_mtem(1,jA,jE) = 2.23816 b_mtem(2,jA,jE) = -3.20847 b_mtem(3,jA,jE) = -4.82853 b_mtem(4,jA,jE) = 20.9192 b_mtem(5,jA,jE) = -27.2819 b_mtem(6,jA,jE) = 11.8655 ! in Na2SO4 jE = jna2so4 b_mtem(1,jA,jE) = 2.56907 b_mtem(2,jA,jE) = 1.13444 b_mtem(3,jA,jE) = -34.6853 b_mtem(4,jA,jE) = 87.9775 b_mtem(5,jA,jE) = -93.2330 b_mtem(6,jA,jE) = 35.9260 ! in HNO3 jE = jhno3 b_mtem(1,jA,jE) = 2.00024 b_mtem(2,jA,jE) = -4.80868 b_mtem(3,jA,jE) = 8.29222 b_mtem(4,jA,jE) = -11.0849 b_mtem(5,jA,jE) = 7.51262 b_mtem(6,jA,jE) = -2.07654 ! in HCl jE = jhcl b_mtem(1,jA,jE) = 2.8009 b_mtem(2,jA,jE) = -6.98416 b_mtem(3,jA,jE) = 14.3146 b_mtem(4,jA,jE) = -22.0068 b_mtem(5,jA,jE) = 17.5557 b_mtem(6,jA,jE) = -5.84917 !---------- ! NH4HSO4 in E jA = jnh4hso4 ! in H2SO4 jE = jh2so4 b_mtem(1,jA,jE) = 0.169160 b_mtem(2,jA,jE) = 2.15094 b_mtem(3,jA,jE) = -9.62904 b_mtem(4,jA,jE) = 18.2631 b_mtem(5,jA,jE) = -17.3333 b_mtem(6,jA,jE) = 6.19835 ! in NH4HSO4 jE = jnh4hso4 b_mtem(1,jA,jE) = -2.34457 b_mtem(2,jA,jE) = 12.8035 b_mtem(3,jA,jE) = -35.2513 b_mtem(4,jA,jE) = 53.6153 b_mtem(5,jA,jE) = -42.7655 b_mtem(6,jA,jE) = 13.7129 ! in (NH4)3H(SO4)2 jE = jlvcite b_mtem(1,jA,jE) = -2.56109 b_mtem(2,jA,jE) = 11.1414 b_mtem(3,jA,jE) = -30.2361 b_mtem(4,jA,jE) = 50.0320 b_mtem(5,jA,jE) = -44.1586 b_mtem(6,jA,jE) = 15.5393 ! in (NH4)2SO4 jE = jnh4so4 b_mtem(1,jA,jE) = -0.97315 b_mtem(2,jA,jE) = 7.06295 b_mtem(3,jA,jE) = -29.3032 b_mtem(4,jA,jE) = 57.6101 b_mtem(5,jA,jE) = -54.9020 b_mtem(6,jA,jE) = 20.2222 ! in NaHSO4 jE = jnahso4 b_mtem(1,jA,jE) = -0.44450 b_mtem(2,jA,jE) = 3.33451 b_mtem(3,jA,jE) = -15.2791 b_mtem(4,jA,jE) = 30.1413 b_mtem(5,jA,jE) = -26.7710 b_mtem(6,jA,jE) = 8.78462 ! in Na3H(SO4)2 jE = jna3hso4 b_mtem(1,jA,jE) = -0.99780 b_mtem(2,jA,jE) = 4.69200 b_mtem(3,jA,jE) = -16.1219 b_mtem(4,jA,jE) = 29.3100 b_mtem(5,jA,jE) = -26.3383 b_mtem(6,jA,jE) = 9.20695 ! in Na2SO4 jE = jna2so4 b_mtem(1,jA,jE) = -0.52694 b_mtem(2,jA,jE) = 7.02684 b_mtem(3,jA,jE) = -33.7508 b_mtem(4,jA,jE) = 70.0565 b_mtem(5,jA,jE) = -68.3226 b_mtem(6,jA,jE) = 25.2692 ! in HNO3 jE = jhno3 b_mtem(1,jA,jE) = 0.572926 b_mtem(2,jA,jE) = -2.04791 b_mtem(3,jA,jE) = 2.1134 b_mtem(4,jA,jE) = 0.246654 b_mtem(5,jA,jE) = -3.06019 b_mtem(6,jA,jE) = 1.98126 ! in HCl jE = jhcl b_mtem(1,jA,jE) = 0.56514 b_mtem(2,jA,jE) = 0.22287 b_mtem(3,jA,jE) = -2.76973 b_mtem(4,jA,jE) = 4.54444 b_mtem(5,jA,jE) = -3.86549 b_mtem(6,jA,jE) = 1.13441 !---------- ! (NH4)3H(SO4)2 in E jA = jlvcite ! in H2SO4 jE = jh2so4 b_mtem(1,jA,jE) = -1.44811 b_mtem(2,jA,jE) = 6.71815 b_mtem(3,jA,jE) = -25.0141 b_mtem(4,jA,jE) = 50.1109 b_mtem(5,jA,jE) = -50.0561 b_mtem(6,jA,jE) = 19.3370 ! in NH4HSO4 jE = jnh4hso4 b_mtem(1,jA,jE) = -3.41707 b_mtem(2,jA,jE) = 13.4496 b_mtem(3,jA,jE) = -34.8018 b_mtem(4,jA,jE) = 55.2987 b_mtem(5,jA,jE) = -48.1839 b_mtem(6,jA,jE) = 17.2444 ! in (NH4)3H(SO4)2 jE = jlvcite b_mtem(1,jA,jE) = -2.54479 b_mtem(2,jA,jE) = 11.8501 b_mtem(3,jA,jE) = -39.7286 b_mtem(4,jA,jE) = 74.2479 b_mtem(5,jA,jE) = -70.4934 b_mtem(6,jA,jE) = 26.2836 ! in (NH4)2SO4 jE = jnh4so4 b_mtem(1,jA,jE) = -2.30561 b_mtem(2,jA,jE) = 14.5806 b_mtem(3,jA,jE) = -55.1238 b_mtem(4,jA,jE) = 103.451 b_mtem(5,jA,jE) = -95.2571 b_mtem(6,jA,jE) = 34.2218 ! in NaHSO4 jE = jnahso4 b_mtem(1,jA,jE) = -2.20809 b_mtem(2,jA,jE) = 13.6391 b_mtem(3,jA,jE) = -57.8246 b_mtem(4,jA,jE) = 117.907 b_mtem(5,jA,jE) = -112.154 b_mtem(6,jA,jE) = 40.3058 ! in Na3H(SO4)2 jE = jna3hso4 b_mtem(1,jA,jE) = -1.15099 b_mtem(2,jA,jE) = 6.32269 b_mtem(3,jA,jE) = -27.3860 b_mtem(4,jA,jE) = 55.4592 b_mtem(5,jA,jE) = -54.0100 b_mtem(6,jA,jE) = 20.3469 ! in Na2SO4 jE = jna2so4 b_mtem(1,jA,jE) = -1.15678 b_mtem(2,jA,jE) = 8.28718 b_mtem(3,jA,jE) = -37.3231 b_mtem(4,jA,jE) = 76.6124 b_mtem(5,jA,jE) = -74.9307 b_mtem(6,jA,jE) = 28.0559 ! in HNO3 jE = jhno3 b_mtem(1,jA,jE) = 0.01502 b_mtem(2,jA,jE) = -3.1197 b_mtem(3,jA,jE) = 3.61104 b_mtem(4,jA,jE) = 3.05196 b_mtem(5,jA,jE) = -9.98957 b_mtem(6,jA,jE) = 6.04155 ! in HCl jE = jhcl b_mtem(1,jA,jE) = -1.06477 b_mtem(2,jA,jE) = 3.38801 b_mtem(3,jA,jE) = -12.5784 b_mtem(4,jA,jE) = 25.2823 b_mtem(5,jA,jE) = -25.4611 b_mtem(6,jA,jE) = 10.0754 !---------- ! NaHSO4 in E jA = jnahso4 ! in H2SO4 jE = jh2so4 b_mtem(1,jA,jE) = 0.68259 b_mtem(2,jA,jE) = 0.71468 b_mtem(3,jA,jE) = -5.59003 b_mtem(4,jA,jE) = 11.0089 b_mtem(5,jA,jE) = -10.7983 b_mtem(6,jA,jE) = 3.82335 ! in NH4HSO4 jE = jnh4hso4 b_mtem(1,jA,jE) = -0.03956 b_mtem(2,jA,jE) = 4.52828 b_mtem(3,jA,jE) = -25.2557 b_mtem(4,jA,jE) = 54.4225 b_mtem(5,jA,jE) = -52.5105 b_mtem(6,jA,jE) = 18.6562 ! in (NH4)3H(SO4)2 jE = jlvcite b_mtem(1,jA,jE) = -1.53503 b_mtem(2,jA,jE) = 8.27608 b_mtem(3,jA,jE) = -28.9539 b_mtem(4,jA,jE) = 55.2876 b_mtem(5,jA,jE) = -51.9563 b_mtem(6,jA,jE) = 18.6576 ! in (NH4)2SO4 jE = jnh4so4 b_mtem(1,jA,jE) = -0.38793 b_mtem(2,jA,jE) = 7.14680 b_mtem(3,jA,jE) = -38.7201 b_mtem(4,jA,jE) = 84.3965 b_mtem(5,jA,jE) = -84.7453 b_mtem(6,jA,jE) = 32.1283 ! in NaHSO4 jE = jnahso4 b_mtem(1,jA,jE) = -0.41982 b_mtem(2,jA,jE) = 4.26491 b_mtem(3,jA,jE) = -20.2351 b_mtem(4,jA,jE) = 42.6764 b_mtem(5,jA,jE) = -40.7503 b_mtem(6,jA,jE) = 14.2868 ! in Na3H(SO4)2 jE = jna3hso4 b_mtem(1,jA,jE) = -0.32912 b_mtem(2,jA,jE) = 1.80808 b_mtem(3,jA,jE) = -8.01286 b_mtem(4,jA,jE) = 15.5791 b_mtem(5,jA,jE) = -14.5494 b_mtem(6,jA,jE) = 5.27052 ! in Na2SO4 jE = jna2so4 b_mtem(1,jA,jE) = 0.10271 b_mtem(2,jA,jE) = 5.09559 b_mtem(3,jA,jE) = -30.3295 b_mtem(4,jA,jE) = 66.2975 b_mtem(5,jA,jE) = -66.3458 b_mtem(6,jA,jE) = 24.9443 ! in HNO3 jE = jhno3 b_mtem(1,jA,jE) = 0.608309 b_mtem(2,jA,jE) = -0.541905 b_mtem(3,jA,jE) = -2.52084 b_mtem(4,jA,jE) = 6.63297 b_mtem(5,jA,jE) = -7.24599 b_mtem(6,jA,jE) = 2.88811 ! in HCl jE = jhcl b_mtem(1,jA,jE) = 1.98399 b_mtem(2,jA,jE) = -4.51562 b_mtem(3,jA,jE) = 8.36059 b_mtem(4,jA,jE) = -12.4948 b_mtem(5,jA,jE) = 9.67514 b_mtem(6,jA,jE) = -3.18004 !---------- ! Na3H(SO4)2 in E jA = jna3hso4 ! in H2SO4 jE = jh2so4 b_mtem(1,jA,jE) = -0.83214 b_mtem(2,jA,jE) = 4.99572 b_mtem(3,jA,jE) = -20.1697 b_mtem(4,jA,jE) = 41.4066 b_mtem(5,jA,jE) = -42.2119 b_mtem(6,jA,jE) = 16.4855 ! in NH4HSO4 jE = jnh4hso4 b_mtem(1,jA,jE) = -0.65139 b_mtem(2,jA,jE) = 3.52300 b_mtem(3,jA,jE) = -22.8220 b_mtem(4,jA,jE) = 56.2956 b_mtem(5,jA,jE) = -59.9028 b_mtem(6,jA,jE) = 23.1844 ! in (NH4)3H(SO4)2 jE = jlvcite b_mtem(1,jA,jE) = -1.31331 b_mtem(2,jA,jE) = 8.40835 b_mtem(3,jA,jE) = -38.1757 b_mtem(4,jA,jE) = 80.5312 b_mtem(5,jA,jE) = -79.8346 b_mtem(6,jA,jE) = 30.0219 ! in (NH4)2SO4 jE = jnh4so4 b_mtem(1,jA,jE) = -1.03054 b_mtem(2,jA,jE) = 8.08155 b_mtem(3,jA,jE) = -38.1046 b_mtem(4,jA,jE) = 78.7168 b_mtem(5,jA,jE) = -77.2263 b_mtem(6,jA,jE) = 29.1521 ! in NaHSO4 jE = jnahso4 b_mtem(1,jA,jE) = -1.90695 b_mtem(2,jA,jE) = 11.6241 b_mtem(3,jA,jE) = -50.3175 b_mtem(4,jA,jE) = 105.884 b_mtem(5,jA,jE) = -103.258 b_mtem(6,jA,jE) = 37.6588 ! in Na3H(SO4)2 jE = jna3hso4 b_mtem(1,jA,jE) = -0.34780 b_mtem(2,jA,jE) = 2.85363 b_mtem(3,jA,jE) = -17.6224 b_mtem(4,jA,jE) = 38.9220 b_mtem(5,jA,jE) = -39.8106 b_mtem(6,jA,jE) = 15.6055 ! in Na2SO4 jE = jna2so4 b_mtem(1,jA,jE) = -0.75230 b_mtem(2,jA,jE) = 10.0140 b_mtem(3,jA,jE) = -50.5677 b_mtem(4,jA,jE) = 106.941 b_mtem(5,jA,jE) = -105.534 b_mtem(6,jA,jE) = 39.5196 ! in HNO3 jE = jhno3 b_mtem(1,jA,jE) = 0.057456 b_mtem(2,jA,jE) = -1.31264 b_mtem(3,jA,jE) = -1.94662 b_mtem(4,jA,jE) = 10.7024 b_mtem(5,jA,jE) = -14.9946 b_mtem(6,jA,jE) = 7.12161 ! in HCl jE = jhcl b_mtem(1,jA,jE) = 0.637894 b_mtem(2,jA,jE) = -2.29719 b_mtem(3,jA,jE) = 0.765361 b_mtem(4,jA,jE) = 4.8748 b_mtem(5,jA,jE) = -9.25978 b_mtem(6,jA,jE) = 4.91773 ! ! ! !---------------------------------------------------------- ! Coefficients for %MDRH(T) = d1 + d2*T + d3*T^2 + d4*T^3 (T in Kelvin) ! valid Temperature Range: 240 - 320 K !---------------------------------------------------------- ! ! SULFATE-POOR SYSTEMS ! AC j_index = 1 d_mdrh(j_index,1) = -58.00268351 d_mdrh(j_index,2) = 2.031077573 d_mdrh(j_index,3) = -0.008281218 d_mdrh(j_index,4) = 1.00447E-05 ! AN j_index = 2 d_mdrh(j_index,1) = 1039.137773 d_mdrh(j_index,2) = -11.47847095 d_mdrh(j_index,3) = 0.047702786 d_mdrh(j_index,4) = -6.77675E-05 ! AS j_index = 3 d_mdrh(j_index,1) = 115.8366357 d_mdrh(j_index,2) = 0.491881663 d_mdrh(j_index,3) = -0.00422807 d_mdrh(j_index,4) = 7.29274E-06 ! SC j_index = 4 d_mdrh(j_index,1) = 253.2424151 d_mdrh(j_index,2) = -1.429957864 d_mdrh(j_index,3) = 0.003727554 d_mdrh(j_index,4) = -3.13037E-06 ! SN j_index = 5 d_mdrh(j_index,1) = -372.4306506 d_mdrh(j_index,2) = 5.3955633 d_mdrh(j_index,3) = -0.019804438 d_mdrh(j_index,4) = 2.25662E-05 ! SS j_index = 6 d_mdrh(j_index,1) = 286.1271416 d_mdrh(j_index,2) = -1.670787758 d_mdrh(j_index,3) = 0.004431373 d_mdrh(j_index,4) = -3.57757E-06 ! CC j_index = 7 d_mdrh(j_index,1) = -1124.07059 d_mdrh(j_index,2) = 14.26364209 d_mdrh(j_index,3) = -0.054816822 d_mdrh(j_index,4) = 6.70107E-05 ! CN j_index = 8 d_mdrh(j_index,1) = 1855.413934 d_mdrh(j_index,2) = -20.29219473 d_mdrh(j_index,3) = 0.07807482 d_mdrh(j_index,4) = -1.017887858e-4 ! AN + AC j_index = 9 d_mdrh(j_index,1) = 1761.176886 d_mdrh(j_index,2) = -19.29811062 d_mdrh(j_index,3) = 0.075676987 d_mdrh(j_index,4) = -1.0116959e-4 ! AS + AC j_index = 10 d_mdrh(j_index,1) = 122.1074303 d_mdrh(j_index,2) = 0.429692122 d_mdrh(j_index,3) = -0.003928277 d_mdrh(j_index,4) = 6.43275E-06 ! AS + AN j_index = 11 d_mdrh(j_index,1) = 2424.634678 d_mdrh(j_index,2) = -26.54031307 d_mdrh(j_index,3) = 0.101625387 d_mdrh(j_index,4) = -1.31544547798e-4 ! AS + AN + AC j_index = 12 d_mdrh(j_index,1) = 2912.082599 d_mdrh(j_index,2) = -31.8894185 d_mdrh(j_index,3) = 0.121185849 d_mdrh(j_index,4) = -1.556534623e-4 ! SC + AC j_index = 13 d_mdrh(j_index,1) = 172.2596493 d_mdrh(j_index,2) = -0.511006195 d_mdrh(j_index,3) = 4.27244597e-4 d_mdrh(j_index,4) = 4.12797E-07 ! SN + AC j_index = 14 d_mdrh(j_index,1) = 1596.184935 d_mdrh(j_index,2) = -16.37945565 d_mdrh(j_index,3) = 0.060281218 d_mdrh(j_index,4) = -7.6161E-05 ! SN + AN j_index = 15 d_mdrh(j_index,1) = 1916.072988 d_mdrh(j_index,2) = -20.85594868 d_mdrh(j_index,3) = 0.081140141 d_mdrh(j_index,4) = -1.07954274796e-4 ! SN + AN + AC j_index = 16 d_mdrh(j_index,1) = 1467.165935 d_mdrh(j_index,2) = -16.01166196 d_mdrh(j_index,3) = 0.063505582 d_mdrh(j_index,4) = -8.66722E-05 ! SN + SC j_index = 17 d_mdrh(j_index,1) = 158.447059 d_mdrh(j_index,2) = -0.628167358 d_mdrh(j_index,3) = 0.002014448 d_mdrh(j_index,4) = -3.13037E-06 ! SN + SC + AC j_index = 18 d_mdrh(j_index,1) = 1115.892468 d_mdrh(j_index,2) = -11.76936534 d_mdrh(j_index,3) = 0.045577399 d_mdrh(j_index,4) = -6.05779E-05 ! SS + AC j_index = 19 d_mdrh(j_index,1) = 269.5432407 d_mdrh(j_index,2) = -1.319963885 d_mdrh(j_index,3) = 0.002592363 d_mdrh(j_index,4) = -1.44479E-06 ! SS + AN j_index = 20 d_mdrh(j_index,1) = 2841.334784 d_mdrh(j_index,2) = -31.1889487 d_mdrh(j_index,3) = 0.118809274 d_mdrh(j_index,4) = -1.53007e-4 ! SS + AN + AC j_index = 21 d_mdrh(j_index,1) = 2199.36914 d_mdrh(j_index,2) = -24.11926569 d_mdrh(j_index,3) = 0.092932361 d_mdrh(j_index,4) = -1.21774e-4 ! SS + AS j_index = 22 d_mdrh(j_index,1) = 395.0051604 d_mdrh(j_index,2) = -2.521101657 d_mdrh(j_index,3) = 0.006139319 d_mdrh(j_index,4) = -4.43756E-06 ! SS + AS + AC j_index = 23 d_mdrh(j_index,1) = 386.5150675 d_mdrh(j_index,2) = -2.4632138 d_mdrh(j_index,3) = 0.006139319 d_mdrh(j_index,4) = -4.98796E-06 ! SS + AS + AN j_index = 24 d_mdrh(j_index,1) = 3101.538491 d_mdrh(j_index,2) = -34.19978105 d_mdrh(j_index,3) = 0.130118605 d_mdrh(j_index,4) = -1.66873e-4 ! SS + AS + AN + AC j_index = 25 d_mdrh(j_index,1) = 2307.579403 d_mdrh(j_index,2) = -25.43136774 d_mdrh(j_index,3) = 0.098064728 d_mdrh(j_index,4) = -1.28301e-4 ! SS + SC j_index = 26 d_mdrh(j_index,1) = 291.8309602 d_mdrh(j_index,2) = -1.828912974 d_mdrh(j_index,3) = 0.005053148 d_mdrh(j_index,4) = -4.57516E-06 ! SS + SC + AC j_index = 27 d_mdrh(j_index,1) = 188.3914345 d_mdrh(j_index,2) = -0.631345031 d_mdrh(j_index,3) = 0.000622807 d_mdrh(j_index,4) = 4.47196E-07 ! SS + SN j_index = 28 d_mdrh(j_index,1) = -167.1252839 d_mdrh(j_index,2) = 2.969828002 d_mdrh(j_index,3) = -0.010637255 d_mdrh(j_index,4) = 1.13175E-05 ! SS + SN + AC j_index = 29 d_mdrh(j_index,1) = 1516.782768 d_mdrh(j_index,2) = -15.7922661 d_mdrh(j_index,3) = 0.058942209 d_mdrh(j_index,4) = -7.5301E-05 ! SS + SN + AN j_index = 30 d_mdrh(j_index,1) = 1739.963163 d_mdrh(j_index,2) = -19.06576022 d_mdrh(j_index,3) = 0.07454963 d_mdrh(j_index,4) = -9.94302E-05 ! SS + SN + AN + AC j_index = 31 d_mdrh(j_index,1) = 2152.104877 d_mdrh(j_index,2) = -23.74998008 d_mdrh(j_index,3) = 0.092256654 d_mdrh(j_index,4) = -1.21953e-4 ! SS + SN + SC j_index = 32 d_mdrh(j_index,1) = 221.9976265 d_mdrh(j_index,2) = -1.311331272 d_mdrh(j_index,3) = 0.004406089 d_mdrh(j_index,4) = -5.88235E-06 ! SS + SN + SC + AC j_index = 33 d_mdrh(j_index,1) = 1205.645615 d_mdrh(j_index,2) = -12.71353459 d_mdrh(j_index,3) = 0.048803922 d_mdrh(j_index,4) = -6.41899E-05 ! CC + AC j_index = 34 d_mdrh(j_index,1) = 506.6737879 d_mdrh(j_index,2) = -3.723520818 d_mdrh(j_index,3) = 0.010814242 d_mdrh(j_index,4) = -1.21087E-05 ! CC + SC j_index = 35 d_mdrh(j_index,1) = -1123.523841 d_mdrh(j_index,2) = 14.08345977 d_mdrh(j_index,3) = -0.053687823 d_mdrh(j_index,4) = 6.52219E-05 ! CC + SC + AC j_index = 36 d_mdrh(j_index,1) = -1159.98607 d_mdrh(j_index,2) = 14.44309169 d_mdrh(j_index,3) = -0.054841073 d_mdrh(j_index,4) = 6.64259E-05 ! CN + AC j_index = 37 d_mdrh(j_index,1) = 756.0747916 d_mdrh(j_index,2) = -8.546826257 d_mdrh(j_index,3) = 0.035798677 d_mdrh(j_index,4) = -5.06629E-05 ! CN + AN j_index = 38 d_mdrh(j_index,1) = 338.668191 d_mdrh(j_index,2) = -2.971223403 d_mdrh(j_index,3) = 0.012294866 d_mdrh(j_index,4) = -1.87558E-05 ! CN + AN + AC j_index = 39 d_mdrh(j_index,1) = -53.18033508 d_mdrh(j_index,2) = 0.663911748 d_mdrh(j_index,3) = 9.16326e-4 d_mdrh(j_index,4) = -6.70354E-06 ! CN + SC j_index = 40 d_mdrh(j_index,1) = 3623.831129 d_mdrh(j_index,2) = -39.27226457 d_mdrh(j_index,3) = 0.144559515 d_mdrh(j_index,4) = -1.78159e-4 ! CN + SC + AC j_index = 41 d_mdrh(j_index,1) = 3436.656743 d_mdrh(j_index,2) = -37.16192684 d_mdrh(j_index,3) = 0.136641377 d_mdrh(j_index,4) = -1.68262e-4 ! CN + SN j_index = 42 d_mdrh(j_index,1) = 768.608476 d_mdrh(j_index,2) = -8.051517149 d_mdrh(j_index,3) = 0.032342332 d_mdrh(j_index,4) = -4.52224E-05 ! CN + SN + AC j_index = 43 d_mdrh(j_index,1) = 33.58027951 d_mdrh(j_index,2) = -0.308772182 d_mdrh(j_index,3) = 0.004713639 d_mdrh(j_index,4) = -1.19658E-05 ! CN + SN + AN j_index = 44 d_mdrh(j_index,1) = 57.80183041 d_mdrh(j_index,2) = 0.215264604 d_mdrh(j_index,3) = 4.11406e-4 d_mdrh(j_index,4) = -4.30702E-06 ! CN + SN + AN + AC j_index = 45 d_mdrh(j_index,1) = -234.368984 d_mdrh(j_index,2) = 2.721045204 d_mdrh(j_index,3) = -0.006688341 d_mdrh(j_index,4) = 2.31729E-06 ! CN + SN + SC j_index = 46 d_mdrh(j_index,1) = 3879.080557 d_mdrh(j_index,2) = -42.13562874 d_mdrh(j_index,3) = 0.155235005 d_mdrh(j_index,4) = -1.91387e-4 ! CN + SN + SC + AC j_index = 47 d_mdrh(j_index,1) = 3600.576985 d_mdrh(j_index,2) = -39.0283489 d_mdrh(j_index,3) = 0.143710316 d_mdrh(j_index,4) = -1.77167e-4 ! CN + CC j_index = 48 d_mdrh(j_index,1) = -1009.729826 d_mdrh(j_index,2) = 12.9145339 d_mdrh(j_index,3) = -0.049811146 d_mdrh(j_index,4) = 6.09563E-05 ! CN + CC + AC j_index = 49 d_mdrh(j_index,1) = -577.0919514 d_mdrh(j_index,2) = 8.020324227 d_mdrh(j_index,3) = -0.031469556 d_mdrh(j_index,4) = 3.82181E-05 ! CN + CC + SC j_index = 50 d_mdrh(j_index,1) = -728.9983499 d_mdrh(j_index,2) = 9.849458215 d_mdrh(j_index,3) = -0.03879257 d_mdrh(j_index,4) = 4.78844E-05 ! CN + CC + SC + AC j_index = 51 d_mdrh(j_index,1) = -803.7026845 d_mdrh(j_index,2) = 10.61881494 d_mdrh(j_index,3) = -0.041402993 d_mdrh(j_index,4) = 5.08084E-05 ! ! SULFATE-RICH SYSTEMS ! AB j_index = 52 d_mdrh(j_index,1) = -493.6190458 d_mdrh(j_index,2) = 6.747053851 d_mdrh(j_index,3) = -0.026955267 d_mdrh(j_index,4) = 3.45118E-05 ! LV j_index = 53 d_mdrh(j_index,1) = 53.37874093 d_mdrh(j_index,2) = 1.01368249 d_mdrh(j_index,3) = -0.005887513 d_mdrh(j_index,4) = 8.94393E-06 ! SB j_index = 54 d_mdrh(j_index,1) = 206.619047 d_mdrh(j_index,2) = -1.342735684 d_mdrh(j_index,3) = 0.003197691 d_mdrh(j_index,4) = -1.93603E-06 ! AB + LV j_index = 55 d_mdrh(j_index,1) = -493.6190458 d_mdrh(j_index,2) = 6.747053851 d_mdrh(j_index,3) = -0.026955267 d_mdrh(j_index,4) = 3.45118E-05 ! AS + LV j_index = 56 d_mdrh(j_index,1) = 53.37874093 d_mdrh(j_index,2) = 1.01368249 d_mdrh(j_index,3) = -0.005887513 d_mdrh(j_index,4) = 8.94393E-06 ! SS + SB j_index = 57 d_mdrh(j_index,1) = 206.619047 d_mdrh(j_index,2) = -1.342735684 d_mdrh(j_index,3) = 0.003197691 d_mdrh(j_index,4) = -1.93603E-06 ! SS + LV j_index = 58 d_mdrh(j_index,1) = 41.7619047 d_mdrh(j_index,2) = 1.303872053 d_mdrh(j_index,3) = -0.007647908 d_mdrh(j_index,4) = 1.17845E-05 ! SS + AS + LV j_index = 59 d_mdrh(j_index,1) = 41.7619047 d_mdrh(j_index,2) = 1.303872053 d_mdrh(j_index,3) = -0.007647908 d_mdrh(j_index,4) = 1.17845E-05 ! SS + AB j_index = 60 d_mdrh(j_index,1) = -369.7142842 d_mdrh(j_index,2) = 5.512878771 d_mdrh(j_index,3) = -0.02301948 d_mdrh(j_index,4) = 3.0303E-05 ! SS + LV + AB j_index = 61 d_mdrh(j_index,1) = -369.7142842 d_mdrh(j_index,2) = 5.512878771 d_mdrh(j_index,3) = -0.02301948 d_mdrh(j_index,4) = 3.0303E-05 ! SB + AB j_index = 62 d_mdrh(j_index,1) = -162.8095232 d_mdrh(j_index,2) = 2.399326592 d_mdrh(j_index,3) = -0.009336219 d_mdrh(j_index,4) = 1.17845E-05 ! SS + SB + AB j_index = 63 d_mdrh(j_index,1) = -735.4285689 d_mdrh(j_index,2) = 8.885521857 d_mdrh(j_index,3) = -0.033488456 d_mdrh(j_index,4) = 4.12458E-05 ! endif ! first if (msoa_flag1 >= 1000) call soa_vbs_load_params( 2 ) return end subroutine load_mosaic_parameters end module module_mosaic_init_aerpar