|
马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?注册
x
IF((NPATCH(1:4).EQ.'WDC3').AND.(INDVAR.EQ.C3)) THEN
L0VAL=L0F(VAL)
L0C3=L0F(C3)
L0C3O=L0F(OLD(C3))
LFDZGN=L0F(DZGNZ)+IZ
DO IX=1,NX
IADD=NY*(IX-1)
DO IY=1,NY
I=IY+IADD
LFVAL=L0VAL+I
LFC3O=L0C3O+I
F(LFVAL)=F(LFC3O)+F(LFC3O)*RG(4)*DT*(1-RG(3))
& /(F(LFDZGN)*(1-RG(4)*DT/F(LFDZGN)))
ENDDO
ENDDO
ENDIF |
|