Télécharger chabar.eso

Retour à la liste

Numérotation des lignes :

  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