潮流程序设计fortran版实用.pdfVIP

$DEBUG $LARGE c for 2002 studens implicit real*8 (a-h,p-z) REAL*8 V(2,50),P(2,50),F(2,50) COMMON NN,NB,NG ,NL,NPV,NGL,NT,ES,AXM,KKK /VPF/V ,P,F 1 /BNDA/ZB(80,5),ZN(50,7),PV A(30,2) /GBDA/GBII(50,4), 1 GBIJ(80,3) /FACT/JF(2,100),DJ(100),UJ(2,2000) OPEN(1,FILE=lf02.res) CALL INPUT(SB,U0,REFA) CALL YM CALL DEVR(U0) DO 18 I=1,NPV K=PVA(I,1) 18 V(1,K)=PV A(I,2) VS=V(1,NN) V(1,NN)=VS*DCOS(REFA/180*3.1415926) V(2,NN)=VS*DSIN(REFA/180*3.1415926) KT=0 27 CALL PQ(1,KT) IF(KT.GT.20) STOP IF(DABS(AXM).LE.ES) GOTO 58 CALL JACO(1) CALL SBW(1,1) KT=KT+1 GOTO 27 58 CONTINUE WRITE(1,*)ITERATING NUMBER: ,KT WRITE(1,7)MAXMUM ERROR,F(1,KKK),F(2,KKK), AT NODE,KKK 7 FORMAT(1X,A,2F18.11,A,I5) CALL PQ(2,1) close(1) STOP END * SUBROUTINE INPUT(SB,U0,REFA) implicit real*8 (a-h,p-z) REAL*8 V(2,50),P(2,50),F(2,50) COMMON NN,NB,NG ,NL,NPV,NGL,NT,ES,AXM,KKK /VPF/V ,P,F 1 /BNDA/ZB(80,5),ZN(50,7),PV A(30,2) /NNOO/NEW(150),NOLD(50,2) COMMON /GBDA/GBII(50,4),GBIJ(80,3) OPEN(7,FILE=lf02.DA T,STATUS=OLD) REWIND(7) READ(7,*)NN,NB,NPV ,NS WRITE(1,*)NN,NB,NPV,NS: WRITE(1,1)NN,NB,NPV ,NS WRITE(1,*)ES,U0,SB,REFA: READ(7,*)ES,U0,SB,REFA WRITE(1,2)ES,U0,SB,REFA DO 8 I=1,NB READ(7,*)(ZB(I,J),J=1,5) 8 CONTINUE DO 18 I=1,NN 18 READ(7,*)(ZN(I,J),J=1,5) DO 38 I=1,NPV 38 READ(7,*)(PV A(I,J),J=1,2) DO 48 I=1,3*NN 48 NEW(I)=0 DO 58 I=1,NB K=DABS(ZB(I,1)) L=DABS(ZB(I,2)) NEW(K)=-1 58 NEW(L)=-1 DO 68 I=1,NPV