subroutine DFLUX(FLUX,SOL,KSTEP,KINC,TIME,NOEL,NPT,COORDS,
     &                 JLTYP,TEMP,PRESS,SNAME)
C
      include 'ABA_PARAM.INC'
C
      dimension FLUX(2), TIME(2), COORDS(3)
      CHARACTER*80 SNAME

      X=COORDS(1)
	  Y=COORDS(2)
	  Z=COORDS(3)
	  
	  TotalTime = 16 !bei der aktuellen Bogenlaenge entspricht das 5mm/s
	  Dist = 0.0
	  R = 100
	  Bogenlsenge = 1/4*3.1415*R
	  Geschw = Bogenlsenge/TotalTime
	  Omega = (2*3.1415/(8*TotalTime))
	  Y_center = 0- R*cos(1*Omega*TIME(2))
	  X_center = 0.0 - R*sin(1*Omega*TIME(2))
	  
	  	  
	  YT=Y_center - Y
	  XT=X_center - X
	  Dist = SQRT((YT*YT)+(XT*XT))	
	  
	  QL=40000.d0                        !LASER POWER [mW/mm2]
      R0=1d0                            !BEAM RADIUS [mm]
      
      q=QL*EXP((XT**2+YT**2)/R0**2)    !Gauss FLUX [mW] ?
!     V=5.d0                              !LASER BEAM VELOCITY [mm/s]
	  
	  if (Dist .le. R0) then
	     FLUX(1) = q
	  else
	     FLUX(1) = 0.0
	  endif

      RETURN
      END