Question

有没有FORTRAN高手帮我弄一下这程序，我这个程序用compaq visual fortran一编译就报错

我是小白，以前没用过fortran，这个程序太长，信得过我的话，加我好友QQ381762478,我将全部程序给你，麻烦帮我改一下，我连那个报错都看不动。。。
还有，我听说compaq visual fortran用的是fortran90，我不知道我那个程序是哪个版本的，有可能是fortran 66的，这有影响吗？下面是部分程序
DIMENSION R(21),H(21),SIGR(2,21),ALF(21),E(21),T(21),AN(2,2,19),
*AM(2,2,20),XM(20),A(20),B(20),C(2,20),CL(2,20),AMN(2,2,20),
*AI(2,2,20),BI(2,20),AIP(2,2,20)BIP(2,20),RE(22),HE(22),DD(2,21),
*SIGE(2,21),SIGP(2,21),HH(21)
DOUBLE PRECISION R,H,T,ALF,E,SIGR,AN,AM,XM,A,B,C,ALFA,EE,CL,F,
*AMN,AI,BI,V,GAMA,NN,DELTA,XO,YO,HH
C CHARACTER*(69)ZS
OPEN(5,FILE='CIN.DATA')
READ(5,*)N
C READ(5,*)ZS
READ(5,*)(R(I),HH(I),E(I),T(I),ALF(I),I=1,N),V,GAMA,NN
CLOSE(5)
WRITE(6,400)
DO 10 I=1,N
10 WRITE(6,420) R(I),HH(I),E(I),T(I),ALF(I)
WRITE(6,430) V,GAMA,NN
WRITE(6,440)
N1=N-1
DO 20 I=1,N1
XM(I)=R(1)/R(I+1)
H(I)=(HH(I)+HH(I+1))/2
A(I)=0.5*(1.0+XM(I)*XM(I))
B(I)=0.5*(1.0-XM(I)*XM(I))
C(1,I)=0.125*(2.0*(1.0+V)*XM(I)*XM(I)+(1.0-V)*XM(I)**4-3.0-V)
C(2,I)=0.125*(2.0*(1.0+V)*XM(I)*XM(I)-(1.0-V)*XM(I)**4-1.0-3.0-V)
CL(1,I)=C(1,I)*GAMA*(2.0*3.1415926*NN/60.0)**2*R(I+1)*R(I+1)
20 CL(2,I)=C(2,I)*GAMA*(2.0*3.1415926*NN/60.0)**2*R(I+1)*R(I+1)
DO 25 I1=1,2
25 BI(I1,1)=CL(I1,1)
AI(1,1,1)=A(1)
AI(1,2,1)=B(1)
AI(2,1,1)=B(1)
AI(2,2,1)=A(1)
I=1
100 AM(1,1,I+1)=A(I+I)
AM(1,2,I+1)=B(I+I)
AM(2,1,I+1)=B(I+I)
AM(2,2,I+1)=A(I+I)
AN(1,1,I)=H(I)/H(I+1)
AN(1,2,I)=0.0
AN(2,1,I)=V*(H(I)/H(I+1)-E(I+1)/E(I))
AN(2,2,I)=E(I+1)/E(I)
DD(2,I)=E(I+1)*(ALF(I)*T(I)-ALF(I+1)*T(T+1))
DD(1,I)=0.0
DO 30 I1=1,2
BIP(I1,I)=0.0
DO 30 I2=1,2
AIP(I1,I2,I)=0.0
DO 30 K=1,2
30 AIP (I1,I2,I)=AIP(I1,I2,I)+AN(I1,K,I)*AI(K,I2,I)
DO 40 I1=1,2
DO 35 I2=1,2
35 BIP(I1,I)=BIP(I1,I)+AN(I1,I2,I)*BI(I2,I)
40 BIP(I1,I)=BIP(I1,I)+DD(I1,I)
DO 42 J1=1,2
BI(J1,I+1)0.0
DO 42 J2=1,2
AI(J1,J2,I+1)=0.0
DO 42 K=1,2
42 AI(J1,J2,I+1)=AI(J1,J2,I+I)+AM(J1,K,I+1)*AIP(K,J2,I)
DO 46 J1=1,2
DO 44 J2=1,2
44 BI(J1,I+1)=BI(J1,I+1)+AM(J1,J2,I+1)*BIP(J2,I)
46 BI(J1,I+1)=BI(J1,I+1)+CL(J1,I+1)
I=I+1
IF(I.LT.N1)GOTO 100
IF(R(1).EQ.0.0) GOTO 300
READ (5,*) XO,YO

Answer 1

AIP(2,2,20)BIP(2,20) 在bip前面加逗号
第43行 最后一个括号里边T+1改成 i+1
第56行 BI(J1,I+1)0.0 改成 BI(J1,I+1)=0.0

没有完整程序，别的就没法改了

现在编译器都是向下兼容的，fortran 77 66什么的都没问题

追问

能加下好友吗，我一编译他有30个错误，我把全部程序跟你，大概200多行