FORTRAN程序改写为VB程序

'在FORTRAN中二维数组PSI(i,j)下标都是从1开始

FOR J=2 TO 4

FOR I=2 TO J+3

DIJ=W*(PSI(I+1,J)+PSI(I-1,J)+PSI(I,J+1)+PSI(I,J-1))/4.0-PSI(I,J)

SDIJ1=SDIJ1+ABS(DIJ)

PSI(I,J)=PSI(I,J)+DIJ

NEXT I,J

IF ABS(SDIJ1-SDIJ0)>EPS THEN

SDIJ0=SDIJ1

GO TO 40

END IF

其实排除了注释和大量的common块之后，真正有效而且需要修改的地方不多。

大部分是个体力活。。。

用module代替common块之后，估计能省10%的代码。。。

goto这种的话，改成if条件判断其实不是很麻烦的~

而且这样一改，你对程序就很熟悉了。时间不会太长，专心改1个多星期我改完了3000行的f77代码。结果与原来的一致，速度也没啥影响。

---