最近在学习写副程式遇到问题
麻烦各位帮忙检查,谢谢
原始码如下:
subroutine forcem(press,th1,th2,nn,n)
implicit real*8 (a-h, o-z)
dimension th1(3),th2(3),n(7)
include "C:\MSC.Software\Marc\2013.1.0\marc2013.1\common\creeps"
b=1.0d0
v=sqrt(2.0d0)
distan=v*(cptim+timinc)
xc=0.5d0*sqrt(2.0d0)*distan
yc=0.5d0*sqrt(2.0d0)*distan
xmin=xc-b/2.0d0
xmax=xc+b/2.0d0
ymin=yc-b/2.0d0
ymax=yc+b/2.0d0
press=0.0d0
if(th1(1).le.xmax.and.th1(1).ge.xmin.and.
th1(2).le.ymax.and.th1(2).ge.ymin)
^
出问题的地方
press=2.0d0
th2(1)=0.0d0
th2(2)=0.0d0
th2(3)=-1.0d0
return
end