Télécharger la plaquette Cast3M


fr en

Télécharger chabar.eso

Retour à la liste

Numérotation des lignes : 

chabar

  1. C CHABAR    SOURCE    CHAT      05/03/07    21:15:01     5041
  2.       SUBROUTINE CHABAR(IUNIT,CENTR,GY,DX,DZ,DENS,XP,YP,ZP,BXP,BYP,BZP)
  3.       IMPLICIT REAL*8 (A-H,O-Z)
  4.       IMPLICIT INTEGER(I-N)
  5. C-----------------------------------------------------------------
  6. CALCUL DU CHAMP D'UNE BARRE RECTANGULAIRE // A L AXE OX
  7. C     IUNIT INDIQUE LES UNITES CHOISIES:
  8. C                             IUNIT=0 OU 1   IUNIT=2         IUNIT=3
  9. C     DIMENSIONS EN                M          MM              MM
  10. C     DENSITE DE COURANT EN       A/M2       A/MM2           A/MM2
  11. C     CHAMP EN                   TESLA       TESLA           GAUSS
  12. C     FMUJ=MU*DENS/4*PI         1E-7*DENS   1E-4*DENS        DENS
  13. C
  14. C
  15. C-----------------------------------------------------------------
  16.       DIMENSION CENTR(3)
  17.       BXP=0.D0
  18.       BYP=0.D0
  19.       BZP=0.D0
  20.       UPSI=0.000001D0
  21.       IF(IUNIT.GT.1) UPSI=0.001D0
  22.       IF(GY.LT.UPSI) RETURN
  23.       X=XP-CENTR(1)
  24.       Y=YP-CENTR(2)
  25.       Z=ZP-CENTR(3)
  26. C
  27.       V=Y
  28.       Y=X
  29.       X=-V
  30. C
  31. CALCUL DE BZP  (ENTREE CHABAY)
  32. C
  33.       FMUJ= 0.0000001D0*DENS
  34.       IF(IUNIT.EQ.2) FMUJ=0.0001D0*DENS
  35.       IF(IUNIT.EQ.3) FMUJ=         DENS
  36.       ICOMP= 1
  37.       X1= X +DX/2.D0
  38.       X2= X -DX/2.D0
  39.       Y1= Y +GY/2.D0
  40.       Y2= Y -GY/2.D0
  41.       Z1= Z +DZ/2.D0
  42.       Z2= Z -DZ/2.D0
  43. C
  44. 10    H= 0.D0
  45.       IF(ABS(X2+X1).LT.UPSI) GOTO 11
  46.       SQ222= SQRT(X2**2+Y2**2+Z2**2)
  47.       SQ221= SQRT(X2**2+Y2**2+Z1**2)
  48.       SQ212= SQRT(X2**2+Y1**2+Z2**2)
  49.       SQ211= SQRT(X2**2+Y1**2+Z1**2)
  50.       SQ122= SQRT(X1**2+Y2**2+Z2**2)
  51.       SQ121= SQRT(X1**2+Y2**2+Z1**2)
  52.       SQ112= SQRT(X1**2+Y1**2+Z2**2)
  53.       SQ111= SQRT(X1**2+Y1**2+Z1**2)
  54.       TAG=FTAN(X2,Y2,Z2,SQ222)-FTAN(X1,Y2,Z2,SQ122)-FTAN(X2,Y1,Z2,SQ212)
  55.      $   -FTAN(X2,Y2,Z1,SQ221)+FTAN(X1,Y2,Z1,SQ121)+FTAN(X2,Y1,Z1,SQ211)
  56.      $   +FTAN(X1,Y1,Z2,SQ112)-FTAN(X1,Y1,Z1,SQ111)
  57.       H=FMUJ*(+FLOZ(UPSI,Z2,Y1,Y2,SQ222,SQ212,SQ122,SQ112)
  58.      $        -FLOZ(UPSI,Z1,Y1,Y2,SQ221,SQ211,SQ121,SQ111)
  59.      $        +FLOZ(UPSI,Y2,Z1,Z2,SQ222,SQ221,SQ122,SQ121)
  60.      $        -FLOZ(UPSI,Y1,Z1,Z2,SQ212,SQ211,SQ112,SQ111)+TAG)
  61. 11    IF(ICOMP.EQ.2) GOTO 12
  62.       BZP= H
  63. C
  64. CALCUL DE BXP  (ENTREE CHABAY)
  65. C
  66.       ICOMP=2
  67.       V = X1
  68.       X1=-Z1
  69.       Z1= V
  70.       V = X2
  71.       X2=-Z2
  72.       Z2= V
  73.       X1=-X1
  74.       X2=-X2
  75.       GOTO10
  76. 12    BXP=-H
  77. 101   V=BYP
  78.       BYP=-BXP
  79.       BXP=V
  80.       END
  81.  
  82.  

© Cast3M 2003 - Tous droits réservés.
Mentions légales