[FLASH-USERS] Using gr_forceDerefInRadius

[Klaus Weide]

On Mon, 1 Jul 2019, vgop at lle.rochester.edu wrote:

> I'm trying to get FLASH to forcibly derefine in some regions of the
> computational domain. It looks like the code used in
> gr_forceDerefInRadius.F90 would fit the bill for this kind of operation -
> however, I don't see where (if?) this routine actually gets called by the
> refinement manager. Based on what Grid_markRefineDerefine looks like, my
> first guess is that I'd add a call to gr_forceDerefInRadius at the end of
> Grid_markRefineDerefine, 

Yes, that's right.

You can see Simulation/SimulationMain/Sedov/Grid_markRefineDerefine.F90
for an example of code that (conditionally) calls gr_forceDerefInRadius.

The call doesn't have to be at the very end - it depends on which priority 
you want to give to this relative to other criteria.

I would recommend putting the call before the

   where (nodetype(:) .NE. LEAF) 
     ...
   end where

lines, at least.

> so I'd need to have a modified version of both
> files. 

You may not need to modify gr_forceDerefInRadius.F90 itself.

Klaus