anchor.diff:     intens = 0.0
      intband = 0.0
      sunit = 0
  +   acoef_ = 0
  +   indexf_ = 0
  +   indexi_ = 0
      !
      Nintens = 0
      !
  
positive.code:        !
         do kitem = 1,maxitems 
           !
           l = len(trim(quantum_numbers(kitem,ilevelf)))
           !
           b_fmt = "(1x,a3)" ; if (l>3) b_fmt = "(1x,a8)"
           !
           write(tunit,b_fmt,advance="no"), trim(quantum_numbers(kitem,ilevelf))
           !
  
negative.code:                 !
                  b_fmt = "(1x,a3)" ; if (l>3) b_fmt = "(1x,a8)"
                  !
                  write(b_fmt,"('(1x,a',i1,')')") nchars_quanta(kitem)
                  !
                  write(out,b_fmt,advance="no"), trim(quantum_numbers(kitem,ilevelf))
                  !
                enddo
  