Replies: 6 comments 9 replies
-
|
Thanks for the comment. I don't immediately know the answer. In order to look into this further, can you let us know which version of OceanBioME you are using and what value of salinity you are using for this test? |
Beta Was this translation helpful? Give feedback.
-
|
Hi @j-ostiguy, This is interesting thank you. Firstly for the simulations with DIC fields, I found a bug in the air/sea flux code that As for the carbon chemistry, it is meant to follow "Dickson, A.G., Sabine, C.L. and Christian, J.R. (2007), Guide to Best Practices for Ocean CO 2 Measurements. PICES Special Publication 3, 191 pp." which looks similar to the Follows paper you linked. Could you share the exact setup you used to calculate them? |
Beta Was this translation helpful? Give feedback.
-
|
Could this be because OceanBioME.jl/src/Models/GasExchange/carbon_dioxide_concentration.jl Lines 31 to 59 in 3c7465a Perhaps when I implemented this I miss understood something about the algorithm. |
Beta Was this translation helpful? Give feedback.
-
|
I checked with |
Beta Was this translation helpful? Give feedback.
-
|
Perhaps I should add a function to output pCO2, or change the default output of cc(; DIC...) to be pCO2 |
Beta Was this translation helpful? Give feedback.
-
|
The difference is because of the units - I was using the same input Alk/DIC for all three algorithms but didn't realize that carbon_chemistry takes in mmol/m^3. It could be helpful to clarify this more explicitly in the documentation (independently of LOBSTER, which I just found lists the units). Thank you for looking into this so thoroughly! |
Beta Was this translation helpful? Give feedback.

Uh oh!
There was an error while loading. Please reload this page.
-
Hello,
I've done a comparison between the carbonate chemistry solver implemented here, the algorithm presented in this paper, and the CO2SYS algorithm (CO2System.jl). For a range of DIC/Alkalinity or temperature, the calculated pCO2 can differ by up to ~5 to 15 uatm which is shown in the figure below.
I've also run Oceananigans simulations with both the OceanBioME carbonate solver and the Follows (2006) algorithm and found that the choice of carbonate solver results in DIC fields that are quite different. I'm wondering if you have any thoughts on this or if this kind of comparison has already been done. I'm also curious to know if the algorithm chosen here has some advantage over the others.
Thank you!
Beta Was this translation helpful? Give feedback.
All reactions