Esope source files

   1 : C APPFIB    SOURCE    CHAT      10/12/16    21:15:13     6819           
   2 : C 
   3 : C====&===1=========2=========3=========4=========5=========6=========7== 
   4 : C Commentaires : Subroutine permettant de définir la relation entre 
   5 : C                le taux de corrosion macroscopique TC et le taux d'éner- 
   6 : C                gie additionel à prendre en compte pour modifier l'endom- 
   7 : C                magement dans le cas du modèle INTIMP 
   8 : C                avec interface imparfaite 
   9 : C Auteurs      : B. Richard 
  10 : C Date         : Février 2008 
  11 : C====&===1=========2=========3=========4=========5=========6=========7== 
  12 :         SUBROUTINE APPFIB(X,Y,NDIM) 
  13 : C 
  14 : C-----DECLARATION GENERALE 
  15 : C 
  16 :         IMPLICIT INTEGER(I-N) 
  17 :       IMPLICIT REAL*8(A-H,O-Z) 
  18 : C 
  19 : C-----DECLARATION DES TYPE DE VARIABLES EN ARGUMENTS 
  20 : C 
  21 :         REAL*8 XTC(20),XYR(20) 
  22 : C 
  23 : C-----VARIABLES LOCALES 
  24 : C 
  25 :         INTEGER I 
  26 : C 
  27 : C-----DEFINITION DES TABLEAUX DE REFERENCE 
  28 : C 
  29 : CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
  30 :         XTC( 1) = 0.00D0 
  31 :         XTC( 2) = 1.62D0 
  32 :         XTC( 3) = 5.61D0 
  33 :         XTC( 4) = 8.58D0 
  34 :         XTC( 5) = 9.45D0 
  35 :         XTC( 6) = 11.58D0 
  36 :         XTC( 7) = 12.35D0 
  37 :         XTC( 8) = 12.66D0 
  38 :         XTC( 9) = 16.54D0 
  39 :         XTC(10) = 17.09D0 
  40 :         XTC(11) = 18.06D0 
  41 :         XTC(12) = 20.45D0 
  42 :         XTC(13) = 22.36D0 
  43 :         XTC(14) = 22.99D0 
  44 :         XTC(15) = 25.07D0 
  45 :  
  46 :         XYR( 1) =  0.00D0 
  47 :         XYR( 2) =  0.00D0 
  48 :         XYR( 3) =  5.00D3 
  49 :         XYR( 4) =  5.30D3 
  50 :         XYR( 5) =  6.50D3 
  51 :         XYR( 6) =  6.00D3 
  52 :         XYR( 7) =  5.00D3 
  53 :         XYR( 8) =  4.80D3 
  54 :         XYR( 9) = -2.00D3 
  55 :         XYR(10) = -4.00D3 
  56 :         XYR(11) = -1.00D4 
  57 :         XYR(12) = -2.50D4 
  58 :         XYR(13) = -7.00D4 
  59 :         XYR(14) = -5.00D5 
  60 :         XYR(15) = -1.00D6 
  61 : CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
  62 : C 
  63 : C-----TEST SI LE DEGRE DE CORROSION MAXI IDENTIFIE  
  64 : C       N EST PAS DEPASSE 
  65 : C 
  66 :         IF (X.GT.XTC(NDIM)) THEN 
  67 :                 Y = XYR(NDIM) 
  68 :                 RETURN 
  69 :         END IF 
  70 : C 
  71 : C-----INTERPOLATION LINEAIRE PAR MORCEAU 
  72 : C 
  73 :         DO I=1,(NDIM-1) 
  74 :                 IF ((XTC(I).LE.X).AND.(X.LE.XTC(I+1))) THEN 
  75 :                         Y=(XYR(I+1)-XYR(I))/(XTC(I+1)-XTC(I))*X+ 
  76 :      &        (XTC(I+1)*XYR(I)-XTC(I)*XYR(I+1))/(XTC(I+1)-XTC(I)) 
  77 :                         RETURN 
  78 :                 END IF 
  79 :         END DO 
  80 :  
  81 :       RETURN 
  82 :         END 
  83 :  
  84 :  
  85 :  
  86 :  
  87 :  
  88 :