def METHOD_NAME(mesh):
    "Return a function space U*H on mesh from the BFDM1 space."
    V = FunctionSpace(mesh, 'BDM', 1)  # Velocity space
    H = FunctionSpace(mesh, 'DG', 1)             # Height space
    W = V * H                                        # Mixed space of both.
    return W