Télécharger coupel.eso

Retour à la liste

Numérotation des lignes :

coupel
  1. C COUPEL SOURCE PV 20/03/24 21:16:14 10554
  2. SUBROUTINE COUPEL(XYZ1,XYZ2,XYZ3,DIS1,DIS2,DIS3,
  3. > NUM0,MELEME,IELEM,JCOLOR)
  4. C Esclave de COUPEK
  5. C Ajout d'un element compose du noeud NUM0
  6. C de l'interpolation entre XYZ1 et XYZ2 avec les poids
  7. C DIS1 et DIS2 et de de l'interpolation entre XYZ1 et XYZ3
  8. C avec les poids DIS1 et DIS3
  9. C
  10. IMPLICIT INTEGER(I-N)
  11. IMPLICIT REAL*8(A-H,O-Z)
  12.  
  13. -INC PPARAM
  14. -INC CCOPTIO
  15. -INC SMCOORD
  16. -INC SMELEME
  17. C
  18. DIMENSION XYZ1(4),XYZ2(4),XYZ3(4)
  19. C
  20. segact mcoord*mod
  21. NBPTI=nbpts
  22. NBPTS=NBPTI+2
  23. SEGADJ,MCOORD
  24. C
  25. COEFF=ABS(DIS1)/(ABS(DIS1)+ABS(DIS2))
  26. NBPTI=NBPTI+1
  27. IREF=(NBPTI-1)*4
  28. DO IE1=1,4
  29. XCOOR(IREF+IE1)=XYZ1(IE1)+COEFF*(XYZ2(IE1)-XYZ1(IE1))
  30. ENDDO
  31. C
  32. COEFF=ABS(DIS1)/(ABS(DIS1)+ABS(DIS3))
  33. NBPTI=NBPTI+1
  34. IREF=(NBPTI-1)*4
  35. DO IE1=1,4
  36. XCOOR(IREF+IE1)=XYZ1(IE1)+COEFF*(XYZ3(IE1)-XYZ1(IE1))
  37. ENDDO
  38. C
  39. NBPTI=NBPTI-2
  40. IELEM=IELEM+1
  41. NUM(1,IELEM)=NUM0
  42. NUM(2,IELEM)=NBPTI+1
  43. NUM(3,IELEM)=NBPTI+2
  44. ICOLOR(IELEM)=JCOLOR
  45. C
  46. RETURN
  47. END
  48.  
  49.  
  50.  
  51.  

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