anchor.diff:     integer(ik) :: imol,iostat_,isotope
      real(rk)    :: gf,gi,mem_t,temp_gamma_n,offset_
      character(55) ch_q,ch_broad
  +   character(len=cl)   :: my_fmt100K,my_fmt0
  +   real(real32) :: nan   
      !
      ln2=log(2.0_rk)
      ln22 = ln2*2.0_rk
  
positive.code:       !
        if (trim(proftype)=='LIFETIME'.or.trim(proftype)=='T-LIFETIME') THEN
          !
          ! allocate the matrix for the sum of the A-coeffs
          !
          allocate(Asum(nlines),stat=info)
          call ArrayStart('Asum',info,size(Asum),kind(Asum))
          !
          Asum = -1.0_rk
          !
        endif
        !
  
negative.code:          write(out,"('Vibrational energies prepared for non-LTE')")
           !
           if(vibpopulation_do) then
              write(print_fmt,'(a,i3,a)') '(i8,f15.6,1x,',ivib2-ivib1+1,'a5,2x,e11.4)'
           else
              write(print_fmt,'(a,i3,a)') '(i8,f15.6,1x,',ivib2-ivib1+1,'a5)'
           endif
           !
           do ivib  = 1,Nvib_states
  