Télécharger la plaquette Cast3M


fr en

Télécharger btrcon.eso

Retour à la liste

Numérotation des lignes : 

btrcon

  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