Skip to content
Snippets Groups Projects
Commit 219a713d authored by tomrink's avatar tomrink
Browse files

snapshot...

parent 20c9e3e0
No related branches found
No related tags found
No related merge requests found
......@@ -284,7 +284,7 @@ def run_best_fit_all(amv_path=None, product_path=None, raob_path=None, gfs_path=
if not full_domain:
for k, file in enumerate(raob_files):
raob_dct, ts = get_raob_dict_cdf(raob_dir+file)
m_d = match_amvs_to_raobs(raob_dct, ts, amv_files=amv_files)
m_d = match_amvs_to_raobs(raob_dct, ts, amv_files)
bf_dct = run_best_fit(m_d, raob_dct, gfs_dir+gfs_files[k])
prd_dct = get_product_at_locs(m_d, ts, prd_files)
out_list.append((bf_dct, prd_dct))
......@@ -1263,7 +1263,7 @@ def make_plot(bin_ranges, bin_values):
#do_plot(x_values, [std_r, std_g, std], ['ICE', 'NOT ICE', 'ALL'], ['blue', 'red', 'black'], title='ACHA - RAOB ', x_axis_label='STD (hPa)', y_axis_label='hPa', invert=True, flip=True)
#do_plot(x_values, [bias_r, bias_g, bias], ['ICE', 'NOT ICE', 'ALL'], ['blue', 'red', 'black'], title='ACHA - RAOB Pres Match', x_axis_label='BIAS (m/s)', y_axis_label='hPa', invert=True, flip=True)
do_plot(x_values, [mad_r, mad_g, mad], ['ICE', 'NOT ICE', 'ALL'], ['blue', 'red', 'black'], title='ACHA - GFS BestFit (ALL)', x_axis_label='MAD (hPa)', y_axis_label='hPa', invert=True, flip=True)
do_plot(x_values, [mad_r, mad_g, mad], ['ICE', 'NOT ICE', 'ALL'], ['blue', 'red', 'black'], title='ACHA - GFS BestFit (ALL)', x_axis_label='MAD (m/s)', y_axis_label='hPa', invert=True, flip=True)
#do_plot(x_values, [mad_r, mad_g, mad], ['ICE', 'NOT ICE', 'ALL'], ['blue', 'red', 'black'], title='ACHA - RAOB Pres Match', x_axis_label='MAD (m/s)', y_axis_label='hPa', invert=True, flip=True)
#do_plot(x_values, [num_vals_r, num_vals_g, num_vals], ['ICE: '+str(num_r), 'NOT ICE: '+str(num_g), 'ALL: '+str(num)], ['blue', 'red', 'black'], title='ACHA - RAOB BestFit', x_axis_label='log(Count)', y_axis_label='hPa', invert=True, flip=True)
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment