Esope source files

   1 : C DEINQT    SOURCE    CHAT      05/01/12    22:41:15     5004           
   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
   8 :       call DerInvariIJJ(sig,p_s,q_s,t_s,n)
   9 :       do i=1,3
  10 :          p_s(i)=-p_s(i)/3.D0
  11 :       enddo
  12 :       aux0=0.D0
  13 :       aux1=0.D0
  14 :       aux2=0.D0
  14 :       call Invari_j2(sig,n,j2)
  15 :       if (j2.ge.1.D-10) then
  15 :         call Invari_j3(sig,n,j3)
  16 :         aux0=dsqrt(3.D0/j2)/2.D0
  17 :         aux1=aux0/j2
  18 :         aux2=1.5D0*j3/j2
  19 :       endif
  20 :       do i=1,n
  21 :          t_s(i)=aux1*(aux2*q_s(i)-t_s(i))
  22 :          q_s(i)=aux0*q_s(i)
  23 :       enddo
  24 :       end
  25 :  
  26 :  
  27 :