import time

    from compas.datastructures import Mesh
    from compas.datastructures import mesh_unify_cycles

    from compas_plotters import MeshPlotter

    from streamlines.streamsystem import Streamsystem
    from streamlines.custom_mesh import StructuralMesh

    # ==========================================================================
    # Import mesh
    # ==========================================================================

    mesh = Mesh()
    mesh.load(HERE)
    mesh_unify_cycles(mesh)

    # ==========================================================================
    # Create PS vector lines
    # ==========================================================================

    vector_tag = 'ps_1_top'
    lines = vector_lines_on_faces(mesh, vector_tag, True, factor=0.05)
    lines = [line for line in map(line_tuple_to_dict, lines)]
    for line in lines:
        line['width'] = 0.60

    # ==========================================================================
    # Instantiate StructuralMesh()
    # ==========================================================================