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

take out lon/lat range check

parent 928d4448
No related branches found
No related tags found
No related merge requests found
...@@ -54,11 +54,11 @@ def compare_aeolus_max_height(aeolus_dict, files_path, grid_value_name='cld_heig ...@@ -54,11 +54,11 @@ def compare_aeolus_max_height(aeolus_dict, files_path, grid_value_name='cld_heig
lat = layers[0, 0] lat = layers[0, 0]
lon = layers[0, 1] lon = layers[0, 1]
# check range # # check range
if lat > G16_lat_range[1] or lat < G16_lat_range[0]: # if lat > G16_lat_range[1] or lat < G16_lat_range[0]:
continue # continue
if lon > G16_lon_range[1] or lon < G16_lon_range[0]: # if lon > G16_lon_range[1] or lon < G16_lon_range[0]:
continue # continue
# ------------------------------------------------ # ------------------------------------------------
a_lons[f_idx].append(lon) a_lons[f_idx].append(lon)
...@@ -123,6 +123,7 @@ def compare_aeolus_max_height(aeolus_dict, files_path, grid_value_name='cld_heig ...@@ -123,6 +123,7 @@ def compare_aeolus_max_height(aeolus_dict, files_path, grid_value_name='cld_heig
a_profs[f_idx] = layers a_profs[f_idx] = layers
total += len(nn_idxs) total += len(nn_idxs)
print(clvrx_files.flist[f_idx], len(nn_idxs))
# loop over Aeolus profile valid in the file # loop over Aeolus profile valid in the file
for k in range(len(nn_idxs)): for k in range(len(nn_idxs)):
j = grd_crds[0][k] j = grd_crds[0][k]
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment