Esope source files

   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 :