Esope source files

   1 : C CVFLCO    SOURCE    CHAT      05/01/12    22:33:14     5004           
   2 : c
   3 :       SUBROUTINE C_VFLCOMACINV (SIGINV,VARINT,OTRVAR,XP,XQ,XT)
   4 :       IMPLICIT INTEGER(I-N)
   5 :       complex*16 siginv(3),varint,otrvar,xp,xq,xt,aux1,aux2,ap,
   6 :      .          c_wwfunc,c_wwfuncder,c_etaconfunc,p,pcap,c_pcapfunc
   7 :       real*8    exm,qaa,pcc,eww,xk1,xk2,etb,exv,xep,paa,cco,
   8 :      .          exl,xnn,aaa,bbb,pc0,exr,cca,phi,alp,pcc2,paa2,qaa2,
   9 :      .          exm2,eww2,cpr,cpm,cfr
  10 :       common /MRScone/ exm,qaa,pcc,eww,xk1,xk2,etb,exv,xep,paa,cco,
  11 :      .                 exl,xnn,aaa,bbb
  12 :       common /MRScapp/ pc0,exr,cca,phi,alp,pcc2,paa2,qaa2,exm2,eww2,
  13 :      .                 cpr,cpm,cfr
  14 :       aux1=c_wwfunc(siginv(3),eww)
  15 :       aux2=(1.D0+siginv(2)/qaa)
  16 :       pcap=c_pcapfunc(otrvar)
  17 :       p=siginv(1)
  18 :       ap=alp*pcap
  19 :       xp=abs(xnn)*c_etaconfunc(varint)*(p-ap)/(p+ap-2.D0*pcc)
  20 :       xq=aux1*aux2**exm+aux1*exm*siginv(2)/qaa*aux2**(exm-1.D0)
  21 :       xt=siginv(2)*aux2**exm*c_wwfuncder(siginv(3),eww)
  22 :       return
  23 :       end
  24 :  
  25 :  
  26 :