Télécharger atrcon.eso

Retour à la liste

Numérotation des lignes :

  1. C ATRCON SOURCE CB215821 17/11/30 21:15:11 9639
  2. SUBROUTINE ATRCON(SANGLE,COUR,RI,RE,Q,TMIN,TMAX,
  3. &X,Z,AX,AY,AZ)
  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 potentiel vecteur induit par un tronc de cone de courant
  11. C COUR, 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 (AX,AY,AZ)
  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 FTRC2
  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/(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. AX=0.
  79. AY=0.
  80. AZ=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=(FTRC2(Z0+H,Q2C1,QC2,C3)-FTRC2(Z0,Q2C1,QC2,C3))
  90. & *DTE(I)
  91. HX=HR*SINT
  92. HY=HR*COST
  93. AX=AX-HX
  94. AY=AY+HY
  95. C
  96. 21 CONTINUE
  97. C
  98. AX=AX*COEF
  99. AY=AY*COEF
  100. C
  101. DTE(IDEB)=DTINI
  102. DTE(IFIN)=DTFIN
  103. TETM(IDEB)=TETINI
  104. TETM(IFIN)=TETFIN
  105. C
  106. RETURN
  107. END
  108.  
  109.  
  110.  
  111.  
  112.  
  113.  

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