SUBROUTINE DFLUX(FLUX,SOL,JSTEP,JINC,TIME,NOEL,NPT,COORDS,JLTYP,

     1                 TEMP,PRESS,SNAME)

C

      INCLUDE 'ABA_PARAM.INC'

      DIMENSION COORDS(3),FLUX(2),TIME(2)

      CHARACTER*80 SNAME

      wu=220.0

      wi=1.0

      effi=0.2

      v=0.004

      r=0.0098

      q=wu*wi*effi

      vr=v/r

      phi=vr*TIME(2)

      x=COORDS(1)

      y=COORDS(2)

      z=COORDS(3)

      x0=r

      y0=0

      z0=0   

      a=0.004

      b=0.002

      c=0.002

      aa=0.001

      f1=1.0

      PI=3.1415926

      beta=0.0

      aacos=a/cos(beta)

      cccos=c*cos(beta)

      aaacos=aa/cos(beta)

      cccos=c*cos(beta)

      heat1=6.0*sqrt(3.0)*q/(aa*b*c*PI*sqrt(PI))*f1

      heat2=6.0*sqrt(3.0)*q/(a*b*c*PI*sqrt(PI))*(2.0-f1)

      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表示哪些角度呢？