- double precision b1,b2,cc,dd,ee,ctilde,dtilde,b2tilde
- integer nloctyp
- common/fourier/ b1(2,maxtor),b2(2,maxtor),cc(2,2,maxtor),
- & dd(2,2,maxtor),ee(2,2,maxtor),ctilde(2,2,maxtor),
- & dtilde(2,2,maxtor),b1tilde(2,maxtor),nloctyp
+ double precision b1,b2,cc,dd,ee,ctilde,dtilde,b2tilde,b1tilde,
+ & b,bnew1,bnew2,eenew,gtb1,gtb2,eeold,gtee
+ integer nloctyp,iloctyp(-ntyp1:ntyp1),itype2loc(-ntyp1:ntyp1)
+ common/fourier/ b1(2,maxres),b2(2,maxres),b(13,-ntyp:ntyp),
+ & bnew1(3,2,-ntyp:ntyp),bnew2(3,2,-ntyp:ntyp),
+ & cc(2,2,-ntyp:ntyp),
+ & dd(2,2,-ntyp:ntyp),eeold(2,2,-ntyp:ntyp),
+ & eenew(2,-ntyp:ntyp),
+ & ee(2,2,maxres),
+ & ctilde(2,2,-ntyp:ntyp),
+ & dtilde(2,2,-ntyp:ntyp),b1tilde(2,maxres),
+ & b2tilde(2,maxres),
+ & gtb1(2,maxres),gtb2(2,maxres),gtEE(2,2,maxres),
+ & nloctyp,iloctyp,itype2loc