Télécharger ffpr15.eso

Retour à la liste

Numérotation des lignes :

  1. C FFPR15 SOURCE CHAT 05/01/12 23:59:06 5004
  2. SUBROUTINE FFPR15(XA)
  3. IMPLICIT INTEGER(I-N)
  4. IMPLICIT REAL*8 (A-H,O-Z)
  5. C************************************************************************
  6. C
  7. C
  8. C
  9. C************************************************************************
  10. DIMENSION FN(15),COOR(3,3),XA(3,*)
  11. C
  12. A=SQRT(2.D0)
  13. AS2=A*0.5D0
  14. C? CA2=2.D0/A/A
  15. CA2=1.D0
  16. C? CA4=4.D0/A/A
  17. CA4=2.D0
  18. A4=4.D0/A
  19.  
  20. COOR(1,1)=AS2
  21. COOR(2,1)=0.D0
  22. COOR(3,1)=0.5D0
  23.  
  24. COOR(1,2)=AS2
  25. COOR(2,2)=AS2
  26. COOR(3,2)=0.5D0
  27.  
  28. COOR(1,3)=0.D0
  29. COOR(2,3)=AS2
  30. COOR(3,3)=0.5D0
  31.  
  32. DO 1 I=1,3
  33. X=COOR(1,I)
  34. Y=COOR(2,I)
  35. Z=COOR(3,I)
  36.  
  37. FN(1)= CA2*(X+Y-A)*(1.D0-Z)*(X+Y+A*Z-AS2)
  38. FN(2)=-CA4*X*(X+Y-A)*(1.D0-Z)
  39. FN(3)=-CA2*X*(1.D0-Z)*(-X+A*Z+AS2)
  40. FN(4)= CA4*X*Y*(1.D0-Z)
  41. FN(5)=-CA2*Y*(1.D0-Z)*(-Y+A*Z+AS2)
  42. FN(6)=-CA4*Y*(X+Y-A)*(1.D0-Z)
  43. FN(7)=-A4*Z*(1.D0-Z)*(X+Y-A)
  44. FN(8)= A4*X*Z*(1.D0-Z)
  45. FN(9) = A4*Y*Z*(1.D0-Z)
  46. FN(10)= CA2*(X+Y-A)*Z*(X+Y+A*(1.D0-Z)-AS2)
  47. FN(11)=-CA4*X*(X+Y-A)*Z
  48. FN(12)=-CA2*X*Z*(-X+A*(1.D0-Z)+AS2)
  49. FN(13)= CA4*X*Y*Z
  50. FN(14)=-CA2*Y*Z*(-Y+A*(1.D0-Z)+AS2)
  51. FN(15)=-CA4*Y*Z*(X+Y-A)
  52.  
  53. DO 2 N=1,3
  54. XA(N,I+15)=
  55. & FN(1)*XA(N,1)+FN(2)*XA(N,2)+FN(3)*XA(N,3)+FN(4)*XA(N,4)
  56. &+ FN(5)*XA(N,5)+FN(6)*XA(N,6)+FN(7)*XA(N,7)+FN(8)*XA(N,8)
  57. &+ FN(9)*XA(N,9)+FN(10)*XA(N,10)+FN(11)*XA(N,11)+FN(12)*XA(N,12)
  58. &+FN(13)*XA(N,13)+FN(14)*XA(N,14)+FN(15)*XA(N,15)
  59. 2 CONTINUE
  60. 1 CONTINUE
  61.  
  62.  
  63.  
  64.  
  65. 1002 format(10(1x,1pe11.4))
  66. RETURN
  67. END
  68.  
  69.  
  70.  

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