Fix bug with restart runs in parallel.

openmc-dev · Aug 19, 2012 · 7fd6176 · 7fd6176
1 parent c1e2a1c
commit 7fd6176
Showing 1 changed file with 12 additions and 10 deletions.
diff --git a/src/state_point.F90 b/src/state_point.F90
@@ -343,22 +343,24 @@ subroutine replay_batch_history
        ! calculate mean keff
        keff = temp(1) / n_realizations
 
-       if (confidence_intervals) then
-          ! Calculate t-value for confidence intervals
-          alpha = ONE - CONFIDENCE_LEVEL
-          t_value = t_percentile(ONE - alpha/TWO, n_realizations)
-       else
-          t_value = ONE
-       end if
+       if (n_realizations > 1) then
+          if (confidence_intervals) then
+             ! Calculate t-value for confidence intervals
+             alpha = ONE - CONFIDENCE_LEVEL
+             t_value = t_percentile(ONE - alpha/TWO, n_realizations)
+          else
+             t_value = ONE
+          end if
 
-       keff_std = t_value * sqrt((temp(2)/n_realizations - keff*keff) &
-            / (n_realizations - 1))
+          keff_std = t_value * sqrt((temp(2)/n_realizations - keff*keff) &
+               / (n_realizations - 1))
+       end if
     else
        keff = k_batch(current_batch)
     end if
 
     ! print out batch keff
-    call print_batch_keff()
+    if (master) call print_batch_keff()
 
     ! Write message at end
     if (current_batch == restart_batch) then