anchor.diff:           ilevelf = ilevelf_ram(i)
            ileveli = ileveli_ram(i)
            gf   = gf_ram(i)
  +         abscoeff   = abscoeff_ram(i)
            energyf = energies(ilevelf)
            energyi = energies(ileveli)
            tranfreq = energyf - energyi
  
positive.code:            !
             if (nchars_>15) cycle
             !
             write(out,b_fmt,advance="no") trim(quantum_numbers(kitem,ilevelf))
             !
           enddo
           !
  
negative.code:            !
             if (nchars_>15) cycle
             !
             write(sunit,b_fmt,advance="no") trim(quantum_numbers(kitem,ileveli))
             !
           enddo
           !
           if ( mod(nint(2*Jf),2)==0 ) then
             !
             write(sunit,'(i7,1x,a1,1x,a1)',advance="no") nint(Jf),trim(quantum_numbers(1,ilevelf)),trim(quantum_numbers(2,ilevelf))
             write(sunit,'(i7,1x,a1,1x,a1)',advance="no") nint(Ji),trim(quantum_numbers(1,ileveli)),trim(quantum_numbers(2,ileveli))
             !
           else
             !
             write(sunit,'(f7.1,1x,a1,1x,a1)',advance="no") Jf,trim(quantum_numbers(1,ilevelf)),trim(quantum_numbers(2,ilevelf))
             write(sunit,'(f7.1,1x,a1,1x,a1)',advance="no") Ji,trim(quantum_numbers(1,ileveli)),trim(quantum_numbers(2,ileveli))
             !
           endif
           !
  