Esope source files

   1 : C APLQDR    SOURCE    CHAT      05/01/12    21:22:53     5004           
   2 :       SUBROUTINE APLQDR(COUR,RI,RE,Y1,Y2,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 plaque rectangulaire
  10 : C  de courant COUR,d'axe Oy:
  11 : C     - rayons interieur et exterieur RI et RE (RI=-RE)
  12 : C     - extension longitudinale Y1 a Y2 par rapport au plan xOz
  13 : C       contenant le point de coordonnees (X,0,Z)
  14 : C       et on calcule le champ (0,AY,0)
  15 : C*********************************************************************
  16 : C
  17 :       EXTERNAL FPLQ4
  18 : C
  19 : C  calcul des grandeurs geometriques reduites
  20 : C
  21 :       AMU=4.E-7*XPI
  22 :       A=(RE-RI)/2.
  23 :       COEF=COUR*AMU/(8.*XPI*A)
  24 :       AX=0.
  25 :       AY=0.
  26 :       AZ=0.
  27 : C
  28 :       Z2=Z**2
  29 :       AMX=A-X
  30 :       APX=A+X
  31 :       XMA2=AMX**2
  32 :       XPA2=APX**2
  33 :       DEXT1=XMA2+Z2
  34 :       DEXT2=XPA2+Z2
  35 :       IF(DEXT1*DEXT2.NE.0) THEN
  36 :          AY=COEF*(FPLQ4(AMX,Y2,Z)-FPLQ4(AMX,Y1,Z)
  37 :      &           -FPLQ4(-APX,Y2,Z)+FPLQ4(-APX,Y1,Z))
  38 :       ELSE
  39 :          AY=XGRAND
  40 :       ENDIF
  41 :       RETURN
  42 :       END
  43 :  
  44 :  
  45 :  
  46 :