c Converts spherical harmonic coefficients from Edmonds to Stacey format
c     Edmonds, A.R. (1960) Angular Momentum and Quantum Mechanics
c     Stacey, F.D. (1977) Physics of the Earth
c
      real c(0:21,0:21),s(0:21,0:21)
      character infile*80,outfile*80
      pi=3.14159275359

      print *,'Enter input file name'
      read (*,'(a)') infile
      open (11,file=infile)

      read (11,*) nmax      
      do 50 i=0,nmax
         do 40 m=0,i
            read (11,*) is1,is2,c(i,m),s(i,m)
c            print *,i,m,is1,is2
            fact=sqrt(1./(2.*pi))*(-1.)**m
            if (m.eq.0) fact=fact/sqrt(2.)
            c(i,m)=c(i,m)*fact
            s(i,m)=-s(i,m)*fact
40       continue
50    continue
51    format (11x,2e15.7)
      close (11)

      print *,'Enter output file name'
      read (*,'(a)') outfile
      open (12,file=outfile)

      write (12,*) nmax
      do 100 i=0,nmax
         do 90 m=0,i
            write (12,101) i,m,c(i,m),s(i,m)
90       continue
100   continue
101   format (2i4,2f10.4)
      close (12)

      stop
      end