[FLASH-USERS] WARNING after gc fillling problem

[Seyit Hocuk]

Hi Gary,

On the first point, are you compiling with -O3. Somehow not using -O3 is 
punishable by Flash3 with zero variable values :).
This is what I noticed.

Seyit



Gary F Forjan wrote:
> 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