anchor.diff:        include 'com5.f'
         include 'com8.f'
         
  -      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)
  +      integer la(nlma),ma(nlma)
  +      real*8 alm(lmax,lmax),blm(lmax,lmax)
  +      real*8 glm(lmax,lmax),hlm(lmax,lmax)
  +      real*8 aalm(lmax,lmax),ablm(lmax,lmax)
  +      real*8 fact1,fact2
   c
   c   constants
   c
  
positive.code: c        aalm(l,m)=glm(l,m)/(anano*escale*fact2*conalm)
  c        ablm(l,m)=hlm(l,m)/(anano*escale*fact2*conblm)
  
  c   wirte a header for cg. file
  
        write(22,2200) nlma,lmax,minc,r(1),r(kc),time/tscale
   2200 format(/, 2x,"nlma=",i3,2x,"lmax=",i3,2x,"minc=",i3,2x,
       $   "r(1)=",f7.4,2x,"r(kc)=",f7.4,2x,"time/tscale=",f9.6)
  
  c      write(22,2206) 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
  
negative.code:       write(21,2100) nlma,lmax,minc,r(1),r(kc),time/tscale
   2100 format(/, 2x,"nlma=",i3,2x,"lmax=",i3,2x,"minc=",i3,2x,
       $   "r(1)=",f7.4,2x,"r(kc)=",f7.4,2x,"time/tscale=",f9.6)
  
  c
  c  write data to cc file
  c
        write(21,2101) (b(lm,1),lm=1,nlma)
  c      write(21,2102) (w(lm,kc),lm=1,nlma)
  