<html><head><meta http-equiv="Content-Type" content="text/html; charset=utf-8"></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class="">Hi all:<div class=""><br class=""></div><div class="">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;</div><div class=""><br class=""></div><div class="">The ECCOv4 python tutorial (at <a href="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" class="">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 "adjusted GM redi parameter (first guess + adjustments)”, which implies that total_diffkr_r009bit11.bin should be the adjusted diffusivity.&nbsp;</div><div class=""><br class=""></div><div class="">On the other hand, the README at&nbsp;<a href="https://ecco.jpl.nasa.gov/drive/files/Version4/Release4/input_init/README" class="">https://ecco.jpl.nasa.gov/drive/files/Version4/Release4/input_init/README</a>&nbsp;says that&nbsp;<span style="caret-color: rgb(0, 0, 0); color: rgb(0, 0, 0);" class="">&nbsp;</span><span style="caret-color: rgb(0, 0, 0); color: rgb(0, 0, 0);" class="">total_diffkr_r009bit11.bin is the&nbsp;</span>"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;<font color="#000000" class=""><span style="caret-color: rgb(0, 0, 0);" class="">total_diffkr_r009bit11 would produce&nbsp;diffusivities far outside the expected range.</span></font></div><div class=""><font color="#000000" class=""><span style="caret-color: rgb(0, 0, 0);" class=""><br class=""></span></font></div><div class=""><font color="#000000" class="">The gencost documentation of the ctrl package (<a href="https://mitgcm.readthedocs.io/en/latest/ocean_state_est/ocean_state_est.html#ctrl-model-parameter-adjustment-capability" style="caret-color: rgb(0, 0, 0);" class="">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;<span style="caret-color: rgb(0, 0, 0);" class="">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></font></div><div class=""><font color="#000000" class=""><span style="caret-color: rgb(0, 0, 0);" class=""><br class=""></span></font></div><div class=""><font color="#000000" class=""><span style="caret-color: rgb(0, 0, 0);" class="">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></font></div><div class=""><font color="#000000" class=""><span style="caret-color: rgb(0, 0, 0);" class=""><br class=""></span></font></div><div class=""><font color="#000000" class=""><span style="caret-color: rgb(0, 0, 0);" class="">Thanks in advance for your advise,</span></font></div><div class=""><font color="#000000" class=""><span style="caret-color: rgb(0, 0, 0);" class="">Christopher</span></font></div><div class=""><font color="#000000" class=""><span style="caret-color: rgb(0, 0, 0);" class=""><br class=""></span></font></div><div class=""><span style="caret-color: rgb(0, 0, 0); color: rgb(0, 0, 0);" class="">----------------------------------------------------------------------------</span><br class="" style="caret-color: rgb(0, 0, 0); color: rgb(0, 0, 0);"><span style="caret-color: rgb(0, 0, 0); color: rgb(0, 0, 0);" class="">Christopher L. Pitt Wolfe</span><br class="" style="caret-color: rgb(0, 0, 0); color: rgb(0, 0, 0);"><span style="caret-color: rgb(0, 0, 0); color: rgb(0, 0, 0);" class="">Associate Professor (Physical Oceanography)</span><br class="" style="caret-color: rgb(0, 0, 0); color: rgb(0, 0, 0);"><span style="caret-color: rgb(0, 0, 0); color: rgb(0, 0, 0);" class="">School of Marine and Atmospheric&nbsp;Sciences</span><br class="" style="caret-color: rgb(0, 0, 0); color: rgb(0, 0, 0);"><span style="caret-color: rgb(0, 0, 0); color: rgb(0, 0, 0);" class="">Stony Brook University</span><br class="" style="caret-color: rgb(0, 0, 0); color: rgb(0, 0, 0);"><a href="mailto:christopher.wolfe@stonybrook.edu" class="" style="caret-color: rgb(0, 0, 0);">christopher.wolfe@stonybrook.edu</a><span style="caret-color: rgb(0, 0, 0); color: rgb(0, 0, 0);" class="">&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;631-632-3152</span><br class="" style="caret-color: rgb(0, 0, 0); color: rgb(0, 0, 0);"><span style="caret-color: rgb(0, 0, 0); color: rgb(0, 0, 0);" class="">----------------------------------------------------------------------------</span></div><div class=""><span style="caret-color: rgb(0, 0, 0); color: rgb(0, 0, 0);" class=""><br class=""></span></div></body></html>