[FLASH-USERS] gr_markInRectangle.F90 not accessing all blocks in lnblocks?

Hansen, Eddie ehansen at pas.rochester.edu
Fri Sep 8 22:30:33 EDT 2023


How many procs are you using? It might be better to just look at a chk file to see if the refinement is giving you the desired result.

******************************************
Eddie Hansen | Research Scientist
Flash Center for Computational Science
Dept. of Physics and Astronomy
University of Rochester, Rochester, NY
Cell 607-341-6126 | flash.rochester.edu
(he / him / his)
******************************************
________________________________
From: flash-users <flash-users-bounces at flash.rochester.edu> on behalf of Lewis,Sean <scl63 at drexel.edu>
Sent: Friday, September 8, 2023 4:35:48 PM
To: flash-users at flash.rochester.edu <flash-users at flash.rochester.edu>
Subject: [FLASH-USERS] gr_markInRectangle.F90 not accessing all blocks in lnblocks?


Hi all,



I am attempting to use gr_markInRectangle.F90 to define a region outside of which blocks are forced to derefine. This is different from the original use of gr_markInRectangle, but even before modification I am not sure if the subroutine has access to all blocks.



gr_markInRectangle loops over all blocks and is called by Grid_markRefineSpecialized.F90 which I call from Grid_markRefineDerefine. To check if gr_markInRectangle is in fact looping over all blocks, I print out “lnblocks” within the subroutine. In doing this, I only ever see block IDs 1-20 each FLASH cycle despite there being >1500 blocks. The initial grid has ~500 blocks and I have set the “contained” parameter in gr_markInRectangle such that any overlap should trigger a refinement flag change. But, I do not see the refinement flag changes reflected in the Rectangle defined, adding to my suspicion that all blocks are not being looped over.



Has anyone successfully used gr_markInRectangle or have some advice on further diagnostics?



-Sean


-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://flash.rochester.edu/pipermail/flash-users/attachments/20230909/598bfee8/attachment.htm>


More information about the flash-users mailing list