zlin1 zp(t)*(1+(Lp+Z)/(Ld-Lp)) Mesh deformation below piston
zlin2 zp(t)*(Ld-Z)/(Ld-Lp) Mesh deformation above piston