Télécharger dvrhin.eso

Retour à la liste

Numérotation des lignes :

  1. C DVRHIN SOURCE CHAT 05/01/12 22:59:15 5004
  2.  
  3. SUBROUTINE DVFLRHMCINV (SIGINV,XPP,XQQ,XQT,XTT)
  4. IMPLICIT INTEGER(I-N)
  5. real*8 siginv(3),xqq,xqt,xtt,xpp,q,t
  6. real*8 kt,ktd,ktdd,aux,aux2,aux3,tat,cst,cst2,
  7. . ktetmc,derktetmc,derderktetmc,aa,bb,snt
  8. real*8 coh,phi,ar,tmax,aap,aan,bbp,bbn
  9. common /rhmcdata/ coh,phi,ar,tmax,aap,aan,bbp,bbn
  10. q=siginv(2)
  11. xqq=0.D0
  12. xqt=0.D0
  13. xtt=0.D0
  14. xpp=0.D0
  15. if (abs(q).lt.1.D-17) return
  16. call Invari_J3_a_t2(q,siginv(3),t)
  17. cst = cos(3.D0*t)
  18. kt = ktetMC(t)
  19. aux = sqrt(q*kt**2+(ar*coh)**2)
  20. if (abs(t).gt.tmax) then
  21. aa = aap
  22. bb = bbp
  23. if (t.lt.0.D0) then
  24. aa = aan
  25. bb = bbn
  26. endif
  27. snt = sin(3.D0*t)
  28. aux2= aa+2.D0*bb*snt
  29. aux3= aa-4.D0*bb*snt
  30. xqq = -(3.D0*bb*snt*aux3/q+(aux2*kt/aux)**2)/(4.D0*aux)
  31. xqt = (3.D0*sqrt(3.D0)*bb)/(4.d0*aux*sqrt(q))*
  32. . (-kt*kt*aux2/aux**2+(4.D0*bb*snt-aa)/q)
  33. xtt = 27.D0*bb**2*(1.D0/q-(kt/aux)**2)/(q*aux*4.D0)
  34. else
  35. ktd = derktetMC(t)
  36. ktdd= derderktetMC(t)
  37. tat = tan(3.D0*t)
  38. cst2= cst**2
  39. xqq =-(kt*ktd*tat-kt*kt)**2/(4.D0*aux**3)+
  40. . (tat)/(2.D0*aux*q)*(tat/2.D0*
  41. . (ktd*ktd+kt*ktdd)+kt*ktd*(1.5D0/cst2-1.D0))
  42. xqt = sqrt(3.D0)/(4.D0*cst*aux*sqrt(q))*
  43. . ((kt*kt*ktd*(kt-ktd*tat))/(aux*aux)+
  44. . (tat*(ktd*ktd+kt*ktdd)+kt*ktd*(3.D0/cst2-2.D0))/q)
  45. xtt = (0.75D0)/(q*aux*cst2)*
  46. . (-(kt*ktd/aux)**2+(ktd*ktd+kt*ktdd+3.D0*kt*ktd*tat)/q)
  47. endif
  48. return
  49. end
  50.  
  51.  
  52.  
  53.  

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