Bugfix in read_table_TOUGH2(): primary variable table for certain TOUGH2 output wasn't being read properly, because of a variant internal header format in the listing
Bugfix in next_table_TOUGH2(): now recognizes (and skips) additional mass flow table in EOS7c listings, so that subsequent tables can be read
Added silent option to mulgrid fit_surface()
mulgrid fit_surface() now raises an exception if columns with more than 4 sides are specified
Added call to use_solver(useUmfpack = False) in fit_surface(), to avoid deprecation warnings about umfpack
Modified t2listing module to handle tables which have blank lines between the header and units lines (e.g. Petrasim EOS3 connection tables)- tidied and modularised some code in the process
Cleanup of t2listing code modifications
Added ability to plot flows in mulgrid layer_plot() method
Renamed mulgrid layer_plot_flow() method to plot_flows()
Added mulgrid plot_colourbar() method, and used it in layer_plot() and slice_plot()
Added ability to plot flows in mulgrid slice_plot() method
Bugfixes in mulgrid slice_plot() and layer_plot()- block names and flow arrows weren't being drawn at the right elevation for some surface blocks on slices; block index for flow arrows wasnt being calculated correctly
Altered calculation of connection node list in mulgrid add_connection(), so that node list is always consistently ordered
Bugfix in slice_plot() and layer_plot()- new plot_flow() method was being called even when no flows are to be plotted
Added t2grid calculate_connection_centres() method
Added flow_arrow_width parameter to mulgrid layer_plot() and slice_plot() methods
Rewrote t2grid flux_matrix() method, using constant or linear flow fitting over each block and pseudo-inverse matrices to assemble into global flux matrix
Rewrote calculation of vertical connection centres and normals in t2grid calculate_connection_centres()
Added flux_matrix parameter to t2listing write_vtk() method, so flux matrices calculated for slice or layer plots can be re-used for write_vtk()
Changed name of mulgrid slice_plot() and layer_plot() flow_key_pos parameter to flow_scale_pos