speedups output files Re: [FLASH-USERS] Details of speedup after restart

[Artur Gawryszczak]

Hi Anshu,

On środa, 6 czerwca 2007, Anshu Dubey wrote:
> In theory, you could change the refinement in flash.par, but it is not
> always as straightforward as that. For instance, if you reduce
> lrefine_max, and in the checkpoint you have blocks that are at the old
> lrefine_max level, you could be in trouble. I don't believe this was ever
> addressed in FLASH2, but there is implementation to restrict the offending
> blocks in FLASH3. Even that I don't think is too well tested.

I did it many times in FLASH2.x and there were no problems. Such situations 
are by default handled safely in mark_grid_refinement, via:

     if ((lrefine(i) > lrefine_max) .and. (nodetype(i) == 1)) then
        refine(i) = .FALSE.
        derefine(i) = .TRUE.
     endif

The derefinement will take place when (mod(nstep, nrefs) == 0) and will strip 
at most one level at a time. The only place in which decrease of lrefine_max 
may lead to confusion is init_mpole from multipole Poisson solver, where 
lrefine_max is used directly for some calculations.

-- 
Cheers,
        Artur