[FLASH-USERS] gr_findAllNeghID

[John ZuHone]

Hi all,

I suspect that this is a question mostly for Klaus, but if anyone here has any thoughts I’d appreciate it. 

I’m developing a new module in which I need to know the grid information about not only the neighbor blocks but also the blocks on the corners. For that, I am attempting to use the gr_findAllNeghID function, which is used currently by the routines which map particle properties to the mesh. I am aware that this function needs to have a sane grid structure, which is typically set up during a call to Grid_fillGuardCells. I believe that calling the function gr_ensureValidNeighborInfo should do the job, if I’m not mistaken.

However, after calling gr_findAllNeghID, I am getting errors like this, which indicate that the mesh is not in the proper state:

 Block handle error for target block:          -1 , proc:          -1 
 . My block is:          31 and proc is:          30 
  and we were trying to find neighbors to guard cell region:           1
           1           1 and my global space is: IAXIS=  -1.00000000000000     
  -1.00000000000000      JAXIS=  -1.00000000000000       -1.00000000000000      
 KAXIS=  -1.00000000000000       -1.00000000000000     
 (1) No block handle.... increase maxblocks_alloc
 DRIVER_ABORT: Damn 1
 Block handle error for target block:          -1 , proc:          -1 
 . My block is:         542 and proc is:         540 
  and we were trying to find neighbors to guard cell region:           1
           1           1 and my global space is: IAXIS=  -1.00000000000000     
  -1.00000000000000      JAXIS=  -1.00000000000000       -1.00000000000000      
 KAXIS=  -1.00000000000000       -1.00000000000000     
 (1) No block handle.... increase maxblocks_alloc
 Block handle error for target block:          -1 , proc:          -1 
 DRIVER_ABORT: Damn 1
 . My block is:        1082 and proc is:        1080 
  and we were trying to find neighbors to guard cell region:           1
           1           1 and my global space is: IAXIS=  -1.00000000000000     
  -1.00000000000000      JAXIS=  -1.00000000000000       -1.00000000000000      
 KAXIS=  -1.00000000000000       -1.00000000000000     
 (1) No block handle.... increase maxblocks_alloc
 DRIVER_ABORT: Damn 1

I’ve tried adding a call to Grid_fillGuardCells right before the call to gr_findAllNeghID (which is overkill for my application, but worth a shot), as well as gr_ensureValidNeighborInfo(10), but neither of these things work.

Does anyone have any idea what I need to do to get this function to work properly? 

Thanks,

John Z