<div dir="ltr">Hi Malia :<div><br></div><div>  When you say the initial EOS data, I presume you mean that you've verified the t = 0 data to be good. Helmholtz is fairly robust on good data, so I suspect if you inspected the data passed to Helmholtz during this call, you will find it to be already corrupted. </div><div><br></div><div>  To add to Dean's comment, in my experience as well, when the code dies in Helmholtz, it is usually the symptom and not the underlying disease. The 'bad' data which manifests itself in this way is usually symptomatic of a more fundamental problem elsewhere in the code. For instance, you may be encountering a numerical instability in the hydrodynamics, or the nuclear burning, but because these modules do not trap NaNs or Infs, they will happily propagate to Helmholtz, where an exception is thrown.   You will need to root out the underlying cause, because the NaNs and Infs will progressively corrupt the simulation.</div><div><br></div><div>  Are you using nuclear burning in this simulation?</div><div><br></div><div>  Best wishes,</div><div><br></div><div>  Bob</div></div><div class="gmail_extra"><br><div class="gmail_quote">On Wed, May 3, 2017 at 5:02 PM, Jenks, Malia T. <span dir="ltr"><<a href="mailto:mjenks@ou.edu" target="_blank">mjenks@ou.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">




<div dir="ltr">
<div id="m_1881886895311193561divtagdefaultwrapper" style="font-size:12pt;color:#000000;font-family:Calibri,Arial,Helvetica,sans-serif" dir="ltr">
<p></p>
<div>
<div>
<div>
<div>I am getting the following error from the helmholtz eos about negative pressure:<br>
<br>
 [eos_helm] Negative total pressure.<br>
   values: dens,temp:   2.321442471007820E-004   1549.01220000000               <br>
            <br>
   values: abar,zbar:                 Infinity                     NaN          <br>
            <br>
   coulomb coupling parameter Gamma:                      NaN                   <br>
            <br>
  Prad   1.45198E-02 Pion  0.00000E+00 Pele           NaN Pcoul  0.00000E+00    <br>
            <br>
  Ptot            NaN df_d    8.31451E+21                   <wbr>                   
<br>
<br>
<br>
</div>
I know that the initial mass fraction supplied to the eos is correct. The initial call of the eos sets abar and zbar correctly. I am not sure what is causing them to error. I also don't know why I am getting a NaN for the electron Pressure. Could anyone help
 point me to what is going wrong and how I can solve it?<br>
<br>
</div>
Thanks,<br>
</div>
<br>
Malia Jenks</div>
<br>
<p></p>
</div>
</div>

</blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr">Dr. Robert Fisher<br>Associate Professor / Graduate Program Director<br>University of Massachusetts/Dartmouth<br>Department of Physics<br>285 Old Westport Road<br>North Dartmouth, Massachusetts 02747<br><a href="mailto:robert.fisher@umassd.edu" target="_blank">robert.fisher@umassd.edu</a><br><a href="http://www.novastella.org" target="_blank">http://www.novastella.org</a></div></div>
</div>