<html><head><meta http-equiv="content-type" content="text/html; charset=utf-8"></head><body style="overflow-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;">Hi Kunal<div><br></div><div>I just looked more carefully and see you have undef SOLVESAPHE</div><div>However I don’t understand your plots. Top left and right panels are the exact same. Are those both DIC for the two runs.</div><div>These are also different times by 9 months which makes it hard to compare - Salt is super different.</div><div><br></div><div><br></div><div>Anyway, here is what I think is happening.</div><div>I added code:</div><div><span style="background-color: rgb(255, 251, 0);"><font color="#002e7a"><br></font></span></div><div><div><span style="background-color: rgb(255, 251, 0);"><font color="#002e7a">cmm - rare, but ice melt can cause prescribed silicate alkalinity to be</font></span></div><div><span style="background-color: rgb(255, 251, 0);"><font color="#002e7a">cmm larger than total alkalinity. Crop here to account for this</font></span></div><div><span style="background-color: rgb(255, 251, 0);"><font color="#002e7a"> IF (siooh3g.GT.0.2 _d 0*ta) siooh3g = 0.2 _d 0*ta</font></span></div><div><span style="background-color: rgb(255, 251, 0);"><font color="#002e7a">cmick - now estimate carbonate alkalinity</font></span></div><div><span style="background-color: rgb(255, 251, 0);"><font color="#002e7a"> cag = ta - bohg - (kwlocal/hguess) + hguess</font></span></div><div><span style="background-color: rgb(255, 251, 0);"><font color="#002e7a"> & - hpo4g - 2.0 _d 0*po4g + h3po4g</font></span></div><div><span style="background-color: rgb(255, 251, 0);"><font color="#002e7a"> & - siooh3g</font></span></div><div><span style="background-color: rgb(255, 251, 0);"><font color="#002e7a">CMM(</font></span></div><div><span style="background-color: rgb(255, 251, 0);"><font color="#002e7a"> IF ( cag .LT. (0.1 _d 0*ta) ) THEN</font></span></div><div><span style="background-color: rgb(255, 251, 0);"><font color="#002e7a"> print *,'CMMcag1', cag,i,j,k,xC(i,j,1,1),yC(i,j,1,1)</font></span></div><div><span style="background-color: rgb(255, 251, 0);"><font color="#002e7a">C myXGlobalLo-1+(bi-1)*sNx+i, myYGlobalLo-1+(bj-1)*sNy+j</font></span></div><div><span style="background-color: rgb(255, 251, 0);"><font color="#002e7a"> print *,'CMMcag2',ta,bohg,(kwlocal/hguess),hguess,hpo4g,</font></span></div><div><span style="background-color: rgb(255, 251, 0);"><font color="#002e7a"> & 2.0 _d 0*po4g,h3po4g, siooh3g</font></span></div><div><span style="background-color: rgb(255, 251, 0);"><font color="#002e7a"> cag = 0.1 _d 0*ta</font></span></div><div><span style="background-color: rgb(255, 251, 0);"><font color="#002e7a"> ENDIF</font></span></div><div><span style="background-color: rgb(255, 251, 0);"><font color="#002e7a">CMM)</font></span></div></div><div><br></div><div>When one inputs fresh water it dilutes alk and should dilute silica but we are prescribing silica and it doesn’t drop. This can make carbonate alkalinity go negative and blow up the carbon system. I put caps here to assume silicate alk is always less than 20% of total and to assume carbonate alkalinity can’t go lower than 10% of total alkalinity. </div><div>This keeps the carbon system stable, but perhaps for your case causes spurious results. Maybe these caps are too strong. </div><div><br></div><div>I notice you give different silica files for DIC and BLING. Maybe you are giving a low enough silica to DIC to keep it stable but high enough to BLING to set off these caps. </div><div><br></div><div>You could try commenting out </div><div><div><span style="background-color: rgb(255, 251, 0);"><font color="#002e7a"> IF (siooh3g.GT.0.2 _d 0*ta) siooh3g = 0.2 _d 0*ta</font></span></div><div><font color="#002e7a"><span style="caret-color: rgb(0, 46, 122);">…</span></font></div><div><span style="color: rgb(0, 46, 122); background-color: rgb(255, 251, 0);">CMM(</span></div><div><span style="background-color: rgb(255, 251, 0);"><font color="#002e7a"> IF ( cag .LT. (0.1 _d 0*ta) ) THEN</font></span></div><div><span style="background-color: rgb(255, 251, 0);"><font color="#002e7a"> print *,'CMMcag1', cag,i,j,k,xC(i,j,1,1),yC(i,j,1,1)</font></span></div><div><span style="background-color: rgb(255, 251, 0);"><font color="#002e7a">C myXGlobalLo-1+(bi-1)*sNx+i, myYGlobalLo-1+(bj-1)*sNy+j</font></span></div><div><span style="background-color: rgb(255, 251, 0);"><font color="#002e7a"> print *,'CMMcag2',ta,bohg,(kwlocal/hguess),hguess,hpo4g,</font></span></div><div><span style="background-color: rgb(255, 251, 0);"><font color="#002e7a"> & 2.0 _d 0*po4g,h3po4g, siooh3g</font></span></div><div><span style="background-color: rgb(255, 251, 0);"><font color="#002e7a"> cag = 0.1 _d 0*ta</font></span></div><div><span style="background-color: rgb(255, 251, 0);"><font color="#002e7a"> ENDIF</font></span></div><div><span style="background-color: rgb(255, 251, 0);"><font color="#002e7a">CMM)</font></span></div></div><div><span style="background-color: rgb(255, 251, 0);"><font color="#002e7a"><br></font></span></div><div>in bling_carbon_chem.F</div><div>and see if it is then the same as DIC</div><div><br></div><div>(I also put caps on the range of pH so this may end up kicking in and keeping the model stable if needed)</div><div><br></div><div>I can try to help further but to do so it would be good to see ALK, Salt, Silica, and pH all for the same timestep, and ideally early in the run when DIC and BLING first diverge.</div><div><br></div><div>Matt</div><div><br></div><div><br></div><div><br></div><div><div><br><blockquote type="cite"><div>On Sep 25, 2023, at 6:58 AM, kunal madkaiker <kunal.madkaiker02@gmail.com> wrote:</div><br class="Apple-interchange-newline"><div><div dir="ltr"><div>Hi all,</div><div><br></div><div>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).</div><div><br></div><div>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).</div><div><br></div><div>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.</div><div><br></div><div>Let me know what you think. I am attaching the required files and plots for your reference.</div><div><br></div><div>Thanks,</div><div>Kunal<br></div></div>
<span id="cid:f_lmyxplyy0"><BLING_OPTIONS.h></span><span id="cid:f_lmyxpwkn1"><DIC_OPTIONS.h></span><span id="cid:f_lmyxrsi92"><data.dic></span><span id="cid:f_lmyxscba3"><data.bling></span><span id="cid:f_lmyxsvcb6"><sterr_file></span><span id="cid:f_lmyxsvc15"><using_bling.gif></span><span id="cid:f_lmyxsvbn4"><using_dic.gif></span>_______________________________________________<br>MITgcm-support mailing list<br>MITgcm-support@mitgcm.org<br>https://urldefense.com/v3/__http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support__;!!Mih3wA!DU6XPS0r_V-iBbpy1DoUX0ePkKRdvMrbR61lSsmKcyG9Zzh1h9aJz2ECIYMYIci2i2Dseq_pUBQUpbrrkaBMwp8yzH-6$ <br></div></blockquote></div><br></div></body></html>