changing code to be u-p formulation by Borja