<div dir="ltr">Hi Ou,<div>It may be possible to close the budgets for the global ocean with single precision output, but this is not the case for regional budgets. For these we really do need to have double precision outputs. Attached are a couple of plots that show this. </div><div><br></div><div>I am looking at mass and heat budgets over the top 1000m. What I did was take the r3/iter3 output, which has double precision, and after loading each file (e.g., UVELMASS, etc.), I converted to single precision (the back to double so the gcmfaces programs work). When I use the double precision output, (tend-ocean convergence - forcing)./ tend < 1e-7 for both mass and heat (see budgetM1000m_errorDouble.eps and budgetHC1000m_errorDouble.eps). However, when the outputs are converted to single precision the errors are much, much larger, (see budgetM1000m_errorSingle.eps and budgetHC1000m_errorSingle.eps). </div><div><br></div><div>One of the touted advantages of ECCO is the ability to do budgets, so we need outputs that allow us to do budget over any arbitrary regions, not just simple ones, like the global average. So can we get the new outputs in double precision please?! </div><div>thanks,</div><div>Martha</div><div><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Thu, Nov 2, 2017 at 1:54 PM, 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 am also able to close the heat budget for the top ~1000m of the global ocean. </div>
<div><br>
</div>
<div>Regards,</div>
<div>Ou</div>
<div><br>
</div>
<span id="m_7778193776502192640OLK_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 class="">
<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><span style="font-weight:bold">Date: </span>Thursday, November 2, 2017 at 10:11 AM<div><div class="h5"><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>Re: [ecco-support] Budgets in ECv4, release 3<br>
</div></div></div><div><div class="h5">
<div><br>
</div>
<div>
<div>
<div dir="ltr">Thanks! Have you done any budgets for the upper-ocean? For vertically integrated budgets the vertical convergences don't play a role.
<div>Martha</div>
<div><br>
</div>
</div>
<div class="gmail_extra"><br>
<div class="gmail_quote">On Thu, Nov 2, 2017 at 1:07 PM, 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 closed the heat budget for the global ocean from top to bottom.</div>
<div><br>
</div>
<div>Regards,</div>
<div>Ou</div>
<div><br>
</div>
<span id="m_7778193776502192640m_5185235751035004130OLK_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><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><span style="font-weight:bold">Date: </span>Thursday, November 2, 2017 at 10:04 AM<span><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><span style="font-weight:bold">Subject: </span>Re: [ecco-support] Budgets in ECv4, release 3<br>
</div>
<div>
<div class="m_7778193776502192640h5">
<div><br>
</div>
<div>
<div>
<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/Versio<wbr>n4/Release3</a>/. The ratio of STD(tendency-convergences-forc<wbr>ing)/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_7778193776502192640m_5185235751035004130m_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="m_7778193776502192640m_5185235751035004130h5">
<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_7778193776502192640m_5185235751035004130m_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/m<wbr>arthabuckley/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" target="_blank">ecco-support@mit.edu</a><br>
<a href="http://mailman.mit.edu/mailman/listinfo/ecco-support" rel="noreferrer" target="_blank">http://mailman.mit.edu/mailman<wbr>/listinfo/ecco-support</a><br>
<br>
</blockquote>
</div>
<br>
<br clear="all">
<div><br>
</div>
-- <br>
<div class="m_7778193776502192640m_5185235751035004130gmail_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/m<wbr>arthabuckley/home</a><br>
</div>
</div>
</div>
</div>
</div>
</div>
</span></div>
<br>
______________________________<wbr>_________________<br>
ecco-support mailing list<br>
<a href="mailto:ecco-support@mit.edu" target="_blank">ecco-support@mit.edu</a><br>
<a href="http://mailman.mit.edu/mailman/listinfo/ecco-support" rel="noreferrer" target="_blank">http://mailman.mit.edu/mailman<wbr>/listinfo/ecco-support</a><br>
<br>
</blockquote>
</div>
<br>
<br clear="all">
<div><br>
</div>
-- <br>
<div class="m_7778193776502192640gmail_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></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>