Suggestion: add option absolute value or not for field_diff_on_grid on data_manip

It could it be interesting to add an option for computed field diff on grid without np.absolute() in order to see if we are bellow or above the reference simulation.

Example used in http://localhost:8888/notebooks/data_manip/extraction/regular_grid.ipynb [...] from data_manip.formats.regular_grid import field_diff_on_grid [...] data_dif = field_diff_on_grid(data1,data2) data_dif = field_diff_on_grid(data1,data2,abs=False)

with modification of the python source as:

def field_diff_on_grid(data1, data2, abs=True): """ Difference between two fields on the same regulare grid @param data1 (list(numpy.array) or numpy.array) list of fields or scalar field @param data2 (list(numpy.array) or numpy.array) list of fields or field @param abs (bool) True (default) if absolute difference, else False @return (list (numpy.array) or numpy.array) difference between two fields """ if abs: if isinstance(data1, list) and isinstance(data1, list): data_i = [] for dat1, dat2 in zip(data1, data2): data_i.append(np.absolute(dat1-dat2)) elif not (isinstance(data1, list) and isinstance(dat1, list)): data_i = np.absolute(data1 - data2) else: raise TelemacException('incompatible type between data1 and data2') else: if isinstance(data1, list) and isinstance(data1, list): data_i = [] for dat1, dat2 in zip(data1, data2): data_i.append(dat1-dat2) elif not (isinstance(data1, list) and isinstance(dat1, list)): data_i = (data1 - data2) else: raise TelemacException('incompatible type between data1 and data2')

return data_i