    div(dev(T(grad(U)))) Gauss linear;
    div(U) Gauss linear;
    div(Usup) Gauss linear;
    div(Ulift0) Gauss linear;
    div(Ulift1) Gauss linear;
    div(dev(T(grad(Ulift0)))) Gauss linear;
    div(dev(T(grad(Ulift1)))) Gauss linear;
    div(dev(T(grad(Usup)))) Gauss linear;
    div((interpolate(Ulift0)&S),U) bounded Gauss linearUpwind grad(U);
    div((interpolate(Ulift0)&S),Usup) bounded Gauss linearUpwind grad(Usup);
    div((interpolate(Ulift0)&S),Ulift0) bounded Gauss linearUpwind grad(Ulift0);

    div((interpolate(Ulift1)&S),U) bounded Gauss linearUpwind grad(U);
    div((interpolate(Ulift1)&S),Usup) bounded Gauss linearUpwind grad(Usup);
    div((interpolate(Ulift1)&S),Ulift1) bounded Gauss linearUpwind grad(Ulift1);

    div((interpolate(Ulift0)&S),Ulift1) bounded Gauss linearUpwind grad(Ulift1);
    div((interpolate(Ulift1)&S),Ulift0) bounded Gauss linearUpwind grad(Ulift0);

    div((interpolate(Usup)&S),Ulift0) bounded Gauss linearUpwind grad(Ulift0);
    div((interpolate(Usup)&S),Ulift1) bounded Gauss linearUpwind grad(Ulift1);
    div((interpolate(U)&S),U) bounded Gauss linearUpwind grad(U);
    div((interpolate(U)&S),Usup) bounded Gauss linearUpwind grad(Usup);
    div((interpolate(U)&S),Ulift0) bounded Gauss linearUpwind grad(Ulift0);
    div((interpolate(U)&S),Ulift1) bounded Gauss linearUpwind grad(Ulift1);
    div((interpolate(Usup)&S),U) bounded Gauss linearUpwind grad(U);
    div((interpolate(Usup)&S),Usup) bounded Gauss linearUpwind grad(Usup);
    div((nut*dev(grad(U).T()))) Gauss linear;
    div((nut*dev(grad(Ulift0).T()))) Gauss linear;
    div((nut*dev(grad(Ulift1).T()))) Gauss linear;
    div((interpolate(Ulift1)&S),Ulift0) bounded Gauss linearUpwind grad(Ulift0);
    div((interpolate(Ulift1)&S),Ulift1) bounded Gauss linearUpwind grad(Ulift1);
    div((nut*dev(grad(Usup).T()))) Gauss linear;
