<div dir="ltr"><div class="gmail_default" style="font-size:large">Hi Everyone,</div><div class="gmail_default" style="font-size:large"><br></div><div class="gmail_default" style="font-size:large">I take this opportunity to ask a question.<br></div><div class="gmail_default" style="font-size:large">I can require conservation of tracers with implicit free surface by setting:</div><div class="gmail_default" style="font-size:large">- linFSConserveTr = .TRUE. in data</div><div class="gmail_default" style="font-size:large">- PTRACERS_linFSConserve(1)= .TRUE. in data.ptracers</div><div class="gmail_default" style="font-size:large"><br></div><div class="gmail_default" style="font-size:large">and, if I use DARWIN instead of DIC package I can also set:</div><div class="gmail_default" style="font-size:large">- darwin_linFSConserve =.TRUE. in data.darwin</div><div class="gmail_default" style="font-size:large"><br></div><div class="gmail_default" style="font-size:large">So my question is: How can I choose which flag to activate? Are these three methods equivalent or are there conflicts if I set them together?</div><div class="gmail_default" style="font-size:large"><br></div><div class="gmail_default" style="font-size:large">Thank you very much for this discussion. I think this topic is really interesting.</div><div class="gmail_default" style="font-size:large"><br></div><div class="gmail_default" style="font-size:large">Francesco</div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">Il giorno mar 5 mar 2024 alle ore 17:13 Gregorio Sandy <<a href="mailto:Sandy.Gregorio@uqar.ca">Sandy.Gregorio@uqar.ca</a>> ha scritto:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div class="msg5870422925765427003">
<div dir="ltr">
<div 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 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><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><span style="font-family:Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)"><br>
</span></div>
<div><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><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><span style="font-family:Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">Good luck,</span></div>
<div><span style="font-family:Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)"><br>
</span></div>
<div><span style="font-family:Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">Sandy</span></div>
<div id="m_3089244907331649317appendonsend"></div>
<hr style="display:inline-block;width:98%">
<div id="m_3089244907331649317divRplyFwdMsg" dir="ltr"><font face="Calibri, sans-serif" style="font-size:11pt" color="#000000"><b>De :</b> MITgcm-support <<a href="mailto:mitgcm-support-bounces@mitgcm.org" target="_blank">mitgcm-support-bounces@mitgcm.org</a>> de la part de Wenrui Jiang <<a href="mailto:wjiang33@jhu.edu" target="_blank">wjiang33@jhu.edu</a>><br>
<b>Envoyé :</b> lundi 4 mars 2024 13:18<br>
<b>À :</b> <a href="mailto:mitgcm-support@mitgcm.org" target="_blank">mitgcm-support@mitgcm.org</a> <<a href="mailto:mitgcm-support@mitgcm.org" target="_blank">mitgcm-support@mitgcm.org</a>><br>
<b>Objet :</b> [MITgcm-support] Closing DIC budget</font>
<div> </div>
</div>
<div dir="ltr">
<table border="0" cellspacing="0" cellpadding="0" align="left" width="100%">
<tbody>
<tr>
<td style="background:rgb(255,0,0);padding:5pt 2pt"></td>
<td width="100%" cellpadding="7px 6px 7px 15px" style="background:rgb(255,224,224);padding:5pt 4pt 5pt 12pt">
<div style="color:rgb(0,0,0)"><span style="color:rgb(0,0,0);font-weight:bold">ATTENTION :</span>
<span style="color:rgb(0,0,0);font-size:10pt">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><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" id="m_3089244907331649317LPlnk840669" target="_blank">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><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><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 style="font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
<br>
</div>
<div 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="m_3089244907331649317x_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: <a href="mailto:wjiang33@jh.edu" target="_blank">wjiang33@jh.edu</a></div>
<div style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
<br>
</div>
</div>
</div>
</div>
</div>
_______________________________________________<br>
MITgcm-support mailing list<br>
<a href="mailto:MITgcm-support@mitgcm.org" target="_blank">MITgcm-support@mitgcm.org</a><br>
<a href="http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support" rel="noreferrer" target="_blank">http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support</a><br>
</div></blockquote></div>