Esope source files

   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     
  62 :       CALL MDCRES(CODU,COD,BETJEF,BETFLU)
  63 :       ENDIF
  64 : C
  65 : C-----------------------------
  66 :       ELSE IF(IFOR.EQ.2) THEN
  67 : C-----------------------------
  68 : C
  69 :       TAU1 = 0.02
  70 :       TAU2 = 2.00
  71 :       DO 8 I = 1,NCOE
  72 :        IF(I.EQ.1)THEN
  73 :        COD(I) = 1.
  74 :         ELSE IF(I.EQ.2)THEN
  75 :         COD(I) = 0.1
  76 :          ELSE IF(I.EQ.3)THEN
  77 :          COD(I) = 0.01
  78 :           ELSE IF(I.EQ.4)THEN
  79 :           COD(I) = 0.001
  80 :            ELSE IF(I.EQ.5)THEN
  81 :            COD(I) = 0.0001
  82 :           ELSE IF(I.EQ.6)THEN
  83 :           COD(I) = 0.00001
  84 :          ELSE IF(I.EQ.7)THEN
  85 :          COD(I) = 0.000001
  86 :         ELSE IF(I.EQ.8)THEN
  87 :         COD(I) = 0.0000001
  88 :        ENDIF          
  89 :     8 CONTINUE 
  90 : C
  91 :       IF (CODL(1,1).EQ.0.D0.AND.CODL(8,8).EQ.0.D0) THEN       
  91 :       CALL MDCREL(CODL,COD,BETJEF,BETFLU)
  92 :       ENDIF      
  93 : C
  94 : C-----------------------------
  95 :       ENDIF
  96 : C-----------------------------
  97 :       RETURN
  98 :       END       
  99 :  
 100 :