Télécharger deinij.eso

Retour à la liste

Numérotation des lignes :

  1. C DEINIJ SOURCE CHAT 05/01/12 22:41:11 5004
  2.  
  3. SUBROUTINE DERINVARIIJJ (SIG,I_S,J2_S,J3_S,N)
  4. IMPLICIT INTEGER(I-N)
  5. integer n,i
  6. real*8 sig(n),I_s(n),J2_s(n),J3_s(n),aux1,J2
  7. call Desviador(sig,J2_s,n)
  8. call Invari_J2(sig,n,J2)
  9. J3_s(1)=J2_s(1)*J2_s(1)
  10. J3_s(2)=J2_s(2)*J2_s(2)
  11. J3_s(3)=J2_s(3)*J2_s(3)
  12. if (n.ge.4) then
  13. J3_s(1)=J3_s(1)+J2_s(4)*J2_s(4)
  14. J3_s(2)=J3_s(2)+J2_s(4)*J2_s(4)
  15. J3_s(4)=2.D0*(J2_s(1)+J2_s(2))*J2_s(4)
  16. if (n.eq.6) then
  17. J3_s(1)=J3_s(1)+J2_s(5)*J2_s(5)
  18. J3_s(2)=J3_s(2)+J2_s(6)*J2_s(6)
  19. J3_s(3)=J3_s(3)+J2_s(5)*J2_s(5)+J2_s(6)*J2_s(6)
  20. J3_s(4)=J3_s(4)+J2_s(5)*J2_s(6)
  21. J3_s(5)=(J2_s(1)+J2_s(3))*J2_s(5)+J2_s(4)*J2_s(6)
  22. J3_s(6)=(J2_s(2)+J2_s(3))*J2_s(6)+J2_s(4)*J2_s(5)
  23. endif
  24. endif
  25. aux1=-J2*2.D0/3.D0
  26. do i=1,3
  27. J3_s(i)=J3_s(i)+aux1
  28. enddo
  29. do i=1,3
  30. I_s(i)=1.D0
  31. enddo
  32. if (n.lt.4) return
  33. do i=4,n
  34. I_s(i)=0.D0
  35. J2_s(i)=2.D0*J2_s(i)
  36. enddo
  37. return
  38. end
  39.  
  40.  
  41.  
  42.  

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