[FLASH-USERS] gr_findAllNeghID

[John ZuHone]

Hi Klaus,

Thanks for your reply, it was helpful.

> 1) Make sure the interface of gr_findAllNeghID is known explicitly where 
> you call it, so you should have something like:
> 
>  use gr_interfaceTypeDecl
>  use gr_interface, ONLY : gr_findAllNeghID

Yes, I already had this in there, but thanks for reminding me to check. 

> 2) You may find the subroutine gr_checkGridState useful
> 
> 3) Subroutine gr_findAllNeghID is also applied in the code that couples 
> FLASH to HYPRE, so maybe you can find a useful hint for use in 
> source/Grid/GridSolvers/HYPRE/paramesh/gr_hypreSetupGrid.F90 .

I had a look at these--the use of gr_findAllNeghID seems pretty identical to what I am doing, so I’m not sure what’s going on here. 

> 4) How many blocks and procs do you have?
> Are the block and proc numbers that occur in the following lines 
> reasonable at all?
> 
>> . My block is:         542 and proc is:         540 
>> . My block is:        1082 and proc is:        1080 
> 
> It looks suspicious that block and proc numbers are so close...

These do look odd--I hadn’t noticed that. The simulation does run fine without my special module turned on. 

> Anyway, does this problem occur even on 1 proc?
> Does it occur always, or only starting at a certain number of blocks, and/or 
> in particular grid configurations?

This did not happen in a slightly different setup where all the blocks were on the same refinement level, and it was also multiple processors--in this case the grid is highly refined. 

> Also, what exactly is the information on neighboring blocks that you want 
> to access? I am asking because there may already be different ways to get 
> at the information, as for example Grid_getBlkNeighLevels for neighboring 
> refinement levels.

Basically, I want to figure out the block number and processor of every leaf block that is touching my leaf block, whether on a face, edge, or corner. It seemed like this was the most straightforward way to do it. 

Best,

John