<div dir="ltr"><div>Dear Christoph,</div><div><br></div><div>Thank you for your email. I'll check the code as you suggested.</div><div><br></div><div>Regards<br></div><div>-Kunal<br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Thu, Apr 1, 2021 at 11:58 PM Christoph Voelker <<a href="mailto:christoph.voelker@awi.de">christoph.voelker@awi.de</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>
    <p>Dear Kunal, <br>
    </p>
    <p>to some extent yes. There is a setting of the model where calcite
      dissolution is calculated as a function of calcite saturation
      state, which is switched on by setting CAR_DISS=.true. in
      DIC_OPTIONS.h. That was originally put into the model by Kasten
      Friis. <br>
    </p>
    <p>In that case, the saturation of CaCO3 is calculated at every
      depth, which implies a calculation of the full carbonate
      chemistry, including the pressure dependence of the equilibrium
      constants. Have a look into the routine calcite_saturation.F if
      you want to understand how. <br>
    </p>
    <p>Note, hoewever, that it is computationally expensive to do so,
      and because of that the standard setting is not to calculate this
      at every time step, but only every now and then. <br>
    </p>
    <p>I hope that helps to get you started. <br>
    </p>
    <p>Cheers, Christoph<br>
    </p>
    <div>Am 01.04.21 um 15:03 schrieb kunal
      madkaiker:<br>
    </div>
    <blockquote type="cite">
      
      <div dir="ltr">
        <div>Dear MITgcm users,</div>
        <div><br>
        </div>
        <div>My question pertains towards the DIC package, which
          provides pH, pCO2, CO2 flux, etc</div>
        <div>as 'surface diagnostics'. As the name suggests, it gives
          output only at surface level.<br>
        </div>
        <div><br>
        </div>
        <div>Is there a way I can simulate these parameters in depth as
          well? Any suggestions?</div>
        <div><br>
        </div>
        <div>-Kunal<br>
        </div>
      </div>
      <br>
      <fieldset></fieldset>
      <pre>_______________________________________________
MITgcm-support mailing list
<a href="mailto:MITgcm-support@mitgcm.org" target="_blank">MITgcm-support@mitgcm.org</a>
<a href="http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support" target="_blank">http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support</a>
</pre>
    </blockquote>
    <pre cols="72">-- 
Christoph Voelker
Alfred Wegener Institute 
Helmholtz Centre for Polar and Marine Research
Am Handelshafen 12
27570 Bremerhaven, Germany
e: <a href="mailto:Christoph.Voelker@awi.de" target="_blank">Christoph.Voelker@awi.de</a>
t: +49 471 4831 1848</pre>
  </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>