<html><head></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; ">Hello, <div><br></div><div>I have been getting an error similar to that described last month here</div><div><a href="http://flash.uchicago.edu/pipermail/flash-users/2010-July/003084.html">http://flash.uchicago.edu/pipermail/flash-users/2010-July/003084.html</a></div><div><br></div><div><br></div><div>The error I get is:</div><div><div><font class="Apple-style-span" face="'Andale Mono'" size="3"><span class="Apple-style-span" style="font-size: 11px;">ERROR After calling Eos, eosData(EOS_EINT) or eosData(EOS_DENS) are zero</span></font></div><div><font class="Apple-style-span" face="'Andale Mono'" size="3"><span class="Apple-style-span" style="font-size: 11px;"> iflag= 1 1 1 1 1 1 1 1</span></font></div><div><font class="Apple-style-span" face="'Andale Mono'" size="3"><span class="Apple-style-span" style="font-size: 11px;"> solnData(EINT_VAR,ib:ie,jb,kb)= 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 </span></font></div><div><font class="Apple-style-span" face="'Andale Mono'" size="3"><span class="Apple-style-span" style="font-size: 11px;"> eosData (eint+1:eint+vecLen) = 1.00000000000000004E-010 1.00000000000000004E-010 1.00000000000000004E-010 1.00000000000000004E-010 1.00000000000000004E-010 1.00000000000000004E-010 1.00000000000000004E-010 1.00000000000000004E-010</span></font></div><div><font class="Apple-style-span" face="'Andale Mono'" size="3"><span class="Apple-style-span" style="font-size: 11px;"> solnData(DENS_VAR,ib:ie,jb,kb)= 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 </span></font></div><div><font class="Apple-style-span" face="'Andale Mono'" size="3"><span class="Apple-style-span" style="font-size: 11px;"> eosData(dens+1:dens+vecLen) = 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 </span></font></div><div><font class="Apple-style-span" face="'Andale Mono'" size="3"><span class="Apple-style-span" style="font-size: 11px;"> Perhaps the initialization routine is wrong..... or</span></font></div><div><font class="Apple-style-span" face="'Andale Mono'" size="3"><span class="Apple-style-span" style="font-size: 11px;"> perhaps the runtime parameter eosMode is wrong.</span></font></div><div><font class="Apple-style-span" face="'Andale Mono'" size="3"><span class="Apple-style-span" style="font-size: 11px;"> Check constants.h to determine value of MODE_DENS_??</span></font></div><div><font class="Apple-style-span" face="'Andale Mono'" size="3"><span class="Apple-style-span" style="font-size: 11px;"> </span></font></div><div><font class="Apple-style-span" face="'Andale Mono'" size="3"><span class="Apple-style-span" style="font-size: 11px;"> Driver_abort called. See log file for details.</span></font></div></div><div><div><font class="Apple-style-span" face="'Andale Mono'" size="3"><span class="Apple-style-span" style="font-size: 11px;">Error message is [Eos_putData] ERROR Density or Internal Energy are zero after a call to EOS!</span></font></div><div><font class="Apple-style-span" face="'Andale Mono'" size="3"><span class="Apple-style-span" style="font-size: 11px;"> Calling MPI_Abort() for shutdown in 2 seconds!</span></font></div></div><div><font class="Apple-style-span" face="'Andale Mono'" size="3"><span class="Apple-style-span" style="font-size: 11px;"><br></span></font></div><div>I've successfully compiled and run a number out of the box simulations, so I don't believe it has anything to do with my libraries or compiler configuration.</div><div><br></div><div>The simulation I'm setting up is a simple version of what will eventually be more complicated. The problem setup is as follows:</div><div>Two gas clouds with the same density and temperature are in pressure equilibrium with their surroundings. They are setup so they're just touching and have some initial velocities so they will collide with each other.</div><div><br></div><div>My flash.par holds the initial density and temperature of the clouds. In the Simulation_init, I call Eos() to get the cloud pressure. Then I set the ambient pressure to be the same as the cloud pressure and the ambient density to be 1/1000 of the cloud density. I then call Eos() again, to get the ambient temperature.</div><div><br></div><div>The system is then setup with these values in Simulation_initBlock with the configuration described above. The only time I call Eos() is in the Simulation_init. I could have just calculated the pressures and temperatures by hand and hardcoded them in. </div><div><br></div><div>If I print out the values before they are put into the simulation with Grid_putPointData, the density and internal energy are not zero, but it looks like the Eos routine doesn't know the correct values.</div><div><br></div><div>Any ideas?</div><div><br></div><div>Thanks for your help.</div></body></html>