潮流程序设计fortran版.docVIP

$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),PVA(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)=PVA(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),PVA(30,2) /NNOO/NEW(150),NOLD(50,2) COMMON /GBDA/GBII(50,4),GBIJ(80,3) OPEN(7,FILE=lf02.DAT,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,*)(PVA(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 K=PVA(I,1) 68 NEW(K)=-2 NEW(NS)=-3 KL=1 DO 338 I=1,3*NN IF(NEW(I).NE.-1) GOTO 338 NEW(I)=KL KL=KL+1 338 CONTINUE DO 348 I=1,3*NN IF(NEW(I).NE.-2) GOTO 348 NEW(I)=KL KL=KL+1 348 CONTINUE NEW(NS)=KL DO 358 I=1,3*NN IF(NEW(I).EQ.0) GOTO 358 K=NEW(I) NOLD(K,1)=I 358 CONTINUE DO 78 I=1,NB K=ABS(ZB(I,1)) L=ABS(ZB(I,2)) ZB(I,1)=ZB(I,1)/K*NEW(K) 78 ZB(I,2)=ZB(I,2)/L*NEW(L) DO 98 I=1,NN K=DABS(ZN(I,1)) 98 ZN(I