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