<div dir="ltr"><div class="gmail_default" style="font-family:arial,helvetica,sans-serif;font-size:small">Hi Klaus,</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">I've decided to give this another go, especially because since the time of this email exchange I have altered IO_writeIntegralQuantities and Simulation_adjustEvolution successfully so I have a better idea how it should work.  </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">I have now tried it and it sort of works.  That is, I can calculate the total explosion energy, which in my case is just the thermal energy above the ambient thermal energy, and a correction factor.  The correction factor should be multiplied by the pressure (or internal energy) of the parcels that have a pressure above the ambient pressure.  However, when I do that it seems like nothing is changed when it returns. I print out the correction factor and it seems right.  It's odd because I have previously been able to use Simulation_adjustEvolution to set off a new explosion.  Would you be interested in taking a look at the code?  It's pretty simple and maybe you could see the problem quickly.</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">Jon</div><div class="gmail_extra"><br><div class="gmail_quote">On Mon, Aug 15, 2016 at 5:14 PM, Klaus Weide <span dir="ltr"><<a href="mailto:klaus@flash.uchicago.edu" target="_blank">klaus@flash.uchicago.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><span>On Wed, 27 Jul 2016, Slavin, Jonathan wrote:<br>
<br>
> An alternative approach that I thought of that could be applied to the<br>
> thermal energy only case is to, after all blocks have been initialized, 1)<br>
> total up the thermal energy, 2) calculate the ratio  sim_expEnergy/E_total<br>
> and 3) adjust the pressure in all the zones in which the pressure is above<br>
> the ambient pressure by multiplying by that factor.<br>
><br>
> While I think I could write the code to calculate the ratio and apply it, I<br>
> don't fully understand the structure of FLASH, so I'm not sure where it<br>
> would need to be called from.  One advantage to this approach is that one<br>
> is pretty much guaranteed that the explosion energy will be what is desired.<br>
<br>
</span>You would need to introduce an additional global loop over all blocks<br>
to do this - or actually, two: a first one for accumulating the current<br>
energy values from all cells of all (leaf) blocks, and a second one for<br>
rescaling them.<br>
<br>
This should be done logically at the end of Driver_initFlash, or the<br>
beginning of Driver_evolveFlash. A convenient place to put something like<br>
this could be your own implementation of Simulation_adjustEvolution; the<br>
call already exists in Driver_evolveFlash.F90. You'd have to make<br>
sure you do the rescaling only in the very first call, for example, by<br>
checking the nstep argument.<br>
<span class="m_-9023627975963660386HOEnZb"><font color="#888888"><br>
Klaus<br>
</font></span></blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="m_-9023627975963660386gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr">______________________________<wbr>__________________________<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:(617)%20496-7981" value="+16174967981" target="_blank">(617) 496-7981</a>       Cambridge, MA 02138-1516<br>cell: <a href="tel:(781)%20363-0035" value="+17813630035" target="_blank">(781) 363-0035</a>             USA<br>______________________________<wbr>__________________________<br><br></div></div></div></div>
</div></div>