[FLASH-USERS] using gr_lrefineMaxRedRadiusFact

[Sean Couch]

Hi Salvatore,

What you are attempting to do can be achieved by adding a few custom calls to Grid_markRefineDerefine.F90.  First, add a copy of Grid_markRefineDerefine.F90 to your Simulation setup directory that you will customize.  Then add to it:

  call gr_enforceMaxRefine(lrefine_max-3) ! Enforce a max level everywhere that is 3 levels lower than lrefine_max
 
  spec(1:3) = 0.0
  spec(4) = RADIUS_OF_REFINED_REGION
  call Grid_markRefineSpecialized(INRADIUS,4,spec,lrefine_max) ! Force refinement to lrefine_max within a sphere at center

The gre_lrefineMaxRedDoByLogR will force decrements as a function of radius which is not, if I understand, exactly what you want.

Hope this helps,
Sean



------------------------------------------------------------------------
Sean M. Couch
Assistant Professor
Department of Physics and Astronomy
Department of Comp. Math, Science, and Engineering
FRIB/National Superconducting Cyclotron Laboratory
Michigan State University
567 Wilson Rd, 3250 BPS
East Lansing, MI  48824
(517) 884-5035
couch at pa.msu.edu <mailto:couch at pa.msu.edu>




> On Oct 14, 2015, at 11:48 AM, Salvatore Cielo <cielo at mpia.de> wrote:
> 
> Hello all,
> 
> I am trying to set up some simulation (in FLASH 4.2.2) with a higher zoom level in the central region 
> (refinement 3 or 4 levels higher than the rest of the box). 
> I would like to keep the AMR structure also in that region, so I though the simplest way would be to use the runtime parameter
> gr_lrefineMaxRedRadiusFact (and the others:  gr_lrefineMaxRedDoByLogR and x[yz]_refine_center).
> 
> But I can't really figure out how it is supposed to work exactly: how does it affect the maximum local refinement level? 
> The description is not really clear... does the level drop repeatedly by one over the value I assign to it? 
> From my early tests it doesn't seem to be so... 
> 
> Does anybody have experience with that?
> 
> Thank you!
> 
> Salvatore
> ____________________________________
> Salvatore Cielo
> 
> Max-Planck-Institute for Astronomy
> Koenigstuhl 17
> D-69117 Heidelberg
> 
> Tel.: (++49|0) 6221 - 528-203
> E-mail:  cielo at mpia.de <mailto:cielo at mpia.de>
> ____________________________________

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://flash.rochester.edu/pipermail/flash-users/attachments/20151014/a3dede09/attachment-0001.htm>