Télécharger aplqci.eso

Retour à la liste

Numérotation des lignes :

  1. C APLQCI SOURCE CHAT 05/01/12 21:22:49 5004
  2. SUBROUTINE APLQCI(SANGLE,COUR,RI,RE,TMIN,TMAX,X,Z,AX,AY,AZ)
  3. IMPLICIT INTEGER(I-N)
  4. IMPLICIT REAL*8(A-H,O-Z)
  5. C
  6. -INC CCREEL
  7. C
  8. C*******************************************************************
  9. C Calcul du potentiel vecteur induit par une couronne circulaire
  10. C de courant COUR, d'axe Oz :
  11. C - rayons interieur et exterieur RI, RE
  12. C - extension angulaire TMIN a TMAX par rapport au plan meridien
  13. C xOz contenant le point de coordonnees (X,0,Z)
  14. C ou on calcule le champ (AX,AY,AZ)
  15. C*******************************************************************
  16. C
  17. SEGMENT SANGLE
  18. REAL*8 TETM(NT4)
  19. REAL*8 TETI(NT4)
  20. REAL*8 DTE(NT4)
  21. ENDSEGMENT
  22. C
  23. C champ magnetique selon Ox, Oy, Or et Oz cree par l'element d'arc de
  24. C longueur DTE situe en TETA, et grandeurs dependant de TETA pour
  25. C l'integration
  26. C
  27. C
  28. C
  29. C calcul des grandeurs geometriques independantes de TETA
  30. C
  31. AMU=4.E-7*XPI
  32. A=(RE-RI)/2.
  33. R=(RE+RI)/2.
  34. COEF=COUR*AMU/(8.*XPI*A)
  35. C
  36. C
  37. C recherche des bornes d'integration
  38. C
  39. IDEB=0
  40. IFIN=0
  41. TTAI=-1.
  42. TTAF=-1.
  43. NT4=TETM(/1)
  44. C
  45. DO 1 I=1,NT4
  46. DMIN=TMIN-TETI(I)
  47. DMAX=TMAX-TETI(I)
  48. TTBI=SIGN(1.D0,DMIN)
  49. TTBF=SIGN(1.D0,DMAX)
  50. TESTI=0.5*(1.-TTAI*TTBI)
  51. TESTF=0.5*(1.-TTAF*TTBF)
  52. IDEB=IDEB+I*INT(TESTI)
  53. IFIN=IFIN+I*INT(TESTF)
  54. TTAI=TTBI
  55. TTAF=TTBF
  56. 1 CONTINUE
  57. C
  58. DTINI=DTE(IDEB)
  59. DTFIN=DTE(IFIN)
  60. TETINI=TETM(IDEB)
  61. TETFIN=TETM(IFIN)
  62. DTE(IDEB)=TETI(IDEB+1)-TMIN
  63. DTE(IFIN)=TMAX-TETI(IFIN)
  64. TETM(IDEB)=0.5*(TETI(IDEB+1)+TMIN)
  65. TETM(IFIN)=0.5*(TETI(IFIN)+TMAX)
  66. C
  67. AX=0.
  68. AY=0.
  69. AZ=0.
  70. DO 21 I=IDEB,IFIN
  71. COST=COS(TETM(I))
  72. SINT=SIN(TETM(I))
  73. DFPL=FPLQ5(X,COST,Z,RE)-FPLQ5(X,COST,Z,RI)
  74. HX=SINT*DFPL
  75. HY=COST*DFPL
  76. AX=AX-HX*DTE(I)
  77. AY=AY+HY*DTE(I)
  78. 21 CONTINUE
  79. AX=AX*COEF
  80. AY=AY*COEF
  81. C
  82. DTE(IDEB)=DTINI
  83. DTE(IFIN)=DTFIN
  84. TETM(IDEB)=TETINI
  85. TETM(IFIN)=TETFIN
  86. C
  87. RETURN
  88. END
  89.  
  90.  
  91.  
  92.  
  93.  

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