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

minor

parent d8ff6113
No related branches found
No related tags found
No related merge requests found
...@@ -32,59 +32,59 @@ num_elems = 5424 ...@@ -32,59 +32,59 @@ num_elems = 5424
num_lines = 5424 num_lines = 5424
def time_dict_to_cld_layers(time_dict): # def time_dict_to_cld_layers(time_dict):
time_dict_layers = {} # time_dict_layers = {}
#
keys = list(time_dict.keys()) # keys = list(time_dict.keys())
for key in keys: # for key in keys:
prof_s = time_dict[key] # prof_s = time_dict[key]
layers = [] # layers = []
prof = prof_s[0] # prof = prof_s[0]
#
if len(prof) == 1: # if len(prof) == 1:
tup = prof[0] # tup = prof[0]
layers.append((tup[0], tup[1], tup[3], tup[4])) # layers.append((tup[0], tup[1], tup[3], tup[4]))
time_dict_layers[key] = layers # time_dict_layers[key] = layers
continue # continue
#
top = -9999.9 # top = -9999.9
last_bot = -9999.9 # last_bot = -9999.9
tup = None # tup = None
for i in range(len(prof)): # for i in range(len(prof)):
tup = prof[i] # tup = prof[i]
#
if i == 0: # if i == 0:
top = tup[3] # top = tup[3]
bot = tup[4] # bot = tup[4]
last_bot = bot # last_bot = bot
else: # else:
if math.fabs(last_bot - tup[3]) > 10.0: # if math.fabs(last_bot - tup[3]) > 10.0:
layers.append((tup[0], tup[1], top, last_bot)) # layers.append((tup[0], tup[1], top, last_bot))
top = tup[3] # top = tup[3]
last_bot = tup[4] # last_bot = tup[4]
#
layers.append((tup[0], tup[1], top, tup[4])) # layers.append((tup[0], tup[1], top, tup[4]))
#
time_dict_layers[key] = layers # time_dict_layers[key] = layers
#
return time_dict_layers # return time_dict_layers
#
#
def time_dict_to_nd_2(time_dict): # def time_dict_to_nd_2(time_dict):
keys = list(time_dict.keys()) # keys = list(time_dict.keys())
for key in keys: # for key in keys:
vals = time_dict[key] # vals = time_dict[key]
if vals is not None: # if vals is not None:
time_dict[key] = np.stack(vals) # time_dict[key] = np.stack(vals)
#
return time_dict # return time_dict
#
#
def get_cloud_layers_dict(filename, lon360=False): # def get_cloud_layers_dict(filename, lon360=False):
a_d = get_aeolus_time_dict(filename, lon360=lon360) # a_d = get_aeolus_time_dict(filename, lon360=lon360)
c_d = time_dict_to_cld_layers(a_d) # c_d = time_dict_to_cld_layers(a_d)
cld_lyr_dct = time_dict_to_nd_2(c_d) # cld_lyr_dct = time_dict_to_nd_2(c_d)
return cld_lyr_dct # return cld_lyr_dct
def compare_aeolus_max_height(aeolus_dict, files_path, grid_value_name='cld_height_acha', one_cld_layer_only=False, highest_layer=False): def compare_aeolus_max_height(aeolus_dict, files_path, grid_value_name='cld_height_acha', one_cld_layer_only=False, highest_layer=False):
......
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