小编naf*_*7me的帖子

这是我解决三维圆柱形翅片问题的方案.但是当我在Fortran 6.2中运行此程序时,会出现以下错误:

错误:数组表达式的形状不符合.[T]"

现在我不明白为什么会这样.我需要快速帮助.有人请帮助我.

        PROGRAM CYLINDRICAL FIN
    DIMENSION T(500,500,500), OLDT(500,500,500),ERR(500,500,500)
    DIMENSION R(500),Y(500),Z(500)

    REAL CC,H,DR,DY,DZ,A,D,RY,YZ,ZR,B,E,TA,RL,YL,ZL
     +P,AC,QF,MF,Q,EF,EFF,QMAX,AS,LC,QT,QF1,QF2
    INTEGER I,J,K,M,N,L,M1,N1,L1,M2,M4,M34,N2,N4,N34,L2,L4,L34

    RL=0.5
    YL=6.283
    ZL=0.04
    M=100
    N=40
    L=20
    M2=((M/2)+1)
    M4=((M/4)+1)
    M34=((3*M/4)+1)
    N2=((N/2)+1)
    N4=((N/4)+1)
    N34=((3*N/4)+1)
    L2=((L/2)+1)
    L4=((L/4)+1)
    L34=((3*L/4)+1)



    DR=RL/M
    DY=YL/N
    DZ=ZL/L

    CC=400.0
    H=10.0
    TA=25
    M1=M-1
    N1=N-1
    L1=L-1

************VARIABLES************
    A=DR*DY*DZ
    D=DR+DY+DZ
    RY=DR*DY
    YZ=DY*DZ
    ZR=DZ*DR
      E=RY+YZ+ZR

************VARIABLES FOR EFFICIENCY AND EFFECTIVENESS (CROSS-SECTION AREA,PERIMETER,M,SURFACE AREA OF FIN)************
    AC=3.1416*DR**2
    P=2*(3.1416*DR+DZ)
    MF=((H*P)/(CC*AC))**(0.5)
    AS=2*3.1416*DR*DZ+3.1416*DR**2
**************************************  distance discritization ******************

    R(1)=0.0
    Y(1)=0.0
    Z(1)=0.0


    R(M+1)=RL
    Y(N+1)=YL
    Z(L+1)=ZL

    DO I=2,M
    R(I)=R(I-1)+DR
    END DO 
    DO J=2,N
    Y(J)=Y(J-1)+DY …