anchor.diff:             !
              if (nchars_>15) cycle
              !
  -           write(out,b_fmt,advance="no") trim(quantum_numbers(kitem,ileveli))
  +           write(sunit,b_fmt,advance="no") trim(quantum_numbers(kitem,ileveli))
              !
            enddo
            !
            if ( mod(nint(2*Jf),2)==0 ) then
              !
  -           write(out,'(i7,1x,a1,1x,a1)',advance="no") nint(Jf),trim(quantum_numbers(1,ilevelf)),trim(quantum_numbers(2,ilevelf))
  -           write(out,'(i7,1x,a1,1x,a1)',advance="no") nint(Ji),trim(quantum_numbers(1,ileveli)),trim(quantum_numbers(2,ileveli))
  +           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(out,'(f7.1,1x,a1,1x,a1)',advance="no") Jf,trim(quantum_numbers(1,ilevelf)),trim(quantum_numbers(2,ilevelf))
  -           write(out,'(f7.1,1x,a1,1x,a1)',advance="no") Ji,trim(quantum_numbers(1,ileveli)),trim(quantum_numbers(2,ileveli))
  +           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
            !
  
positive.code:          endif
           ilog_gamma2 = int(log(gamma2)*(1.0_rk/(log(10.0_rk)*0.001_rk))) + 2**14
           !
           write(gfunit,*) &
                       ielion,wl2ind,iener,iloggf,ilog_gamma1,ilog_gamma2,ilog_n
       enddo
       !
       !ielion, wl, gf-value, xi, igr, igs, igw
  
negative.code:             !
           enddo
           !
           write(sunit,'(1x,i1)',advance="no") ierror_nu
           write(sunit,'(i1)',advance="no") ierror_S
           write(sunit,'(i1)',advance="no") HITRAN_Air%ierr
           write(sunit,'(i1)',advance="no") HITRAN_Self%ierr
           write(sunit,'(i1)',advance="no") HITRAN_n%ierr
           write(sunit,'(i1)',advance="no") HITRAN_Delta%ierr
           !
           write(sunit,"(' 0 0 0 0 0 0 ',f7.1,f7.1)",advance="yes") real(gtot(ilevelf)),real(gtot(ileveli))
           !
         else
            !
            ! Not HITRAN: the main detailed stick-output is written to out and 
            ! also a simple two-column-version into .stick 
            !
            ! write to .stick-file
            write(sunit,my_fmt) tranfreq,abscoef
  