<p><b>ringler@lanl.gov</b> 2010-04-27 16:57:02 -0600 (Tue, 27 Apr 2010)</p><p><br>
in the computation of pv_edge at boundary points, initial thickness fields to zero before using<br>
</p><hr noshade><pre><font color="gray">Modified: branches/ocean_projects/z_level_mrp/mpas/src/core_ocean/module_time_integration.F
===================================================================
--- branches/ocean_projects/z_level_mrp/mpas/src/core_ocean/module_time_integration.F        2010-04-27 22:18:36 UTC (rev 216)
+++ branches/ocean_projects/z_level_mrp/mpas/src/core_ocean/module_time_integration.F        2010-04-27 22:57:02 UTC (rev 217)
@@ -1016,6 +1016,8 @@
          cell2 = cellsOnEdge(2,iEdge)
          do k = 1,nVertLevels
            if(uBC(k,iEdge).eq.1) then
+             h1 = 0.0
+             h2 = 0.0
              if(cell1.gt.0) h1 = h(k,cell1)
              if(cell2.gt.0) h2 = h(k,cell2)
              pv_edge(k,iEdge) = fEdge(iEdge) / ( max(h1,h2) )

</font>
</pre>