diff --git a/modules/amv/intercompare.py b/modules/amv/intercompare.py index ff7308644a7357484b6c98b7c96a235008e4d777..eefd67bea495d783c82ddbe3b0d7287678be44ee 100644 --- a/modules/amv/intercompare.py +++ b/modules/amv/intercompare.py @@ -356,14 +356,6 @@ def amv_raob_diff(amv_lat, amv_lon, amv_spd, amv_dir, amv_press, raob_dct): rdir = raob[lev_idx, raob_dir_idx] rspd = raob[lev_idx, raob_spd_idx] - # ddiff = np.abs(adir - rdir) - # if ddiff > 180: - # ddiff = 360 - ddiff - # if rdir < adir: - # ddiff = -1 * ddiff - # else: - # if rdir > adir: - # ddiff = -1 * ddiff ddiff = direction_difference(adir, rdir) dir_diffs.append(ddiff) @@ -581,15 +573,6 @@ def raob_gfs_diff(raob_dct, gfs_filename): gdir = g_wdir[j, g_lev_idx] print(rspd, gspd, ':', rdir, gdir) - - # ddiff = np.abs(adir - rdir) - # if ddiff > 180: - # ddiff = 360 - ddiff - # if rdir < adir: - # ddiff = -1 * ddiff - # else: - # if rdir > adir: - # ddiff = -1 * ddiff ddiff = direction_difference(gdir, rdir) dir_diffs.append(ddiff) @@ -673,14 +656,6 @@ def amv_gfs_diff(amvs, gfs_filename): rspd = wspd[j, lev_idx] rdir = wdir[j, lev_idx] - # ddiff = np.abs(adir - rdir) - # if ddiff > 180: - # ddiff = 360 - ddiff - # if rdir < adir: - # ddiff = -1 * ddiff - # else: - # if rdir > adir: - # ddiff = -1 * ddiff ddiff = direction_difference(adir, rdir) dir_diffs.append(ddiff) @@ -735,6 +710,14 @@ def amv_gfs_diff(amvs, gfs_filename): return None +# ddiff = np.abs(adir - rdir) +# if ddiff > 180: +# ddiff = 360 - ddiff +# if rdir < adir: +# ddiff = -1 * ddiff +# else: +# if rdir > adir: +# ddiff = -1 * ddiff def direction_difference(dir_a, dir_b): diff = np.abs(dir_a - dir_b)