def METHOD_NAME(self, cases, ID: int, nnodes: int, nelements: int,
                   surfaces: np.ndarray, is_cs: np.ndarray):
    #results_form = []
    geometry_form = [
        ('NodeID', 0, []),
        ('ElementID', 1, []),
        ('SurfaceID', 2, []),
        ('isControlSurface', 3, []),
    ]
    nids = np.arange(1, nnodes + 1)
    eids = np.arange(1, nelements + 1)
    assert len(eids) == nelements, len(eids)
    nid_res = GuiResult(ID, 'NodeID', 'NodeID', 'node', nids)
    eid_res = GuiResult(ID, 'ElementID', 'ElementID', 'centroid', eids)
    surface_res = GuiResult(ID, 'SurfaceID', 'SurfaceID', 'centroid', surfaces)
    cs_res = GuiResult(ID, 'isControlSurface', 'isControlSurface', 'centroid', is_cs)
    i = 0
    cases[i] = (nid_res, (0, 'NodeID'))
    cases[i + 1] = (eid_res, (0, 'ElementID'))
    cases[i + 2] = (surface_res, (0, 'SurfaceID'))
    cases[i + 3] = (cs_res, (0, 'isControlSurface'))
    form = [
        ('Geometry', None, geometry_form),
    ]
    return form, cases, nids, eids