Sources Esope | Cast3M

vfsgmh
C VFSGMH    SOURCE    CHAT      05/01/13    04:06:56     5004
 
      SUBROUTINE VFLSIGMIEHE (SIG,VECM)
      IMPLICIT INTEGER(I-N)
      integer m,i
      real*8  sig(3),vecm(3)
      real*8  desv(3),aux1,aux2,alpha1,alpha2,alpha3,vecmaux(3)
      real*8  valor
      real*8            sigy0,kiso,siginf,velo,cpar,mpar
      common /miehdata/ sigy0,kiso,siginf,velo,cpar,mpar
      m=nint(mpar)
      call yieldmieheb(sig,valor)
      aux1=valor**(1-2*m)/2.D0*sqrt(2.D0/3.D0)
      aux2=aux1*cpar*3.D0**(2*m)/(2.D0**(2*m-1)+1.D0)
      aux1=aux1*(1.D0-cpar)
      call Desviador(sig,desv,3)
      alpha1=(desv(1)-desv(2))**(2*m-1)
      alpha2=(desv(2)-desv(3))**(2*m-1)
      alpha3=(desv(3)-desv(1))**(2*m-1)
      vecmaux(1)=alpha1-alpha3
      vecmaux(2)=alpha2-alpha1
      vecmaux(3)=alpha3-alpha2
      do i=1,3
       vecmaux(i)=aux1*vecmaux(i)+aux2*desv(i)**(2*m-1)
      enddo
      call Desviador(vecmaux,vecm,3)
      return
      end