Movement Heat Flux using DFLUX subroutine

Dear All,

I am currently working on Heat movement using DFLUX subroutine.
I wrote a subroutine using double ellipsoidal Gaussian distribution.
Unfortunately, results are false.

Attached my model and my subroutine.
Anyone can tell me what is the possible error?

Many Thanks;
Rami


-- 
*Cordialement;*

*Rami GHORBEL*

*Material Engineer*

*Phd Student*

*Phone :* (00216) 22 667 771

  ----------

      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
      Real t, x, y, z, x1, y1, z1, Lweld, v
      Real Q, a, b, c, pi, x0, y0, z0, heat, shape1
C time
      t=TIME(2)
      a=0.3
      b=0.5
      c=0.2
      cf=0.003
	pi=3.141592
	
      FLUX(1)=0
      FLUX(2)=0
      Lweld=60
      v=5
	x=COORDS(1)
	y=COORDS(2)
	z=COORDS(3)
c         try
      if (JSTEP. eq. 1) then
      
	x0=0
      y0=0
      z0=2
      Q=0.5*10*900
      x1=x-x0
	y1=y-v*t
	z1=z-z0

      heat=6*sqrt(3.0)*Q/(a*b*c*pi*sqrt(pi))
      shape1=exp(-3*(x1)**2/a**2-3*(y1)**2/b**2-3*(z1)**2/c**2)
      
      FLUX(1)=heat*shape1
      
      end if
      
      RETURN
      END


[Non-text portions of this message have been removed]