/BATCH,LIST
/title, plt-94s,vz, PLTRACE,EMAG with SOLID122 (E normal to B)
!/sho,jpeg
/PREP7
/COM
/COM This problem is created from plt-54s to add Electrostatic field
/COM in consideration of the particle tracing in magnetostatics.
/COM The Electrostatic problem is solved at the beginning with el122
/COM The second stage is solving the Magnetostatics with el96.
/COM To plot particle traces the new label is used in PLTRACE,EMAG.
/COM
!/NOPR
!/OUT,scratch
pi=3.141592653589793
h=0.         ! magnetic field for 2 Tesla flux density
charge=0.1            ! charge of particle
mass=0.01             ! mass of particle
x=1.5                   ! initial location of particle
y=1.5
z=2.0
vx=2                  ! velocity of particle
vy=10
vz=.1
!
/PREP7
ET,1,122
MP,perx,1,10
BLOCK,,3,,3,,3
ESIZE,0.2
VMESH,ALL
fini

/SOLU
ANTYPE,STAT
NSEL,S,LOC,Y,0
NSEL,R,LOC,X,0
D,ALL,VOLT,0.3e1

NSEL,S,LOC,Y,0
NSEL,R,LOC,X,3
D,ALL,VOLT,0.6e1

NSEL,S,LOC,Y,3
D,ALL,VOLT,0
NSEL,ALL
SOLVE
FINISH

/post1
set,1,1
nsel,all
esel,all
rsys,0
etable,efx,ef,x
etable,efy,ef,y
etable,efz,ef,z
/out,ef122-n
pretab
/out
finish

!
/PREP7
ET,1,96
MP,MURX,1,1
alls
DDELE,all,all
FINISH
!
/SOLU
ANTYPE,STAT
NSLE,S,CORNER
NSEL,R,LOC,X,0
D,ALL,MAG,0
NSLE,S,CORNER
NSEL,R,LOC,X,3
D,ALL,MAG,-h*3
NSEL,ALL
ldread,EF,1,1,,,,rth
SOLVE
FINISH

/out
/POST1
SET
/VIEW, 1, 0.2562    , 0.7335    , 0.6296
/ANG,  1, -14.00

!TRPOINT,x-.2,y,z,vx,vy,vz,charge,mass
TRPOINT,x,y,z,vx,vy,vz,charge,mass

PLTRACE,EMAG
/USER,  1   
/VIEW,  1,  0.562574149555    ,  0.622783779240    ,  0.543737703831
/ANG,   1,   2.39841772833  
/REPLO  
/VIEW,  1,  0.408832962040    ,  0.768944886893    ,  0.491507243151
/ANG,   1,  -2.11428610090  
/REPLO  
/VIEW,  1,  0.244073773275    ,  0.905013564007    ,  0.348394090309
/ANG,   1,  -6.10640119876  
/REPLO  
/VIEW,  1,  0.179669309142E-01,  0.995820447933    ,  0.895478915007E-01
/ANG,   1,  -32.2675758855  
/REPLO  
/VIEW,  1,  0.215656483206    ,  0.937252451279    ,  0.273952776630
/ANG,   1,  -5.06904730261  
/REPLO  
TRTIME,0,0, ,0.02,0
pltr,emag,,,,,,,,-3.8,1.5
/RGB,INDEX,100,100,100, 0
/RGB,INDEX, 80, 80, 80,13
/RGB,INDEX, 60, 60, 60,14
/RGB,INDEX, 0, 0, 0,15
/REPLOT
/PLOPTS,INFO,0
/PLOPTS,LEG1,0
/PLOPTS,LEG2,0
/PLOPTS,LEG3,0
/PLOPTS,FRAME,0
/PLOPTS,TITLE,0
/PLOPTS,MINM,0
/PLOPTS,FILE,0
/PLOPTS,LOGO,0
/PLOPTS,WINS,1
/PLOPTS,WP,0
/PLOPTS,DATE,0
/TRIAD,OFF
/REPLOT
!*
/DIST,1,0.729,1
/REP,FAST


/title, PLTRACE,EMAG with B=(3,0,0) and -4.5 < Escale < 4.5
/device,anim,2
/seg,delete
/seg,mult,,0.1
*do,i,1,30
pltr,emag,,,,,,,,-4.5+i*0.3,1.5
*enddo
*do,i,1,30
pltr,emag,,,,,,,,4.2-i*0.3,1.5
*enddo
/seg,off
anim
/anfile,save,,avi
fini