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

      subroutine sccfc(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 1300 mhz.
c-----------------------------------------------------------------------
c
      include 'ucom.h'
c
      dimension fc(14),beta(6),fct(6,14),wb(6,14),a(6),b(6),c(6),iop(2)
      dimension tab(3)
c-----------------------------------------------------------------------
      data nte/6/
      data iflagd/0/
c     fourier coefficients for side-coupled cavities
      data beta(1),(fct(1,i),i=1,14)
     1 / .50,  1.7528799,   .6026982,  .0956489,  -.0007633,
     2    -.0047612,  -.0012945,  -.0001396,   .0000224,   .0000119,
     3     .0000021,   .0000001,  -.0000001,  -.0000000,   .0000000/
      data beta(2),(fct(2,i),i=1,14)
     1 / .60,  1.7640248,   .6204240,   .0513374,  -.0325247,
     2    -.0134447,  -.0014003,   .0006422,   .0002947,   .0000346,
     3    -.0000129,  -.0000063,  -.0000008,   .0000002,   .0000001/
      data beta(3),(fct(3,i),i=1,14)
     1 / .70,  1.7541771,   .5444473,  -.0529861,  -.0708168,
     2    -.0098096,   .0054686,   .0022115,  -.0002226,  -.0002832,
     3     .0000049,   .0000307,  -.0000010,  -.0000048,   .0000005/
      data beta(4),(fct(4,i),i=1,14)
     1 / .80,  1.7389319,   .4496130,  -.1549319,  -.0900213,
     2     .0094835,   .0149095,   .0006300,  -.0022191,  -.0003900,
     3     .0002888,   .0000932,  -.0000322,  -.0000169,   .0000029/
      data beta(5),(fct(5,i),i=1,14)
     1 / .90,  1.7255494,   .3671034,  -.2342443,  -.0869619,
     2     .0341585,   .0178781,  -.0047510,  -.0034920,   .0006180,
     3     .0006624,  -.0000661,  -.0001196,   .0000064,   .0000224/
      data beta(6),(fct(6,i),i=1,14)
     1 /1.00,  1.6879765,   .1966623,  -.3083845,  -.0259680,
     2     .0667456,   .0009515,  -.0142383,   .0009093,   .0029194,
     3    -.0004541,  -.0005850,   .0001495,   .0001138,  -.0000412/
      if(iflagd.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(1),fct(1,i),wb(1,i),iop(1),1,a(1),b(1),c(1))
   10 continue
   20 pi=4.*atan(1.)
      do 30 i=1,14
      call spl1d2 (nte,beta(1),fct(1,i),wb(1,i),1,betac,tab(1))
   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++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++*