1 | subroutine genfc(betac,fc) |
---|
2 | c---return 14 fourier coefficients for chalk river cavity with beta=betac |
---|
3 | c---coefficients determined from superfish runs at 9000 mhz. |
---|
4 | c-------------------------------------------------------------------------- |
---|
5 | save |
---|
6 | c |
---|
7 | include 'constcom.h' |
---|
8 | include 'ucom.h' |
---|
9 | c |
---|
10 | dimension fc(14),bpb(8),fcpb(8,14),bfc(7),fcfc(7,14) |
---|
11 | dimension fcpherm(14) |
---|
12 | c-------------------------------------------------------------------------- |
---|
13 | c* |
---|
14 | c--fourieicients for croac pillbox |
---|
15 | data bpb(1),(fcpb(1,i),i=1,14) |
---|
16 | 1 / .30, 1.6981021, .3662608, -.0374986, -.0155366, |
---|
17 | 2 .0010105, .0005714, .0000118, -.0000016, .0000002, |
---|
18 | 3 -.0000009, .0000000, .0000002, -.0000001, -.0000001/ |
---|
19 | data bpb(2),(fcpb(2,i),i=1,14) |
---|
20 | 1 / .40, 1.6549106, .1754135, -.1274388, .0110193, |
---|
21 | 2 .0104464, -.0037673, -.0005845, .0006565, -.0000161, |
---|
22 | 3 -.0001009, .0000136, .0000132, -.0000026, -.0000014/ |
---|
23 | data bpb(3),(fcpb(3,i),i=1,14) |
---|
24 | 1 / .50, 1.6120163, .0178196, -.1418520, .0527615, |
---|
25 | 2 -.0014920, -.0064592, .0027495, -.0002412, -.0002598, |
---|
26 | 3 .0001155, -.0000057, -.0000090, -.0000009, .0000030/ |
---|
27 | data bpb(4),(fcpb(4,i),i=1,14) |
---|
28 | 1 / .60, 1.5735265, -.1010958, -.1124353, .0820944, |
---|
29 | 2 -.0261573, -.0007972, .0054602, -.0028485, .0005526, |
---|
30 | 3 .0002084, -.0001970, .0000638, .0000036, -.0000134/ |
---|
31 | data bpb(5),(fcpb(5,i),i=1,14) |
---|
32 | 1 / .70, 1.5411520, -.1843384, -.0670372, .0864748, |
---|
33 | 2 -.0476563, .0143338, .0009418, -.0040842, .0028027, |
---|
34 | 3 -.0010348, .0001308, .0001436, -.0001164, .0000520/ |
---|
35 | data bpb(6),(fcpb(6,i),i=1,14) |
---|
36 | 1 / .80, 1.5142367, -.2442793, -.0209533, .0760618, |
---|
37 | 2 -.0599447, .0309120, -.0103484, .0002884, .0029203, |
---|
38 | 3 -.0029874, .0020227, -.0010015, .0002743, .0001098/ |
---|
39 | data bpb(7),(fcpb(7,i),i=1,14) |
---|
40 | 1 / .90, 1.4917677, -.2864900, .0199212, .0576660, |
---|
41 | 2 -.0610854, .0408767, -.0200853, .0065650, -.0000880, |
---|
42 | 3 -.0019109, .0018311, -.0012107, .0006312, -.0003039/ |
---|
43 | data bpb(8),(fcpb(8,i),i=1,14) |
---|
44 | 1 / 1.00, 1.4746846, -.3155308, .0512825, .0397525, |
---|
45 | 2 -.0583078, .0470460, -.0283883, .0128171, -.0032004, |
---|
46 | 3 -.0012525, .0024120, -.0020447, .0012211, -.0005597/ |
---|
47 | c fourier coefficients for croac |
---|
48 | data bfc(1),(fcfc(1,i),i=1,14) |
---|
49 | 1 / .40, 1.7892894, .6857973, .0430116, -.0098477, |
---|
50 | 2 -.0026227, -.0043393, -.0010201, -.0022380, -.0011799, |
---|
51 | 3 .0004360, .0003532, .0000438, .0000046, -.0000096/ |
---|
52 | data bfc(2),(fcfc(2,i),i=1,14) |
---|
53 | 1 / .50, 1.7160123, .3352655, -.1762127, -.0164347, |
---|
54 | 2 -.0183020, -.0079387, .0102261 , .0001315, .0005291, |
---|
55 | 3 .0007131, -.0011001 , -.0004295 , .0003561, .0001671/ |
---|
56 | data bfc(3),(fcfc(3,i),i=1,14) |
---|
57 | 1 / .60, 1.6692360, .1337539 , -.2745378, .0105648, |
---|
58 | 2 .0273667, .0067394, .0002206 , -.0023561 , -.0027576, |
---|
59 | 3 -.0001529, .0031047, -.0008543, -.0011834, .0006437/ |
---|
60 | data bfc(4),(fcfc(4,i),i=1,14) |
---|
61 | 1 / .70, 1.6172791, -.0585426, -.2732258, .1142691, |
---|
62 | 2 .0282834, -.0302966, -.0028203, .0103538, -.0034405, |
---|
63 | 3 -.0008613, .0011889, -.0008872, .0007548, -.0003886/ |
---|
64 | data bfc(5),(fcfc(5,i),i=1,14) |
---|
65 | 1 / .80, 1.5757450 , -.1905119 , -.2102027, .1857369, |
---|
66 | 2 -.0428980, -.0295982, .0277137, -.0076888 , -.0024615, |
---|
67 | 3 .0033408, -.0018291, .0008091, -.0002476, -.0001563/ |
---|
68 | data bfc(6),(fcfc(6,i),i=1,14) |
---|
69 | 1 / .90, 1.5511584, -.2490323, -.1520595, .1927862 , |
---|
70 | 2 -.0789607, -.0074175, .0215380 , -.0071498, .0021416, |
---|
71 | 3 -.0064930, .0070933, -.0013774 , -.0040154, .0042958/ |
---|
72 | data bfc(7),(fcfc(7,i),i=1,14) |
---|
73 | 1 /1.00, 1.5203916, -.3152512, -.0761705, .1846429, |
---|
74 | 2 -.1254612, .0316493, .0201999, -.0251474, .0110114, |
---|
75 | 3 .0000393, -.0026879, .0010176, .0005092, -.0006987/ |
---|
76 | c---fourier coefficients for phermex |
---|
77 | data (fcpherm(i),i=1,14) |
---|
78 | 1 / 1.4392398, -.3917328, .1409807, -.0177985, |
---|
79 | 2 -.0502258, .0836854, -.0922970, .0834457, -.0633808, |
---|
80 | 3 .0381197, -.0126771, -.0087312, .0234238, -.0300528/ |
---|
81 | if(ltype.eq.4)call interp(8,bpb,fcpb,betac,fc) |
---|
82 | if(ltype.eq.3)call interp(7,bfc,fcfc,betac,fc) |
---|
83 | if(ltype.ne.5)return |
---|
84 | do 10 i=1,14 |
---|
85 | 10 fc(i)=fcpherm(i) |
---|
86 | return |
---|
87 | end |
---|
88 | c++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++* |
---|