snapshot...

modules/aeolus/aeolus_amv.py

@@ -1074,6 +1074,8 @@ def make_plot():
     bin_ranges = tup_r[0]  # same for all
     bin_pres_r = tup_r[1]
     bin_pres_g = tup_g[1]
+    bin_spd_r = tup_r[2]
+    bin_spd_g = tup_g[2]
 
     x_values = []
     num_pres_r = []
@@ -1084,6 +1086,10 @@ def make_plot():
     pres_bias_r = []
     pres_mad_g = []
     pres_bias_g = []
+    spd_mad_r = []
+    spd_bias_r = []
+    spd_mad_g = []
+    spd_bias_g = []
 
     num_r = 0
     num_g = 0
@@ -1095,15 +1101,21 @@ def make_plot():
         x_values.append(np.average(bin_ranges[i]))
         num_pres_r.append((bin_pres_r[i].shape[0])/num_r)
         num_pres_g.append((bin_pres_g[i].shape[0])/num_g)
-        #num_spd.append(bin_spd[i].shape[0])
-        #num_dir.append(bin_dir[i].shape[0])
+
         pres_mad_r.append(np.average(np.abs(bin_pres_r[i])))
         pres_bias_r.append(np.average(bin_pres_r[i]))
         pres_mad_g.append(np.average(np.abs(bin_pres_g[i])))
         pres_bias_g.append(np.average(bin_pres_g[i]))
 
+        spd_mad_r.append(np.average(np.abs(bin_spd_r[i])))
+        spd_bias_r.append(np.average(bin_spd_r[i]))
+        spd_mad_g.append(np.average(np.abs(bin_spd_g[i])))
+        spd_bias_g.append(np.average(bin_spd_g[i]))
+
     #do_plot(x_values, [pres_mad_r, pres_mad_g], ['RAOB', 'GFS'], ['blue', 'red'], title='ACHA - BestFit', x_axis_label='MAD', y_axis_label='hPa', invert=True, flip=True)
-    do_plot(x_values, [pres_bias_r, pres_bias_g], ['RAOB', 'GFS'], ['blue', 'red'], title='ACHA - BestFit', x_axis_label='BIAS', y_axis_label='hPa', invert=True, flip=True)
+    #do_plot(x_values, [spd_mad_r, spd_mad_g], ['RAOB', 'GFS'], ['blue', 'red'], title='ACHA - BestFit', x_axis_label='MAE (m/s)', y_axis_label='hPa', invert=True, flip=True)
+    #do_plot(x_values, [pres_bias_r, pres_bias_g], ['RAOB', 'GFS'], ['blue', 'red'], title='ACHA - BestFit', x_axis_label='BIAS', y_axis_label='hPa', invert=True, flip=True)
+    #do_plot(x_values, [spd_bias_r, spd_bias_g], ['RAOB', 'GFS'], ['blue', 'red'], title='ACHA - BestFit', x_axis_label='BIAS (m/s)', y_axis_label='hPa', invert=True, flip=True)
     #do_plot(x_values, [num_pres_r, num_pres_g], ['RAOB:'+str(num_r), 'GFS:'+str(num_g)], ['blue', 'red'], x_axis_label='Normalized Count', y_axis_label='hPa', invert=True, flip=True)