[FLASH-USERS] Calculate Gravity around Init Phase

[Klaus Weide]

On Mon, 14 Mar 2016, Stephen Ro wrote:

> Hello,
> 
> I'm an amateur FLASH user and could use some guidance. Is there a way to
> re-iterate the initial conditions after computing the (multipole) gravity
> from the first iteration? 

You would have to write at least some code for this, but it may be 
minimal.

You should write a routine that calls Simulation_initBlock (or a modified 
version thereof, making sure it does not overwrite GPOT_VAR or any other 
variables you want to keep from the Gravity solve) in a loop over blocks.
This routine needs to be called one after the initial Gravity calls
(see Driver_initFlash).
A convenient place to put this could be Simulation_adjustEvolution.F90, 
which gets called at the beginning of each time step; you can check the
nstep argument to make sure it only does something when it is called for 
the first time in a simulation.

Klaus