<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, &
+ 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>