关于fidap的subroutine的编写 请教

frank1981

以下是本人写的一个subroutine 是想计算一个依存于表面流体流速而变化的htransfer值 但是下面这个subr应该有很多的错误  请高手指教 谢谢 流速和htra的关系式是现成的htra=VEL**0.8 还有就是怎么把这个SUBR代回到主程序里呢 比如主程序的
写的是htra(set=1，subroutine) 好像也不对吧  很急
SUBROUTINE USRGAP (NELT,NE,NG,COEF,VARI,DVARI,NDFCD,LDOFU,SHP,DSDX,XYZL,PROP,TIME,NPTS,NDP,MNDP,IERR,ZIGN,ZIGT,GAP,TEMP2)
include "IMPLCT.COM"
include "PARUSR.COM"
      DIMENSION COEF(NPTS), GAP(NPTS), TEMP2(NPTS)
      DIMENSION SHP(NPTS,MNDP), DSDX(NPTS,NDFCD,MNDP), XYZL(NPTS,NDFCD)
      DIMENSION PROP(*), VARI(NPTS,*), DVARI(NPTS,NDFCD,*),LDOFU(*)
      DO 100 I=1, NPTS
UX=VARI(I,LDOFU(KDU))
UY=VARI(I,LDOFU(KDV))
UZ=VARI(I,LDOFU(KDW))
VEL=SQRT(UX*UX+UY*UY+UZ*UZ)
IF(VEL.GT.1)
htra=VEL**0.8
COEF(N)=100
IF(VEL.LE.1)
COEF(N)=10
100 CONTINUE
RETURN
END

frank1981

flunt公司的人给了我些建议 说是htra=VEL**0.8这个语句不能在subr里面用 而且主程序里stra(set=1,subr,reft=,,,,)的reft要做设定  我现在就是不知道htra=VEL**0.8这个语句到底怎么表达出来了  谁能给个建议阿。。。

chaobo1982

估计这里用fidap的人不多