anchor.diff:  
         character*40 ttl
         character*2  titre
  -      common /char/ titre(30), ttl
  +      common /char2/ titre(30), ttl
   
         character*2  ttire(9)
         data ttire /'s ', 'p*', 'p ', 'd*', 'd ', 'f*', 'f ','g*', 'g '/
  
positive.code:          open(15,file=trim(header)//'chi.dat')
           do
              read(15,'(A)',end=311,err=311) buf
              dict_out = trim(dict_out)//trim(buf)//NEW_LINE('A')
           end do
   311     close(15)
        else
  
negative.code: c titre = identification of the wave functions  s,p*,p,........
        character*40 ttl
        character*2  titre
        common /char2/ titre(30), ttl
  
  c  -- This read commented out to make input easier, not used for
  c     PHASE calculations
  