From 4cfe686750f4f1c924ec81b95a46384bc1e65dd5 Mon Sep 17 00:00:00 2001 From: tomrink <rink@ssec.wisc.edu> Date: Thu, 27 May 2021 22:13:09 -0500 Subject: [PATCH] minor --- modules/amv/aeolus.py | 106 +++++++++++++++++++++--------------------- 1 file changed, 53 insertions(+), 53 deletions(-) diff --git a/modules/amv/aeolus.py b/modules/amv/aeolus.py index 6c86e503..25d344b0 100644 --- a/modules/amv/aeolus.py +++ b/modules/amv/aeolus.py @@ -32,59 +32,59 @@ num_elems = 5424 num_lines = 5424 -def time_dict_to_cld_layers(time_dict): - time_dict_layers = {} - - keys = list(time_dict.keys()) - for key in keys: - prof_s = time_dict[key] - layers = [] - prof = prof_s[0] - - if len(prof) == 1: - tup = prof[0] - layers.append((tup[0], tup[1], tup[3], tup[4])) - time_dict_layers[key] = layers - continue - - top = -9999.9 - last_bot = -9999.9 - tup = None - for i in range(len(prof)): - tup = prof[i] - - if i == 0: - top = tup[3] - bot = tup[4] - last_bot = bot - else: - if math.fabs(last_bot - tup[3]) > 10.0: - layers.append((tup[0], tup[1], top, last_bot)) - top = tup[3] - last_bot = tup[4] - - layers.append((tup[0], tup[1], top, tup[4])) - - time_dict_layers[key] = layers - - return time_dict_layers - - -def time_dict_to_nd_2(time_dict): - keys = list(time_dict.keys()) - for key in keys: - vals = time_dict[key] - if vals is not None: - time_dict[key] = np.stack(vals) - - return time_dict - - -def get_cloud_layers_dict(filename, lon360=False): - a_d = get_aeolus_time_dict(filename, lon360=lon360) - c_d = time_dict_to_cld_layers(a_d) - cld_lyr_dct = time_dict_to_nd_2(c_d) - return cld_lyr_dct +# def time_dict_to_cld_layers(time_dict): +# time_dict_layers = {} +# +# keys = list(time_dict.keys()) +# for key in keys: +# prof_s = time_dict[key] +# layers = [] +# prof = prof_s[0] +# +# if len(prof) == 1: +# tup = prof[0] +# layers.append((tup[0], tup[1], tup[3], tup[4])) +# time_dict_layers[key] = layers +# continue +# +# top = -9999.9 +# last_bot = -9999.9 +# tup = None +# for i in range(len(prof)): +# tup = prof[i] +# +# if i == 0: +# top = tup[3] +# bot = tup[4] +# last_bot = bot +# else: +# if math.fabs(last_bot - tup[3]) > 10.0: +# layers.append((tup[0], tup[1], top, last_bot)) +# top = tup[3] +# last_bot = tup[4] +# +# layers.append((tup[0], tup[1], top, tup[4])) +# +# time_dict_layers[key] = layers +# +# return time_dict_layers +# +# +# def time_dict_to_nd_2(time_dict): +# keys = list(time_dict.keys()) +# for key in keys: +# vals = time_dict[key] +# if vals is not None: +# time_dict[key] = np.stack(vals) +# +# return time_dict +# +# +# def get_cloud_layers_dict(filename, lon360=False): +# a_d = get_aeolus_time_dict(filename, lon360=lon360) +# c_d = time_dict_to_cld_layers(a_d) +# cld_lyr_dct = time_dict_to_nd_2(c_d) +# 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): -- GitLab