[FLASH-USERS] WARNING after gc filling problem

[Gary F Forjan]

Hi Seyit

Thanks for the suggestion but compiling the code with the -O3 flag did not resolve the problem of zero values for the variables.  I have to think Klaus is right that something is wrong in Paramesh but I am not sure what in my Config and flash.par file is causing it.  Perhaps my composition.  I include only the line

REQUIRES Simulation/SimulationComposition/Ionize

in my Config file to implement a proton+elec plasma but the Flash.h file shows that NSPECIES is defined as 5.  I am not sure how FLASH determines that.

Gary

Gary Forjan
Department of Computational and Data Sciences
George Mason University, Fairfax, Va
gforjan at gmu.edu

----- Original Message -----
From: Seyit Hocuk <seyit at astro.rug.nl>
Date: Monday, September 14, 2009 6:29 am
Subject: Re: [FLASH-USERS] WARNING after gc fillling problem

> 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
> 
>