anchor.diff:        do 31 lc=mc,nlaf
          lm=lm+1
          mclm(lm)=mc
  -c      write(22,2201) lm, mclm(lm)
  -c 2201 format(/, 2x,i3,2x,i4)     
  +c      write(21,2104) lm, mclm(lm)
  +c 2104 format(/, 2x,i3,2x,i4)     
      31 continue
      35 continue
   
  
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
  