<div dir="ltr"><div>Dear FLASH Users,<br><br></div>Since I can't submit to the flash-bugs list, and there's no public bug tracking system, I'm submitting my bug report here. If there's a better forum for such communications, by all means let me know.<br>
<br>I've come across a bug in the Ionize unit 
which prevents an update to the internal energy due to evolution of the 
population fractions. Specifically, the file in question is 
[FLASH4.2.2]:<br>
<br>source/physics/sourceTerms/Ionize/IonizeMain/Ionize.F90 .<br><div><br></div><div>The
 lines in question are lines 128 - 140; an update to the internal energy
 is calculated and stored in the local variable 'ei', but never 
propagated to the 'solnData' array. When the equation of state is 
invoked later on (in MODE_DENS_EI mode, no less), it uses the old internal 
energy as input. The updated internal energy from the Ionize unit is thus thrown away.<br>
<br>To fix this bug, the internal energy should be updated in the 
solnData array, along with the total energy (the following 2 lines 
replace line 137 in IonizeMain/Ionize.F90):<br><br>solnData(EINT_VAR,i,j,k) = ei<br>solnData(ENER_VAR,i,j,k) = ei + ek<br>
</div><br>I do not have any data on the effects of this bug, but 
it does seem to have existed since at least FLASH 2.5. Please let me 
know if you need more information.<br><br>Best Regards,<br><br>Jason Galyardt, PhD<br>
Dept. of Physics and Astronomy<br>University of Georgia<br></div>