def METHOD_NAME(interpolator):
    model = gempy.create_data(extent=[0, 2500, 0, 1000, 0, 1000],
                              resolution=[50, 20, 20],
                              path_o=input_path2 + "jan_models/fixture_model_orientations.csv",
                              path_i=input_path2 + "jan_models/fixture_model_surfaces.csv")
    # Assigning series to surface as well as their order (timewise)
    gp.map_stack_to_surfaces(model,
                             {"Fault_Series": ('fault'), "Strat_Series1": ('rock3'),
                              "Strat_Series2": ('rock2', 'rock1'),
                              "Basement_Series": ('basement')})
    model.set_is_fault(['Fault_Series'])
    model.set_aesara_function(interpolator)
    return model