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

snapshot...

parent 89a8b4fc
No related merge requests found
...@@ -295,14 +295,20 @@ def get_grid_values_all(h5f, grid_name, scale_factor_name='scale_factor', add_of ...@@ -295,14 +295,20 @@ def get_grid_values_all(h5f, grid_name, scale_factor_name='scale_factor', add_of
attr = attrs.get(scale_factor_name) attr = attrs.get(scale_factor_name)
if attr is None: if attr is None:
raise GenericException('Attribute: '+scale_factor_name+' not found for variable: '+grid_name) raise GenericException('Attribute: '+scale_factor_name+' not found for variable: '+grid_name)
scale_factor = attr[0] if np.isscalar(attr):
scale_factor = attr
else:
scale_factor = attr[0]
grd_vals = grd_vals * scale_factor grd_vals = grd_vals * scale_factor
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 None:
raise GenericException('Attribute: '+add_offset_name+' not found for variable: '+grid_name) raise GenericException('Attribute: '+add_offset_name+' not found for variable: '+grid_name)
add_offset = attr[0] if np.isscalar(attr):
add_offset = attr
else:
add_offset = attr[0]
grd_vals = grd_vals + add_offset grd_vals = grd_vals + add_offset
if range_name is not None: if range_name is not None:
...@@ -317,7 +323,10 @@ def get_grid_values_all(h5f, grid_name, scale_factor_name='scale_factor', add_of ...@@ -317,7 +323,10 @@ def get_grid_values_all(h5f, grid_name, scale_factor_name='scale_factor', add_of
attr = attrs.get(fill_value_name) attr = attrs.get(fill_value_name)
if attr is None: if attr is None:
raise GenericException('Attribute: '+fill_value_name+' not found for variable: '+grid_name) raise GenericException('Attribute: '+fill_value_name+' not found for variable: '+grid_name)
fill_value = attr[0] if np.isscalar(attr):
fill_value = attr
else:
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)
return grd_vals return grd_vals
......
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