Esope source files

   1 : C BYRETO    SOURCE    CHAT      05/01/12    21:44:37     5004           
   2 :       SUBROUTINE BYRETO (XX,YX,NX,YT,XT,NT,KASN,ZZ)
   3 :       IMPLICIT INTEGER(I-N)
   4 :       IMPLICIT REAL*8 (A-H,O-Z)
   5 :       DIMENSION XT(NT),YT(NT),KASN(NT),ZZ(NT)
   6 :       DIMENSION XX(NX),YX(NX)
   7 : C
   8 : C     operateur FUITE
   9 : C>>>  interpolation de YX defini sur XX dans XT : resultat ZZ
  10 : C     XT et ZZ ordonne
  11 : C     ZZ-> YT selon la numerotation de la ligne obj.maillage                                                                   BTR00090
  12 : C            via le tableau d indirection KASN
  13 : 
  14 :       I=1
  15 :       ZZ(1)=YX(1)
  16 :       ZZ(NT)=YX(NX)
  17 : 
  18 :       DO 10 K=2,NT-1
  19 : 
  20 :  20      CONTINUE
  21 :          IF (XT(K).GE.XX(I)) THEN
  22 :             I=I+1
  23 :             GOTO 20
  24 :          ENDIF
  25 : 
  26 :          DX1=XT(K)-XX(I-1)
  27 :          DX2=XX(I)-XT(K)
  28 :          ZZ(K)=(DX1*YX(I)+DX2*YX(I-1))/(DX1+DX2)
  29 : 
  30 :   10  CONTINUE
  31 : 
  32 :       DO 11 K=1,NT
  33 :       J=KASN(K)
  34 :       YT(J)=ZZ(K)
  35 :   11  CONTINUE
  36 : 
  37 :       RETURN
  38 :       END
  39 :  
  40 :