added TSL velocity gradient calculation, fixed TSL flux jacobian