热源子程序编译
浏览:120338 回答:1
为什么在fortran编译时出现下列情况,我的子程序是按照例子上面输的,格式也没错。
fortcom: Error: Illegal character in statement label field [颹
fortcom: Error: Illegal character in statement label field [籡
fortcom: Error: Illegal character in statement label field [縘。
下面是我输入的子程序。
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.005
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
fortcom: Error: Illegal character in statement label field [颹
fortcom: Error: Illegal character in statement label field [籡
fortcom: Error: Illegal character in statement label field [縘。
下面是我输入的子程序。
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.005
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