Index | Rundeck: E6TomaF40 | Created: Sat May 18 03:30:40 EDT 2024 |
get_sulfate | Module: | File: TRACERS_AEROSOLS_Koch_e4.f |
Author : Dorothy Koch
Version:
Declaration:
subroutine get_sulfate(pl, temp_in, fcloud, wa_vol, wmxtr, sulfin, sulfinom, sulfinc, sulfout, tr_left, tmg, tmd, airm, lhx, dt_sulf, fcld0)
- bygasc : used from constant
- mair : used from constant
- teeny : used from constant
- mb2kg : used from constant
- gasc : used from constant
- lhe : used from constant
- trname : used from oldtracer_mod
- mass2vol : used from oldtracer_mod
- tr_mm : used from oldtracer_mod
- tr_rkd : used from oldtracer_mod
- tr_dhd : used from oldtracer_mod
- n_h2o2_s : used from tracer_com
- n_so2 : used from tracer_com
- ntm : used from tracer_com
- lm : used from tracer_com
- n_so4 : used from tracer_com
- n_h2o2 : used from tracer_com
- coupled_chem : used from tracer_com
- aqchem_count : used from tracer_com
- aqchem_list : used from tracer_com
- ntx : used from clouds
- dxypij : used from clouds
- dtsrc : used from model_com
- sulfin :
real*8, dimension(aqchem_count), intent(out)
- amount of precursor used to make product from the gas phase (kg)
- sulfinom :
real*8, dimension(aqchem_count), intent(out)
- sulfin=sulfinom*tmg, for avoiding divisions sulfinom=sulfin/tmg
- sulfinc :
real*8, dimension(aqchem_count), intent(out)
- amount of precursor used to make product from the condensate (kg)
- sulfout :
real*8, dimension(aqchem_count), intent(out)
- total amount of product generated (kg)
- tr_left :
real*8, dimension(aqchem_count), intent(out)
- is the amount of precursor left after product is made
- fcloud :
real*8, intent(in)
- cloud fraction available for tracer condensation. fcloud=fplume
- fcld0 :
real*8, intent(in)
- updated cloud fraction, given the current state of large-scale
- lhx :
real*8, intent(in)
- latent heat of evaporation or sublimation (J/Kg). When equal to lhe
- finc :
real*8
- XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX
- airm :
real*8, intent(in)
- layer pressure depth (mb). Multiply by mb2kg to convert to air mass
- amass :
real*8
- airmass in kg, calculated by airm*mb2kg*dxypij
- pl :
real*8, intent(in)
- pressure at current altitude (mbar)
- press :
real*8
- pressure at current altitude (Pa)
- temp :
real*8
- temperature to be used (K), always greater than 230K
- tfac :
real*8
- exponent factor for temperature dependence calculations (mol/J).
- clwc :
real*8
- cloud liquid water content (volume water/volume air)
- temp_in :
real*8, intent(in)
- temperature at current altitude (K)
- wmxtr :
real*8, intent(in)
- cloud water mixing ratio (kg water/kg air)
- k1so2dissoc0 :
real*8, parameter
- first dissociation rate of dissolved SO2 at 298K.
Initial Value= 1.3d-2
- dh1so2dissoc :
real*8, parameter
- enthalpy of dissociation for the first dissociation rate
Initial Value= -1.6736d4
- k1so2dissoc :
real*8
- first dissociation rate of dissolved SO2 at current temp.
- hplus :
real*8, parameter
- concentration of H+ (molar) for pH=4.5
Initial Value= 10.d0**(-4.5d0)
- henry :
real*8
- modified henry constant of the current tracer at current
- kso2h2o20 :
real*8, parameter
- reaction rate of SO2 + H2O2 at 298K. SO2 + H2O2 --> SO3 + H2O
Initial Value= 6.357d14
- dhso2h2o2 :
real*8, parameter
- enthalpy of reaction of SO2 + H2O2
Initial Value= 3.95d4
- kso2h2o2 :
real*8
- reaction rate of SO2 + H2O2 at current temp.
- ix :
integer
- index of current species
- is :
integer
- index of current species in ntx array
- ih :
integer
- index of current species in ntx array
- isx :
integer
- index of SO2 species in aqchem_list array
- ihx :
integer
- index of H2O2 species in aqchem_list array
- tmg :
real*8, dimension(ntx), intent(in)
- amount of tracer in the gas phase in the cloudy area (kg).
- tmd :
real*8, dimension(ntx), intent(in)
- amount of tracer in the aqueous phase (kg).
- tmgmol :
real*8, dimension(aqchem_count)
- amount of gas phase tracer in cloudy area (moles)
- tmdmol :
real*8, dimension(aqchem_count)
- amount of gas phase tracer in cloudy area (moles)
- tmgrate :
real*8, dimension(aqchem_count)
- is the new concentration of species that resulted from the
- tmdrate :
real*8, dimension(aqchem_count)
- is the new concentration of species that resulted from the
- wa_vol :
real*8, intent(in)
- cloud water volume (liters)
- dso4g :
real*8
- amount of sulfate produced from the gas phase (moles/kg/kg)
- dso4d :
real*8
- amount of sulfate produced from the condensate phase (moles/kg/kg)
- dso4gt :
real*8
- amount of sulfate produced from the gas phase (moles)
- dso4dt :
real*8
- amount of sulfate produced from the condensate phase (moles)
- n :
integer
- index for tracer number loop
- dt_sulf :
real*8, dimension(ntx), intent(inout)
- accumulated diagnostic of sulfate chemistry changes