<p><b>laura@ucar.edu</b> 2010-12-09 13:23:56 -0700 (Thu, 09 Dec 2010)</p><p>corrected the calculation of surface_pressure in subroutine compute_large_step variables<br>
</p><hr noshade><pre><font color="gray">Modified: branches/atmos_physics/src/core_nhyd_atmos/module_time_integration.F
===================================================================
--- branches/atmos_physics/src/core_nhyd_atmos/module_time_integration.F        2010-12-09 20:23:01 UTC (rev 638)
+++ branches/atmos_physics/src/core_nhyd_atmos/module_time_integration.F        2010-12-09 20:23:56 UTC (rev 639)
@@ -1129,8 +1129,8 @@
 
          !calculation of the surface pressure:
          surface_pressure(iCell) = 0.5*gravity/rdzw(1)                                              &amp;
-                       * (1.25*(rtheta_p(1,iCell) + rtheta_base(1,iCell)) * (1. + qvapor(1,iCell))  &amp;
-                       -  0.25*(rtheta_p(2,iCell) + rtheta_base(2,iCell)) * (1. + qvapor(2,iCell)))
+                       * (1.25* rho(1,iCell) * (1. + qvapor(1,iCell))  &amp;
+                       -  0.25* rho(2,iCell) * (1. + qvapor(2,iCell)))
          surface_pressure(iCell) = surface_pressure(iCell) + pressure_p(1,iCell) + pressure_b(1,iCell)
 
 

</font>
</pre>