<div dir="ltr">You are right about the patterns being non-random. I am glad that you are convinced! So when will the double precision output be available online? :)<div>thanks,</div><div>Martha</div><div><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Mon, Nov 6, 2017 at 12:08 PM, Menemenlis, Dimitris (329C) <span dir="ltr"><<a href="mailto:Dimitris.Menemenlis@jpl.nasa.gov" target="_blank">Dimitris.Menemenlis@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">
Hi Martha, the non-random patterns in your figures
<div>convince me that double precision is needed.
<div>
<div><br>
</div>
<div>Cheers, Dimitris<span class=""><br>
<div><br>
</div>
<div>
<div>
<blockquote type="cite">
<div>On Nov 6, 2017, at 8:56 AM, Martha Buckley <<a href="mailto:marthabuckley@gmail.com" target="_blank">marthabuckley@gmail.com</a>> wrote:</div>
<br class="m_7171254137721343140Apple-interchange-newline">
<div><span style="font-family:Helvetica;font-size:12px;font-style:normal;font-variant-caps:normal;font-weight:normal;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;float:none;display:inline!important">I
didn't chose the most instructive colorbar to show the magnitude of the errors. Attached are plots with colorbars to indicate the magnitude of the errors. In particular, note that the errors due to single precision in the mass budget are on the order of 10%
(saturated on the 1% colorbar). Note these are the errors for a particular time step based on monthly output. </span>
<div style="font-family:Helvetica;font-size:12px;font-style:normal;font-variant-caps:normal;font-weight:normal;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px">
<br>
</div>
<div style="font-family:Helvetica;font-size:12px;font-style:normal;font-variant-caps:normal;font-weight:normal;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px">
The issue with budgets with single precision output is that the tendancy is always a small residual of larger terms. While of course the uncertainty of ECCO estimates is larger than even these 10% errors, I think this is not particularly relevant to this discussion.
ECCO in the end is also a forward model and MITgcm developers have take great pains to insure that budgets close. I don't know why we'd sacrifice the ability to close budgets rigorously in order to save a factor of 2 in space by storing as single precision.
Also, for the user it is a lot easier if the budgets close to machine precision rather than having a user wonder if an error is "acceptable" or not. </div>
<div style="font-family:Helvetica;font-size:12px;font-style:normal;font-variant-caps:normal;font-weight:normal;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px">
<br>
</div>
<div style="font-family:Helvetica;font-size:12px;font-style:normal;font-variant-caps:normal;font-weight:normal;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px">
Martha</div>
</div>
</blockquote>
<br>
</div>
<br>
</div>
</span></div>
</div>
</div>
</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>