<p><b>laura@ucar.edu</b> 2012-07-12 11:44:27 -0600 (Thu, 12 Jul 2012)</p><p>point the local surface pressure psfc_p to the surface pressure extrapolated to the surface using the actual pressure sounding<br>
</p><hr noshade><pre><font color="gray">Modified: branches/atmos_physics/src/core_atmos_physics/mpas_atmphys_interface_nhyd.F
===================================================================
--- branches/atmos_physics/src/core_atmos_physics/mpas_atmphys_interface_nhyd.F        2012-07-12 17:40:40 UTC (rev 2019)
+++ branches/atmos_physics/src/core_atmos_physics/mpas_atmphys_interface_nhyd.F        2012-07-12 17:44:27 UTC (rev 2020)
@@ -163,7 +163,7 @@
  u =&gt; diag  % uReconstructZonal % array
  v =&gt; diag  % uReconstructMeridional % array
 
-!calculation of the surface pressure:
+!calculation of the surface pressure using hydrostatic assumption down to the surface::
  do i = its,ite
     tem1 = zgrid(2,i)-zgrid(1,i)
     tem2 = zgrid(3,i)-zgrid(2,i)
@@ -173,14 +173,13 @@
                     * (rho1 + 0.5*(rho2-rho1)*tem1/(tem1+tem2))
     sfc_pressure(i) = sfc_pressure(i) + pressure_p(1,i) + pressure_b(1,i)
  enddo
-!enddo
 
 !arrays located at theta points:
- do j = jts, jte
- do i = its, ite
-    psfc_p(i,j) = diag % surface_pressure % array(i)
- enddo
- enddo
+!do j = jts, jte
+!do i = its, ite
+!   psfc_p(i,j) = diag % surface_pressure % array(i)
+!enddo
+!enddo
 
  do j = jts, jte
  do k = kts, kte
@@ -274,7 +273,8 @@
  enddo
  enddo
 
-!interpolation of pressure and temperature from theta points to the surface:
+!ldf (2012-06-22): recalculates the pressure at the surface as an extrapolation of the
+!pressures in the 2 layers above the surface, as was originally done:
  k = kts
  do j = jts,jte
  do i = its,ite
@@ -284,9 +284,8 @@
     w1 = (z0-z2)/(z1-z2)
     w2 = 1.-w1
     t2_p(i,k,j)    = w1*t_p(i,k,j)+w2*t_p(i,k+1,j)
-    !use surface pressure calculated in subroutine recover_large_step_variables.
-    !pres2_p(i,k,j) = w1*pres_p(i,k,j)+w2*pres_p(i,k+1,j)
-    pres2_p(i,k,j) = psfc_p(i,j)
+    pres2_p(i,k,j) = w1*pres_p(i,k,j)+w2*pres_p(i,k+1,j)
+    psfc_p(i,j) = pres2_p(i,k,j)
  enddo
  enddo
 

</font>
</pre>