求大神帮我看一下焊接的子程序，编译时老是出错。
浏览：14694 回答：2
SUBROUTINE DFLUX(FLUX,SOL,JSTEP,JINC,TIME,NOEL,NPT,COORDS,JLTYP,
1 TEMP,PRESS,SNAME)
C
INCLUDE 'ABA_PARAM.INC'
parameter(one=1.d0)
DIMENSION COORDS(3),FLUX(2),TIME(2)
CHARACTER*80 SNAME
q=800
v=0.05
d=v*TIME(2)
x=COORDS(1)
y=COORDS(2)
z=COORDS(3)
x0=0
y0=0
z0=0.0025
a=0.0028
b=0.0034
c=0.0032
PI=3.1415
heat=6*sqrt(3.0）*q/(a*b*c*PI*sqrt(PI))
shape=exp(-3*(x-x0)**2/c**2-3*(y-y0-d)**2/a**2-3*(z-z0)**2/b**2)
c JLTYP=1,表示为体热源
JLTYP=1
if (JSTEP.eq.one) then
FLUX(1)=heat*shape
endif
RETURN
END
![J{W[@V0JEGJ_@GUVN`0)]DL.png](http://img.jishulink.com/caenet/forums/upload/2015/05/25/106/198628637507665.png "J{W[@V0JEGJ_@GUVN`0)]DL.png")
1 TEMP,PRESS,SNAME)
C
INCLUDE 'ABA_PARAM.INC'
parameter(one=1.d0)
DIMENSION COORDS(3),FLUX(2),TIME(2)
CHARACTER*80 SNAME
q=800
v=0.05
d=v*TIME(2)
x=COORDS(1)
y=COORDS(2)
z=COORDS(3)
x0=0
y0=0
z0=0.0025
a=0.0028
b=0.0034
c=0.0032
PI=3.1415
heat=6*sqrt(3.0）*q/(a*b*c*PI*sqrt(PI))
shape=exp(-3*(x-x0)**2/c**2-3*(y-y0-d)**2/a**2-3*(z-z0)**2/b**2)
c JLTYP=1,表示为体热源
JLTYP=1
if (JSTEP.eq.one) then
FLUX(1)=heat*shape
endif
RETURN
END
adiou 
