Télécharger deinqt.eso

Retour à la liste

Numérotation des lignes :

  1. C DEINQT SOURCE CB215821 16/04/21 21:16:22 8920
  2.  
  3. SUBROUTINE DERINVARIPQT (SIG,P_S,Q_S,T_S,N)
  4. c la derivada t_s falta dividirla por sin(3t)
  5. c no se hace pq luego eso cancela con las derivadas de g(t), MODELO MRSLade
  6. IMPLICIT INTEGER(I-N)
  7. integer n,i
  8. real*8 sig(n),p_s(n),q_s(n),t_s(n),j2,j3,aux0,aux1,aux2
  9. call DerInvariIJJ(sig,p_s,q_s,t_s,n)
  10. do i=1,3
  11. p_s(i)=-p_s(i)/3.D0
  12. enddo
  13. aux0=0.D0
  14. aux1=0.D0
  15. aux2=0.D0
  16. call Invari_j2(sig,n,j2)
  17. if (j2.ge.1.D-10) then
  18. call Invari_j3(sig,n,j3)
  19. aux0=SQRT(3.D0/j2)/2.D0
  20. aux1=aux0/j2
  21. aux2=1.5D0*j3/j2
  22. endif
  23. do i=1,n
  24. t_s(i)=aux1*(aux2*q_s(i)-t_s(i))
  25. q_s(i)=aux0*q_s(i)
  26. enddo
  27. end
  28.  
  29.  
  30.  
  31.  
  32.  
  33.  

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