[FLASH-USERS] WARNING after gc fillling problem

[Gary F Forjan]

Hello FLASH Users

I have a question similar to Latif's in her email of 9/7/09, only my problem is more serious.  I am getting zero values for my ENER_VAR, DENS_VAR, and other variables and my simulation is crashing after a call to EOS.  I have a 1-dimensional domain and am using the following boundary condition types

xl_boundary_type = "user"
xr_boundary_type = "reflect"

To implement the user defined boundary condition, I have followed the WindTunnel example and just modified the Grid_applyBCEdge.F90 routine to use my own variable values which appear to be correct.  I also make a call to EOS in my Simulation_initBlock.F90 routine using temperature and density to obtain the pressure and energy.  I get the following error at the start of the flash run after several refinements have taken place (Nblockx=1)

WARNING after gc filling: min. unk(DENS_VAR)=0.000000000000000000000          PE=1     block=1                                type=1                               
  1  2.32E-15 2.24E-15 2.16E-15 2.08E-15 2.01E-15 1.94E-15 1.87E-15 1.80E-15 1.74E-15 1.68E-15 1.62E-15 1.56E-15  0.0      0.0      0.0      0.0    
WARNING after gc filling: min. unk(ENER_VAR)=0.000000000000000                PE=1     block=1                                type=1                               
  1  2.24E+14 2.24E+14 2.24E+14 2.24E+14 2.24E+14 2.24E+14 2.24E+14 2.24E+14 2.24E+14 2.24E+14 2.24E+14 2.24E+14  0.0      0.0      0.0      0.0    
WARNING after gc filling: min. unk(EINT_VAR)=0.000000000000000                PE=1     block=1                                type=1                               
  1  2.24E+14 2.24E+14 2.24E+14 2.24E+14 2.24E+14 2.24E+14 2.24E+14 2.24E+14 2.24E+14 2.24E+14 2.24E+14 2.24E+14  0.0      0.0      0.0      0.0    
 DRIVER_ABORT: 
 [Eos_putData] ERROR Density or Internal Energy are zero after a call to EOS!


Would anyone know what is causing this problem and or where I should start looking in the code or my setup files to track it down?  Thanks for any information.  I have attached my Config, Simulation_initBlock.F90, Grid_applyBCEdge.F90, and error output files as well in case they are of any help.    

One other point: I want to use just a simple fully ionized proton+electron plasma so that I can use radiative cooling in my simulation.  This was easy in FLASH2 since you could just specify in the Config file

REQUIRES materials/composition/prot+elec

For FLASH3, I assume I need to just specify

REQUIRES Simulation/SimulationComposition/Ionize

in the Config file. 

Thanks

Gary

Gary Forjan
Department of Computational and Data Sciences
George Mason University, Fairfax, Va
gforjan at gmu.edu
-------------- next part --------------
An embedded and charset-unspecified text was scrubbed...
Name: Grid_applyBCEdge.F90
URL: <http://flash.rochester.edu/pipermail/flash-users/attachments/20090914/27251bfb/attachment.ksh>
-------------- next part --------------
An embedded and charset-unspecified text was scrubbed...
Name: Simulation_initBlock.F90
URL: <http://flash.rochester.edu/pipermail/flash-users/attachments/20090914/27251bfb/attachment-0001.ksh>
-------------- next part --------------
An embedded and charset-unspecified text was scrubbed...
Name: Config
URL: <http://flash.rochester.edu/pipermail/flash-users/attachments/20090914/27251bfb/attachment-0002.ksh>
-------------- next part --------------
An embedded and charset-unspecified text was scrubbed...
Name: flash_output
URL: <http://flash.rochester.edu/pipermail/flash-users/attachments/20090914/27251bfb/attachment-0003.ksh>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: flash.log
Type: application/octet-stream
Size: 101882 bytes
Desc: not available
URL: <http://flash.rochester.edu/pipermail/flash-users/attachments/20090914/27251bfb/attachment.obj>