Esope source files

   1 : C DVCAIN    SOURCE    CHAT      05/01/12    22:58:54     5004           
   2 : c
   3 :       SUBROUTINE DVFLCAMACINV (SIGINV,VARINT,OTRVAR,XPP,XQQ,XQT,XTT)
   4 :       IMPLICIT INTEGER(I-N)
   5 :       real*8    siginv(3),xpp,xqq,xqt,xtt,q,qqa,fr,pcap,gtet,
   6 :      .          gtetd,gtetdd,pcapfunc,etaconfunc,
   7 :      .          wwfunc,wwfuncder,wwfuncderder,varint,otrvar
   8 :       real*8    pc0,exr,cca,phi,alp,pcc,paa,qaa,exm,eww,
   9 :      .          cpr,cpm,cfr
  10 :       common /MRScapp/ pc0,exr,cca,phi,alp,pcc,paa,qaa,exm,eww,
  11 :      .                 cpr,cpm,cfr
  12 :       pcap  =pcapfunc(varint)
  13 :       gtet  =wwfunc(siginv(3),eww)
  14 :       gtetd =wwfuncder(siginv(3),eww)
  15 :       gtetdd=wwfuncderder(siginv(3),eww)
  16 :       fr    =(cfr*pcap-pcc)*etaconfunc(otrvar)
  17 :       q     =siginv(2)
  18 :       qqa   =(1.D0+q/qaa)
  19 :       xpp=2.D0/(cpr*pcap)**2
  20 :       xqq=(2.D0*gtet*gtet/fr/fr)*(qqa**(2.D0*exm)+
  21 :      .     exm*(2.D0*exm-1.D0)*q*q*qqa**(2.D0*exm-2.D0)/qaa/qaa+
  22 :      .     4.D0*exm*q*qqa**(2.D0*exm-1.D0)/qaa)
  23 :       xqt=(4.D0*gtet*gtetd*q/fr/fr)*
  24 :      .    (qqa**(2.D0*exm)+exm*q*qqa**(2.D0*exm-1.D0)/qaa)
  25 :       xtt=(2.D0*q*q*qqa**(2.D0*exm))/fr/fr*(gtetd**2+gtet*gtetdd)
  26 :       return
  27 :       end
  28 :  
  29 :  
  30 :