anchor.diff:        dimension la(nlma),ma(nlma)
         dimension alm(lmax,lmax),blm(lmax,lmax)
         dimension glm(lmax,lmax),hlm(lmax,lmax)
  +      dimension aalm(lmax,lmax),ablm(lmax,lmax)
   c
   c   constants
   c
  
positive.code:        tl3=-ma(lm)*(lmax+1)/minc
         la(lm)=tl1+tl2+tl3
  c PRINT lm, la(lm), ma(lm) HERE
  c      write(22,2202) lm, la(lm), ma(lm)
  c 2202 format(/,2x,i3,2x,i4,2x,i4)	  
     36 continue
  c
  
  
negative.code: c for the poloidal magnetic potential b(lm,1) into real and
  c imaginary part and assign two new indices
  c  
        do 39 lm=1,nlma
         c1=real(b(lm,1))
         c2=aimag(b(lm,1))
  c assign new indices
  