snapshot...

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)