PROFILE report

[kate]

I used to get profile reports, but haven't in quite some time now. I went into the debugger to find out why it's not printing. There's this test:

Code: Select all

IF (thread_count.eq.numthreads) THEN

thread_count is 2 at this point, being initialized to 0, then incremented both during initialization and in the wrap-up. numthreads is the number of processes, here 4. Therefore the if condition fails and nothing gets printed beyond this:

Code: Select all

 Elapsed CPU time (seconds):

 Node   #  3 CPU:      83.468
 Node   #  1 CPU:      83.468
 Node   #  2 CPU:      83.468
 Node   #  0 CPU:      83.468

What's the best way to fix this?

[hugobastos]

Hi Kate, did you find a correction for this? I'm missing the profile reports too...

[kate]

No, I got no reply. I'm guessing one possible fix for those of us who only use MPI would be to comment out the IF test.

[kate]

Not safe for OpenMP, but I got desperate for a profile report and made this change:

Code: Select all

diff --git a/ROMS/Utility/timers.F b/ROMS/Utility/timers.F
index a943335..548e09e 100644
--- a/ROMS/Utility/timers.F
+++ b/ROMS/Utility/timers.F
@@ -176,7 +176,7 @@
             END DO
           END IF
         END DO
-        IF (thread_count.eq.numthreads) THEN
+!        IF (thread_count.eq.numthreads) THEN
           thread_count=0
 #ifdef DISTRIBUTE
           op_handle(0:Nregion)='SUM'
@@ -269,7 +269,7 @@
   60      FORMAT (/,' All percentages are with respect to total time =',&
      &            5x,f12.3)
 #endif-
-        END IF
+!        END IF
 !$OMP END CRITICAL (FINALIZE_WCLOCK)
       END IF
       RETURN

Sure enough, the report came back.

[arango]

This is in my to do list. I had noticed it too. The logic is delicate because of shared-memory. I am busy with other codes now.