<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=us-ascii">
<style type="text/css" style="display:none;"> P {margin-top:0;margin-bottom:0;} </style>
</head>
<body dir="ltr">
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0); background-color: rgb(255, 255, 255);" class="elementToProof">
Hi Jean-Michel,</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0); background-color: rgb(255, 255, 255);" class="elementToProof">
<br>
</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0); background-color: rgb(255, 255, 255);" class="elementToProof">
Many thanks for your quick and clear response.</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0); background-color: rgb(255, 255, 255);" class="elementToProof">
<br>
</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0); background-color: rgb(255, 255, 255);" class="elementToProof">
I agree that I should try option 2), which appears to calculate both the viscosity and diffusivity for the tracers (simply by using a constant Prandtl number)? Yes, the Smagorinsky model essentially parameterises subgrid mixing by using the resolved shear
and grid resolution, and so does not explicitly account for stratification. I hope that my very fine vertical grid resolves most buoyancy-driven (and mechanical) mixing. </div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0); background-color: rgb(255, 255, 255);" class="elementToProof">
<br>
</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0); background-color: rgb(255, 255, 255);" class="elementToProof">
MITgcm does have a TKE 1.5-order closure scheme (i.e., package ggl90), however, it seems only suitable for hydrostatic simulations, when the mixing length in the scheme doesn't need to be limited to the grid resolution. This at least is my understanding from
the paper by Blanke and Delecluse (1993).</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0); background-color: rgb(255, 255, 255);" class="elementToProof">
<br>
</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0); background-color: rgb(255, 255, 255);" class="elementToProof">
My grid is certainly not isotropic, but I have no bottom topography so the grid cells are at least aligned with the horizontal. I am not sure what errors my grid may introduce when using the Smagorinsky scheme? I guess I can test the results using the Smagorinsky
scheme against simply specifying a profile of viscosity and diffusivity using the relevant Reynolds number (i.e. setting <span style="font-size:12pt;background-color:rgb(255, 255, 255);display:inline !important" class="ContentPasted0">viscArNr).</span></div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0); background-color: rgb(255, 255, 255);" class="elementToProof">
<br>
</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0); background-color: rgb(255, 255, 255);" class="elementToProof">
<br>
</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0); background-color: rgb(255, 255, 255);" class="elementToProof">
Cheers,</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0); background-color: rgb(255, 255, 255);" class="elementToProof">
Paul</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0); background-color: rgb(255, 255, 255);" class="elementToProof">
<br>
</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0); background-color: rgb(255, 255, 255);" class="elementToProof">
<br>
</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0); background-color: rgb(255, 255, 255);" class="elementToProof">
<br>
</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0); background-color: rgb(255, 255, 255);" class="elementToProof">
<br>
</div>
<div class="elementToProof">
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
<br>
</div>
<div id="Signature">
<div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
<div style="margin:0px;font-size:15px;font-family:"Segoe UI", "Segoe UI Web (West European)", "Segoe UI", -apple-system, system-ui, Roboto, "Helvetica Neue", sans-serif;color:rgb(32, 31, 30);text-align:start;background-color:rgb(255, 255, 255)">
<p style="font-size:11pt;font-family:Calibri, sans-serif;margin:0px"><span style="color: rgb(12, 100, 192);"><u> </u></span></p>
<p style="font-size:11pt;font-family:Calibri, sans-serif;margin:0px"><span style="margin:0px;color:rgb(31, 78, 121)">Dr Paul Burns MInstP FRMetS</span></p>
<p style="font-size:11pt;font-family:Calibri, sans-serif;margin:0px"><span style="margin:0px;color:rgb(31, 78, 121)">Research Fellow</span></p>
<p style="font-size:11pt;font-family:Calibri, sans-serif;margin:0px"><span style="margin:0px;color:rgb(31, 78, 121)">Geophysical and Astrophysical Fluid Dynamics Group</span></p>
<p style="font-size:11pt;font-family:Calibri, sans-serif;margin:0px"><span style="margin:0px;color:rgb(31, 78, 121)">Department of Mathematics</span></p>
<p style="font-size:11pt;font-family:Calibri, sans-serif;margin:0px"><span style="margin:0px;color:rgb(31, 78, 121)">Harrison Building</span></p>
<p style="font-size:11pt;font-family:Calibri, sans-serif;margin:0px"><span style="margin:0px;color:rgb(31, 78, 121)">University of Exeter</span></p>
</div>
<br>
</div>
</div>
</div>
</div>
<div id="appendonsend"></div>
<hr style="display:inline-block;width:98%" tabindex="-1">
<div id="divRplyFwdMsg" dir="ltr"><font face="Calibri, sans-serif" style="font-size:11pt" color="#000000"><b>From:</b> MITgcm-support <mitgcm-support-bounces@mitgcm.org> on behalf of Jean-Michel Campin <jmc@mit.edu><br>
<b>Sent:</b> 01 February 2023 15:25<br>
<b>To:</b> mitgcm-support@mitgcm.org <mitgcm-support@mitgcm.org><br>
<b>Subject:</b> Re: [MITgcm-support] Vertical diffusion for non-hydrostatic ocean model with a stretched grid</font>
<div> </div>
</div>
<div class="BodyFragment"><font size="2"><span style="font-size:11pt;">
<div class="PlainText">CAUTION: This email originated from outside of the organisation. Do not click links or open attachments unless you recognise the sender and know the content is safe.<br>
<br>
<br>
Hi Paul,<br>
<br>
There is some confusion about Smagorinsky in MITgcm since we have 2 implementations:<br>
1) a Smagorinsky like horizontal viscosity that is function of horizontal shear,<br>
with the choice of harmonic (viscC2smag) or bi-harmonic (viscC4smag) viscosity.<br>
2) a simple 3-D Smagorinsky harmonic viscosity that work best when the grid is<br>
near isotropic (mostly for Non-Hydrostatic simulation) and does not acount for<br>
stratification.<br>
Need ALLOW_SMAG_3D_DIFFUSIVITY to be defined and uses run-time parameter:<br>
useSmag3D = .TRUE., with smag3D_coeff,<br>
see experiment: tutorial_deep_convection, with input.smag3d/data<br>
<br>
You may want to try this second implementation.<br>
Alternatively, you could always set a vertical profile of (background)<br>
vertical viscosity using Nr values to set vector "viscArNr" in main parameter<br>
input file "data" (note that the fist value corresponds to the very top,<br>
i.e., surface, vertical index k=1, and has no effect).<br>
<br>
Cheers,<br>
Jean-Michel<br>
<br>
On Wed, Feb 01, 2023 at 12:55:22PM +0000, Burns, Paul wrote:<br>
> Hi,<br>
><br>
> I was hoping that someone could give me some advice about the best way to define the vertical viscosity (and diffusivity) in MITgcm for a non-hydrostatic model using a stretched vertical grid. I've used a stretched vertical grid to resolve very small features
of interest (~0.2 m) in the upper ocean and to keep the simulation practical.<br>
><br>
> The viscosity should vary with the grid resolution (e.g., Smagorinsky, 1993), however, the MITgcm documentation says that the Smagorinsky scheme isn't implemented in the vertical. So would a reasonable approach be to compute a vertical profile of viscosities
using a vertical Reynolds number and use that to set viscArNr?<br>
><br>
><br>
> Many thanks,<br>
> Paul<br>
><br>
><br>
><br>
><br>
> Dr Paul Burns MInstP FRMetS<br>
><br>
> Research Fellow<br>
><br>
> Geophysical and Astrophysical Fluid Dynamics Group<br>
><br>
> Department of Mathematics<br>
><br>
> Harrison Building<br>
><br>
> University of Exeter<br>
><br>
<br>
> _______________________________________________<br>
> MITgcm-support mailing list<br>
> MITgcm-support@mitgcm.org<br>
> <a href="https://eur03.safelinks.protection.outlook.com/?url=http%3A%2F%2Fmailman.mitgcm.org%2Fmailman%2Flistinfo%2Fmitgcm-support&data=05%7C01%7Cp.burns2%40exeter.ac.uk%7Cc5f922888565484cca2d08db04689db2%7C912a5d77fb984eeeaf321334d8f04a53%7C0%7C0%7C638108619586622927%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000%7C%7C%7C&sdata=lIiEnoR5jjVDjN6nbhJqucfgF1b8zQyyxo0%2F9ryH%2F%2Bw%3D&reserved=0">
https://eur03.safelinks.protection.outlook.com/?url=http%3A%2F%2Fmailman.mitgcm.org%2Fmailman%2Flistinfo%2Fmitgcm-support&data=05%7C01%7Cp.burns2%40exeter.ac.uk%7Cc5f922888565484cca2d08db04689db2%7C912a5d77fb984eeeaf321334d8f04a53%7C0%7C0%7C638108619586622927%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000%7C%7C%7C&sdata=lIiEnoR5jjVDjN6nbhJqucfgF1b8zQyyxo0%2F9ryH%2F%2Bw%3D&reserved=0</a><br>
<br>
_______________________________________________<br>
MITgcm-support mailing list<br>
MITgcm-support@mitgcm.org<br>
<a href="https://eur03.safelinks.protection.outlook.com/?url=http%3A%2F%2Fmailman.mitgcm.org%2Fmailman%2Flistinfo%2Fmitgcm-support&data=05%7C01%7Cp.burns2%40exeter.ac.uk%7Cc5f922888565484cca2d08db04689db2%7C912a5d77fb984eeeaf321334d8f04a53%7C0%7C0%7C638108619586622927%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000%7C%7C%7C&sdata=lIiEnoR5jjVDjN6nbhJqucfgF1b8zQyyxo0%2F9ryH%2F%2Bw%3D&reserved=0">https://eur03.safelinks.protection.outlook.com/?url=http%3A%2F%2Fmailman.mitgcm.org%2Fmailman%2Flistinfo%2Fmitgcm-support&data=05%7C01%7Cp.burns2%40exeter.ac.uk%7Cc5f922888565484cca2d08db04689db2%7C912a5d77fb984eeeaf321334d8f04a53%7C0%7C0%7C638108619586622927%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000%7C%7C%7C&sdata=lIiEnoR5jjVDjN6nbhJqucfgF1b8zQyyxo0%2F9ryH%2F%2Bw%3D&reserved=0</a><br>
</div>
</span></font></div>
</body>
</html>