<p><b>duda</b> 2010-04-19 18:37:15 -0600 (Mon, 19 Apr 2010)</p><p>BRANCH COMMIT<br>
<br>
Decide whether an edge or vertex is owned by a block by looking<br>
at the first non-zero entry of cellsOn{Edge,Vertex} rather than<br>
just the first entry. The 0 cell represents a null/missing cell,<br>
and it should be possible for this cell to appear first in the<br>
cellsOn{Edge,Vertex} lists.<br>
<br>
M    module_block_decomp.F<br>
</p><hr noshade><pre><font color="gray">Modified: branches/lateral_boundary_conditions/src/framework/module_block_decomp.F
===================================================================
--- branches/lateral_boundary_conditions/src/framework/module_block_decomp.F        2010-04-16 00:07:00 UTC (rev 190)
+++ branches/lateral_boundary_conditions/src/framework/module_block_decomp.F        2010-04-20 00:37:15 UTC (rev 191)
@@ -142,8 +142,11 @@
       edgeIDListLocal(:) = edgeIDList(:)
 
       do i=1,nEdges
-if (cellsOnEdge(1,i) == 0) write(0,*) 'block_decomp_partitioned_edge_list: cellsOnEdge(1,iEdge) is never expected to be 0'
-         if (hash_search(h, cellsOnEdge(1,i))) then
+         do j=1,maxCells
+            if (cellsOnEdge(j,i) /= 0) exit
+         end do
+if (j &gt; maxCells) write(0,*) 'Error in block_decomp_partitioned_edge_list: edge/vertex is not adjacent to any valid cells'
+         if (hash_search(h, cellsOnEdge(j,i))) then
             lastEdge = lastEdge + 1
             edgeIDList(lastEdge) = edgeIDListLocal(i)
          else

</font>
</pre>