PPT Slide
########################################################################
########################################################################
operator D_FW(f,t) := (ə>f[0][0] - ɘ>f[0][0])/dt
operator D_BW(f,t) := (ɘ>f[0][0] - <-1>f[0][0])/dt
operator D_LF(f,t,t) := (ə>f[0][0] - 2*ɘ>f[0][0] + <-1>f[0][0])/(dt^2)
operator D_LF(f,rho,rho):=(ɘ>f[1][0]-2*ɘ>f[0][0]+ɘ>f[-1][0])/(drho^2)
operator D_LF(f,z,z) := (ɘ>f[0][1] - 2*ɘ>f[0][0] + ɘ>f[0][-1])/(dz^2)
.......................................................................
########################################################################
# Residual definitions (equations of motion)
########################################################################
[1:1] [1:Nz] := QFIT(Phi,rho);
ɘ>r[0][0]*D_BW2(Phi,t)+rho*D_BW2(Phi,rho)+z*D_CADV(Phi,z);
D_LF(Phi,t,t)=D_LF(Phi,z,z)+D_LF(Phi,rho,rho)+D_LF(Phi,rho)/rho;
ɘ>r[0][0]*D_BW2(Phi,t)+rho*D_CADV(Phi,rho)+z*D_FW2(Phi,z);
ɘ>r[0][0]*D_BW2(Phi,t)+rho*D_CADV(Phi,rho)+z*D_BW2(Phi,z);