def METHOD_NAME(self, avl_filename: str,
                      name: str='main', plot: bool=True) -> None:
    model_name = name
    #key = self.case_keys[self.icase]
    #case = self.result_cases[key]
    #skip_reading = self._remove_old_adb_geometry(adb_filename)
    #if skip_reading:
        #return
    log = self.gui.log
    model = read_avl(avl_filename, log=log, debug=False)
    self.gui.model_type = 'avl'
    nodes, elements, line_elements, surfaces, is_cs = model.get_nodes_elements()
    #self.model_type = model.model_type
    nxyz_nodes = nodes.shape[0]
    nquad_elements = elements.shape[0]
    nline_elements = 0
    if line_elements:
        nline_elements = line_elements.shape[0]
        assert nline_elements > 0, nline_elements
    nnodes = nxyz_nodes
    nelements = nquad_elements + nline_elements
    self.gui.nnodes = nnodes
    self.gui.nelements = nelements
    grid = self.gui.grid
    grid.Allocate(self.gui.nelements, 1000)
    points = vtkPoints()
    points.SetNumberOfPoints(self.gui.nnodes)
    #vectorReselt.SetNumberOfComponents(3)
    self.gui.nid_map = {}
    assert nodes is not None
    nid = 0
    #print('nxyz_nodes=%s' % nxyz_nodes)
    mmax = np.amax(nodes, axis=0)
    mmin = np.amin(nodes, axis=0)
    dim_max = (mmax - mmin).max()
    self.gui.create_global_axes(dim_max)
    for i in range(nxyz_nodes):
        points.InsertPoint(nid, nodes[i, :])
        nid += 1
    log.info('nnodes=%s nquad_elements=%s nline_elements=%s' % (
        nxyz_nodes, nquad_elements, nline_elements))
    #elements -= 1
    for eid in range(nquad_elements):
        elem = vtkQuad()
        node_ids = elements[eid, :]
        point_ids = elem.GetPointIds()
        point_ids.SetId(0, node_ids[0])
        point_ids.SetId(1, node_ids[1])
        point_ids.SetId(2, node_ids[2])
        point_ids.SetId(3, node_ids[3])
        grid.InsertNextCell(elem.GetCellType(), point_ids)
    for eid in range(nline_elements):
        elem = vtkLine()
        node_ids = line_elements[eid, :]
        point_ids = elem.GetPointIds()
        point_ids.SetId(0, node_ids[0])
        point_ids.SetId(1, node_ids[1])
        grid.InsertNextCell(elem.GetCellType(), point_ids)
    grid.SetPoints(points)
    grid.Modified()
    # load results - regions/loads
    self.gui.scalar_bar_actor.VisibilityOn()
    self.gui.scalar_bar_actor.Modified()
    note = ''
    self.gui.isubcase_name_map = {1: ['AVL%s' % note, '']}
    cases = {}
    ID = 1
    form, cases, node_ids, element_ids = self._fill_avl_case(
        cases, ID, nnodes, nelements, surfaces, is_cs)
    self.gui.node_ids = node_ids
    self.gui.element_ids = element_ids
    self.gui._finish_results_io2(model_name, form, cases)