[MITgcm-support] Query regarding BLING

[Matthew Mazloff]

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/4ff81e78/attachment-0001.html>