|
马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?注册
x
C For SOLVE(TEM1)-------- phase-1 specific heat
C Property name: PRPT01
IF(ISTEP.GE.1 .AND.ISTEP.LE.LSTEP ) THEN
IF(IZ.GE.1 .AND.IZ.LE.NZ ) THEN
LFISPH =L0F(ISPH1 )
LFTEM1=L0F(INAME('TEM1 '))
DO 91402 IX=IXF ,IXL
IADD=NY*(IX-1)
DO 91402 IY=IYF ,IYL
I=IY+IADD
L0ISPH=LFISPH+I
L0TEM1=LFTEM1+I
IF(F(L0TEM1).LE.973)THEN
F(L0ISPH )=513.76-0.33504*F(L0TEM1)+6.8929E-04*F(L0TEM1)**2
ELSE IF(F(L0TEM1).GE.973.AND.F(L0TEM1).LE.1023)THEN
F(L0ISPH )=-10539+11.7*F(L0TEM1)
ELSE IF(F(L0TEM1).GE.1023.AND.F(L0TEM1).LE.1100)THEN
F(L0ISPH )=11873-10.208*F(L0TEM1)
ELSE IF(F(L0TEM1).GE.1100.AND.F(L0TEM1).LE.1379)THEN
F(L0ISPH )=644
ELSE
F(L0ISPH )=354.34+0.21*F(L0TEM1)
ENDIF
91402 CONTINUE
ENDIF
ENDIF
RETURN
905 CONTINUE
这是GROUP9里的一个比热的程序
LFISPH =L0F(ISPH1 )
LFTEM1=L0F(INAME('TEM1 '))
这个代表啥?
还有关于F()是怎么回事??
急!!
谢谢!!!1 |
|