DFLUX子程序环形路径坐标变换

wulan

      vr=v/r
      phi=vr*TIME(2)

      x=COORDS(1)
      y=COORDS(2)
      z=COORDS(3)

      IF(z .GE. 0.0) THEN
      theta=acos(x/sqrt(x*x+z*z))
      ELSE
      theta=2.0*PI-acos(x/sqrt(x*x+z*z))
      ENDIF

      deta=theta-phi
      x1=sqrt(x*x+z*z)*cos(deta)
      y1=y
      z1=sqrt(x*x+z*z)*sin(deta)
我的理解是PHI为热源走过弧长所对应的圆心角，想请教下theta和deta表示哪些角度呢？