snapshot...

modules/amv/intercompare.py

@@ -1161,9 +1161,6 @@ def analyze_best_fit_all(amvs_dct, bfs_dct, bin_size=200):
 def analyze_best_fit_single(amvs, bfs, bin_size=200):
     bfs = np.array(bfs)
 
-    # keep = bfs[:, 4] == 0
-    # keep_idxs = bfs[keep, 0]
-    # keep_idxs = keep_idxs.astype(np.int32)
     keep_idxs = bfs[:, 0]
     keep_idxs = keep_idxs.astype(np.int32)
 
@@ -1172,6 +1169,7 @@ def analyze_best_fit_single(amvs, bfs, bin_size=200):
     diff = amv_p - bf_p
     mad = np.average(np.abs(diff))
     bias = np.average(diff)
+    print('********************************************************')
     print('num of best fits: ', bf_p.shape[0])
     print('press, MAD: ', mad)
     print('press, bias: ', bias)
@@ -1214,7 +1212,6 @@ def analyze_best_fit_single(amvs, bfs, bin_size=200):
     vd_std = np.std(vd)
     print('VD bias/rms: ', vd_mean, np.sqrt(vd_mean**2 + vd_std**2))
     print('------------------------------------------')
-    print('------------------------------------------')
 
     x_values = []
     num_pres = []