<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) > 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>