<p><b>mpetersen@lanl.gov</b> 2010-09-07 16:11:37 -0600 (Tue, 07 Sep 2010)</p><p>Correcting error in z-level pressure computation.<br>
</p><hr noshade><pre><font color="gray">Modified: trunk/mpas/src/core_ocean/module_time_integration.F
===================================================================
--- trunk/mpas/src/core_ocean/module_time_integration.F        2010-09-07 22:10:11 UTC (rev 494)
+++ trunk/mpas/src/core_ocean/module_time_integration.F        2010-09-07 22:11:37 UTC (rev 495)
@@ -1388,8 +1388,9 @@
          pZLevel(1,iCell) = rho(1,iCell)*gravity &amp;
             * (h(1,iCell)-0.5*hZLevel(1)) 
          do k=2,nVertLevels
-            delta_p = rho(k,iCell)*gravity*hZLevel(k)
-            pZLevel(k,iCell) = pZLevel(k-1,iCell) + 0.5*delta_p
+            pZLevel(k,iCell) = pZLevel(k-1,iCell)  &amp;
+              + 0.5*gravity*(  rho(k-1,iCell)*hZLevel(k-1) &amp;
+                             + rho(k  ,iCell)*hZLevel(k  ))
          end do
 
       end do

</font>
</pre>