anchor.diff:         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)	  
  +c      write(21,2105) lm, la(lm), ma(lm)
  +c 2105 format(/,2x,i3,2x,i4,2x,i4)	  
      36 continue
   c
   
  
positive.code:       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)     
     31 continue
     35 continue
  
  
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
  