Question: How do I switch an unknown which was previously set?

In the final line, my output is in the form of tau-c*z. How do I switch the tau-c*z back to 
zeta, because I want to integrate T with respect to zeta.
Thank you. 
> restart;
> with(plots);
> U := V(zeta)*exp(I*(k*tau-w*z));
> zeta := tau-c*z;
> Uz := diff(U, z);
> Utau := diff(U, tau);
> Utautau := diff(Utau, tau);
> P := convert(U, trig);
>
> A := abs(P)^2;
> B := evalc(A);
> Q := simplify(B, trig);

> I*Uz+alpha*U*Q+beta*Utautau = 0;
> R := simplify(%);
> S := R/exp(I*(k*tau-w*z));
> k := c/(2*beta);
> w := c^2/(4*beta)-theta;
> S;
> T := simplify(S);
>
Please Wait...