source: PSPA/parmelaPSPA/trunk/genfc.f @ 315

Last change on this file since 315 was 12, checked in by lemeur, 12 years ago

parmela pspa initial

File size: 4.5 KB
Line 
1      subroutine genfc(betac,fc)
2c---return 14 fourier coefficients for chalk river cavity with beta=betac
3c---coefficients determined from superfish runs at 9000 mhz.
4c--------------------------------------------------------------------------
5      save
6c
7      include 'constcom.h'
8      include 'ucom.h'
9c
10      dimension fc(14),bpb(8),fcpb(8,14),bfc(7),fcfc(7,14)
11      dimension fcpherm(14)
12c--------------------------------------------------------------------------
13c*
14c--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/
47c  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/
76c---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
88c++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++*
Note: See TracBrowser for help on using the repository browser.