<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)) &
+ diag % surface_pressure % array(iCell) = 0.5*gravity/rdzw(1) &
* (1.25* rho_zz(1,iCell) * (1. + scalars(state % index_qv, 1, iCell)) &
- 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>