[FLASH-USERS] ERROR After calling Eos

[Sean Lynch]

Hello, 

I have been getting an error similar to that described last month here
http://flash.uchicago.edu/pipermail/flash-users/2010-July/003084.html


The error I get is:
ERROR After calling Eos, eosData(EOS_EINT) or eosData(EOS_DENS) are zero
 iflag=           1           1           1           1           1           1           1           1
 solnData(EINT_VAR,ib:ie,jb,kb)=   0.0000000000000000        0.0000000000000000        0.0000000000000000        0.0000000000000000        0.0000000000000000        0.0000000000000000        0.0000000000000000        0.0000000000000000     
 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
 solnData(DENS_VAR,ib:ie,jb,kb)=   0.0000000000000000        0.0000000000000000        0.0000000000000000        0.0000000000000000        0.0000000000000000        0.0000000000000000        0.0000000000000000        0.0000000000000000     
 eosData(dens+1:dens+vecLen)   =   0.0000000000000000        0.0000000000000000        0.0000000000000000        0.0000000000000000        0.0000000000000000        0.0000000000000000        0.0000000000000000        0.0000000000000000     
   Perhaps the initialization routine is wrong..... or
   perhaps the runtime parameter eosMode is wrong.
      Check constants.h to determine value of MODE_DENS_??
  
 Driver_abort called. See log file for details.
Error message is [Eos_putData] ERROR Density or Internal Energy are zero after a call to EOS!
 Calling MPI_Abort() for shutdown in   2 seconds!

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.

The simulation I'm setting up is a simple version of what will eventually be more complicated. The problem setup is as follows:
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.

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.

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. 

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.

Any ideas?

Thanks for your help.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://flash.rochester.edu/pipermail/flash-users/attachments/20100824/5aa1f102/attachment.htm>