<html><head><meta http-equiv="content-type" content="text/html; charset=utf-8"></head><body dir="auto">Thanks!<br><br><div dir="ltr">Christopher<div><br></div></div><div dir="ltr"><br><blockquote type="cite">On Jun 4, 2021, at 2:41 PM, Wang, Ou (US 329B) <ou.wang@jpl.nasa.gov> wrote:<br><br></blockquote></div><blockquote type="cite"><div dir="ltr">
<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>
<div class="WordSection1">
<p class="MsoNormal">Hi Christopher,<o:p></o:p></p>
<p class="MsoNormal"><o:p> </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> </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> </o:p></span></p>
<p class="MsoNormal"><span style="color:black">Ou</span><o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"><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"><ecco-support-bounces@mit.edu> on behalf of "Wang, Ou (US 329B)" <ou.wang@jpl.nasa.gov><br>
<b>Reply-To: </b>"ECCO support list, wider membership" <ecco-support@mit.edu><br>
<b>Date: </b>Friday, June 4, 2021 at 9:17 AM<br>
<b>To: </b>"ECCO support list, wider membership" <ecco-support@mit.edu><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> </o:p></p>
</div>
<p class="MsoNormal">Hi Christopher,<o:p></o:p></p>
<p class="MsoNormal"> <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"> </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"> <o:p></o:p></p>
<p class="MsoNormal">Ou<o:p></o:p></p>
<p class="MsoNormal"> <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"><ecco-support-bounces@mit.edu> on behalf of Christopher Wolfe <christopher.wolfe@stonybrook.edu><br>
<b>Reply-To: </b>"ECCO support list, wider membership" <ecco-support@mit.edu><br>
<b>Date: </b>Friday, June 4, 2021 at 6:44 AM<br>
<b>To: </b>"ecco-support@mit.edu" <ecco-support@mit.edu><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"> <o:p></o:p></p>
</div>
<p class="MsoNormal">Hi all: <o:p></o:p></p>
<div>
<p class="MsoNormal"> <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. <o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> <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 total_kapredi_r009bit11.bin is the "adjusted GM redi parameter (first guess
+ adjustments)”, which implies that total_diffkr_r009bit11.bin should be the adjusted diffusivity. <o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> <o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">On the other hand, the README at <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> says that <span style="color:black"> total_diffkr_r009bit11.bin
is the </span>"vert. diff. of release 1 (this field plus xx is the total)”, which means that I should add xx_diffkr.0000000129.data to this to get the full diffusivity. However, the values in this file are order 1, so adding it to <span style="color:black">total_diffkr_r009bit11
would produce 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 non-dimensional and need to be scaled by one over the square root of the weights, which are in 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 Sciences<br>
Stony Brook University<br>
</span><a href="mailto:christopher.wolfe@stonybrook.edu">christopher.wolfe@stonybrook.edu</a><span style="color:black"> 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>
<span>_______________________________________________</span><br><span>ecco-support mailing list</span><br><span>ecco-support@mit.edu</span><br><span>http://mailman.mit.edu/mailman/listinfo/ecco-support</span><br></div></blockquote></body></html>