<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 => domain % blocklist
do while (associated(block))
call compute_injection_rate(block % state % time_levs(1) % state, block % mesh, gVor, localArea)
block => block % next
enddo
- ! write(6,*) ' calling dmpar_sum_real'
+ write(6,*) ' calling dmpar_sum_real', gVor, localArea
block => domain % blocklist
do while (associated(block))
call dmpar_sum_real(domain % dminfo, gVor, globalgVor)
call dmpar_sum_real(domain % dminfo, localArea, globalArea)
block => block % next
end do
- ! write(6,*) ' filling pvForcingAmp'
block => 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 => 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>