diff --git a/modules/amv/aeolus.py b/modules/amv/aeolus.py index bc72a3584fda93632632e5c95dca166ad82ac19a..7c68b246598c17fd5d4d150b0298597106ef7c89 100644 --- a/modules/amv/aeolus.py +++ b/modules/amv/aeolus.py @@ -545,71 +545,6 @@ def get_grid_values(h5f, grid_name, scale_factor_name='scale_factor', add_offset return grd_vals -# assume profile height obs monotonic decreasing -# def get_aeolus_prof_dict(filename, min_num_levs=10, rcvch=0, ll2bct=1): -# prof_dict = {} -# -# with open(filename) as file: -# last_hgt = None -# last_sgn = None -# prof = [] -# idx = 0 -# key = None -# for line in file: -# toks = line.split() -# -# lat = float(toks[1]) -# lon = float(toks[2]) -# yr = int(float(toks[3])) -# mon = int(float(toks[4])) -# dy = int(float(toks[5])) -# hr = int(float(toks[6])) -# mn = int(float(toks[7])) -# az = float(toks[8]) -# ws = float(toks[9]) -# ws_err = float(toks[10]) -# hgt = float(toks[11]) -# conf = float(toks[13]) -# ch = float(toks[15]) -# clr = float(toks[16]) -# -# #if not (conf == 0 and ch == rcvch and clr == ll2bct and hgt > 0): -# if not (hgt > 0): -# continue -# -# tup = (lat, lon, az, ws, ws_err, hgt, ch, clr) -# -# # start: sign of delta hgt changes -# if last_hgt is not None: -# sgn = 1 if hgt - last_hgt >= 0 else -1 -# -# if (last_hgt is None) or (sgn != last_sgn and sgn == 1): -# if len(prof) >= min_num_levs: -# nda = np.array(prof) -# nda[:, 2] = azimuth_to_met(nda[:, 2]) -# prof_dict[key] = nda -# -# # ---- reset key make a new profile obj -# dto = datetime.datetime(year=yr, month=mon, day=dy, hour=hr, minute=mn) -# dto = dto.replace(tzinfo=timezone.utc) -# timestamp = dto.timestamp() -# key = (idx, timestamp) -# idx += 1 -# prof = [] -# -# if last_hgt is not None: -# last_sgn = sgn -# -# last_hgt = hgt -# -# # if conf == 0 and ch == rcvch and clr == ll2bct: -# prof.append(tup) -# -# file.close() -# -# return prof_dict - - def analyze_aelolus_cld_layers(cld_layer_dct): keys = list(cld_layer_dct.keys())