[MITgcm-support] Query regarding BLING

[Matthew Mazloff]

Ah, I see PTRACERS_ref can be used for freshwater value depending on setup. You must have this on for your BLING setup. 

I think the issue is that you are driving ALK close to zero but allowing silica to stay large, violating the assumptions in the model. 

Either you could make the runoff have high alk (perhaps with PTRACERS_ref) or prescribe silica low there when there is significant runoff. Or perhaps try using SIBLING.  Or perhaps use my caps but make them less stringent. Or, if DIC package is stable, just remove my caps and see what happens.

Matt




> On Sep 25, 2023, at 8:19 AM, Matthew Mazloff <mmazloff at ucsd.edu> wrote:
> 
> Are you giving ICs in data.ptracers? If so PTRACERS_ref shouldn’t be being used, I think. 
> 
> Matt
> 
> 
>> On Sep 25, 2023, at 8:43 AM, kunal madkaiker <kunal.madkaiker02 at gmail.com> wrote:
>> 
>> Thanks Matt,
>> 
>> I tried with #undef CARBONCHEM_SOLVESAPHE but the issue persists.
>> 
>> 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)
>> 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).
>> 
>> Kunal
>> 
>> On Mon, Sep 25, 2023 at 7:35 PM Matthew Mazloff <mmazloff at ucsd.edu <mailto:mmazloff at ucsd.edu>> wrote:
>>> Hi Kunal
>>> 
>>> The default bling chemistry is the same as the DIC package. Is your run with the funny behavior using 
>>> #undef CARBONCHEM_SOLVESAPHE
>>> If so, can you check that you get the desired outcome with this undefined.
>>> 
>>> There may be a bug in the BLING SOLVESAPHE code.
>>> 
>>> Matt
>>> 
>>> 
>>> 
>>> > On Sep 25, 2023, at 6:58 AM, kunal madkaiker <kunal.madkaiker02 at gmail.com <mailto:kunal.madkaiker02 at gmail.com>> wrote:
>>> > 
>>> > Hi all,
>>> > 
>>> > 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).
>>> > 
>>> > 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).
>>> > 
>>> > 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.
>>> > 
>>> > Let me know what you think. I am attaching the required files and plots for your reference.
>>> > 
>>> > Thanks,
>>> > Kunal
>>> > <BLING_OPTIONS.h><DIC_OPTIONS.h><data.dic><data.bling><sterr_file><using_bling.gif><using_dic.gif>_______________________________________________
>>> > MITgcm-support mailing list
>>> > MITgcm-support at mitgcm.org <mailto:MITgcm-support at mitgcm.org>
>>> > https://urldefense.com/v3/__http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support__;!!Mih3wA!DU6XPS0r_V-iBbpy1DoUX0ePkKRdvMrbR61lSsmKcyG9Zzh1h9aJz2ECIYMYIci2i2Dseq_pUBQUpbrrkaBMwp8yzH-6$ 
>>> 
>>> _______________________________________________
>>> MITgcm-support mailing list
>>> MITgcm-support at mitgcm.org <mailto:MITgcm-support at mitgcm.org>
>>> http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support <https://urldefense.com/v3/__http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support__;!!Mih3wA!G6DkxmOYF85FN3jomwPi4vliPg8Z4L_Z5zOGMMAAOaOrgFn_wRfAMOu2L7iHgrz6wMrNqKEULF8HTwv5ENaqyznlGaGr$>
>> <using_dic.gif><using_bling.gif>_______________________________________________
>> MITgcm-support mailing list
>> MITgcm-support at mitgcm.org
>> https://urldefense.com/v3/__http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support__;!!Mih3wA!G6DkxmOYF85FN3jomwPi4vliPg8Z4L_Z5zOGMMAAOaOrgFn_wRfAMOu2L7iHgrz6wMrNqKEULF8HTwv5ENaqyznlGaGr$ 
> 

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.mitgcm.org/pipermail/mitgcm-support/attachments/20230925/699d0367/attachment.html>