<html xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40">
<head>
<meta http-equiv="Content-Type" content="text/html; charset=utf-8">
<meta name="Generator" content="Microsoft Word 15 (filtered medium)">
<style><!--
/* Font Definitions */
@font-face
        {font-family:"Cambria Math";
        panose-1:2 4 5 3 5 4 6 3 2 4;}
@font-face
        {font-family:Calibri;
        panose-1:2 15 5 2 2 2 4 3 2 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
        {margin:0in;
        font-size:11.0pt;
        font-family:"Calibri",sans-serif;}
a:link, span.MsoHyperlink
        {mso-style-priority:99;
        color:blue;
        text-decoration:underline;}
span.EmailStyle19
        {mso-style-type:personal-reply;
        font-family:"Calibri",sans-serif;
        color:windowtext;}
.MsoChpDefault
        {mso-style-type:export-only;
        font-size:10.0pt;}
@page WordSection1
        {size:8.5in 11.0in;
        margin:1.0in 1.0in 1.0in 1.0in;}
div.WordSection1
        {page:WordSection1;}
--></style>
</head>
<body lang="EN-US" link="blue" vlink="purple" style="word-wrap:break-word">
<div class="WordSection1">
<p class="MsoNormal">Hi Christopher,<o:p></o:p></p>
<p class="MsoNormal"><o:p>&nbsp;</o:p></p>
<p class="MsoNormal">One correction to the alternative way I mentioned;<o:p></o:p></p>
<p class="MsoNormal">Alternatively, you can calculate it by summing up <b>the first guess of a constant 1e-5 m^2/s (not
<span style="color:black">total_diffkr_r009bit11.bin)</span></b><span style="color:black"> and xx_diffkr.dim.0000000129.data.<o:p></o:p></span></p>
<p class="MsoNormal"><span style="color:black"><o:p>&nbsp;</o:p></span></p>
<p class="MsoNormal"><span style="color:black">See Forget et al. (2015; doi:10.5194/gmd-8-3071-2015), Table 1 for the first guess of the mixing parameters.<o:p></o:p></span></p>
<p class="MsoNormal"><span style="color:black"><o:p>&nbsp;</o:p></span></p>
<p class="MsoNormal"><span style="color:black">Ou</span><o:p></o:p></p>
<p class="MsoNormal"><o:p>&nbsp;</o:p></p>
<p class="MsoNormal"><o:p>&nbsp;</o:p></p>
<div style="border:none;border-top:solid #B5C4DF 1.0pt;padding:3.0pt 0in 0in 0in">
<p class="MsoNormal"><b><span style="font-size:12.0pt;color:black">From: </span></b><span style="font-size:12.0pt;color:black">&lt;ecco-support-bounces@mit.edu&gt; on behalf of &quot;Wang, Ou (US 329B)&quot; &lt;ou.wang@jpl.nasa.gov&gt;<br>
<b>Reply-To: </b>&quot;ECCO support list, wider membership&quot; &lt;ecco-support@mit.edu&gt;<br>
<b>Date: </b>Friday, June 4, 2021 at 9:17 AM<br>
<b>To: </b>&quot;ECCO support list, wider membership&quot; &lt;ecco-support@mit.edu&gt;<br>
<b>Subject: </b>Re: [ecco-support] [EXTERNAL] How to get the total diffusivity in ECCOv4<o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal"><o:p>&nbsp;</o:p></p>
</div>
<p class="MsoNormal">Hi Christopher,<o:p></o:p></p>
<p class="MsoNormal">&nbsp;<o:p></o:p></p>
<p class="MsoNormal">PO.DAAC is now hosting V4r4’s files. You can get V4r4’s dimensional vertical diffusivity field (first guess plus the control adjustment) from PO.DAAC at
<a href="https://urldefense.us/v3/__https:/archive.podaac.earthdata.nasa.gov/podaac-ops-cumulus-protected/ECCO_L4_OCEAN_3D_MIX_COEFFS_LLC0090GRID_V4R4/OCEAN_3D_MIXING_COEFFS_ECCO_V4r4_native_llc0090.nc__;!!PvBDto6Hs4WbVuu7!bO3-ZfkkeOEujVIhL2nn3AEb9q1Zl86CHb1ZTrYDdah4HPKuSZ_F2v6ZPAv3JA9Giy4$">
https://archive.podaac.earthdata.nasa.gov/podaac-ops-cumulus-protected/ECCO_L4_OCEAN_3D_MIX_COEFFS_LLC0090GRID_V4R4/OCEAN_3D_MIXING_COEFFS_ECCO_V4r4_native_llc0090.nc</a>. The variable is called “DIFFKR” in the NetCDF file. Alternatively, you can calculate
 it by summing up <span style="color:black">total_diffkr_r009bit11.bin and xx_diffkr.dim.0000000129.data that are available on ECCO Drive at https://ecco.jpl.nasa.gov/drive/files/Version4/Release4/input_init/total_diffkr_r009bit11.bin and
<a href="https://ecco.jpl.nasa.gov/drive/files/Version4/Release4/other/adjustments/xx_diffkr.dim.0000000129.data">
https://ecco.jpl.nasa.gov/drive/files/Version4/Release4/other/adjustments/xx_diffkr.dim.0000000129.data</a>.
</span><o:p></o:p></p>
<p class="MsoNormal"><span style="color:black">&nbsp;</span><o:p></o:p></p>
<p class="MsoNormal"><span style="color:black">The description about </span>total_kapredi_r009bit11.bin in the ECCOv4 Python tutorial is incorrect and will be corrected. The field in total_kapredi_r009bit11.bin is the Redi parameter from Release
<span style="color:red">1</span>, not the adjusted Redi parameter from Release 4.
<o:p></o:p></p>
<p class="MsoNormal">&nbsp;<o:p></o:p></p>
<p class="MsoNormal">Ou<o:p></o:p></p>
<p class="MsoNormal">&nbsp;<o:p></o:p></p>
<div style="border:none;border-top:solid #B5C4DF 1.0pt;padding:3.0pt 0in 0in 0in">
<p class="MsoNormal"><b><span style="font-size:12.0pt;color:black">From: </span></b><span style="font-size:12.0pt;color:black">&lt;ecco-support-bounces@mit.edu&gt; on behalf of Christopher Wolfe &lt;christopher.wolfe@stonybrook.edu&gt;<br>
<b>Reply-To: </b>&quot;ECCO support list, wider membership&quot; &lt;ecco-support@mit.edu&gt;<br>
<b>Date: </b>Friday, June 4, 2021 at 6:44 AM<br>
<b>To: </b>&quot;ecco-support@mit.edu&quot; &lt;ecco-support@mit.edu&gt;<br>
<b>Subject: </b>[EXTERNAL] [ecco-support] How to get the total diffusivity in ECCOv4</span><o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">&nbsp;<o:p></o:p></p>
</div>
<p class="MsoNormal">Hi all: <o:p></o:p></p>
<div>
<p class="MsoNormal">&nbsp;<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">I’m trying to figure out what the vertical diffusivity field in ECCOv4r4 is, but I’m a little confused about which input file actually contains this data.&nbsp;<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">&nbsp;<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">The ECCOv4 python tutorial (at <a href="https://urldefense.us/v3/__https:/ecco-v4-python-tutorial.readthedocs.io/ECCO_v4_Loading_LLC_compact_binary_files.html*Example-3:-Load-a-3D-'compact'-llc-binary-file-with-3rd-dimension-=-Depth__;Iw!!PvBDto6Hs4WbVuu7!eGOlWbFIVdKaXoIQQEcSSqRPb2UW2kabBqJXw8Hdg5AizGFceCu2pSW1rdXfRM0vvIA$">
https://ecco-v4-python-tutorial.readthedocs.io/ECCO_v4_Loading_LLC_compact_binary_files.html#Example-3:-Load-a-3D-'compact'-llc-binary-file-with-3rd-dimension-=-Depth</a>) says that&nbsp;total_kapredi_r009bit11.bin is the &quot;adjusted GM redi parameter (first guess
 &#43; adjustments)”, which implies that total_diffkr_r009bit11.bin should be the adjusted diffusivity.&nbsp;<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">&nbsp;<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">On the other hand, the README at&nbsp;<a href="https://ecco.jpl.nasa.gov/drive/files/Version4/Release4/input_init/README">https://ecco.jpl.nasa.gov/drive/files/Version4/Release4/input_init/README</a>&nbsp;says that&nbsp;<span style="color:black">&nbsp;total_diffkr_r009bit11.bin
 is the&nbsp;</span>&quot;vert. diff. of release 1 (this field plus xx is the total)”, which means that I should add&nbsp;xx_diffkr.0000000129.data to this to get the full diffusivity. However, the values in this file are order 1, so adding it to&nbsp;<span style="color:black">total_diffkr_r009bit11
 would produce&nbsp;diffusivities far outside the expected range.</span><o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><span style="color:black"><br>
<br>
<br>
</span><o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><span style="color:black">The gencost documentation of the ctrl package (<a href="https://urldefense.us/v3/__https:/mitgcm.readthedocs.io/en/latest/ocean_state_est/ocean_state_est.html*ctrl-model-parameter-adjustment-capability__;Iw!!PvBDto6Hs4WbVuu7!eGOlWbFIVdKaXoIQQEcSSqRPb2UW2kabBqJXw8Hdg5AizGFceCu2pSW1rdXfrjEs71w$">https://mitgcm.readthedocs.io/en/latest/ocean_state_est/ocean_state_est.html#ctrl-model-parameter-adjustment-capability</a>)
 says that the xx adjustments are&nbsp;non-dimensional and need to be scaled by one over the square root of the weights, which are in&nbsp;r2.wdiffkrFldv2.data. Using this produces adjustments on the order of 1e-6 to 1e-5, which seems about right. However, the data.ctrl
 namelist also specifies WC01 smoothing and some thresholding. Do these need to be applied as well?</span><o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><span style="color:black"><br>
<br>
<br>
</span><o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><span style="color:black">Is the best way to get this data simply to run the model for a bit and have it dump diffkr? It would be convenient to be able to calculate these quantities without having to run the model, however.</span><o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><span style="color:black"><br>
<br>
<br>
</span><o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><span style="color:black">Thanks in advance for your advise,</span><o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><span style="color:black">Christopher</span><o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><span style="color:black"><br>
<br>
<br>
</span><o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><span style="color:black">----------------------------------------------------------------------------<br>
Christopher L. Pitt Wolfe<br>
Associate Professor (Physical Oceanography)<br>
School of Marine and Atmospheric&nbsp;Sciences<br>
Stony Brook University<br>
</span><a href="mailto:christopher.wolfe@stonybrook.edu">christopher.wolfe@stonybrook.edu</a><span style="color:black">&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;631-632-3152<br>
----------------------------------------------------------------------------</span><o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><span style="color:black"><br>
<br>
<br>
</span><o:p></o:p></p>
</div>
</div>
</body>
</html>