<div dir="ltr"><div dir="ltr"><div>A more direct way of doing this may be to simply 
disallow changes to the internal energy in the hydro update. If for 
example you're using the split PPM solver, you can comment out<br><br>              solnData(EINT_VAR,i,j,k) = einternal<br><br></div>in source/physics/Hydro/HydroMain/split/PPM/hy_ppm_updateSoln.F90.<br></div><br clear="all"><div class="gmail_extra"><div><div class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr">Max Katz<br>Ph.D. Candidate, Department of Physics and Astronomy<br>Stony Brook University<br><a href="http://astro.sunysb.edu/mkatz/" target="_blank">http://astro.sunysb.edu/mkatz/</a><br></div></div></div></div></div></div></div>
<br><div class="gmail_quote">On Wed, Jun 10, 2015 at 1:20 PM, Andrea Derdzinski <span dir="ltr"><<a href="mailto:aderdzinski@astro.columbia.edu" target="_blank">aderdzinski@astro.columbia.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>Hi FLASH users,<br><br> I am running a simulation in which I want to hold the internal energy fixed. The internal energy has a radial dependence, so I need to read in the grid coordinate information at every time step and set the internal energy to its initial profile eint(r).   My idea was to modify the EOS unit (using the ideal gamma EOS, so changing eos_idealGamma.F90) such that it reads in the grid coordinate information and resets the internal energy to it's initial state. However, this module does not receive the block information as an input, so I am unsure of the best way to retrieve the coordinate information. <br>Does anyone have any suggestions for how to do this?<br><br><br></div><div>Thanks,<br></div><div>-Andrea<br></div><br><br clear="all"><div><div><br></div></div></div>
</blockquote></div><br></div></div>