|
马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
以下是本人写的一个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 |
|