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 :
© Cast3M 2003 - All rights reserved.
Disclaimer