anchor.diff:              !
               !if (info>0) cycle
               !
  -            if (nint(Ji-QN%Jref)==0) then 
  +            if (nint(Ji-QN%Jref)<=0) then 
                 Nvib_states = Nvib_states + 1
               endif 
               !  info = 1
  
positive.code:            !
             ! check of Jref is appropriate, integer/half integer
             !
             if (nint(Ji-QN%Jref)==0) then
               !
               ivib = ivib + 1 
               !           
  
negative.code:                   !read(quantum_numbers(QN%dens_col-4,ileveli),*) ndensity
                    !
                    abscoef=cmcoef*acoef*gtot(ilevelf)*exp(-c2/temp*ener_rot)*ndensity*&
                            (1.0_rk-exp(-c2/temp*tranfreq))/(tranfreq**2*partfunc)
                    !
                    ! split into a product of vib and rot parts 
                 else
  