Esope source files

   1 : C BTRCON    SOURCE    CHAT      05/01/12    21:44:27     5004           
   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 :