Télécharger la plaquette Cast3M


fr en

Télécharger trjral.eso

Retour à la liste

Numérotation des lignes : 

trjral

  1. C TRJRAL    SOURCE    CHAT      05/01/13    03:50:51     5004
  2.       SUBROUTINE TRJRAL(IPARPO,IZLAR,MELEME,IZSH)
  3. C
  4. C
  5. C
  6. CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
  7. C
  8. C     PASSE DES COORDONNEES DE REFERENCES AUX COORDONNEES REELLES
  9. C
  10. C     IPARPO SEGMENT DES COORDONNEES DE REFERENCES
  11. C     IZLAR SEGMENT DES COORDONNEES REELLES
  12. C     MELEME MAILLAGE
  13. C     IZSH SEGMENT DE TRAVAIL POUR LE CALCUL DES FONCTIONS DE FORME
  14. C
  15. C     SHP FONCTION DE FORME
  16. C     XYZL COORDONNEES REELLES DES NOEUDS DE L ELEMENT
  17. C     XR COORDONNEE DE LA PARTICULE DANS L ELEMENT DE REFERENCE
  18. C
  19. CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
  20. C
  21. C
  22.       IMPLICIT INTEGER(I-N)
  23.       IMPLICIT REAL*8 (A-H,O-Z)
  24. C
  25. C
  26.  
  27. -INC PPARAM
  28. -INC CCOPTIO
  29. -INC SMCOORD
  30. -INC SMELEME
  31. C
  32.       SEGMENT  IPARPO
  33.         INTEGER NAPAR(NPOS),NUMP(NPOS)
  34.         REAL*8 CREF(NDIM,NPOS),TPAR(NPOS)
  35.       ENDSEGMENT
  36.       POINTEUR IZLAR.IPARPO
  37. C
  38.       SEGMENT IZSH
  39.             REAL*8 SHP(6,MNO9),SHY(12,MNO9),XYZL(3,MNO9)
  40.       ENDSEGMENT
  41. C
  42. C
  43. C
  44. C
  45. C
  46. C
  47.       DIMENSION XR(3)
  48. C
  49. C
  50. C
  51. C     WRITE(6,*) ' TRJRAL '
  52.       NDIM=CREF(/1)
  53.       NPOS=NAPAR(/1)
  54.       SEGINI IZLAR
  55. C     write(6,*)' SEGINI IZLAR',IZLAR
  56.       DO 1 IPOS=1,NPOS
  57.       IEL=NAPAR(IPOS)
  58. C     write(6,*) iel,nump(ipos),cref(1,ipos),cref(2,ipos),tpar(ipos)
  59.       IZLAR.NAPAR(IPOS)=IEL
  60.       IZLAR.TPAR(IPOS)=TPAR(IPOS)
  61.       IZLAR.NUMP(IPOS)=NUMP(IPOS)
  62.       CALL MELNEL(IEL,MELEME,IPT1,NEL0,1)
  63.       NOEL=IPT1.NUM(/1)
  64.       IELL=IEL-NEL0
  65.       CALL DOXE(XCOOR,IDIM,NOEL,IPT1.NUM,IELL,XYZL)
  66.       ITY=IPT1.ITYPEL
  67.       XR(3)=0.D0
  68.       DO 4 ID1=1,NDIM
  69.       XR(ID1)=CREF(ID1,IPOS)
  70. C     WRITE(6,*)' CREFIJ ',ID1,IPOS,CREF(ID1,IPOS)
  71.    4  CONTINUE
  72.       CALL SHAPE(XR(1),XR(2),XR(3),ITY,SHP,IRET)
  73.       DO 5 ID1=1,NDIM
  74.       IZLAR.CREF(ID1,IPOS)=0.D0
  75.       DO 6 K=1,NOEL
  76.       IZLAR.CREF(ID1,IPOS)=IZLAR.CREF(ID1,IPOS)+SHP(1,K)*XYZL(ID1,K)
  77.    6  CONTINUE
  78.    5  CONTINUE
  79. C     WRITE(6,*)IPOS,' CLEFIJ ',(IZLAR.CREF(ID1,IPOS),ID1=1,NDIM)
  80.    1  CONTINUE
  81.       RETURN
  82.       END
  83.  
  84.  
  85.  

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