Esope source files

   1 : C DVCOIN    SOURCE    CHAT      05/01/12    22:58:58     5004           
   2 : c
   3 :       SUBROUTINE DVFLCOMACINV (SIGINV,VARINT,OTRVAR,XPP,XQQ,XQT,XTT)
   4 :       IMPLICIT INTEGER(I-N)
   5 :       real*8    siginv(3),xpp,xqq,xqt,xtt,auxq,wwd0,wwd1,wwd2,q,
   6 :      .          wwfunc,wwfuncder,wwfuncderder,pcap,pcapfunc,p,ap,
   7 :      .          varint,otrvar,etaconfunc
   8 :       real*8    exm,qaa,pcc,eww,xk1,xk2,etb,exv,xep,paa,cco,
   9 :      .          exl,xnn,aaa,bbb,pc0,exr,cca,phi,alp,pcc2,paa2,qaa2,
  10 :      .          exm2,eww2,cpr,cpm,cfr
  11 :       common /MRScone/ exm,qaa,pcc,eww,xk1,xk2,etb,exv,xep,paa,cco,
  12 :      .                 exl,xnn,aaa,bbb
  13 :       common /MRScapp/ pc0,exr,cca,phi,alp,pcc2,paa2,qaa2,exm2,eww2,
  14 :      .                 cpr,cpm,cfr
  15 :       wwd0=wwfunc(siginv(3),eww)
  16 :       wwd1=wwfuncder(siginv(3),eww)
  17 :       wwd2=wwfuncderder(siginv(3),eww)
  18 :       pcap=pcapfunc(otrvar)
  19 :       p=siginv(1)
  20 :       ap=alp*pcap
  21 :       q=siginv(2)
  22 :       auxq=(1.D0+q/qaa)
  23 :       xpp=abs(xnn)*etaconfunc(varint)*(2.D0*ap-2.D0*pcc)/
  24 :      .    (p+ap-2.D0*pcc)**2
  25 :       xqq=(wwd0*(exm-1.D0)*exm*q*auxq**(exm-2.D0))/(qaa*qaa)+
  26 :      .    (2.D0*wwd0*exm*auxq**(exm-1.D0))/qaa
  27 :       xqt=wwd1*auxq**exm+(wwd1*exm*q*auxq**(exm-1.D0))/qaa
  28 :       xtt=q*auxq**exm*wwd2
  29 :       return
  30 :       end
  31 :  
  32 :  
  33 :