Télécharger byreto.eso

Retour à la liste

Numérotation des lignes :

  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.  
  41.  

© Cast3M 2003 - Tous droits réservés.
Mentions légales