Télécharger cvflco.eso

Retour à la liste

Numérotation des lignes :

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

© Cast3M 2003 - Tous droits réservés.
Mentions légales