def METHOD_NAME(units, uy_north):
    """Test that solution is independent of units.
    Parameters:
        units (dict): Dictionary with keys as those in
            :class:`~pp.models.material_constants.MaterialConstants`.
        uy_north (float): Value of displacement on the north boundary.
    """
    params = {
        "suppress_export": True,  # Suppress output for tests
        "fracture_indices": [0, 1],
        "cartesian": True,
        "uy_north": uy_north,
        "max_iterations": 10,
    }
    reference_params = copy.deepcopy(params)
    # Create model and run simulation
    setup_0 = LinearModel(reference_params)
    pp.run_time_dependent_model(setup_0, reference_params)
    params["units"] = pp.Units(**units)
    setup_1 = LinearModel(params)
    pp.run_time_dependent_model(setup_1, params)
    variables = [
        setup_1.displacement_variable,
        setup_1.interface_displacement_variable,
        setup_1.contact_traction_variable,
    ]
    variable_units = ["m", "m", "Pa"]
    compare_scaled_primary_variables(setup_0, setup_1, variables, variable_units)
    secondary_variables = ["stress", "displacement_jump"]
    secondary_units = ["Pa * m", "m"]
    domain_dimensions = [2, 1]
    compare_scaled_model_quantities(
        setup_0, setup_1, secondary_variables, secondary_units, domain_dimensions
    )