Télécharger la plaquette Cast3M


fr en

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