<p><b>laura@ucar.edu</b> 2013-04-18 15:52:55 -0600 (Thu, 18 Apr 2013)</p><p>Save the model-top pressure to output files.<br>
</p><hr noshade><pre><font color="gray">Modified: branches/atmos_physics/src/core_atmos_physics/mpas_atmphys_driver.F
===================================================================
--- branches/atmos_physics/src/core_atmos_physics/mpas_atmphys_driver.F        2013-04-18 21:50:27 UTC (rev 2776)
+++ branches/atmos_physics/src/core_atmos_physics/mpas_atmphys_driver.F        2013-04-18 21:52:55 UTC (rev 2777)
@@ -60,7 +60,8 @@
 
     !physics prep step:
 #ifdef non_hydrostatic_core
-    call MPAS_to_physics(block%mesh,block%state%time_levs(1)%state,block%diag)
+    call MPAS_to_physics(block%mesh,block%state%time_levs(1)%state,block%diag, &amp;
+                         block%diag_physics)
 #elif hydrostatic_core
     call MPAS_to_physics(block%state%time_levs(1)%state,block%diag)
 #endif

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        2013-04-18 21:50:27 UTC (rev 2776)
+++ branches/atmos_physics/src/core_atmos_physics/mpas_atmphys_interface_nhyd.F        2013-04-18 21:52:55 UTC (rev 2777)
@@ -111,7 +111,7 @@
  end subroutine deallocate_forall_physics
 
 !=============================================================================================
- subroutine MPAS_to_physics(mesh,state,diag)
+ subroutine MPAS_to_physics(mesh,state,diag,diag_physics)
 !=============================================================================================
 
 !input variables:
@@ -119,6 +119,9 @@
  type(state_type),intent(in):: state
  type(diag_type) ,intent(in):: diag
 
+!inout variables:
+ type(diag_physics_type),intent(inout):: diag_physics
+
 !local variables:
  integer:: i,k,j
  real(kind=RKIND):: z0,z1,z2,w1,w2
@@ -322,6 +325,13 @@
  enddo
  enddo
 
+!save the model-top pressure:
+ do j = jts,jte
+ do i = its,ite
+    diag_physics % plrad % array(i) = pres2_p(i,kte+1,j) 
+ enddo
+ enddo
+
 !formats: 
  201 format(3i8,10(1x,e15.8))
  202 format(2i6,10(1x,e15.8))

</font>
</pre>