<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 &
* (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) &
+ + 0.5*gravity*( rho(k-1,iCell)*hZLevel(k-1) &
+ + rho(k ,iCell)*hZLevel(k ))
end do
end do
</font>
</pre>