From a1d57905623ce56de278d204378494cdceef4970 Mon Sep 17 00:00:00 2001
From: Thomas White <taw@physics.org>
Date: Thu, 19 Nov 2015 13:08:11 +0100
Subject: Monitor average B

---
 src/scaling.c | 8 ++++++++
 1 file changed, 8 insertions(+)

diff --git a/src/scaling.c b/src/scaling.c
index dfe76d69..05e741ca 100644
--- a/src/scaling.c
+++ b/src/scaling.c
@@ -438,6 +438,14 @@ void scale_all(Crystal **crystals, int n_crystals, int nthreads,
 		STATUS("Log residual went from %e to %e, %i crystals\n",
 		       bef_res, new_res, ninc);
 
+		int i;
+		double meanB = 0.0;
+		for ( i=0; i<n_crystals; i++ ) {
+			meanB += crystal_get_Bfac(crystals[i]);
+		}
+		meanB /= n_crystals;
+		STATUS("Mean B = %e\n", meanB);
+
 		reflist_free(full);
 
 	} while ( fabs(new_res-old_res) >= 0.01*old_res );
-- 
cgit v1.2.3