Télécharger druco1.eso

Retour à la liste

Numérotation des lignes :

  1. C DRUCO1 SOURCE CHAT 05/01/12 22:55:14 5004
  2. C DRUCO1 SOURCE JFG 02/07/04
  3. SUBROUTINE DRUCO1(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,IGAU1
  13. ENDSEGMENT
  14. C
  15. BETA = 1.2D0
  16. ALC=(BETA-1.)/(2.*BETA-1.)
  17. BEC=BETA/(2.*BETA-1.)
  18. C
  19. SX=S1(1)
  20. SY=S1(2)
  21. SZ=S1(3)
  22. SXY=S1(4)
  23. SXZ=S1(5)
  24. SYZ=S1(6)
  25. A = SX-SY
  26. B = SY-SZ
  27. C = SX-SZ
  28. DJ2=(A*A+B*B+C*C+SXY*SXY+SYZ*SYZ+SXZ*SXZ)
  29. DJ2=DJ2 * (1.D0/6.D0)
  30. DJ2=3.D0*DJ2
  31. DJ2=SQRT(DJ2)
  32. C
  33. DF(1) = ((2.*SX) - SY -SZ)/2./DJ2 + ALC
  34. DF(2) = (-SX +(2.*SY) -SZ)/2./DJ2 + ALC
  35. DF(3) = (-SX -SY +(2.*SZ))/2./DJ2 + ALC
  36. DF(4) = 6. * SXY/2./DJ2
  37. DF(5) = 6. * SXZ/2./DJ2
  38. DF(6) = 6. * SYZ/2./DJ2
  39. C
  40. RETURN
  41. END
  42.  
  43.  
  44.  
  45.  
  46.  

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