semipublic.f 838 B

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