Numi extraction channel at Recycler ring 20 inch quad (for PANDIRA) Original version from the 1987 User's Guide, Chapter 10.5 6-29-99 simple panowski quad 08-23-00 put in tables for Ferrite and Sm-Cobalt 04-14-03 update to PANDIRA 7.0 04-15-03 build model of magnets in extraction channel by 20 inch quad 04-22-03 add in MI dipole 05-07-03 add sheet steel between quad and EPB ; Copyright 1998, by the University of California. ; Unauthorized commercial use is prohibited. ® kprob=0, ; Defines a POISSON or PANDIRA problem ;dx=0.08, ;dy=0.08, ; Mesh size xreg = -6.0, -4.5, 4.5, 11.0, 17.0, ; make smaller mesh near pole tip kreg = 0, 25, 30, 160, 180, 240, kmax = 245, ; mesh points yreg = -25.0, -13.0 -4.5, 4.5, 10.0, lreg = 0, 100, 155, 235, 275, 285, lmax = 290, ; mesh points conv = 2.54, ; inches nbslo = 0, ; Dirchelt boundary condtion at bottom lmin = 10 ,ltop = 30, ; limits on where to calculate fields for outpan ktype = 2, ; symmetry midplane for this nterm = 11, ; number of harmonics terms nptc = 172, ; number of arc points rint = 0.39370, ; radius to calculate harmonics rnorm = 0.39370, ; normilization radius angle = 360.0, ; angle to calculate harmonics mode=0 & ; Materials have variable permeability ; entire universe &po x= -30.000000, y= -68.000000 & &po x= 30.000000, y= -68.000000 & &po x= 30.000000, y= 20.000000 & &po x= -30.000000, y= 20.000000 & &po x= -30.000000, y= -68.000000 & ® mat=2, & ;FLUX RETURN Quadrant 1 &po x= 3.625000, y= 0.000000 & &po x= 4.125000, y= 0.000000 & &po x= 4.125000, y= 4.125000 & &po x= 0.000000, y= 4.125000 & &po x= 0.000000, y= 3.625000 & &po x= 3.625000, y= 3.625000 & &po x= 3.625000, y= 0.000000 & ® mat=2, & ;FLUX RETURN Quadrant 2 &po x= -3.625000, y= 0.000000 & &po x= -4.125000, y= 0.000000 & &po x= -4.125000, y= 4.125000 & &po x= -0.000000, y= 4.125000 & &po x= -0.000000, y= 3.625000 & &po x= -3.625000, y= 3.625000 & &po x= -3.625000, y= 0.000000 & ® mat=2, & ;FLUX RETURN Quadrant 3 &po x= -3.625000, y= -0.000000 & &po x= -4.125000, y= -0.000000 & &po x= -4.125000, y= -4.125000 & &po x= -0.000000, y= -4.125000 & &po x= -0.000000, y= -3.625000 & &po x= -3.625000, y= -3.625000 & &po x= -3.625000, y= -0.000000 & ® mat=2, & ;FLUX RETURN Quadrant 4 &po x= 3.625000, y= -0.000000 & &po x= 4.125000, y= -0.000000 & &po x= 4.125000, y= -4.125000 & &po x= 0.000000, y= -4.125000 & &po x= 0.000000, y= -3.625000 & &po x= 3.625000, y= -3.625000 & &po x= 3.625000, y= -0.000000 & ® mat=2, & ;POLETIP quadrant 1 &po x= 2.625000, y= 2.625000 & &po x= 0.488870, y= 2.625000 & &po x= 0.524735, y= 2.478689 & &po x= 0.566746, y= 2.332379 & &po x= 0.612703, y= 2.186069 & &po x= 0.662147, y= 2.039758 & &po x= 0.716002, y= 1.893448 & &po x= 0.776327, y= 1.747137 & &po x= 0.846187, y= 1.600827 & &po x= 0.929703, y= 1.454516 & &po x= 1.032382, y= 1.308206 & &po x= 1.161895, y= 1.161895 & &po x= 1.308206, y= 1.032382 & &po x= 1.454516, y= 0.929703 & &po x= 1.600827, y= 0.846187 & &po x= 1.747137, y= 0.776327 & &po x= 1.893448, y= 0.716002 & &po x= 2.039758, y= 0.662147 & &po x= 2.186069, y= 0.612703 & &po x= 2.332379, y= 0.566746 & &po x= 2.478689, y= 0.524735 & &po x= 2.625000, y= 0.488870 & &po x= 2.625000, y= 2.625000 & ® mat=2, & ;POLETIP quadrant 2 &po x= -2.625000, y= 2.625000 & &po x= -0.488870, y= 2.625000 & &po x= -0.524735, y= 2.478689 & &po x= -0.566746, y= 2.332379 & &po x= -0.612703, y= 2.186069 & &po x= -0.662147, y= 2.039758 & &po x= -0.716002, y= 1.893448 & &po x= -0.776327, y= 1.747137 & &po x= -0.846187, y= 1.600827 & &po x= -0.929703, y= 1.454516 & &po x= -1.032382, y= 1.308206 & &po x= -1.161895, y= 1.161895 & &po x= -1.308206, y= 1.032382 & &po x= -1.454516, y= 0.929703 & &po x= -1.600827, y= 0.846187 & &po x= -1.747137, y= 0.776327 & &po x= -1.893448, y= 0.716002 & &po x= -2.039758, y= 0.662147 & &po x= -2.186069, y= 0.612703 & &po x= -2.332379, y= 0.566746 & &po x= -2.478689, y= 0.524735 & &po x= -2.625000, y= 0.488870 & &po x= -2.625000, y= 2.625000 & ® mat=2, & ;POLETIP quadrant 3 &po x= -2.625000, y= -2.625000 & &po x= -0.488870, y= -2.625000 & &po x= -0.524735, y= -2.478689 & &po x= -0.566746, y= -2.332379 & &po x= -0.612703, y= -2.186069 & &po x= -0.662147, y= -2.039758 & &po x= -0.716002, y= -1.893448 & &po x= -0.776327, y= -1.747137 & &po x= -0.846187, y= -1.600827 & &po x= -0.929703, y= -1.454516 & &po x= -1.032382, y= -1.308206 & &po x= -1.161895, y= -1.161895 & &po x= -1.308206, y= -1.032382 & &po x= -1.454516, y= -0.929703 & &po x= -1.600827, y= -0.846187 & &po x= -1.747137, y= -0.776327 & &po x= -1.893448, y= -0.716002 & &po x= -2.039758, y= -0.662147 & &po x= -2.186069, y= -0.612703 & &po x= -2.332379, y= -0.566746 & &po x= -2.478689, y= -0.524735 & &po x= -2.625000, y= -0.488870 & &po x= -2.625000, y= -2.625000 & ® mat=2, & ;POLETIP quadrant 4 &po x= 2.625000, y= -2.625000 & &po x= 0.488870, y= -2.625000 & &po x= 0.524735, y= -2.478689 & &po x= 0.566746, y= -2.332379 & &po x= 0.612703, y= -2.186069 & &po x= 0.662147, y= -2.039758 & &po x= 0.716002, y= -1.893448 & &po x= 0.776327, y= -1.747137 & &po x= 0.846187, y= -1.600827 & &po x= 0.929703, y= -1.454516 & &po x= 1.032382, y= -1.308206 & &po x= 1.161895, y= -1.161895 & &po x= 1.308206, y= -1.032382 & &po x= 1.454516, y= -0.929703 & &po x= 1.600827, y= -0.846187 & &po x= 1.747137, y= -0.776327 & &po x= 1.893448, y= -0.716002 & &po x= 2.039758, y= -0.662147 & &po x= 2.186069, y= -0.612703 & &po x= 2.332379, y= -0.566746 & &po x= 2.478689, y= -0.524735 & &po x= 2.625000, y= -0.488870 & &po x= 2.625000, y= -2.625000 & ® mat=8, mshape=1, mtid=2 & ;BRICK for Quadrant 1 RHS &po x= 2.625000, y= 2.625000 & &po x= 2.625000, y= 0.625000 & &po x= 3.625000, y= 0.625000 & &po x= 3.625000, y= 2.625000 & &po x= 2.625000, y= 2.625000 & ® mat=9, mshape=1, mtid=3 & ;BRICK for Quadrant 1 top &po x= 2.625000, y= 2.625000 & &po x= 0.625000, y= 2.625000 & &po x= 0.625000, y= 3.625000 & &po x= 2.625000, y= 3.625000 & &po x= 2.625000, y= 2.625000 & ® mat=10, mshape=1, mtid=4 & ; BRICK for Quadrant 2 LHS &po x= -2.625000, y= 2.625000 & &po x= -2.625000, y= 0.625000 & &po x= -3.625000, y= 0.625000 & &po x= -3.625000, y= 2.625000 & &po x= -2.625000, y= 2.625000 & ® mat=11, mshape=1, mtid=5 & ;BRICK for Quadrant 2 top &po x= -2.625000, y= 2.625000 & &po x= -0.625000, y= 2.625000 & &po x= -0.625000, y= 3.625000 & &po x= -2.625000, y= 3.625000 & &po x= -2.625000, y= 2.625000 & ® mat=12, mshape=1, mtid=6 & ;BRICK for Quadrant 3 LHS &po x= -2.625000, y= -2.625000 & &po x= -2.625000, y= -0.625000 & &po x= -3.625000, y= -0.625000 & &po x= -3.625000, y= -2.625000 & &po x= -2.625000, y= -2.625000 & ® mat=13, mshape=1, mtid=7 & ;BRICK for Quadrant 3 Bottom &po x= -2.625000, y= -2.625000 & &po x= -0.625000, y= -2.625000 & &po x= -0.625000, y= -3.625000 & &po x= -2.625000, y= -3.625000 & &po x= -2.625000, y= -2.625000 & ® mat=14, mshape=1, mtid=8 & ;BRICK for Quadrant 4RHS &po x= 2.625000, y= -2.625000 & &po x= 2.625000, y= -0.625000 & &po x= 3.625000, y= -0.625000 & &po x= 3.625000, y= -2.625000 & &po x= 2.625000, y= -2.625000 & ® mat=15, mshape=1, mtid=9 & ;BRICK for Quadrant 4 bottom &po x= 2.625000, y= -2.625000 & &po x= 0.625000, y= -2.625000 & &po x= 0.625000, y= -3.625000 & &po x= 2.625000, y= -3.625000 & &po x= 2.625000, y= -2.625000 & ;****************************************************************************** ® mat=2, & ; adde 1 inch thick sheet steel between quad and EPB &po x = -5.500, y = -5.000 & &po x = 15.500, y = -5.000 & &po x = 15.500, y = -6.000 & &po x = -5.500, y = -6.000 & &po x = -5.500, y = -5.000 & ;!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! ;Iron for EPB ® mat=2, & &PO X= 2.521, Y= -28.357 & &PO X= 14.640, Y= -17.911 & &PO X= 6.479, Y= -8.443 & &PO X= -5.640, Y= -18.889 & &PO X= 2.521, Y= -28.357 & ;LHS upper coil ® mat=1,cur 29120.0 & &PO X= 0.223, Y= -17.268 & &PO X= -1.973, Y= -19.161 & &PO X= -0.080, Y= -21.358 & &PO X= 2.117, Y= -19.464 & &PO X= 0.223, Y= -17.268 & ;RHS Lower Coil ® mat=1,cur -29120.0 & &PO X= 8.777, Y= -19.532 & &PO X= 10.973, Y= -17.639 & &PO X= 9.080, Y= -15.442 & &PO X= 6.883, Y= -17.336 & &PO X= 8.777, Y= -19.532 & ;LHS lower coil ® mat=1,cur 29120.0 & &PO X= 4.989, Y= -22.797 & &PO X= 2.793, Y= -24.690 & &PO X= 0.899, Y= -22.494 & &PO X= 3.096, Y= -20.600 & &PO X= 4.989, Y= -22.797 & ;RHS upper coil ® mat=1,cur -29120.0 & &PO X= 4.011, Y= -14.003 & &PO X= 6.207, Y= -12.110 & &PO X= 8.101, Y= -14.306 & &PO X= 5.904, Y= -16.200 & &PO X= 4.011, Y= -14.003 & ;Gap for beam ® mat=1 & &PO X= 0.899, Y= -22.494 & &PO X= 9.080, Y= -15.442 & &PO X= 8.101, Y= -14.306 & &PO X= -0.080, Y= -21.358 & &PO X= 0.899, Y= -22.494 & ;########################################################## ; Main injector dipole ;Upper ® mat = 2, & &PO X= -13.000, Y= -44.750 & &PO X= 13.000, Y= -44.750 & &PO X= 16.500, Y= -47.750 & &PO X= 16.500, Y= -56.000 & &PO X= 10.375, Y= -56.000 & &PO X= 10.375, Y= -50.500 & &PO X= 5.375, Y= -50.500 & &PO X= 3.250, Y= -55.000 & &PO X= -3.250, Y= -55.000 & &PO X= -5.375, Y= -50.500 & &PO X= -10.375, Y= -50.500 & &PO X= -10.375, Y= -56.000 & &PO X= -16.500, Y= -56.000 & &PO X= -16.500, Y= -47.750 & &PO X= -13.000, Y= -44.750 & ;Lower ® mat = 2, & &PO X= -13.000, Y= -67.250 & &PO X= 13.000, Y= -67.250 & &PO X= 16.500, Y= -64.250 & &PO X= 16.500, Y= -56.000 & &PO X= 10.375, Y= -56.000 & &PO X= 10.375, Y= -61.500 & &PO X= 5.375, Y= -61.500 & &PO X= 3.250, Y= -57.000 & &PO X= -3.250, Y= -57.000 & &PO X= -5.375, Y= -61.500 & &PO X= -10.375, Y= -61.500 & &PO X= -10.375, Y= -56.000 & &PO X= -16.500, Y= -56.000 & &PO X= -16.500, Y= -64.250 & &PO X= -13.000, Y= -67.250 & ;RHS Lower coil ® mat = 1, cur = -18488 & &PO X= 5.875, Y= -57.000 & &PO X= 5.875, Y= -61.000 & &PO X= 9.875, Y= -61.000 & &PO X= 9.875, Y= -57.000 & &PO X= 5.875, Y= -57.000 & ;RHS upper coil ® mat = 1, cur = -18488 & &PO X= 5.875, Y= -55.000 & &PO X= 5.875, Y= -51.000 & &PO X= 9.875, Y= -51.000 & &PO X= 9.875, Y= -55.000 & &PO X= 5.875, Y= -55.000 & ;LHS lower coil ® mat = 1, cur = 18488 & &PO X= -5.875, Y= -57.000 & &PO X= -5.875, Y= -61.000 & &PO X= -9.875, Y= -61.000 & &PO X= -9.875, Y= -57.000 & &PO X= -5.875, Y= -57.000 & ;LHS upper coil ® mat = 1, cur = 18488 & &PO X= -5.875, Y= -55.000 & &PO X= -5.875, Y= -51.000 & &PO X= -9.875, Y= -51.000 & &PO X= -9.875, Y= -55.000 & &PO X= -5.875, Y= -55.000 & &mt mtid=1,aeasy=90,gamper=1 ; Material ID, Easy axis direction, Gamma along hard axis hcept=-3500,bcept=3800. & ; H and B intercepts &mt mtid=2,aeasy= 180,gamper=1,hcept=-3500,bcept=3800. & ; Sr Ferrite no comp &mt mtid=3,aeasy= 270,gamper=1,hcept=-3500,bcept=3800. & ; Sr Ferrite no comp &mt mtid=4,aeasy= 180,gamper=1,hcept=-3500,bcept=3800. & ; Sr Ferrite no comp &mt mtid=5,aeasy= 90,gamper=1,hcept=-3500,bcept=3800. & ; Sr Ferrite no comp &mt mtid=6,aeasy= 0,gamper=1,hcept=-3500,bcept=3800. & ; Sr Ferrite no comp &mt mtid=7,aeasy= 90,gamper=1,hcept=-3500,bcept=3800. & ; Sr Ferrite no comp &mt mtid=8,aeasy= 0,gamper=1,hcept=-3500,bcept=3800. & ; Sr Ferrite no comp &mt mtid=9,aeasy= 270,gamper=1,hcept=-3500,bcept=3800. & ; Sr Ferrite no comp