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

snapshot...

parent e7427171
Branches
No related tags found
No related merge requests found
...@@ -499,46 +499,52 @@ def get_grid_values(h5f, grid_name, j_c, i_c, half_width, num_j=None, num_i=None ...@@ -499,46 +499,52 @@ def get_grid_values(h5f, grid_name, j_c, i_c, half_width, num_j=None, num_i=None
i_r = i_c + num_i + 1 i_r = i_c + num_i + 1
grd_vals = hfds[j_l:j_r, i_l:i_r] grd_vals = hfds[j_l:j_r, i_l:i_r]
if fill_value is not None: if fill_value is not None:
grd_vals = np.where(grd_vals == fill_value, np.nan, grd_vals) grd_vals = np.where(grd_vals == fill_value, np.nan, grd_vals)
if scale_factor_name is not None: if scale_factor_name is not None:
attr = attrs.get(scale_factor_name) attr = attrs.get(scale_factor_name)
if attr is None: if attr is not None:
raise GenericException('Attribute: '+scale_factor_name+' not found for dataset: '+grid_name)
if np.isscalar(attr): if np.isscalar(attr):
scale_factor = attr scale_factor = attr
else: else:
scale_factor = attr[0] scale_factor = attr[0]
grd_vals = grd_vals * scale_factor grd_vals = grd_vals * scale_factor
else:
print('Attribute: ' + scale_factor_name + ' not found for dataset: ' + grid_name)
if add_offset_name is not None: if add_offset_name is not None:
attr = attrs.get(add_offset_name) attr = attrs.get(add_offset_name)
if attr is None: if attr is not None:
raise GenericException('Attribute: '+add_offset_name+' not found for dataset: '+grid_name)
if np.isscalar(attr): if np.isscalar(attr):
add_offset = attr add_offset = attr
else: else:
add_offset = attr[0] add_offset = attr[0]
grd_vals = grd_vals + add_offset grd_vals = grd_vals + add_offset
else:
print('Attribute: '+add_offset_name+' not found for dataset: '+grid_name)
if range_name is not None: if range_name is not None:
attr = attrs.get(range_name) attr = attrs.get(range_name)
if attr is None: if attr is not None:
raise GenericException('Attribute: '+range_name+' not found for dataset: '+grid_name)
low = attr[0] low = attr[0]
high = attr[1] high = attr[1]
grd_vals = np.where(grd_vals < low, np.nan, grd_vals) grd_vals = np.where(grd_vals < low, np.nan, grd_vals)
grd_vals = np.where(grd_vals > high, np.nan, grd_vals) grd_vals = np.where(grd_vals > high, np.nan, grd_vals)
else:
print('Attribute: '+range_name+' not found for dataset: '+grid_name)
elif fill_value_name is not None: elif fill_value_name is not None:
attr = attrs.get(fill_value_name) attr = attrs.get(fill_value_name)
if attr is None: if attr is not None:
raise GenericException('Attribute: '+fill_value_name+' not found for dataset: '+grid_name)
if np.isscalar(attr): if np.isscalar(attr):
fill_value = attr fill_value = attr
else: else:
fill_value = attr[0] fill_value = attr[0]
grd_vals = np.where(grd_vals == fill_value, np.nan, grd_vals) grd_vals = np.where(grd_vals == fill_value, np.nan, grd_vals)
else:
print('Attribute: '+fill_value_name+' not found for dataset: '+grid_name)
return grd_vals return grd_vals
... ...
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please to comment