<p><b>laura@ucar.edu</b> 2010-07-23 16:58:15 -0600 (Fri, 23 Jul 2010)</p><p>merged with revision 395 from trunk<br>
</p><hr noshade><pre><font color="gray">Modified: branches/atmos_physics/src/core_sw/mpas_interface.F
===================================================================
--- branches/atmos_physics/src/core_sw/mpas_interface.F        2010-07-23 22:57:41 UTC (rev 433)
+++ branches/atmos_physics/src/core_sw/mpas_interface.F        2010-07-23 22:58:15 UTC (rev 434)
@@ -52,15 +52,32 @@
use grid_types
use time_integration
+ use timer
+ use global_diagnostics
implicit none
type (domain_type), intent(inout) :: domain
integer, intent(in) :: itimestep
real (kind=RKIND), intent(in) :: dt
+ type (block_type), pointer :: block_ptr
call timestep(domain, dt)
+ if (mod(itimestep, config_stats_interval) == 0) then
+ block_ptr => domain % blocklist
+ if(associated(block_ptr % next)) then
+ write(0,*) 'Error: computeGlobalDiagnostics assumes ',&
+ 'that there is only one block per processor.'
+ end if
+
+ call timer_start("global_diagnostics")
+ call computeGlobalDiagnostics(domain % dminfo, &
+ block_ptr % time_levs(2) % state, block_ptr % mesh, &
+ itimestep, dt)
+ call timer_stop("global_diagnostics")
+ end if
+
end subroutine mpas_timestep
</font>
</pre>