Télécharger la plaquette Cast3M


fr en

Télécharger typflu.eso

Retour à la liste

Numérotation des lignes : 

typflu

  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