<p><b>mpetersen@lanl.gov</b> 2010-11-15 13:47:44 -0700 (Mon, 15 Nov 2010)</p><p>BRANCH COMMIT<br>
<br>
This small change causes output to not match.  The difference is at<br>
vertices on the boundary.<br>
<br>
old version: Cycles the vertex loop for boundary vertices, so<br>
pv_vertex is zero.<br>
<br>
new version: pv_vertex computed and generally nonzero for boundary<br>
vertices.<br>
</p><hr noshade><pre><font color="gray">Modified: branches/ocean_projects/topography2_mrp/src/core_ocean/module_time_integration.F
===================================================================
--- branches/ocean_projects/topography2_mrp/src/core_ocean/module_time_integration.F        2010-11-15 19:50:10 UTC (rev 614)
+++ branches/ocean_projects/topography2_mrp/src/core_ocean/module_time_integration.F        2010-11-15 20:47:44 UTC (rev 615)
@@ -1388,10 +1388,7 @@
       ! Compute height at vertices, pv at vertices, and average pv to edge locations
       !  ( this computes pv_vertex at all vertices bounding real cells and distance-1 ghost cells )
       !
-      VTX_LOOP: do iVertex = 1,nVertices
-         do i=1,vertexDegree
-            if (cellsOnVertex(i,iVertex) &gt; nCells) cycle VTX_LOOP
-         end do
+      do iVertex = 1,nVertices
          do k=1,maxLevelVertexBot(iVertex)
             h_vertex = 0.0
             do i=1,vertexDegree
@@ -1401,7 +1398,7 @@
 
             pv_vertex(k,iVertex) = (fVertex(iVertex) + vorticity(k,iVertex)) / h_vertex
          end do
-      end do VTX_LOOP
+      end do
 
       !
       ! Compute pv at cell centers

</font>
</pre>