File tree 2 files changed +2
-2
lines changed
2 files changed +2
-2
lines changed Original file line number Diff line number Diff line change @@ -29,7 +29,7 @@ sca_T=log((298.15-TEMP_DOXY)/(273.15+TEMP_DOXY))
29
29
30
30
# O2 solubility separated in Scorr and Tcorr (SCOR WG 142 coefficients: Garcia and Gordon 1992, Benson and Krause 1984 refit )
31
31
TCorr <- 44.6596 * (exp(2.00907 + 3.22014 * sca_T + 4.05010 * sca_T ^ 2 + 4.94457 * sca_T ^ 3 - 2.56847e-1 * sca_T ^ 4 + 3.88767 * sca_T ^ 5 ))
32
- Scorr <- exp(PSAL_DOXY * (- 6.24523e-3 - 7.37614e-3 * sca_T - 1.03410e-3 * sca_T ^ 2 - 8.17083e-3 * sca_T ^ 3 )- 4.88682e-7 * PSAL_DOXY ^ 2 )
32
+ Scorr <- exp(PSAL_DOXY * (- 6.24523e-3 - 7.37614e-3 * sca_T - 1.03410e-2 * sca_T ^ 2 - 8.17083e-3 * sca_T ^ 3 )- 4.88682e-7 * PSAL_DOXY ^ 2 )
33
33
34
34
35
35
PPOX = 1013.25 * DOXY * (xO2 * (atm_press - pH2Osatsal ))/ (TCorr * Scorr )* exp(Vm * PRES_DOXY / (R * (TEMP_DOXY + 273.15 ))) # en mbar
Original file line number Diff line number Diff line change @@ -32,7 +32,7 @@ sca_T=log((298.15-TEMP_DOXY)/(273.15+TEMP_DOXY))
32
32
33
33
# O2 solubility separated in Scorr and Tcorr (SCOR WG 142 coefficients: Garcia and Gordon 1992, Benson and Krause 1984 refit )
34
34
TCorr <- 44.6596 * (exp(2.00907 + 3.22014 * sca_T + 4.05010 * sca_T ^ 2 + 4.94457 * sca_T ^ 3 - 2.56847e-1 * sca_T ^ 4 + 3.88767 * sca_T ^ 5 ))
35
- Scorr <- exp(PSAL_DOXY * (- 6.24523e-3 - 7.37614e-3 * sca_T - 1.03410e-3 * sca_T ^ 2 - 8.17083e-3 * sca_T ^ 3 )- 4.88682e-7 * PSAL_DOXY ^ 2 )
35
+ Scorr <- exp(PSAL_DOXY * (- 6.24523e-3 - 7.37614e-3 * sca_T - 1.03410e-2 * sca_T ^ 2 - 8.17083e-3 * sca_T ^ 3 )- 4.88682e-7 * PSAL_DOXY ^ 2 )
36
36
37
37
DOXY = PPOX * (TCorr * Scorr )/ (exp(Vm * PRES_DOXY / (R * (TEMP_DOXY + 273.15 )))* (xO2 * (atm_press - pH2Osatsal ))) # (DOXY en micromol / l)
38
38
You can’t perform that action at this time.
0 commit comments