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

snapshot...

parent e2d0bb63
No related branches found
No related tags found
No related merge requests found
...@@ -488,6 +488,44 @@ def run_best_fit(raob_to_amv_dct, raob_dct, gfs_filename=None): ...@@ -488,6 +488,44 @@ def run_best_fit(raob_to_amv_dct, raob_dct, gfs_filename=None):
return out_dct return out_dct
def run_compare():
raob_dir = '/home/rink/data/raob/'
raob_files = ['raob_soundings20191117_0000.cdf',
'raob_soundings20191118_0000.cdf',
'raob_soundings20191119_0000.cdf',
'raob_soundings20191120_0000.cdf',
'raob_soundings20191121_0000.cdf',
'raob_soundings20191122_0000.cdf',
'raob_soundings20191123_0000.cdf',
'raob_soundings20191124_0000.cdf',
'raob_soundings20191125_0000.cdf',
'raob_soundings20191126_0000.cdf',
'raob_soundings20191127_0000.cdf',
'raob_soundings20191128_0000.cdf',
'raob_soundings20191129_0000.cdf',
'raob_soundings20191130_0000.cdf']
gfs_dir = '/home/rink/data/gfs/'
gfs_files = ['gfs.19111612_F012.h5',
'gfs.19111712_F012.h5',
'gfs.19111812_F012.h5',
'gfs.19111912_F012.h5',
'gfs.19112012_F012.h5',
'gfs.19112112_F012.h5',
'gfs.19112212_F012.h5',
'gfs.19112312_F012.h5',
'gfs.19112412_F012.h5',
'gfs.19112512_F012.h5',
'gfs.19112612_F012.h5',
'gfs.19112712_F012.h5',
'gfs.19112812_F012.h5',
'gfs.19112912_F012.h5']
for k, file in enumerate(raob_files):
raob_dct, ts = get_raob_dict_cdf(raob_dir+file)
compare_gfs_raob(raob_dct, gfs_dir+gfs_files[k])
def compare_gfs_raob(raob_dct, gfs_filename): def compare_gfs_raob(raob_dct, gfs_filename):
keys = list(raob_dct.keys()) keys = list(raob_dct.keys())
...@@ -545,7 +583,6 @@ def compare_gfs_raob(raob_dct, gfs_filename): ...@@ -545,7 +583,6 @@ def compare_gfs_raob(raob_dct, gfs_filename):
r_4 = np.where(raob_prs == 100)[0] r_4 = np.where(raob_prs == 100)[0]
spd_intrp = np.interp(gfs_press, raob_prs, raob_v) spd_intrp = np.interp(gfs_press, raob_prs, raob_v)
#print(raob_v[r_0] - gfs_v[key_idx, g_0], raob_v[r_1] - gfs_v[key_idx, g_1], raob_v[r_3] - gfs_v[key_idx, g_3])
if len(r_0) > 0: if len(r_0) > 0:
r_0 = r_0[0] r_0 = r_0[0]
...@@ -580,6 +617,7 @@ def compare_gfs_raob(raob_dct, gfs_filename): ...@@ -580,6 +617,7 @@ def compare_gfs_raob(raob_dct, gfs_filename):
print(np.average(np.abs(lvl_2))) print(np.average(np.abs(lvl_2)))
print(np.average(np.abs(lvl_3))) print(np.average(np.abs(lvl_3)))
print(np.average(np.abs(lvl_4))) print(np.average(np.abs(lvl_4)))
print('**************************************')
return None return None
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please to comment