Esope source files

   1 : C AARCTR    SOURCE    CHAT      05/01/12    21:15:09     5004           
   2 :       SUBROUTINE AARCTR(SANGLE,COUR,HI,HE,RI,RE,QB,QH,TMIN,TMAX,
   3 :      &XP,ZP,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  Calcul du potentiel vecteur induit par un arc de cercle de courant COUR,
  10 : C  d'axe Oz, de section trapezoidale asymetrique par rapport a xOy:
  11 : C     - rayons interieur et exterieur RI et RE
  12 : C     - hauteurs interieure et exterieure HI et HE
  13 : C     - pentes basse et haute QB et QH
  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,0)
  17 : C*******************************************************************
  18 : C
  19 :       SEGMENT SANGLE
  20 :          REAL*8 TETM(NT4)
  21 :          REAL*8 TETI(NT4)
  22 :          REAL*8 DTE(NT4)
  23 :       ENDSEGMENT
  24 : C
  25 : C  champ magnetique selon Ox, Oy, Or et Oz cree par l'element d'arc de
  26 : C  longueur DTE situe en TETA, et grandeurs dependant de TETA pour
  27 : C  l'integration
  28 : C
  29 : C
  30 :       EXTERNAL FARCA
  31 : C
  32 : C  calcul des grandeurs geometriques independantes de TETA
  33 : C
  34 :       X=XP
  35 :       Z=ZP
  36 :       AMU=4.E-7*XPI
  37 :       R=(RI+RE)/2.
  38 :       A=(RE-RI)/2.
  39 :       B=(HI+HE)/4.
  40 :       COEF=COUR*AMU/(16.*XPI*A*B)
  41 :       RE2=RE**2
  42 :       RI2=RI**2
  43 : C
  44 :       QH2=QH**2
  45 :       AH=1.+QH2
  46 :       QHR=QH*R
  47 :       QHRPZ=QHR+Z
  48 :       QHZ=QH*Z
  49 :       CH=X**2+B**2+QHR**2+Z**2-2.*B*QHRPZ+2.*QHR*Z
  50 :       DH=QH
  51 :       EH=B-QHRPZ
  52 : C
  53 :       QB2=QB**2
  54 :       AB=1.+QB2
  55 :       QBR=QB*R
  56 :       QBRPZ=QBR+Z
  57 :       QBZ=QB*Z
  58 :       CB=X**2+B**2+QBR**2+Z**2+2.*B*QBRPZ+2.*QBR*Z
  59 :       DB=QB
  60 :       EB=-B-QBRPZ
  61 : C
  62 : C  recherche des bornes d'integration
  63 : C
  64 :       IDEB=0 
  65 :       IFIN=0
  66 :       TTAI=-1.     
  67 :       TTAF=-1.
  68 :       NT4=TETM(/1)
  69 : C
  70 :       DO 1 I=1,NT4
  71 :          DMIN=TMIN-TETI(I)      
  72 :          DMAX=TMAX-TETI(I)
  73 :          TTBI=SIGN(1.D0,DMIN)
  74 :          TTBF=SIGN(1.D0,DMAX)
  75 :          TESTI=0.5*(1.-TTAI*TTBI) 
  76 :          TESTF=0.5*(1.-TTAF*TTBF) 
  77 :          IDEB=IDEB+I*INT(TESTI)
  78 :          IFIN=IFIN+I*INT(TESTF)
  79 :          TTAI=TTBI
  80 :          TTAF=TTBF
  81 :  1    CONTINUE   
  82 : C
  83 :       DTINI=DTE(IDEB)
  84 :       DTFIN=DTE(IFIN)
  85 :       TETINI=TETM(IDEB)
  86 :       TETFIN=TETM(IFIN)
  87 :       DTE(IDEB)=TETI(IDEB+1)-TMIN
  88 :       DTE(IFIN)=TMAX-TETI(IFIN)
  89 :       TETM(IDEB)=0.5*(TETI(IDEB+1)+TMIN)
  90 :       TETM(IFIN)=0.5*(TETI(IFIN)+TMAX)
  91 : C
  92 :       AX=0.
  93 :       AY=0.
  94 :       AZ=0.
  95 : C
  96 :       DO 21 I=IDEB,IFIN
  97 :          COST=COS(TETM(I))
  98 :          SINT=SIN(TETM(I))
  99 :          BH=-2.*(-B*QH+QH*QHRPZ+X*COST)
 100 :          BB=-2.*(B*QB+QB*QBRPZ+X*COST)
 101 : C
 102 :          HR11=RE2*LOG(SQRT(AH*RE2+BH*RE+CH)+DH*RE+EH)
 103 :      &       -RI2*LOG(SQRT(AH*RI2+BH*RI+CH)+DH*RI+EH)
 104 :          HR12=RE2*LOG(SQRT(AB*RE2+BB*RE+CB)+DB*RE+EB)
 105 :      &       -RI2*LOG(SQRT(AB*RI2+BB*RI+CB)+DB*RI+EB)
 106 :          HR1=HR11-HR12
 107 :          HR1=0.5*HR1
 108 :          HR211=FARCA(RE,AH,BH,CH,DH,EH,X,COST,SINT)
 109 :          HR212=FARCA(RI,AH,BH,CH,DH,EH,X,COST,SINT)
 110 :          HR21=HR211-HR212
 111 :          HR221=FARCA(RE,AB,BB,CB,DB,EB,X,COST,SINT)
 112 :          HR222=FARCA(RI,AB,BB,CB,DB,EB,X,COST,SINT)
 113 :          HR22=HR221-HR222
 114 :          HR2=HR21-HR22
 115 :          HR=HR1-HR2
 116 :          HR=HR*DTE(I)
 117 : C
 118 :          HX=HR*SINT
 119 :          HY=HR*COST
 120 : C
 121 :          AX=AX-HX
 122 :          AY=AY+HY
 123 : C
 124 :  21   CONTINUE
 125 : C
 126 :       AX=AX*COEF
 127 :       AY=AY*COEF
 128 : C
 129 :       DTE(IDEB)=DTINI
 130 :       DTE(IFIN)=DTFIN
 131 :       TETM(IDEB)=TETINI
 132 :       TETM(IFIN)=TETFIN
 133 : C
 134 :       RETURN
 135 :       END
 136 :  
 137 :  
 138 :  
 139 :