<p><b>ringler@lanl.gov</b> 2010-04-06 09:35:46 -0600 (Tue, 06 Apr 2010)</p><p><br>
code modified: module_io_output.F<br>
<br>
reason: test required for domain % blocklist % mesh % nEdgesOnEdge % array(i) .eq. 0 in order<br>
to prevent addressing array out of bounds<br>
</p><hr noshade><pre><font color="gray">Modified: branches/lateral_boundary_conditions/src/framework/module_io_output.F
===================================================================
--- branches/lateral_boundary_conditions/src/framework/module_io_output.F        2010-04-01 18:47:27 UTC (rev 177)
+++ branches/lateral_boundary_conditions/src/framework/module_io_output.F        2010-04-06 15:35:46 UTC (rev 178)
@@ -187,8 +187,12 @@
                                                                            domain % blocklist % mesh % edgesOnEdge % array(j,i))
          end do
          do j=domain % blocklist % mesh % nEdgesOnEdge % array(i)+1,2*domain % blocklist % mesh % maxEdges
-            edgesOnEdge(j,i) = domain % blocklist % mesh % indexToEdgeID % array( &amp;
+            if(domain % blocklist % mesh % nEdgesOnEdge % array(i) .eq. 0) then
+               edgesOnEdge(j,i) = domain % blocklist % mesh % nEdgesSolve + 1
+            else
+               edgesOnEdge(j,i) = domain % blocklist % mesh % indexToEdgeID % array( &amp;
                                                                            domain % blocklist % mesh % nEdgesOnEdge % array(i))
+            endif
          end do
       end do
       do i=1,domain % blocklist % mesh % nVerticesSolve

</font>
</pre>