<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 =&gt; domain % blocklist
+       if(associated(block_ptr % next)) then
+           write(0,*) 'Error: computeGlobalDiagnostics assumes ',&amp;
+                      'that there is only one block per processor.'
+       end if
+
+       call timer_start(&quot;global_diagnostics&quot;)
+       call computeGlobalDiagnostics(domain % dminfo, &amp;
+                block_ptr % time_levs(2) % state, block_ptr % mesh, &amp;
+                itimestep, dt)
+       call timer_stop(&quot;global_diagnostics&quot;)
+   end if
+
 end subroutine mpas_timestep
 
 

</font>
</pre>