Changes between Version 7 and Version 8 of CarbonCycleUnits
- Timestamp:
- 11/29/17 11:21:39 (6 years ago)
Legend:
- Unmodified
- Added
- Removed
- Modified
-
CarbonCycleUnits
v7 v8 26 26 27 27 == Inorganic Carbon Equilibrium for Climber, subroutine incche == 28 28 29 '''incche.f'''' calculates the carbonate system equilibrium in the same way than in CLIMBER code. It takes temperature, salinity, DIC and alkalinity as input values and returns sCO2, xpCO2, xCO2, xHCO3, xCO3 as output, where xCO2, xHCO3 and xCO3 are the aqueous components of the carbonate system (carbonic acid, bicarbonate and carbonate), xpCO2 is the partial pressure of the CO,,2,, gas in water according to Henry's Law and sCO2 is the solubility of CO,,2,, in seawater. 30 [[BR]] 29 31 32 This subroutine uses basically the equations described in: F. Millero, Geochimica et Cosmochimica Acta, 59, 661-677, 1995, using some original equations determined by Weiss 1974, Mar. Chem., 2, 203-215. 33 [[BR]] 34 35 incche subroutine calls several other subroutines to calculate different specific constants and variables, such as total boric acid in seawater, first and second dissociation constants of carbonic acid in seawater, Henry's law constant for seawater to calculate solubility of CO,,2,, and so on. 36 [[BR]] 37 38 According to the equations and parameters used in the intermediate calculations, the units for these constants are mol×kg^-1^ or mol×kg^-1^×atm^-1^ 39 [[BR]] 40 41 Therefore, everything indicates that units of xCO2, xHCO3 and xCO3 in the output and DIC and alkalinity in the input must be mol×kg^-1^ too. These are the units provided to MEDUSA in '''mod_iloveclim_o2s.f''' in lines 416 to 419 and converted there to the units required by MEDUSA, mol×m^-3^, multiplying by 1.0e+03 42 [[BR]] 43 44 In lines 258 to 260 of '''ocn_bio.f''' incche subroutine is called using directly variables TM, SM, ODIC, and OALK. Following the previous analysis, OALK and ODIC units must be mol×kg^-1^ or mol×L^-1^ 30 45 31 46 == Scale factors for unit conversion in OCYCC ==