Télécharger btrcon.eso

Retour à la liste

Numérotation des lignes :

  1. C BTRCON SOURCE CB215821 17/11/30 21:15:17 9639
  2. SUBROUTINE BTRCON(SANGLE,COUR,RI,RE,Q,TMIN,TMAX,
  3. &X,Z,BX,BY,BZ)
  4. IMPLICIT INTEGER(I-N)
  5. IMPLICIT REAL*8(A-H,O-Z)
  6. C
  7. -INC CCREEL
  8. C
  9. C*******************************************************************
  10. C Calcul du champ induit par un tronc de cone de courant COUR,
  11. C d'axe Oz :
  12. C - rayons interieur et exterieur RI et RE
  13. C - pente Q
  14. C - extension angulaire TMIN a TMAX par rapport au plan meridien
  15. C xOz contenant le point de coordonnees (X,0,Z)
  16. C ou on calcule le champ (BX,BY,BZ)
  17. C*******************************************************************
  18. C
  19. C
  20. SEGMENT SANGLE
  21. REAL*8 TETM(NT4)
  22. REAL*8 TETI(NT4)
  23. REAL*8 DTE(NT4)
  24. ENDSEGMENT
  25. C
  26. C champ magnetique selon Ox, Oy, Or et Oz cree par l'element d'arc de
  27. C longueur DTE situe en TETA, et grandeurs dependant de TETA pour
  28. C l'integration
  29. C
  30. C
  31. EXTERNAL FTRC
  32. C
  33. C calcul des grandeurs geometriques independantes de TETA
  34. C
  35. AMU=4.E-7*XPI
  36. R=(RI+RE)/2.
  37. A=(RE-RI)/2.
  38. Z0=Q*RI
  39. C changement de repere
  40. ZR=Z+Q*R
  41. H=2.*Q*A
  42. Q2=Q**2
  43. COEF=AMU*COUR*Q/(4.*XPI*H)
  44. C1=X**2+ZR**2
  45. Q2C1=Q2*C1
  46. C3=1.+Q2
  47. C
  48. C recherche des bornes d'integration
  49. C
  50. IDEB=0
  51. IFIN=0
  52. TTAI=-1.
  53. TTAF=-1.
  54. NT4=TETM(/1)
  55. C
  56. DO 1 I=1,NT4
  57. DMIN=TMIN-TETI(I)
  58. DMAX=TMAX-TETI(I)
  59. TTBI=SIGN(1.D0,DMIN)
  60. TTBF=SIGN(1.D0,DMAX)
  61. TESTI=0.5*(1.-TTAI*TTBI)
  62. TESTF=0.5*(1.-TTAF*TTBF)
  63. IDEB=IDEB+I*INT(TESTI)
  64. IFIN=IFIN+I*INT(TESTF)
  65. TTAI=TTBI
  66. TTAF=TTBF
  67. 1 CONTINUE
  68. C
  69. DTINI=DTE(IDEB)
  70. DTFIN=DTE(IFIN)
  71. TETINI=TETM(IDEB)
  72. TETFIN=TETM(IFIN)
  73. DTE(IDEB)=TETI(IDEB+1)-TMIN
  74. DTE(IFIN)=TMAX-TETI(IFIN)
  75. TETM(IDEB)=0.5*(TETI(IDEB+1)+TMIN)
  76. TETM(IFIN)=0.5*(TETI(IFIN)+TMAX)
  77. C
  78. BX=0.
  79. BY=0.
  80. BZ=0.
  81. C
  82. DO 21 I=IDEB,IFIN
  83. C
  84. COST=COS(TETM(I))
  85. SINT=SIN(TETM(I))
  86. XC=X*COST
  87. C2=-2.*(XC+Q*ZR)
  88. QC2=Q*C2
  89. HR=Q*(-FTRC(Z0+H,Q2C1,QC2,C3,ZR)+FTRC(Z0,Q2C1,QC2,C3,ZR))
  90. & *DTE(I)
  91. C
  92. HX=HR*COST
  93. HY=HR*SINT
  94. HZ=(FTRC((Z0+H)/Q,C1,C2,C3,XC)-FTRC(Z0/Q,C1,C2,C3,XC))*DTE(I)
  95. C
  96. BX=BX+HX
  97. BY=BY+HY
  98. BZ=BZ+HZ
  99. C
  100. 21 CONTINUE
  101. C
  102. BX=BX*COEF
  103. BY=BY*COEF
  104. BZ=BZ*COEF
  105. C
  106. DTE(IDEB)=DTINI
  107. DTE(IFIN)=DTFIN
  108. TETM(IDEB)=TETINI
  109. TETM(IFIN)=TETFIN
  110. C
  111. RETURN
  112. END
  113.  
  114.  
  115.  
  116.  
  117.  
  118.  

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