plot a Serafin results file using Python API about pyteltools HOT 4 CLOSED

import matplotlib.pyplot as plt
from matplotlib import cm
import matplotlib.tri as mtri
from mpl_toolkits.axes_grid1 import make_axes_locatable
import numpy as np

from pyteltools.slf import Serafin


with Serafin.Read('../PyTelTools_validation/data/Yen/fis_yen-exp.slf', 'en') as resin:
    # Read header (SerafinHeader is stored in `header` attribute of `Serafin` class)
    resin.read_header()

    # Display a summary
    print(resin.header.summary())

    # Get time (in seconds) and display it
    resin.get_time()
    print(resin.time)

    values = resin.read_var_in_frame(10, 'B')

    fig, axes = plt.subplots(1)
    fig.set_size_inches(10.0, 6.0)
    levels = np.linspace(np.nanmin(values), np.nanmax(values), 100)

    triang = mtri.Triangulation(resin.header.x, resin.header.y, resin.header.ikle_2d -1)
    axes.tricontourf(triang, values, cmap='terrain', levels=levels)
    axes.set_aspect('equal')

    divider = make_axes_locatable(axes)
    cax = divider.append_axes('right', size='5%', pad=0.2)
    cmap = cm.ScalarMappable(cmap='terrain')
    cmap.set_array(levels)
    fig.colorbar(cmap, cax=cax)

    plt.show()