Télécharger typflu.eso

Retour à la liste

Numérotation des lignes :

  1. C TYPFLU SOURCE CHAT 05/01/13 03:57:44 5004
  2. C TYPFLU SOURCE AMAND 00/03/27
  3. SUBROUTINE TYPFLU(CODU,CODL,COD,BETJEF,BETFLU)
  4. IMPLICIT INTEGER(I-N)
  5. IMPLICIT REAL*8(A-H,O-Z)
  6. DIMENSION CODU(9,9),CODL(8,8),COD(8)
  7. C
  8. SEGMENT BETJEF
  9. REAL*8 AA,BETA,RB,ALPHA,EX,XNU,GFC,GFT,CAR,ETA,TDEF,
  10. & TCON,DPSTF1,DPSTF2,TETA,PDT,TP00
  11. INTEGER ICT,ICC,IMOD,IVIS,ITER,
  12. & ISIM,IBB,IGAU,IZON
  13. ENDSEGMENT
  14. SEGMENT BETFLU
  15. REAL*8 DATCOU,DATCUR,DATSEC,E28,PGTZO,PGDUR,TAU1,TAU2,
  16. & TP0,TZER
  17. INTEGER ITYPE,IMD,NBRC,NCOE,NTZERO,NTPS,IFOR
  18. ENDSEGMENT
  19. C
  20. C
  21. C VALEURS MAX DES PAS DE TEMPS
  22. C------------------------------
  23. MC = NBRC+1
  24. NTPS = 199
  25. NTZERO = 199
  26. NCOE = 8
  27. NC = NCOE+1
  28. C NTZERO1 = NTZERO+1
  29. C NTPS1 = NTPS+1
  30. TZER = 1.D0
  31. PGTZO = 1.05
  32. PGDUR = 1.05
  33. C
  34. C
  35. C DEFINITION DES COEF DE RELAX DES COURBES DE MODULES DE RAIDEUR
  36. C
  37. IF(IFOR.EQ.1) THEN
  38. TAU1 = 2.
  39. TAU2 = 20.
  40.  
  41. C----------------------------
  42. DO 7 I=1,NCOE
  43. IF(I.EQ.1)THEN
  44. COD(I)=1.
  45. ELSE IF(I.EQ.2)THEN
  46. COD(I)=0.1
  47. ELSE IF(I.EQ.3)THEN
  48. COD(I)=0.01
  49. ELSE IF(I.EQ.4)THEN
  50. COD(I)=0.001
  51. ELSE IF(I.EQ.5)THEN
  52. COD(I)=0.0001
  53. ELSE IF(I.EQ.6)THEN
  54. COD(I)=0.00001
  55. ELSE IF(I.EQ.7)THEN
  56. COD(I)=0.000001
  57. ELSE IF(I.EQ.8)THEN
  58. COD(I)=0.0000001
  59. ENDIF
  60. 7 CONTINUE
  61. C
  62. IF (CODU(1,1).EQ.0.D0.AND.CODU(9,9).EQ.0.D0) THEN
  63. CALL MDCRES(CODU,COD,BETJEF,BETFLU)
  64. ENDIF
  65. C
  66. C-----------------------------
  67. ELSE IF(IFOR.EQ.2) THEN
  68. C-----------------------------
  69. C
  70. TAU1 = 0.02
  71. TAU2 = 2.00
  72. DO 8 I = 1,NCOE
  73. IF(I.EQ.1)THEN
  74. COD(I) = 1.
  75. ELSE IF(I.EQ.2)THEN
  76. COD(I) = 0.1
  77. ELSE IF(I.EQ.3)THEN
  78. COD(I) = 0.01
  79. ELSE IF(I.EQ.4)THEN
  80. COD(I) = 0.001
  81. ELSE IF(I.EQ.5)THEN
  82. COD(I) = 0.0001
  83. ELSE IF(I.EQ.6)THEN
  84. COD(I) = 0.00001
  85. ELSE IF(I.EQ.7)THEN
  86. COD(I) = 0.000001
  87. ELSE IF(I.EQ.8)THEN
  88. COD(I) = 0.0000001
  89. ENDIF
  90. 8 CONTINUE
  91. C
  92. IF (CODL(1,1).EQ.0.D0.AND.CODL(8,8).EQ.0.D0) THEN
  93. CALL MDCREL(CODL,COD,BETJEF,BETFLU)
  94. ENDIF
  95. C
  96. C-----------------------------
  97. ENDIF
  98. C-----------------------------
  99. RETURN
  100. END
  101.  
  102.  
  103.  

© Cast3M 2003 - Tous droits réservés.
Mentions légales