Sources Esope | Cast3M

byreto
C BYRETO    SOURCE    CHAT      05/01/12    21:44:37     5004
      SUBROUTINE BYRETO (XX,YX,NX,YT,XT,NT,KASN,ZZ)
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8 (A-H,O-Z)
      DIMENSION XT(NT),YT(NT),KASN(NT),ZZ(NT)
      DIMENSION XX(NX),YX(NX)
C
C     operateur FUITE
C>>>  interpolation de YX defini sur XX dans XT : resultat ZZ
C     XT et ZZ ordonne
C     ZZ-> YT selon la numerotation de la ligne obj.maillage                                                                   BTR00090
C            via le tableau d indirection KASN
 
      I=1
      ZZ(1)=YX(1)
      ZZ(NT)=YX(NX)
 
      DO 10 K=2,NT-1
 
 20      CONTINUE
         IF (XT(K).GE.XX(I)) THEN
            I=I+1
            GOTO 20
         ENDIF
 
         DX1=XT(K)-XX(I-1)
         DX2=XX(I)-XT(K)
         ZZ(K)=(DX1*YX(I)+DX2*YX(I-1))/(DX1+DX2)
 
  10  CONTINUE
 
      DO 11 K=1,NT
      J=KASN(K)
      YT(J)=ZZ(K)
  11  CONTINUE
 
      RETURN
      END