anchor.diff:  c for the poloidal magnetic potential b(lm,1) into real and
   c imaginary part and assign two new indices
   c  
  -      do 37 lm=1,nlma
  +      do 39 lm=1,nlma
          c1=real(b(lm,1))
          c2=aimag(b(lm,1))
   c assign new indices
  
positive.code:        alm(l,m)=c1
         blm(l,m)=c2
  
  c   Convertion starts here
  c   the following code converts between Gauss coefficients (glm, hlm)
  c   and the spherical harmonic coefficients (alm,blm)
  c   of the magnetic potential at harmonic degree l
  
negative.code:       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
  