<div dir="ltr"><div class="gmail_default" style="font-family:arial,helvetica,sans-serif;font-size:small">Hi Jonathan,<br><br></div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif;font-size:small">I have also noticed this problem.<br><br></div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif;font-size:small">The energy one sees in the .dat file is typically less than the input energy for sim_expEnergy. Also noticed that this is geometry dependent.<br><br></div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif;font-size:small">I input on the order of 1e10 ergs, and in cartesian I get out 9.2e9; for cylindrical I get even less! 3.2e8... which is huge. Note I am going based off of memory here. I cannot get into my HPC system now.<br><br></div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif;font-size:small">If I am remembering things correctly: seems that in the source code's volume function in Simulation_init.f90 doesn't always choose the volume you're choosing to distribute the energy over.<br><br></div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif;font-size:small">I think it is the section where you'll find a comment that notes the code will calculate the initial volume and interior pressure. It is assuming pi r^2 for to get the energy density... where as you may want something different for other chosen geometries.<br><br></div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif;font-size:small">If I remember correctly the output energy also changes depending where you put the hotspot.<br><br></div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif;font-size:small">Sorry for the haphazard e-mail, but I am experiencing similar issues!<br></div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif;font-size:small"><br></div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif;font-size:small">Best,<br></div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif;font-size:small">Marissa<br></div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif;font-size:small"><br><br> <br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Tue, Jul 26, 2016 at 10:30 AM, Slavin, Jonathan <span dir="ltr"><<a href="mailto:jslavin@cfa.harvard.edu" target="_blank">jslavin@cfa.harvard.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 class="gmail_default" style="font-family:arial,helvetica,sans-serif">Hi all,</div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif"><br></div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif">I've been running supernova remnant evolution simulations. The initialization is done with code taken from the supplied Sedov simulation with certain small changes. One issue that I've run into is that the total energy is not very accurate. That is, though the energy is conserved, the total energy, initially input as thermal energy, is below that which is given in the sim_expEnergy parameter by ~5%. Has anyone else run into this? I don't think that it has to do with any of the changes that I've made to Simulation_initBlock or Simulation_init. I'm wondering if increasing the sim_nsubzones parameter might help. Any advice would be appreciated.</div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif"><br></div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif">Thanks,</div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif">Jon</div><span class="HOEnZb"><font color="#888888"><div><br></div>-- <br><div data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr">________________________________________________________<br>Jonathan D. Slavin Harvard-Smithsonian CfA<br><a href="mailto:jslavin@cfa.harvard.edu" target="_blank">jslavin@cfa.harvard.edu</a> 60 Garden Street, MS 83<br>phone: <a href="tel:%28617%29%20496-7981" value="+16174967981" target="_blank">(617) 496-7981</a> Cambridge, MA 02138-1516<br>cell: <a href="tel:%28781%29%20363-0035" value="+17813630035" target="_blank">(781) 363-0035</a> USA<br>________________________________________________________<br><br></div></div></div></div>
</font></span></div>
</blockquote></div><br><br clear="all"><br>-- <br><div class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><span><font size="1"><span style="color:rgb(0,0,0)"><span style="font-family:monospace,monospace"><div>Marissa Adams <span><font size="1"><span style="color:rgb(0,0,0)"><span style="font-family:monospace,monospace">E-mail:</span></span><span style="color:rgb(0,0,0)"><span style="font-family:monospace,monospace"><a> </a><a href="mailto:madams@pas.rochester.edu" target="_blank">madams@pas.rochester.edu</a></span></span></font></span><br>Graduate Student <span><font size="1"><span style="color:rgb(0,0,0)"><span style="font-family:monospace,monospace">Ph: (585) 402-5779</span></span></font></span></div><div>Department of Physics & Astronomy <span><font size="1"><span style="color:rgb(0,0,0)"><span style="font-family:monospace,monospace">Website: <a href="http://www.pas.rochester.edu/%7Emadams/" target="_blank">http://www.pas.rochester.edu/~madams/</a></span></span></font></span></div>University of Rochester<br>478 Bausch & Lomb Hall<br>P.O. Box 270171</span></span></font></span></div></div></div></div></div></div></div></div>
</div>