<div dir="ltr">Hi Kaitlin, <div><br></div><div>I believe, and maybe this is what Julian was also alluding to is, that you can have only Redi and no GM if you make a choice of parameters like: </div><div>'</div><div>GM_background_K = 0.,<br>GM_isopycK = 200.,'<br></div><div><br></div><div>GM_background_K is K_gm and GM_isopyckK is K_redi. In default if you set only one the other is taken to be equal, otherwise you set them individually. This is probably what Martin meant by turning off GM and only keeping Redi, please correct me if I am wrong. They are 2 different parameterization, which cleanly combine in a single tensor. </div><div><br></div><div>Best,</div><div>Dhruv</div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Wed, Feb 20, 2019 at 7:05 AM Naughten, Kaitlin A. <<a href="mailto:kaight@bas.ac.uk">kaight@bas.ac.uk</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 id="gmail-m_-461945393649954043divtagdefaultwrapper" style="font-size:12pt;color:rgb(0,0,0);font-family:Calibri,Helvetica,sans-serif" dir="ltr">
<p style="margin-top:0px;margin-bottom:0px">Hi Julian,</p>
<p style="margin-top:0px;margin-bottom:0px"><br>
</p>
<p style="margin-top:0px;margin-bottom:0px">Thanks for your reply, and for pointing me towards the streamfunction calculation. It looks like since I have GM_AdvForm=.false. (the default) I'm using the skew-flux form even though <span>GM_BOLUS_ADVEC is defined.
This means the streamfunction doesn't actually get calculated. I guess a lot of the variables which are zero in my output only relate to the bolus advective form. </span></p>
<p style="margin-top:0px;margin-bottom:0px"><span><br>
</span></p>
<p style="margin-top:0px;margin-bottom:0px"><span>According to the documentation (<a href="https://urldefense.proofpoint.com/v2/url?u=https-3A__mitgcm.readthedocs.io_en_latest_phys-5Fpkgs_gmredi.html&d=DwMFAg&c=slrrB7dE8n7gBJbeO0g-IQ&r=3huS5z96YHlUbTqcqC5pEg&m=qMZixtwdhH-XWepZS7WdMJSmyBQjqdFvnY1NLjtajCo&s=3X0aqIiRYHqFu8bdd2ajIGPQlkZvpPiQmndEjZWqEVQ&e=" class="gmail-m_-461945393649954043OWAAutoLink" target="_blank">https://mitgcm.readthedocs.io/en/latest/phys_pkgs/gmredi.html</a>), the skew-flux form
can achieve the "GM" part of GM-Redi just by modifying two elements of the Redi tensor, and doesn't need to calculate streamfunctions or bolus velocities or anything like that. So it appears I do have both parts of the scheme turned on, even though many variables
are zero.</span></p>
<p style="margin-top:0px;margin-bottom:0px"><span><br>
</span></p>
<p style="margin-top:0px;margin-bottom:0px"><span>I don't see how it's possible to only have the Redi part, given the MITgcm options, but it sounds like Martin has done that? Maybe it was an older version of the code.</span></p>
<p style="margin-top:0px;margin-bottom:0px"><span><br>
</span></p>
<p style="margin-top:0px;margin-bottom:0px"><span>All the best,</span></p>
<p style="margin-top:0px;margin-bottom:0px"><span>Kaitlin</span></p>
<p style="margin-top:0px;margin-bottom:0px"><br>
</p>
<div id="gmail-m_-461945393649954043Signature">
<div id="gmail-m_-461945393649954043divtagdefaultwrapper" dir="ltr" style="font-size:12pt;color:rgb(0,0,0);font-family:Calibri,Helvetica,sans-serif,EmojiFont,"Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols">
<p style="margin-top:0px;margin-bottom:0px"></p>
<p style="font-family:Calibri,Helvetica,sans-serif,EmojiFont,"Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols;font-size:16px">
</p>
<div><span style="font-size:11pt;color:rgb(117,123,128)">Dr Kaitlin Naughten</span></div>
<span style="font-size:11pt;color:rgb(117,123,128)"></span>
<div><span style="font-size:11pt;color:rgb(117,123,128)">British Antarctic Survey</span></div>
<p></p>
</div>
</div>
</div>
<hr style="display:inline-block;width:98%">
<div id="gmail-m_-461945393649954043divRplyFwdMsg" dir="ltr"><font face="Calibri, sans-serif" style="font-size:11pt" color="#000000"><b>From:</b> MITgcm-support <<a href="mailto:mitgcm-support-bounces@mitgcm.org" target="_blank">mitgcm-support-bounces@mitgcm.org</a>> on behalf of Julian Mak <<a href="mailto:julian.mak@physics.ox.ac.uk" target="_blank">julian.mak@physics.ox.ac.uk</a>><br>
<b>Sent:</b> 20 February 2019 11:24:59<br>
<b>To:</b> MITgcm Support<br>
<b>Subject:</b> Re: [MITgcm-support] GM-Redi question</font>
<div> </div>
</div>
<div style="background-color:rgb(255,255,255)">
<p>I think it's "GM_isopycK" which is the Redi part? By default if that is not specified I think it takes on the value of GM (although I don't know if this would be "GM_background_K" or whatever the kGM value that the Visbeck scheme spits out in this case).</p>
<p><br>
</p>
<p>A data.gmredi I used is given below; ignore the MMB stuff which is our GEOMETRIC implementation where we hijacked the existing Visbeck scheme in a rather disgusting fashion (I haven't tried very hard tidying it up...)<br>
</p>
<p><br>
</p>
<p>Best wishes,<br>
Julian<br>
</p>
<p><br>
</p>
<p>%%%%%%%%%%%%%%%%%%%%<br>
</p>
<p><br>
</p>
<p>&GM_PARM01 <br>
GM_AdvForm = .FALSE.<br>
GM_background_K = 0.,<br>
GM_isopycK = 200.,<br>
GM_taper_scheme = 'gkw91',<br>
GM_maxSlope = 5.e-3,<br>
GM_Kmin_horiz = 0.,<br>
GM_Scrit = 4.e-3,<br>
GM_Sd = 1.e-3,<br>
GM_Visbeck_alpha = 0.016,<br>
GM_Visbeck_minval_K= 50.,<br>
GM_Visbeck_maxval_K= 1.e+4,<br>
# MMB params<br>
MMB_lambda = 1.e-7,<br>
MMB_alpha = 0.04,<br>
MMB_minval_K = 50.,<br>
MMB_maxval_K = 1.e+4,<br>
# MMB_minval_K = 1155,<br>
# MMB_maxval_K = 1155,<br>
energy_init = 1.e-2,<br>
energy_kappa = 0,<br>
energy_local = .false.,</p>
<p><br>
</p>
<p>%%%%%%%%%%%%%%%%%%%%</p>
<p><br>
</p>
<div class="gmail-m_-461945393649954043x_moz-cite-prefix">On 20/02/2019 10:34, Naughten, Kaitlin A. wrote:<br>
</div>
<blockquote type="cite">
<div id="gmail-m_-461945393649954043x_divtagdefaultwrapper" dir="ltr" style="font-size:12pt;color:rgb(0,0,0);font-family:Calibri,Helvetica,sans-serif">
<p style="margin-top:0px;margin-bottom:0px">Hello everyone,</p>
<p style="margin-top:0px;margin-bottom:0px"><br>
</p>
<p style="margin-top:0px;margin-bottom:0px">I am a bit confused about the GM-Redi scheme, and how to turn on the "GM" part (bolus or skew fluxes) versus just the "Redi" part (tracer diffusion along isentropes). An old message from Martin (<a href="https://urldefense.proofpoint.com/v2/url?u=http-3A__mailman.mitgcm.org_pipermail_mitgcm-2Dsupport_2016-2DMarch_010355.html&d=DwMFAg&c=slrrB7dE8n7gBJbeO0g-IQ&r=3huS5z96YHlUbTqcqC5pEg&m=qMZixtwdhH-XWepZS7WdMJSmyBQjqdFvnY1NLjtajCo&s=MRoT3cqh1pVjTcv9E0n5wsOVplrjIe6CK9QGiWF59mE&e=" class="gmail-m_-461945393649954043x_OWAAutoLink" target="_blank">http://mailman.mitgcm.org/pipermail/mitgcm-support/2016-March/010355.html</a>)
said that he only used the Redi part, but I'm not sure what options control this.</p>
<p style="margin-top:0px;margin-bottom:0px"><br>
</p>
<p style="margin-top:0px;margin-bottom:0px">I switched on all the GM-Redi diagnostics in my configuration and discovered that many variables are zero everywhere (GM_PsiX, GM_PsiY, GM_ubT, GM_vbT, GM_hTrsl, GM_baseS, GM_rLamb, plus two components of the Redi tensor
which I think are supposed to be zero: GM_Kuz and GM_Kvz). Does this mean that I only have the Redi part switched on? And if so, why is the GM part not on?</p>
<p style="margin-top:0px;margin-bottom:0px"><br>
</p>
<p style="margin-top:0px;margin-bottom:0px">Here are the active lines of <span style="font-size:12pt">my data.gmredi:</span></p>
<p style="margin-top:0px;margin-bottom:0px"><br>
</p>
<div> &GM_PARM01</div>
<div> GM_background_K=10.,</div>
<div> GM_taper_scheme='ldd97',</div>
<div> GM_Kmin_horiz=50.,</div>
<div> GM_Visbeck_alpha=0.07, </div>
<div> GM_Visbeck_length=50.E3,</div>
<div> GM_Visbeck_depth=2000.,</div>
<div> GM_Visbeck_maxVal_K=300.,</div>
<br>
<p style="margin-top:0px;margin-bottom:0px">And here are the active lines of my GM_REDI_OPTIONS.h:</p>
<p style="margin-top:0px;margin-bottom:0px"><span style="font-size:12pt"><br>
</span></p>
<p style="margin-top:0px;margin-bottom:0px"><span style="font-size:12pt">#ifndef GMREDI_OPTIONS_H</span><br>
</p>
<div>#define GMREDI_OPTIONS_H</div>
<div>#include "PACKAGES_CONFIG.h"</div>
<div>#include "CPP_OPTIONS.h"</div>
<div>#ifdef ALLOW_GMREDI</div>
<div>#define GM_VISBECK_VARIABLE_K</div>
<div>#undef OLD_VISBECK_CALC</div>
<div>#undef GM_K3D</div>
<div>#undef GM_K3D_PASSIVE</div>
<div>#define GM_NON_UNITY_DIAGONAL</div>
<div>#define GM_EXTRA_DIAGONAL</div>
<div>#define GM_BOLUS_ADVEC</div>
<div>#undef GM_BOLUS_BVP</div>
<div>#endif /* ALLOW_GMREDI */</div>
<div>#endif /* GMREDI_OPTIONS_H */</div>
<div><br>
</div>
<div>This is probably a stupid question but hopefully that means it has a quick and easy answer!</div>
<div><br>
</div>
<div>Many thanks,</div>
<div>Kaitlin</div>
<p style="margin-top:0px;margin-bottom:0px"><br>
</p>
<div id="gmail-m_-461945393649954043x_Signature">
<div id="gmail-m_-461945393649954043x_divtagdefaultwrapper" dir="ltr">
<p></p>
<div><span style="font-size:11pt;color:rgb(117,123,128)">Dr Kaitlin Naughten</span></div>
<span style="font-size:11pt;color:rgb(117,123,128)"></span>
<div><span style="font-size:11pt;color:rgb(117,123,128)">British Antarctic Survey</span></div>
</div>
</div>
</div>
<br>
<br>
This email and any attachments are intended solely for the use of the named recipients. If you are not the intended recipient you must not use, disclose, copy or distribute this email or any of its attachments and should notify the sender immediately and delete
this email from your system. <br>
UK Research and Innovation has taken every reasonable precaution to minimise risk of this email or any attachments containing viruses or malware but the recipient should carry out its own virus and malware checks before opening the attachments. UK Research
and Innovation does not accept any liability for any losses or damages which the recipient may sustain due to presence of any viruses.
<br>
Opinions, conclusions or other information in this message and attachments that are not related directly to UK Research and Innovation business are solely those of the author and do not represent the views of UK Research and Innovation.<br>
<br>
<br>
<fieldset class="gmail-m_-461945393649954043x_mimeAttachmentHeader"></fieldset>
<pre class="gmail-m_-461945393649954043x_moz-quote-pre">_______________________________________________
MITgcm-support mailing list
<a class="gmail-m_-461945393649954043x_moz-txt-link-abbreviated" href="mailto:MITgcm-support@mitgcm.org" target="_blank">MITgcm-support@mitgcm.org</a>
<a class="gmail-m_-461945393649954043x_moz-txt-link-freetext" href="https://urldefense.proofpoint.com/v2/url?u=http-3A__mailman.mitgcm.org_mailman_listinfo_mitgcm-2Dsupport&d=DwMFAg&c=slrrB7dE8n7gBJbeO0g-IQ&r=3huS5z96YHlUbTqcqC5pEg&m=qMZixtwdhH-XWepZS7WdMJSmyBQjqdFvnY1NLjtajCo&s=rJAyt1xstYmTTDfN02t-tWC90c-srw-K0EGi6DtIb34&e=" target="_blank">http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support</a>
</pre>
</blockquote>
<pre class="gmail-m_-461945393649954043x_moz-signature" cols="72">--
Julian Mak,
Department of Physics,
University of Oxford,
Clarendon Laboratory, Parks Road
Oxford, OX1 3PU
e-mail: <a class="gmail-m_-461945393649954043x_moz-txt-link-abbreviated" href="mailto:julian.mak@physics.ox.ac.uk" target="_blank">julian.mak@physics.ox.ac.uk</a>
www: <a class="gmail-m_-461945393649954043x_moz-txt-link-freetext" href="https://urldefense.proofpoint.com/v2/url?u=https-3A__sites.google.com_site_julianclmak_home&d=DwMFAg&c=slrrB7dE8n7gBJbeO0g-IQ&r=3huS5z96YHlUbTqcqC5pEg&m=qMZixtwdhH-XWepZS7WdMJSmyBQjqdFvnY1NLjtajCo&s=E1OXDc28ZD97AQY00ipa6TKTD-N-vyxoskTTdNi4yAE&e=" target="_blank">https://sites.google.com/site/julianclmak/home</a></pre>
</div>
<br>
<br>
This email and any attachments are intended solely for the use of the named recipients. If you are not the intended recipient you must not use, disclose, copy or distribute this email or any of its attachments and should notify the sender immediately and delete
this email from your system. <br>
UK Research and Innovation has taken every reasonable precaution to minimise risk of this email or any attachments containing viruses or malware but the recipient should carry out its own virus and malware checks before opening the attachments. UK Research
and Innovation does not accept any liability for any losses or damages which the recipient may sustain due to presence of any viruses.
<br>
Opinions, conclusions or other information in this message and attachments that are not related directly to UK Research and Innovation business are solely those of the author and do not represent the views of UK Research and Innovation.<br>
<br>
</div>
_______________________________________________<br>
MITgcm-support mailing list<br>
<a href="mailto:MITgcm-support@mitgcm.org" target="_blank">MITgcm-support@mitgcm.org</a><br>
<a href="https://urldefense.proofpoint.com/v2/url?u=http-3A__mailman.mitgcm.org_mailman_listinfo_mitgcm-2Dsupport&d=DwICAg&c=slrrB7dE8n7gBJbeO0g-IQ&r=3huS5z96YHlUbTqcqC5pEg&m=qMZixtwdhH-XWepZS7WdMJSmyBQjqdFvnY1NLjtajCo&s=rJAyt1xstYmTTDfN02t-tWC90c-srw-K0EGi6DtIb34&e=" rel="noreferrer" target="_blank">https://urldefense.proofpoint.com/v2/url?u=http-3A__mailman.mitgcm.org_mailman_listinfo_mitgcm-2Dsupport&d=DwICAg&c=slrrB7dE8n7gBJbeO0g-IQ&r=3huS5z96YHlUbTqcqC5pEg&m=qMZixtwdhH-XWepZS7WdMJSmyBQjqdFvnY1NLjtajCo&s=rJAyt1xstYmTTDfN02t-tWC90c-srw-K0EGi6DtIb34&e=</a><br>
</blockquote></div>