Télécharger drutr1.eso

Retour à la liste

Numérotation des lignes :

drutr1
  1. C DRUTR1 SOURCE CHAT 05/01/12 22:55:22 5004
  2. C DRUTR1 SOURCE JFG 02/07/04
  3. SUBROUTINE DRUTR1(S1,DF,BETINSA)
  4. C
  5. IMPLICIT INTEGER(I-N)
  6. IMPLICIT REAL*8(A-H,O-Z)
  7. DIMENSION S1(6),DF(6)
  8. C
  9. SEGMENT BETINSA
  10. REAL*8 RT,RC,YOUN,XNU,GFT,GFC,CAR
  11. REAL*8 DKT,DKC,SEQT,SEQC,ENDT,ENDC
  12. INTEGER IFIS,IPLA,IBB,IGAU
  13. ENDSEGMENT
  14. C
  15. ALT=1.D0
  16. BET=2.D0
  17. C
  18. SX=S1(1)
  19. SY=S1(2)
  20. SZ=S1(3)
  21. SXY=S1(4)
  22. SXZ=S1(5)
  23. SYZ=S1(6)
  24. A = SX-SY
  25. B = SY-SZ
  26. C = SX-SZ
  27. DJ2=(A*A+B*B+C*C+SXY*SXY+SYZ*SYZ+SXZ*SXZ)
  28. DJ2=DJ2 * (1.D0/6.D0)
  29. DJ2=3.D0*DJ2
  30. DJ2=SQRT(DJ2)
  31. C
  32. DF(1) = ((2.*SX) - SY -SZ)/2./DJ2 + ALT
  33. DF(2) = (-SX +(2.*SY) -SZ)/2./DJ2 + ALT
  34. DF(3) = (-SX -SY +(2.*SZ))/2./DJ2 + ALT
  35. DF(4) = 6. * SXY/2./DJ2
  36. DF(5) = 6. * SXZ/2./DJ2
  37. DF(6) = 6. * SYZ/2./DJ2
  38. C
  39. RETURN
  40. END
  41.  
  42.  
  43.  
  44.  
  45.  

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