123456789101112131415161718192021222324252627282930313233343536373839 |
- SUBROUTINE SEMIPUBLIC (PRFT,PRFR,TC,USCT,UCET,USET,FSC,FCE,AA,
- 1 PSC,PCE,BB,TAUSC,TAUCE,TAUSE)
- IMPLICIT REAL*8 (A-H,O-Z)
- DIMENSION W(3000), IACT(51)
- INCLUDE 'param.inc'
- NMAX = 50
- EPS = 1.0D-6
- DO 10, I=1,NMAX
-
- USCT0 = USCT
- UCET0 = UCET
- USET0 = USET
- C PRINT '(9F8.3)', USCT, UCET,USET, TAUSC, TAUCE, TAUSE,PSC,PCE,BB
- CALL PRFROAD (PRFR,TC,USCT,UCET,USET,FSC,FCE,AA,
- 1 PSC,PCE,BB,TAUSC,TAUCE,TAUSE)
- ID = 10
- CALL DECENTRALIZE (FSC,FCE,AA,PSC,PCE,BB,
- 1 TAUSC,TAUCE,TAUSE,USCT,UCET,USET,TC,ID)
- DU = (USCT-USCT0)**2+(UCET-UCET0)**2+(USET-USET0)**2
-
- IF (DU.LT.EPS) GOTO 15
- 10 CONTINUE
- PRINT *, 'No convergence after ', NMAX, ' iterations'
- 15 CONTINUE
- RETURN
- END
|