diff --git a/modules/amv/intercompare.py b/modules/amv/intercompare.py index 0f74928c96d0b03dbd8148eca2fcdb53f8dfe70b..294e7b6c180f40e641bd9044a5632b5887f39ba2 100644 --- a/modules/amv/intercompare.py +++ b/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 = []