<div dir="ltr"><div>Hi Matt,</div><div><br></div><div>I tried commenting the lines you suggested in bling_carbon_chem.F as well as increasing PTRACERS_ref of DIC and ALK. Also define/undef USE_SIBLING isn't changing anything.</div><div>The values of DIC/ALK at runoff points keep on diluting as the model evolves. However, pH is decreasing there instead of increasing (refer Fig1).</div><div><br></div><div>Besides DIC and ALK, the PTRACERS_ref values for other tracers (O2, PO4, NO3, SiO4, etc) is Zero. Still the conc. are increasing at the runoff points (Fig2).</div><div><br></div><div>Are the PTRACERS_ref values being overwritten somehow? Is there a way I can completely switch off tracers at runoff points?</div><div><br></div><div>Thanks,</div><div>Kunal<br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Mon, Sep 25, 2023 at 10:06 PM kunal madkaiker <<a href="mailto:kunal.madkaiker02@gmail.com">kunal.madkaiker02@gmail.com</a>> wrote:<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 dir="ltr"><div>Thanks Matt, I'll try these fixes and let you know.</div><div><br></div><div>Hi Francesco, as of now I am prescribing a constant pCO2 value of 391.E-6 atm.</div><div><br></div><div>Kunal<br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Mon, Sep 25, 2023 at 9:13 PM Dattilo, Francesco <<a href="mailto:fdattilo@ogs.it" target="_blank">fdattilo@ogs.it</a>> wrote:<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 dir="auto">I'm not really expert, but is pCO2 value corrected?<div dir="auto">CO2 partial pressure could also affect the values of dic and alkalinity.</div><div dir="auto"><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">Il lun 25 set 2023, 17:14 kunal madkaiker <<a href="mailto:kunal.madkaiker02@gmail.com" target="_blank">kunal.madkaiker02@gmail.com</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 dir="ltr"><div>Thanks Matt,</div><div><br></div><div>I tried with 
#undef CARBONCHEM_SOLVESAPHE but the issue persists.</div><div><br></div><div>Just an observation, when I keep PTRACERS_ref=0 in data.dic, it doesn't impact the DIC and ALK at those runoff points (ref. using_dic.gif)</div><div>However, whether I keep PTRACERS_ref=0 or any other value, it impacts the conc. of those points (which keep on drifting as the run progresses, ref. using_bling.gif).</div><div><br></div><div>Kunal<br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Mon, Sep 25, 2023 at 7:35 PM Matthew Mazloff <<a href="mailto:mmazloff@ucsd.edu" rel="noreferrer" target="_blank">mmazloff@ucsd.edu</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Hi Kunal<br>
<br>
The default bling chemistry is the same as the DIC package. Is your run with the funny behavior using <br>
#undef CARBONCHEM_SOLVESAPHE<br>
If so, can you check that you get the desired outcome with this undefined.<br>
<br>
There may be a bug in the BLING SOLVESAPHE code.<br>
<br>
Matt<br>
<br>
<br>
<br>
> On Sep 25, 2023, at 6:58 AM, kunal madkaiker <<a href="mailto:kunal.madkaiker02@gmail.com" rel="noreferrer" target="_blank">kunal.madkaiker02@gmail.com</a>> wrote:<br>
> <br>
> Hi all,<br>
> <br>
> I have some issues resolving the carbon chemistry using the BLING package. My domain has river runoff included as a point source. When using the 'DIC' package, DIC and ALK conc. get diluted causing pH to increase at those points (which is correct).<br>
> <br>
> However, when I use the same setup with the BLING package, the DIC and ALK conc. are getting diluted but the pH decreases rapidly (5-6 pH).<br>
> <br>
> Secondly, has the SOLVESAPHE algorithm been fully implemented in BLING or not yet? I tried defining 'CARBONCHEM_SOLVESAPHE' in BLING_OPTIONS.h and 'selectpHsolver=1' in data.bling but then I got the segmentation error. Otherwise, it works fine.<br>
> <br>
> Let me know what you think. I am attaching the required files and plots for your reference.<br>
> <br>
> Thanks,<br>
> Kunal<br>
> <BLING_OPTIONS.h><DIC_OPTIONS.h><data.dic><data.bling><sterr_file><using_bling.gif><using_dic.gif>_______________________________________________<br>
> MITgcm-support mailing list<br>
> <a href="mailto:MITgcm-support@mitgcm.org" rel="noreferrer" target="_blank">MITgcm-support@mitgcm.org</a><br>
> <a href="https://urldefense.com/v3/__http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support__;!!Mih3wA!DU6XPS0r_V-iBbpy1DoUX0ePkKRdvMrbR61lSsmKcyG9Zzh1h9aJz2ECIYMYIci2i2Dseq_pUBQUpbrrkaBMwp8yzH-6$" rel="noreferrer noreferrer" target="_blank">https://urldefense.com/v3/__http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support__;!!Mih3wA!DU6XPS0r_V-iBbpy1DoUX0ePkKRdvMrbR61lSsmKcyG9Zzh1h9aJz2ECIYMYIci2i2Dseq_pUBQUpbrrkaBMwp8yzH-6$</a> <br>
<br>
_______________________________________________<br>
MITgcm-support mailing list<br>
<a href="mailto:MITgcm-support@mitgcm.org" rel="noreferrer" target="_blank">MITgcm-support@mitgcm.org</a><br>
<a href="http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support" rel="noreferrer noreferrer" target="_blank">http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support</a><br>
</blockquote></div>
_______________________________________________<br>
MITgcm-support mailing list<br>
<a href="mailto:MITgcm-support@mitgcm.org" rel="noreferrer" target="_blank">MITgcm-support@mitgcm.org</a><br>
<a href="http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support" rel="noreferrer noreferrer" target="_blank">http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support</a><br>
</blockquote></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>
</blockquote></div>
</blockquote></div>