def METHOD_NAME(interpolator):
    METHOD_NAME = gp.create_model('Test_uncomformities')
    # Importing the data from CSV-files and setting extent and resolution
    gp.init_data(METHOD_NAME, [0, 10., 0, 2., 0, 5.], [100, 3, 100],
                 path_o=input_path + '/05_toy_fold_unconformity_orientations.csv',
                 path_i=input_path + '/05_toy_fold_unconformity_interfaces.csv',
                 default_values=True)
    gp.map_stack_to_surfaces(METHOD_NAME,
                             {"Flat_Series": 'Flat',
                               "Inclined_Series": 'Inclined',
                               "Fold_Series": ('Basefold', 'Topfold', 'basement')})
    # Create the aesara model
    METHOD_NAME.set_aesara_function(interpolator)
    return METHOD_NAME