<div dir="ltr">Hi Ou,<div>Thanks for the reply. Can you tell me more about the type of budget that you closed, e.g., the diagnostic (mass, heat, salt, etc) and the layer/geographical area. The reason that I am asking is that a budget over a single grid cell may not need 64-bit precision, but when integrating over depth, etc, the errors may accumulate. I am currently looking at budgets over the top 1000m for mass, heat, salt. </div><div><br></div><div>A very long time ago, I did a budget calculation using MITgcm output that did not close due to 32-bit output, and I changed to 64-bit output and it closed (since then I have always used 64-bit output so I haven't verified that this is necessary). </div><div><br></div><div>thanks,</div><div>Martha</div><div><br></div><div><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Thu, Nov 2, 2017 at 12:03 AM, Wang, Ou (329B) <span dir="ltr"><<a href="mailto:Ou.Wang@jpl.nasa.gov" target="_blank">Ou.Wang@jpl.nasa.gov</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div style="word-wrap:break-word;color:rgb(0,0,0);font-size:14px;font-family:Calibri,sans-serif">
<div>Hi Martha,</div>
<div><br>
</div>
<div>I have closed the heat budget for the global ocean using the nctiles output available on
<a href="ftp://ecco.jpl.nasa.gov/Version4/Release3" target="_blank">ftp://ecco.jpl.nasa.gov/<wbr>Version4/Release3</a>/. The ratio of STD(tendency-convergences-<wbr>forcing)/STD(tendency) is ~1e-7. </div>
<div><br>
</div>
<div>Although the netCDF files are all in 64-bit precision, some of the original binary diagnostic outputs that were converted to the tiled netCDF files are in 32-bit precision. So apparently it is not necessary for the files to be in 64-bit precision to close
budgets. </div>
<div><br>
</div>
<div>Release 3 has geothermal flux. See Section 4.2.1 (page 10) of Chris Piecuch’s
<a href="ftp://ecco.jpl.nasa.gov/Version4/Release3/doc/evaluating_budgets_in_eccov4r3.pdf" target="_blank">
note</a> on budget calculation. </div>
<div><br>
</div>
<div>Regards,</div>
<div>Ou Wang</div>
<div><br>
</div>
<span id="m_3189574810513820701OLK_SRC_BODY_SECTION">
<div style="font-family:Calibri;font-size:11pt;text-align:left;color:black;BORDER-BOTTOM:medium none;BORDER-LEFT:medium none;PADDING-BOTTOM:0in;PADDING-LEFT:0in;PADDING-RIGHT:0in;BORDER-TOP:#b5c4df 1pt solid;BORDER-RIGHT:medium none;PADDING-TOP:3pt">
<span style="font-weight:bold">From: </span><<a href="mailto:ecco-support-bounces@mit.edu" target="_blank">ecco-support-bounces@mit.edu</a>> on behalf of Martha Buckley <<a href="mailto:marthabuckley@gmail.com" target="_blank">marthabuckley@gmail.com</a>><br>
<span style="font-weight:bold">Reply-To: </span>"ECCO support list, wider membership" <<a href="mailto:ecco-support@mit.edu" target="_blank">ecco-support@mit.edu</a>><br>
<span style="font-weight:bold">Date: </span>Tuesday, October 31, 2017 at 7:35 AM<br>
<span style="font-weight:bold">To: </span>"ECCO support list, wider membership" <<a href="mailto:ecco-support@mit.edu" target="_blank">ecco-support@mit.edu</a>><br>
<span style="font-weight:bold">Subject: </span>[ecco-support] Budgets in ECv4, release 3<br>
</div><div><div class="h5">
<div><br>
</div>
<div>
<div>
<div dir="ltr">Hi all,
<div>I am trying to close mass, heat, and salt budgets in ECCo v4 release 3. I closed the budgets to machine precision ((tendancy - ocean convergences - surfaceFluxes)/tendancy ~ 1e-7)) for an older iteration. However, I unable to close the budgets to machine
precision for the files on the website for release3 (the budgets visually look fine though). </div>
<div><br>
</div>
<div>I'm trying to figure out if I have a bug in my code or if something changed with the new output. A few questions:</div>
<div>(1) Are the files 64-bit precision, which is needed to close budgets? This seems to be the case since ncdump says the files are "double". </div>
<div>(2) Did anything change between "r2/iter3" and release3 that might explain this? </div>
<div>(3) Has anyone closed budgets for release3 using the nctiles output available on the ECCO website?</div>
<div>thanks,</div>
<div>Martha</div>
<div><br>
</div>
<div><br>
</div>
<div><br clear="all">
<div><br>
</div>
-- <br>
<div class="m_3189574810513820701gmail_signature" data-smartmail="gmail_signature">Martha W. Buckley<br>
<a href="mailto:marthabuckley@gmail.com" target="_blank">marthabuckley@gmail.com</a><br>
<a href="http://sites.google.com/site/marthabuckley/home" target="_blank">http://sites.google.com/site/<wbr>marthabuckley/home</a><br>
</div>
</div>
</div>
</div>
</div>
</div></div></span>
</div>
<br>______________________________<wbr>_________________<br>
ecco-support mailing list<br>
<a href="mailto:ecco-support@mit.edu">ecco-support@mit.edu</a><br>
<a href="http://mailman.mit.edu/mailman/listinfo/ecco-support" rel="noreferrer" target="_blank">http://mailman.mit.edu/<wbr>mailman/listinfo/ecco-support</a><br>
<br></blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature" data-smartmail="gmail_signature">Martha W. Buckley<br><a href="mailto:marthabuckley@gmail.com" target="_blank">marthabuckley@gmail.com</a><br><a href="http://sites.google.com/site/marthabuckley/home" target="_blank">http://sites.google.com/site/marthabuckley/home</a><br></div>
</div>