Put in fourth-order derivative for computing gradient of velocity field
fix merge conflicts