Sources Esope | Cast3M

vycavr
C VYCAVR    SOURCE    CHAT      05/01/13    04:10:32     5004
 
      SUBROUTINE VYICAMACVAR (SIG,NDIMS,VARINT,OTRVAR,VALOR)
      IMPLICIT INTEGER(I-N)
      integer ndims
      real*8    siginv(3),varint,otrvar,valor,pcapfunc,pcapfuncder,
     .          pcap,pcapder,ecf,aux1,aux2,aux3,wwfunc,etaconfunc,
     .          sig(ndims)
      real*8    pc0,exr,cca,phi,alp,pcc,paa,qaa,exm,eww,
     .          cpr,cpm,cfr
      common /MRScapp/ pc0,exr,cca,phi,alp,pcc,paa,qaa,exm,eww,
     .                 cpr,cpm,cfr
      call InvariantesPQT(sig,ndims,siginv)
      pcap   =pcapfunc(varint)
      pcapder=pcapfuncder(varint)
      ecf    =etaconfunc(otrvar)
      aux1   =(siginv(1)-cpm*pcap)
      aux2   =siginv(2)*(1.D0+siginv(2)/qaa)**exm*wwfunc(siginv(3),eww)
      aux3   =ecf*(cfr*pcap-pcc)
      valor  =-2.D0*(aux1/(cpr*pcap)**3*(cpm*cpr*pcap+cpr*aux1)+
     .               ecf*aux2**2*cfr/aux3**3)*pcapder
      return
      end