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

      subroutine rtmfc(betac,fc)
c     return the 14 fourier coefficients for the side-coupled cavity
c     with beta=betac.  the coefficients were determined from superfish
c     runs at 2380 mhz.
c-----------------------------------------------------------------------
      dimension fc(14),beta(10),fct(10,14),wb(10,14),a(10),b(10),c(10),
     1  iop(2),tab(3)
      data nte/10/
      data iflagc/0/
c-----------------------------------------------------------------------
c     fourier coefficients for side-coupled cavities
      data beta(1),(fct(1,i),i=1,14)
     1 /  .55,  1.8162297,   .8202260,   .1140768,  -.0584496,
     2   -.0336300,  -.0049842,   .0021335,   .0012629,   .0001983,
     3   -.0000658,  -.0000385,  -.0000048,   .0000021,   .0000007/
      data beta(2),(fct(2,i),i=1,14)
     1 /  .60,  1.8124060,   .7848132,   .0516917,  -.0955583,
     2   -.0379814,   .0010221,   .0055377,   .0015881,  -.0002348,
     3   -.0002880,  -.0000617,   .0000150,   .0000116,   .0000019/
      data beta(3),(fct(3,i),i=1,14)
     1 /  .65,  1.8022308,   .7233542,  -.0174431,  -.1191904,
     2   -.0328106,   .0065496,   .0058709,   .0007324,  -.0004280,
     3   -.0001649,   .0000036,   .0000115,   .0000009,  -.0000006/
      data beta(4),(fct(4,i),i=1,14)
     1 /  .70,  1.7905805,   .6480495,  -.1110423,  -.1537068,
     2   -.0203962,   .0213467,   .0089263,  -.0013231,  -.0017684,
     3   -.0002145,   .0002390,   .0000967,  -.0000121,  -.0000180/
      data beta(5),(fct(5,i),i=1,14)
     1 /  .75,  1.7714580,   .5402790,  -.2114199,  -.1572789,
     2    .0175976,   .0346638,   .0010967,  -.0071227,  -.0012653,
     3    .0012362,   .0004182,  -.0001846,  -.0001021,   .0000184/
      data beta(6),(fct(6,i),i=1,14)
     1 /  .80,  1.7587262,   .4692263,  -.2725046,  -.1512699,
     2    .0417247,   .0363340,  -.0061000,  -.0079966,   .0010337,
     3    .0017934,  -.0001993,  -.0004341,   .0000382,   .0001169/
      data beta(7),(fct(7,i),i=1,14)
     1 /  .85,  1.7418385,   .3826707,  -.3312878,  -.1309669,
     2    .0694703,   .0355463,  -.0137761,  -.0090094,   .0025033,
     3    .0021663,  -.0003746,  -.0004995,   .0000317,   .0001075/
      data beta(8),(fct(8,i),i=1,14)
     1 /  .90,  1.7210430,   .2811841,  -.3830460,  -.0831754,
     2    .1023592,   .0178458,  -.0266090,  -.0024255,   .0069314,
     3   -.0001750,  -.0018239,   .0002990,   .0004796,  -.0001707/
      data beta(9),(fct(9,i),i=1,14)
     1 /  .95,  1.7047580,   .2064316,  -.4116058,  -.0459696,
     2    .1180010,   .0047165,  -.0315829,   .0012915,   .0079611,
     3   -.0009923,  -.0018505,   .0003863,   .0004041,  -.0001152/
      data beta(10),(fct(10,i),i=1,14)
     1 / 1.00,  1.6841249,   .1129266,  -.4366304,   .0183279,
     2    .1320116,  -.0293359,  -.0344015,   .0159707,   .0070763,
     3   -.0062417,  -.0007853,   .0019354,  -.0001402,  -.0004452/
      if(iflagc.ne.0)go to 20
c     fit splines to fct vs beta for each term
      do 10 i=1,14
      iop(1)=1
      iop(2)=1
      wb(1,i)=0.
      wb(nte,i)=0.
      call spl1d1(nte,beta,fct(1,i),wb(1,i),iop,1,a,b,c)
   10 continue
   20 pi=4.*atan(1.)
      do 30 i=1,14
      call spl1d2 (nte,beta,fct(1,i),wb(1,i),1,betac,tab)
   30 fc(i)=tab(1)
c     normalize coefficients
      sum=0.
      t=1.
      sgn=1.
      do 40 l=1,14
      sum=sum+fc(l)*sgn/t
      t=t+2.
   40 sgn=-sgn
      fac=2.*sum/pi
      do 50 l=1,14
   50 fc(l)=fc(l)/fac
      return
      end
c++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++*