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

snapshot...

parent 43332cb1
No related branches found
No related tags found
No related merge requests found
...@@ -11,15 +11,16 @@ clavrx_dir = '/apollo/cloud/scratch/ICING/' ...@@ -11,15 +11,16 @@ clavrx_dir = '/apollo/cloud/scratch/ICING/'
dir_fmt = '%Y_%m_%d_%j' dir_fmt = '%Y_%m_%d_%j'
# dir_list = [f.path for f in os.scandir('.') if f.is_dir()] # dir_list = [f.path for f in os.scandir('.') if f.is_dir()]
ds_dct = {} ds_dct = {}
ice_dict = None
no_ice_dict = None
time_keys = None
pirep_file = '/home/rink/data/pireps/pireps_2019010000_2019063023.csv' pirep_file = '/home/rink/data/pireps/pireps_2019010000_2019063023.csv'
ice_dict, no_ice_dict = pirep_icing(pirep_file)
print('num obs: ice, no ice', len(ice_dict), len(no_ice_dict))
time_keys = list(ice_dict.keys())
hist = np.zeros(20) def setup():
ice_dict, no_ice_dict = pirep_icing(pirep_file)
print('num obs: ice, no ice', len(ice_dict), len(no_ice_dict))
time_keys = list(ice_dict.keys())
def get_clavrx_datasource(timestamp): def get_clavrx_datasource(timestamp):
...@@ -41,7 +42,7 @@ def get_grid_values(h5f, grid_name, j_c, i_c, half_width, scale_factor_name='sca ...@@ -41,7 +42,7 @@ def get_grid_values(h5f, grid_name, j_c, i_c, half_width, scale_factor_name='sca
if j_l < 0 or i_l < 0: if j_l < 0 or i_l < 0:
return None return None
j_r = j_c+half_width+1 j_r = j_c+half_width+1
i_r = j_c+half_width+1 i_r = i_c+half_width+1
if j_r >= ylen or i_r >= xlen: if j_r >= ylen or i_r >= xlen:
return None return None
......
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