<p><b>ringler@lanl.gov</b> 2011-10-18 16:32:04 -0600 (Tue, 18 Oct 2011)</p><p><br>
added diagnostic print statements in order to figure out why global reduction is not correct<br>
</p><hr noshade><pre><font color="gray">Modified: branches/pv_based_swm/mpas/src/core_pvsw/module_time_integration.F
===================================================================
--- branches/pv_based_swm/mpas/src/core_pvsw/module_time_integration.F        2011-10-18 18:03:46 UTC (rev 1107)
+++ branches/pv_based_swm/mpas/src/core_pvsw/module_time_integration.F        2011-10-18 22:32:04 UTC (rev 1108)
@@ -77,25 +77,23 @@
 
     ! find amplitude of forcing based on user-specified enstrophy injection rate
       if(config_forcing) then
-    !   write(6,*) ' calling compute_injection_rate '
         block =&gt; domain % blocklist
         do while (associated(block))
            call compute_injection_rate(block % state % time_levs(1) % state, block % mesh, gVor, localArea)
            block =&gt; block % next
         enddo
-    !   write(6,*) ' calling dmpar_sum_real'
+        write(6,*) ' calling dmpar_sum_real', gVor, localArea
         block =&gt; domain % blocklist
         do while (associated(block))
            call dmpar_sum_real(domain % dminfo, gVor, globalgVor)
            call dmpar_sum_real(domain % dminfo, localArea, globalArea)
            block =&gt; block % next
         end do
-    !   write(6,*) ' filling pvForcingAmp'
         block =&gt; domain % blocklist
         do while (associated(block))
           block % state % time_levs(1) % state % pvForcingAmp % scalar = config_enstrophy_injection_rate / (globalgVor / globalArea)  ! units of s^-2
-          write(6,10) 'pvForcingAmp :',block % state % time_levs(1) % state % xtime % scalar, block % state % time_levs(1) % state % pvForcingAmp % scalar, globalgVor
-          10 format(a15,3e20.10)
+          write(6,10) 'pvForcingAmp :',block % state % time_levs(1) % state % xtime % scalar, block % state % time_levs(1) % state % pvForcingAmp % scalar, globalgVor, globalArea
+          10 format(a15,4e20.10)
           block =&gt; block % next
         enddo
       endif
@@ -948,7 +946,7 @@
        enddo
       enddo
 
-      write(6,*) ' exiting compute_injection_rate ', gVor
+      write(6,*) ' exiting compute_injection_rate ', gVor, localArea
 
    end subroutine compute_injection_rate
 

</font>
</pre>