solitary wave benchmark updates
comments
fixes in the melt_material_4 test
use the phase shift to calculate the shape error of the solitary wave
remove the boundary integral from the LHS and use the boundary integral on the RHS for a boundary condition grad p_f = rho_s g