Sources Esope | Cast3M

vflcov
C VFLCOV    SOURCE    CHAT      05/01/13    04:06:25     5004
 
      SUBROUTINE VFLCOMACINV (SIGINV,VARINT,OTRVAR,XP,XQ,XT)
      IMPLICIT INTEGER(I-N)
      real*8    siginv(3),varint,otrvar,xp,xq,xt,aux1,aux2,
     .          wwfunc,etaconfunc,wwfuncder,pcap,pcapfunc,p,q,ap
      real*8    exm,qaa,pcc,eww,xk1,xk2,etb,exv,xep,paa,cco,
     .          exl,xnn,aaa,bbb,pc0,exr,cca,phi,alp,pcc2,paa2,qaa2,
     .          exm2,eww2,cpr,cpm,cfr
      common /MRScone/ exm,qaa,pcc,eww,xk1,xk2,etb,exv,xep,paa,cco,
     .                 exl,xnn,aaa,bbb
      common /MRScapp/ pc0,exr,cca,phi,alp,pcc2,paa2,qaa2,exm2,eww2,
     .                 cpr,cpm,cfr
      p=siginv(1)
      q=siginv(2)
      aux1=wwfunc(siginv(3),eww)
      aux2=(1.D0+q/qaa)
      pcap=pcapfunc(otrvar)
      ap=alp*pcap
      xp=abs(xnn)*etaconfunc(varint)*(p-ap)/(p+ap-2.D0*pcc)
      xq=aux1*aux2**exm+aux1*exm/qaa*aux2**(exm-1.D0)*q
      xt=aux2**exm*wwfuncder(siginv(3),eww)*q
      return
      end