<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
<style type="text/css" style="display:none;"> P {margin-top:0;margin-bottom:0;} </style>
</head>
<body dir="ltr">
<div class="elementToProof" style="font-family: Aptos, Aptos_EmbeddedFont, Aptos_MSFontService, Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
Hi Sandy,</div>
<div class="elementToProof" style="font-family: Aptos, Aptos_EmbeddedFont, Aptos_MSFontService, Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
<br>
</div>
<div class="elementToProof" style="font-family: Aptos, Aptos_EmbeddedFont, Aptos_MSFontService, Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
Thank you very much for the suggestion. Although I have not figured out everything yet, but it is now much closer to a closure. Thank you!<br>
</div>
<div class="elementToProof" style="font-family: Aptos, Aptos_EmbeddedFont, Aptos_MSFontService, Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
<br>
</div>
<div class="elementToProof" style="font-family: Aptos, Aptos_EmbeddedFont, Aptos_MSFontService, Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
Best,
<div>Wenrui</div>
<br>
</div>
<div id="appendonsend"></div>
<hr style="display:inline-block;width:98%" tabindex="-1">
<div id="divRplyFwdMsg" dir="ltr"><font face="Calibri, sans-serif" style="font-size:11pt" color="#000000"><b>From:</b> MITgcm-support <mitgcm-support-bounces@mitgcm.org> on behalf of Gregorio Sandy <Sandy.Gregorio@uqar.ca><br>
<b>Sent:</b> Tuesday, March 5, 2024 11:13 AM<br>
<b>To:</b> mitgcm-support@mitgcm.org <mitgcm-support@mitgcm.org><br>
<b>Subject:</b> Re: [MITgcm-support] Closing DIC budget</font>
<div> </div>
</div>
<style type="text/css" style="display:none">
<!--
p
{margin-top:0;
margin-bottom:0}
-->
</style>
<div dir="ltr">
<p><strong><br>
</strong><strong><span style="background:red; color:white; font-family:"Times New Roman",serif; font-size:12pt"> External Email - Use Caution </span></strong></p>
<p><strong><span style="background:red; color:white; font-family:"Times New Roman",serif; font-size:12pt"><br>
</span></strong></p>
<p><strong><span style="background:red; color:white; font-family:"Times New Roman",serif; font-size:12pt"><br>
</span></strong></p>
<div>
<div class="x_elementToProof" style="font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif; font-size:12pt; color:rgb(0,0,0)">
<span style="font-family:Arial,Helvetica,sans-serif">Hi Wenrui,</span></div>
<div class="x_elementToProof" style="font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif; font-size:12pt; color:rgb(0,0,0)">
<span style="font-family:Arial,Helvetica,sans-serif"><br>
</span></div>
<div class="x_elementToProof"><span style="font-family:Arial,Helvetica,sans-serif; font-size:12pt; color:rgb(0,0,0)">I had a problem as well in closing my tracers budget (by just looking at timeseries of my tracers integrated over my domain, I could see a drift).
I don't know if it can help you but I added the following line in data.ptracers for all my tracers (1 to 5):</span></div>
<div class="x_elementToProof"><span style="font-family:Arial,Helvetica,sans-serif; font-size:12pt; color:rgb(0,0,0)"><br>
</span></div>
<div class="x_elementToProof"><span style="font-family:Arial,Helvetica,sans-serif; font-size:12pt; color:rgb(0,0,0)"><b>PTRACERS_linFSConserve(1)=.TRUE.,</b></span></div>
<pre><div class="x_elementToProof"><span style="font-family:Arial,Helvetica,sans-serif; font-size:12pt; color:rgb(0,0,0)">PTRACERS_linFSConserve applies mean Free-Surf source/sink at surface.</span></div></pre>
<div class="x_elementToProof"><span style="font-family:Arial,Helvetica,sans-serif; font-size:12pt; color:rgb(0,0,0)">Good luck,</span></div>
<div class="x_elementToProof"><span style="font-family:Arial,Helvetica,sans-serif; font-size:12pt; color:rgb(0,0,0)"><br>
</span></div>
<div class="x_elementToProof"><span style="font-family:Arial,Helvetica,sans-serif; font-size:12pt; color:rgb(0,0,0)">Sandy</span></div>
<div id="x_appendonsend"></div>
<hr tabindex="-1" style="display:inline-block; width:98%">
<div id="x_divRplyFwdMsg" dir="ltr"><font face="Calibri, sans-serif" color="#000000" style="font-size:11pt"><b>De :</b> MITgcm-support <mitgcm-support-bounces@mitgcm.org> de la part de Wenrui Jiang <wjiang33@jhu.edu><br>
<b>Envoyé :</b> lundi 4 mars 2024 13:18<br>
<b>À :</b> mitgcm-support@mitgcm.org <mitgcm-support@mitgcm.org><br>
<b>Objet :</b> [MITgcm-support] Closing DIC budget</font>
<div> </div>
</div>
<style type="text/css" style="display:none">
<!--
p
{margin-top:0;
margin-bottom:0}
-->
</style>
<div dir="ltr">
<table border="0" cellspacing="0" cellpadding="0" align="left" width="100%">
<tbody>
<tr>
<td style="background:#ff0000; padding:5pt 2pt 5pt 2pt"></td>
<td width="100%" cellpadding="7px 6px 7px 15px" style="background:#ffe0e0; padding:5pt 4pt 5pt 12pt; word-wrap:break-word">
<div style="color:#000000"><span style="color:#000000; font-weight:bold">ATTENTION :</span>
<span style="color:#000000; font-size:10.0pt">Ce courriel provient d'une source externe. Assurez-vous que cette dernière est sûre avant d'ouvrir une pièce jointe ou de cliquer sur un lien. Dans le doute, signalez le message.
</span></div>
</td>
</tr>
</tbody>
</table>
<br>
<div>
<div class="x_x_elementToProof"><span style="font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif; font-size:11pt; color:rgb(0,0,0)">Dear MITgcm support group,<br>
<br>
Thanks for reading this. My name is Wenrui Jiang, and I am a PhD student from Johns Hopkins. I want to apologize in advance for my inexperience.<br>
<br>
I am trying to close the DIC budget on the biogeochemical tutorial example (<a href="https://mitgcm.readthedocs.io/en/latest/examples/global_oce_biogeo/global_oce_biogeo.html" originalsrc="https://mitgcm.readthedocs.io/en/latest/examples/global_oce_biogeo/global_oce_biogeo.html" shash="LJLkw2k7DT3JyUOzDRbN0rEq87HhroTztBJekzUviVCKkRAN8+ZI2S7Wv8n0gQydYdKwaoxeos5Mb7Q4tjII7tUc3CtSXKMyOfumix21mWQfDIFrFp4zQXoaCpHsUFwbEKtbbw8pXb6ecegGrkpOd8Bmts8i/gRCEJneDQXHiBc=" id="LPlnk840669">https://mitgcm.readthedocs.io/en/latest/examples/global_oce_biogeo/global_oce_biogeo.html</a>).
The diagnostics that I include are:<br>
<br>
TRAC01 |SMR MR|mol/m^3 |Dissolved Inorganic Carbon (DIC) [mol C/m^3] concentration<br>
ForcTr01|SMR MR|mol/m^3/s |DIC forcing tendency<br>
AB_gTr01|SMR MR|mol/m^3/s |DIC tendency from Adams-Bashforth<br>
Tp_gTr01|SMR MR|mol/m^3/s |DIC total transport tendency (before gchem_forcing_sep)<br>
ADVrTr01|WM LR|mol/m^3.m^3/s |Vertical Advective Flux of DIC<br>
ADVxTr01|UU MR|mol/m^3.m^3/s |Zonal Advective Flux of DIC<br>
ADVyTr01|VV MR|mol/m^3.m^3/s |Meridional Advective Flux of DIC<br>
DFrETr01|WM LR|mol/m^3.m^3/s |Vertical Diffusive Flux of DIC (Explicit part)<br>
DFxETr01|UU MR|mol/m^3.m^3/s |Zonal Diffusive Flux of DIC<br>
DFyETr01|VV MR|mol/m^3.m^3/s |Meridional Diffusive Flux of DIC<br>
DFrITr01|WM LR|mol/m^3.m^3/s |Vertical Diffusive Flux of DIC (Implicit part)<br>
<br>
DICBIOA | 15 |SM P MR |mol/m3/sec |Biological Productivity (mol/m3/s)<br>
DICCARB | 15 |SM P MR |mol eq/m3/sec |Carbonate chg-biol prod and remin (mol eq/m3/s)<br>
DICTFLX | 1 |SM P L1 |mol/m3/sec |Tendency of DIC due to air-sea exch (mol/m3/s<br>
<br>
So far, I managed to close the Tp_gTr01 term with the advective and diffusive flux in the interior to machine precision. However, at the surface there is a large residual, which I am not able to close with ForcTR01. I am wondering if this is a dilution effect
from PTRACER_ref.<br>
<br>
I am also having a hard time matching Tp_gTr01 and the tendency diagnosed from DIC and Eta snapshots using the GCHEM diagnostics. The GCHEM diagnostics seems to have very small value (10^-10 mol/m^3/s) which is much smalller than the difference.<br>
<br>
I would really appreciate it if you could tell me whether:<br>
<br>
1. I am using the correct diagnostics in the correct way;<br>
2.there are any model setups that prevent a tracer closure.<br>
</span></div>
<div class="x_x_elementToProof"><span style="font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif; font-size:11pt; color:rgb(0,0,0)"><br>
</span></div>
<div class="x_x_elementToProof"><span style="font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif; font-size:11pt; color:rgb(0,0,0)">Thank you very much!<br>
</span></div>
<div class="x_x_elementToProof" style="font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif; font-size:12pt; color:rgb(0,0,0)">
<br>
</div>
<div class="x_x_elementToProof" style="font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif; font-size:12pt; color:rgb(0,0,0)">
Sincerely,
<div>Wenrui</div>
<br>
</div>
<div id="x_x_Signature">
<div style="font-family:Calibri,Arial,Helvetica,sans-serif; font-size:12pt; color:rgb(0,0,0)">
<span style="font-size:16pt; color:rgb(12,136,42)">Wenrui Jiang</span></div>
<div style="font-family:Calibri,Arial,Helvetica,sans-serif; font-size:12pt; color:rgb(0,0,0)">
<br>
</div>
<div style="font-family:Calibri,Arial,Helvetica,sans-serif; font-size:12pt; color:rgb(0,0,0)">
Ph.D. Student,</div>
<div style="font-family:Calibri,Arial,Helvetica,sans-serif; font-size:12pt; color:rgb(0,0,0)">
Department of Earth & Planetary Sciences,</div>
<div style="font-family:Calibri,Arial,Helvetica,sans-serif; font-size:12pt; color:rgb(0,0,0)">
Johns Hopkins University.</div>
<div style="font-family:Calibri,Arial,Helvetica,sans-serif; font-size:12pt; color:rgb(0,0,0)">
<br>
</div>
<div style="font-family:Calibri,Arial,Helvetica,sans-serif; font-size:12pt; color:rgb(0,0,0)">
Email: wjiang33@jh.edu</div>
<div style="font-family:Calibri,Arial,Helvetica,sans-serif; font-size:12pt; color:rgb(0,0,0)">
<br>
</div>
</div>
</div>
</div>
</div>
</div>
</body>
</html>