Télécharger la plaquette Cast3M


fr en

Télécharger ellp23.eso

Retour à la liste

Numérotation des lignes : 

ellp23

  1. C ELLP23    SOURCE    KK2000    14/04/09    21:15:22     8027
  2.         SUBROUTINE ELLP23(CARACT,COOR,GAMA,ZXX,XCOR,VALDE1,VALDE2,
  3.      *                     ZS,NP,NBELEM,XPI)
  4. C
  5.       IMPLICIT INTEGER(I-N)
  6.         IMPLICIT REAL*8     (A-H,O-Y)
  7.         IMPLICIT COMPLEX*16 (Z)
  8. C
  9. C CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
  10. C
  11. C   OPERATEUR ELFE LAPLACE POUTRE
  12. C
  13. C   CALCUL POUR LA POUTRE N› INP LES VALEURS DE Ua ET DE Ub EN LOCAL
  14. C
  15. C   PARAMETRES :
  16. C
  17. C   CARACT : TABLEAU DES CARACTERISTIQUE DES POUTRES (9,NP)
  18. C   COOR   : TABLEAU DES COORDONNEES DES NOEUDS (3,2*NP)
  19. C   GAMA   : TABLEAU DES VECTEUR DEFINISSANT LE PLAN LOCALE OXY (3,2*NP)
  20. C   ZXX    : VECTEUR (24*NP) DONNANT LE VECTEUR SOLUTION DU TREILLIS
  21. C            DANS LE REPERE GLOBAL
  22. C   XCOR   : TABLEAU DONNANT LES COORD. DES NOEUDS DE L'OBJET MAILLAGE
  23. C            DANS LE REPERE GLOBAL
  24. C   ZS     : VALEUR DE S = S0 + I W
  25. C   NP     : NOMBRE DE POUTRES DU TREILLIS
  26. C   NBELEM : NOMBRE D'ELEMENTS DU SOUS-MAILLAGE
  27. C   XPI    : VALEUR PRECISE DE PI
  28. C
  29. C
  30. C   SORTIES :
  31. C
  32. C   VALDE1 : MODULE DE UX, UY, UZ POUR LES 2 POINTS DE CHAQUE ELEMENT
  33. C            DU SOUS-MAILLAGE
  34. C   VALDE2 : PHASE DE UX, UY, UZ POUR LES 2 POINTS DE CHAQUE ELEMENT
  35. C            DU SOUS-MAILLAGE
  36. C
  37. C
  38. C   AUTEUR : SAINT-DIZIER
  39. C   DATE   : 19 JUILLET 1990
  40. C
  41. C CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
  42. C
  43.         DIMENSION CARACT(12,*),COOR(3,*),GAMA(3,*)
  44.         DIMENSION ZXX(*),XCOR(2,3,NBELEM)
  45.         DIMENSION VALDE1(2,NBELEM,3),VALDE2(2,NBELEM,3)
  46. C
  47.         DIMENSION ZU(24),ZU1(24),ZU2(24)
  48.         DIMENSION R1(3,3),R2(3,3),ZD(12),ZE(12,12)
  49. C
  50.         DO 100 INP = 1 , NP
  51. C
  52.           DO 101 II = 1 , 24
  53.             ZU(II) = ZXX(24*(INP-1)+II)
  54. 101       CONTINUE
  55. C
  56. C
  57. C -------------- ACQUISITION DES CARACTERISTIQUES DE LA POUTRE
  58. C                ---------------------------------------------
  59.         CE    = CARACT( 1,INP)
  60.         CNU   = CARACT( 2,INP)
  61.         CRO   = CARACT( 3,INP)
  62.         CSE   = CARACT( 4,INP)
  63.         CC    = CARACT( 5,INP)
  64.         CIP   = CARACT( 6,INP)
  65.         CIY   = CARACT( 7,INP)
  66.         CIZ   = CARACT( 8,INP)
  67.         CKCY  = CARACT( 9,INP)
  68.         CKCZ  = CARACT(10,INP)
  69.         CAM   = CARACT(11,INP)
  70.         CETA  = CARACT(12,INP)
  71. C
  72.         ZI = (0.D0 , 1.D0)
  73. C
  74.         ZCE = CE*(CMPLX(1.D0)+CETA*ZI)
  75.         ZCP   = SQRT(ZCE/CRO)
  76.         ZCMU  = ZCE/(CMPLX(2.D0)*(CMPLX(1.D0)+CNU))
  77.         ZCR   = SQRT(ZCMU/CRO)
  78. C
  79.         ZALFA2= ZS*ZS*CRO*CSE + CAM*ZS
  80.         ZALFA2= ZALFA2 / (ZCE*CSE)
  81.         ZALFAX= SQRT(ZALFA2)
  82. C
  83.         ZALFA2= ZS*ZS*CRO*CIP + CAM*ZS
  84.         ZALFA2= ZALFA2 / (ZCMU*CC)
  85.         ZALFAR= SQRT(ZALFA2)
  86. C
  87.         ZALFA4= CRO*CRO*CIZ*ZS*ZS*ZS*ZS/ZCMU/CKCY
  88.      *        + CAM*CRO*CIZ*ZS*ZS*ZS/ZCMU/CKCY/CSE
  89.      *        + CRO*CSE*ZS*ZS
  90.      *        + CAM*ZS
  91.         ZALFA4= ZALFA4/(ZCE*CIZ)
  92.         ZALFA2= SQRT(ZALFA4)
  93.         ZALFAY= SQRT(ZALFA2)
  94. C
  95.         ZZ3   = CRO*CIZ*(CMPLX(1.D0)+ZCE/(ZCMU*CKCY))*ZS*ZS
  96.      *        + CAM*ZCE*CIZ*ZS/(ZCMU*CKCY*CSE)
  97.         ZZ3   = ZZ3 / (ZCE*CIZ*ZALFA2)
  98. C
  99.         ZALFA4= CRO*CRO*CIY*ZS*ZS*ZS*ZS/ZCMU/CKCZ
  100.      *        + CAM*CRO*CIY*ZS*ZS*ZS/ZCMU/CKCZ/CSE
  101.      *        + CRO*CSE*ZS*ZS
  102.      *        + CAM*ZS
  103.         ZALFA4= ZALFA4/(ZCE*CIY)
  104.         ZALFA2= SQRT(ZALFA4)
  105.         ZALFAZ= SQRT(ZALFA2)
  106. C
  107.         ZZ4   = CRO*CIY*(CMPLX(1.D0)+ZCE/(ZCMU*CKCZ))*ZS*ZS
  108.      *        + CAM*ZCE*CIY*ZS/(ZCMU*CKCZ*CSE)
  109.         ZZ4   = ZZ4 / (ZCE*CIY*ZALFA2)
  110. C
  111.         N1  = 2*INP-1
  112.         N2  = 2*INP
  113.         XXA  = COOR (1,N1)
  114.         XYA  = COOR (2,N1)
  115.         XZA  = COOR (3,N1)
  116.         XXB  = COOR (1,N2)
  117.         XYB  = COOR (2,N2)
  118.         XZB  = COOR (3,N2)
  119. C
  120.         XL = SQRT((XXB-XXA)**2 + (XYB-XYA)**2 + (XZB-XZA)**2)
  121. C
  122. C ------------------------------ VECTEUR UNITAIRE OX REPERE LOCALE
  123. C                                ---------------------------------
  124.         XI1 = (XXB-XXA)/XL
  125.         XI2 = (XYB-XYA)/XL
  126.         XI3 = (XZB-XZA)/XL
  127. C
  128. C ------------------------- VECTEUR UNITAIRE OY REPERE LOCALE
  129. C                           ---------------------------------
  130.         GX = GAMA(1,INP)
  131.         GY = GAMA(2,INP)
  132.         GZ = GAMA(3,INP)
  133.         GG = SQRT(GX*GX + GY*GY + GZ*GZ)
  134.         GX = GX/GG
  135.         GY = GY/GG
  136.         GZ = GZ/GG
  137. C
  138.         DELTA =  SQRT (1.D0 - (XI1*GX + XI2*GY + XI3*GZ)**2)
  139. C
  140.         DET=-(GX*XI3-GZ*XI1)**2-(GY*XI1-GX*XI2)**2-(GY*XI3-GZ*XI2)**2
  141. C
  142.         XJ1 = (XI2*(GY*XI1-GX*XI2)-XI3*(GX*XI3-GZ*XI1))*DELTA/DET
  143.         XJ2 = (XI3*(GZ*XI2-GY*XI3)-XI1*(GY*XI1-GX*XI2))*DELTA/DET
  144.         XJ3 = (XI1*(GX*XI3-GZ*XI1)-XI2*(GZ*XI2-GY*XI3))*DELTA/DET
  145. C
  146. C ---------------------------- VECTEUR UNITAIRE OZ REPERE LOCALE
  147. C                              ---------------------------------
  148.         XK1 = XI2*XJ3 - XI3*XJ2
  149.         XK2 = XI3*XJ1 - XI1*XJ3
  150.         XK3 = XI1*XJ2 - XI2*XJ1
  151. C
  152.         R1(1,1) = XI1
  153.         R1(1,2) = XJ1
  154.         R1(1,3) = XK1
  155.         R1(2,1) = XI2
  156.         R1(2,2) = XJ2
  157.         R1(2,3) = XK2
  158.         R1(3,1) = XI3
  159.         R1(3,2) = XJ3
  160.         R1(3,3) = XK3
  161. C
  162. C
  163.         R2(1,1) = XJ2*XK3 - XJ3*XK2
  164.         R2(1,2) = XJ3*XK1 - XJ1*XK3
  165.         R2(1,3) = XJ1*XK2 - XJ2*XK1
  166.         R2(2,1) = XI3*XK2 - XI2*XK3
  167.         R2(2,2) = XI1*XK3 - XI3*XK1
  168.         R2(2,3) = XI2*XK1 - XI1*XK2
  169.         R2(3,1) = XI2*XJ3 - XI3*XJ2
  170.         R2(3,2) = XI3*XJ1 - XI1*XJ3
  171.         R2(3,3) = XI1*XJ2 - XI2*XJ1
  172. C
  173. C -----------------------------  MATRICE DIMENSIONNEMENT ZD
  174. C                                --------------------------
  175.         ZD( 1) = XL**2/ZCP
  176.         ZD( 2) = XL**2/ZCP
  177.         ZD( 3) = XL**2/ZCP
  178.         ZD( 4) = XL/ZCR
  179.         ZD( 5) = XL**2/ZCP*ZALFAZ
  180.         ZD( 6) = XL**2/ZCP*ZALFAY
  181.         ZD( 7) = XL**2/ZCP*ZALFAX
  182.         ZD( 8) = XL**2/ZCP*(ZALFAY**3)
  183.         ZD( 9) = XL**2/ZCP*(ZALFAZ**3)
  184.         ZD(10) = XL/ZCR*ZALFAR
  185.         ZD(11) = XL**2/ZCP*(ZALFAZ**2)
  186.         ZD(12) = XL**2/ZCP*(ZALFAY**2)
  187. C
  188. C
  189. C ------------------------------- MATRICE DES EFFORTS ZE
  190. C                                 ----------------------
  191.         DO 20 I = 1 , 12
  192.             DO 21 J = 1 , 12
  193.                 ZE(I,J) = 0.D0
  194. 21          CONTINUE
  195. 20      CONTINUE
  196. C
  197. C
  198.         ZA5 = CMPLX(1.D0)
  199.         ZA6 = CMPLX(1.D0)
  200.         ZB6 = CMPLX(1.D0) / (ZCMU*CKCY*CSE)
  201.         ZB5 = CMPLX(1.D0) / (ZCMU*CKCZ*CSE)
  202. C
  203.         ZA8 = (CRO*CSE*ZS*ZS+CAM*ZS)/(ZCMU*ZCMU*CKCY*CKCY*CSE*CSE)
  204.      *      -  CMPLX(1.D0)/(ZCE*CIZ)
  205.         ZA9 = (CRO*CSE*ZS*ZS+CAM*ZS)/(ZCMU*ZCMU*CKCZ*CKCZ*CSE*CSE)
  206.      *      -  CMPLX(1.D0)/(ZCE*CIY)
  207.         ZB8 = CRO*ZS*ZS*(CMPLX(1.D0)/(ZCMU*CKCY)+CMPLX(1.D0)/ZCE)
  208.      *      + CAM*ZS/(ZCMU*CKCY*CSE)
  209.         ZB9 = CRO*ZS*ZS*(CMPLX(1.D0)/(ZCMU*CKCZ)+CMPLX(1.D0)/ZCE)
  210.      *      + CAM*ZS/(ZCMU*CKCZ*CSE)
  211. C
  212.         ZA12= CMPLX(1.D0)/(ZCE*CIZ)
  213.         ZA11= CMPLX(1.D0)/(ZCE*CIY)
  214.         ZB12= (CRO*CSE*ZS*ZS+CAM*ZS)/(ZCMU*CKCY*CSE)
  215.         ZB11= (CRO*CSE*ZS*ZS+CAM*ZS)/(ZCMU*CKCZ*CSE)
  216. C
  217. C
  218.         ZE( 1, 1) = CMPLX(1.D0)
  219.         ZE( 2, 2) = CMPLX(1.D0)
  220.         ZE( 3, 3) = CMPLX(1.D0)
  221.         ZE( 4, 4) = CMPLX(1.D0)
  222.         ZE( 5, 5) = ZA5
  223.         ZE( 5, 9) = ZB5
  224.         ZE( 6, 6) = ZA6
  225.         ZE( 6, 8) = ZB6
  226.         ZE( 7, 7) = CMPLX(1.D0)/(ZCE*CSE)
  227.         ZE( 8, 6) = ZB8
  228.         ZE( 8, 8) = ZA8
  229.         ZE( 9, 5) = ZB9
  230.         ZE( 9, 9) = ZA9
  231.         ZE(10,10) =-CMPLX(1.D0)/(ZCMU*CC)
  232.         ZE(11, 3) = ZB11
  233.         ZE(11,11) = ZA11
  234.         ZE(12, 2) = ZB12
  235.         ZE(12,12) = ZA12
  236. C
  237. C
  238. C -- CACUL DE ZU1 = ZR * ZU
  239. C
  240.         DO 27 I = 1 , 24
  241.           ZU1(I) = CMPLX(0.D0)
  242.           ZU2(I) = CMPLX(0.D0)
  243. 27      CONTINUE
  244. C
  245.         DO 26 I = 1 , 8
  246.             DO 25 J = 1 , 3
  247.                 DO 24 K = 1 , 3
  248.                   ZU1((I-1)*3+J)=ZU1((I-1)*3+J)+R2(J,K)*ZU((I-1)*3+K)
  249. 24              CONTINUE
  250. 25          CONTINUE
  251. 26      CONTINUE
  252. C
  253. C -- CACUL DE ZU2 = ZE * ZU1
  254. C
  255.         DO 30 I = 1,12
  256.             DO 31 J = 1,12
  257.               ZU2(I) = ZU2(I) + ZE(I,J)*ZU1(J)
  258. 31          CONTINUE
  259. 30      CONTINUE
  260. C
  261.         DO 32 I = 1,12
  262.             DO 33 J = 1,6
  263.               ZU2(I+12) = ZU2(I+12) + ZE(I,J  )*ZU1(J+12)
  264.               ZU2(I+12) = ZU2(I+12) - ZE(I,J+6)*ZU1(J+18)
  265. 33          CONTINUE
  266. 32      CONTINUE
  267. C
  268. C
  269. C -- CACUL DE ZU2 = ZD * ZU2
  270. C
  271. C
  272.         DO 34 I = 1,12
  273.             ZU2(I   ) = ZU2(I   )/ZD(I)
  274.             ZU2(I+12) = ZU2(I+12)/ZD(I)
  275. 34      CONTINUE
  276. C
  277. C
  278. C ---------------  NBELEM : NOMBRE D'ELEMENTS DANS LE SOUS MAILLAGE
  279. C
  280.         DO 110 IN = 1 , NBELEM
  281. C
  282.             XX1 = XCOR ( 1,1,IN)
  283.             XY1 = XCOR ( 1,2,IN)
  284.             XZ1 = XCOR ( 1,3,IN)
  285. C
  286.             XX2 = XCOR ( 2,1,IN)
  287.             XY2 = XCOR ( 2,2,IN)
  288.             XZ2 = XCOR ( 2,3,IN)
  289. C
  290.             CALL ELLP07(XXA,XYA,XZA,XXB,XYB,XZB,
  291.      *                  XX1,XY1,XZ1,XX2,XY2,XZ2,IOUI)
  292. C
  293.             IF (IOUI.EQ.0) THEN
  294. C
  295.              XXD1 = SQRT((XXA-XX1)**2+(XYA-XY1)**2+(XZA-XZ1)**2)
  296.              XXD2 = SQRT((XXA-XX2)**2+(XYA-XY2)**2+(XZA-XZ2)**2)
  297.              XXL1= XL -XXD1
  298.              XXL2= XL -XXD2
  299. C
  300. C              XXDi  :  DISTANCE DE A AU PT M DE CALCUL DU DEPLACEMENT
  301. C
  302.              CALL ELLP31(XXD1,INP,ZS,ZALFAX,ZALFAR,ZALFAY,ZALFAZ,
  303.      *       Z10X ,Z11X ,Z1,Z2,Z30X ,Z31X ,Z32X ,Z33X ,Z34X ,ZZ3,
  304.      *       Z40X ,Z41X ,Z42X ,Z43X ,Z44X ,ZZ4)
  305. C
  306.              CALL ELLP31(XXL1,INP,ZS,ZALFAX,ZALFAR,ZALFAY,ZALFAZ,
  307.      *       Z10XL,Z11XL,Z1,Z2,Z30XL,Z31XL,Z32XL,Z33XL,Z34XL,ZZ3,
  308.      *       Z40XL,Z41XL,Z42XL,Z43XL,Z44XL,ZZ4)
  309. C
  310.              ZUX1 = ZU2(19)*Z10XL-ZU2(13)*Z11XL-ZU2(7)*Z10X-ZU2(1)*Z11X
  311. C
  312.              ZUY1 =-(ZZ3*Z31X -Z33X )*ZU2(2) - (ZZ3*Z30X -Z32X )*ZU2(6)
  313.      *             + Z31X *ZU2(12) + Z30X *ZU2(8)
  314.      *             -(ZZ3*Z31XL-Z33XL)*ZU2(14)+(ZZ3*Z30XL-Z32XL)*ZU2(18)
  315.      *             + Z31XL*ZU2(24) - Z30XL*ZU2(20)
  316. C
  317.              ZUZ1 =-(ZZ4*Z41X -Z43X )*ZU2(3) - (ZZ4*Z40X -Z42X )*ZU2(5)
  318.      *             + Z41X *ZU2(11) + Z40X *ZU2(9)
  319.      *             -(ZZ4*Z41XL-Z43XL)*ZU2(15)+(ZZ4*Z40XL-Z42XL)*ZU2(17)
  320.      *             + Z41XL*ZU2(23) - Z40XL*ZU2(21)
  321. C
  322. C
  323.              ZUX1 = ZUX1 * ZD(1)
  324.              ZUY1 = ZUY1 * ZD(2)
  325.              ZUZ1 = ZUZ1 * ZD(3)
  326. C
  327.              ZX1 = R1(1,1)*ZUX1 + R1(1,2)*ZUY1 + R1(1,3)*ZUZ1
  328.              ZY1 = R1(2,1)*ZUX1 + R1(2,2)*ZUY1 + R1(2,3)*ZUZ1
  329.              ZZ1 = R1(3,1)*ZUX1 + R1(3,2)*ZUY1 + R1(3,3)*ZUZ1
  330. C
  331.              CALL ELLP31(XXD2,INP,ZS,ZALFAX,ZALFAR,ZALFAY,ZALFAZ,
  332.      *       Z10X ,Z11X ,Z1,Z2,Z30X ,Z31X ,Z32X ,Z33X ,Z34X ,ZZ3,
  333.      *       Z40X ,Z41X ,Z42X ,Z43X ,Z44X ,ZZ4)
  334. C
  335.              CALL ELLP31(XXL2,INP,ZS,ZALFAX,ZALFAR,ZALFAY,ZALFAZ,
  336.      *       Z10XL,Z11XL,Z1,Z2,Z30XL,Z31XL,Z32XL,Z33XL,Z34XL,ZZ3,
  337.      *       Z40XL,Z41XL,Z42XL,Z43XL,Z44XL,ZZ4)
  338. C
  339. C
  340.             ZUX2 = ZU2(19)*Z10XL-ZU2(13)*Z11XL-ZU2(7)*Z10X-ZU2(1)*Z11X
  341. C
  342.             ZUY2 =-(ZZ3*Z31X -Z33X )*ZU2(2) - (ZZ3*Z30X -Z32X )*ZU2(6)
  343.      *            + Z31X *ZU2(12) + Z30X *ZU2(8)
  344.      *            -(ZZ3*Z31XL-Z33XL)*ZU2(14)+(ZZ3*Z30XL-Z32XL)*ZU2(18)
  345.      *            + Z31XL*ZU2(24) - Z30XL*ZU2(20)
  346. C
  347.             ZUZ2 =-(ZZ4*Z41X -Z43X )*ZU2(3) - (ZZ4*Z40X -Z42X )*ZU2(5)
  348.      *            + Z41X *ZU2(11) + Z40X *ZU2(9)
  349.      *            -(ZZ4*Z41XL-Z43XL)*ZU2(15)+(ZZ4*Z40XL-Z42XL)*ZU2(17)
  350.      *            + Z41XL*ZU2(23) - Z40XL*ZU2(21)
  351. C
  352. C
  353.             ZUX2 = ZUX2 * ZD(1)
  354.             ZUY2 = ZUY2 * ZD(2)
  355.             ZUZ2 = ZUZ2 * ZD(3)
  356. C
  357. C
  358.             ZX2 = R1(1,1)*ZUX2 + R1(1,2)*ZUY2 + R1(1,3)*ZUZ2
  359.             ZY2 = R1(2,1)*ZUX2 + R1(2,2)*ZUY2 + R1(2,3)*ZUZ2
  360.             ZZ2 = R1(3,1)*ZUX2 + R1(3,2)*ZUY2 + R1(3,3)*ZUZ2
  361. C
  362. C
  363.             PRZX1 = ZX1
  364.             PIZX1 = ZX1*(0.D0,-1.D0)
  365.             PRZY1 = ZY1
  366.             PIZY1 = ZY1*(0.D0,-1.D0)
  367.             PRZZ1 = ZZ1
  368.             PIZZ1 = ZZ1*(0.D0,-1.D0)
  369.             PRZX2 = ZX2
  370.             PIZX2 = ZX2*(0.D0,-1.D0)
  371.             PRZY2 = ZY2
  372.             PIZY2 = ZY2*(0.D0,-1.D0)
  373.             PRZZ2 = ZZ2
  374.             PIZZ2 = ZZ2*(0.D0,-1.D0)
  375. C
  376.             IF (ABS(PRZX1).LT.1.E-38.AND.ABS(PIZX1).LT.1.E-38) THEN
  377.                PPX1 = 0.D0
  378.             ELSE
  379.                PPX1 = ATAN2(PIZX1,PRZX1)*180.D0/XPI
  380.             END IF
  381. C
  382.             IF (ABS(PRZY1).LT.1.E-38.AND.ABS(PIZY1).LT.1.E-38) THEN
  383.                PPY1 = 0.D0
  384.             ELSE
  385.                PPY1 = ATAN2(PIZY1,PRZY1)*180.D0/XPI
  386.             END IF
  387. C
  388.             IF (ABS(PRZZ1).LT.1.E-38.AND.ABS(PIZZ1).LT.1.E-38) THEN
  389.                PPZ1 = 0.D0
  390.             ELSE
  391.                PPZ1 = ATAN2(PIZZ1,PRZZ1)*180.D0/XPI
  392.             END IF
  393. C
  394. C
  395.             IF (ABS(PRZX2).LT.1.E-38.AND.ABS(PIZX2).LT.1.E-38) THEN
  396.                PPX2 = 0.D0
  397.             ELSE
  398.                PPX2 = ATAN2(PIZX2,PRZX2)*180.D0/XPI
  399.             END IF
  400. C
  401.             IF (ABS(PRZY2).LT.1.E-38.AND.ABS(PIZY2).LT.1.E-38) THEN
  402.                PPY2 = 0.D0
  403.             ELSE
  404.                PPY2 = ATAN2(PIZY2,PRZY2)*180.D0/XPI
  405.             END IF
  406. C
  407.             IF (ABS(PRZZ2).LT.1.E-38.AND.ABS(PIZZ2).LT.1.E-38) THEN
  408.                PPZ2 = 0.D0
  409.             ELSE
  410.                PPZ2 = ATAN2(PIZZ2,PRZZ2)*180.D0/XPI
  411.             END IF
  412. C
  413.             VALDE1(1,IN,1) = ABS(ZX1)
  414.             VALDE1(1,IN,2) = ABS(ZY1)
  415.             VALDE1(1,IN,3) = ABS(ZZ1)
  416.             VALDE1(2,IN,1) = ABS(ZX2)
  417.             VALDE1(2,IN,2) = ABS(ZY2)
  418.             VALDE1(2,IN,3) = ABS(ZZ2)
  419. C
  420.             VALDE2(1,IN,1) = PPX1
  421.             VALDE2(1,IN,2) = PPY1
  422.             VALDE2(1,IN,3) = PPZ1
  423.             VALDE2(2,IN,1) = PPX2
  424.             VALDE2(2,IN,2) = PPY2
  425.             VALDE2(2,IN,3) = PPZ2
  426. C
  427. C
  428.         END IF
  429. C
  430. 110     CONTINUE
  431. C
  432. 100     CONTINUE
  433. C
  434.         END
  435.  
  436.  
  437.  
  438.  

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