<p><b>laura@ucar.edu</b> 2012-01-13 11:31:48 -0700 (Fri, 13 Jan 2012)</p><p>clean-up of the surface pressure<br>
</p><hr noshade><pre><font color="gray">Modified: branches/atmos_physics/src/core_init_nhyd_atmos/module_test_cases.F
===================================================================
--- branches/atmos_physics/src/core_init_nhyd_atmos/module_test_cases.F        2012-01-13 18:20:04 UTC (rev 1361)
+++ branches/atmos_physics/src/core_init_nhyd_atmos/module_test_cases.F        2012-01-13 18:31:48 UTC (rev 1362)
@@ -4130,9 +4130,10 @@
       end if    ! config_met_interp
 
      
-      ! Calculate surface pressure
+      ! Calculate surface pressure (This is an ad-hoc calculation. The actual surface pressure is actually re-calculated at
+      !the top of the subroutine MPAS_to_physics in ./core_physics/module_physics_interface_nhyd.F
       do iCell=1,grid%nCells
-         diag % surface_pressure % array(iCell) = 0.5*gravity*(zgrid(2,i)-zgrid(1,i)) &amp;
+         diag % surface_pressure % array(iCell) = 0.5*gravity/rdzw(1)                                              &amp;
                                                 * (1.25* rho_zz(1,iCell) * (1. + scalars(state % index_qv, 1, iCell))  &amp;
                                                 -  0.25* rho_zz(2,iCell) * (1. + scalars(state % index_qv, 2, iCell)))
          diag % surface_pressure % array(iCell) = diag % surface_pressure % array(iCell) + pp(1,iCell) + ppb(1,iCell)

</font>
</pre>