anchor.diff:        do 41 j=0,i
            write(22,2201) glm(i,j)
            write(22,2202) hlm(i,j)
  - 2201 format(2x,f15.5)
  - 2202 format(4x,f15.5)
  + 2201 format(2x,f10.5)
  + 2202 format(4x,f10.5)
      41 continue
      43 continue
   
  
positive.code: c   define MAG harmonic-Gauss harmonic conversion factors
  c   conalm,conblm
  c
        conalm=fact1
        conblm=-fact1
  
  c      if (l .le. 0 .or. m .lt. 0 .or. m .gt. l) then
  c       write(6,'(''bad l or m in getgauss'')')
  
negative.code: 
  c           return
  c      else ! form dimensionless fully normalized potential coeffs
  c        aalm(l,m)=glm(l,m)/(anano*escale*fact2*conalm)
  c        ablm(l,m)=hlm(l,m)/(anano*escale*fact2*conblm)
  
  c      write(21,2106) l,m,aalm(l,m),ablm(l,m),glm(l,m),hlm(l,m)
  c 2206 format(/,2x,2i3,2x,2(f15.5),2x,2(f15.5))
  
  c      return
  c      endif
  