Sources Esope | Cast3M

vysgp2
C VYSGP2    SOURCE    CB215821  16/04/21    21:18:40     8920
 
      SUBROUTINE VYISIGPOWDER2(SIG,VECM)
      IMPLICIT INTEGER(I-N)
      integer i
      real*8  sig(3),vecm(3),aux2,aux3,xi1,xj2
      real*8               nn1,nn2,eta0,aa1,aa2,bb1,bb2,
     .                     pia,pib,sigy0,pcc,rrr
      common /poder2_data/ nn1,nn2,eta0,aa1,aa2,bb1,bb2,
     .                     pia,pib,sigy0,pcc,rrr
      call Invari_I1(sig,3,xi1)
      call desviador(sig,vecm,3)
      if (xi1.lt.pia) then
         aux2 = 2.D0*(xi1-pia)/rrr**2
         do i=1,3
           vecm(i)= vecm(i)+aux2
         enddo
      else
         aux2 = bb1/3.D0
         call Invari_J2(sig,3,xj2)
         aux3 = 1.D0 / SQRT(2.D0*xj2)
         do i=1,3
           vecm(i)= aux3*vecm(i)+aux2
         enddo
      endif
      return
      end