anchor.diff:                    !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*ener_rot*tranfreq))/(tranfreq**2*partfunc)
  +                          (1.0_rk-exp(-c2/temp*tranfreq))/(tranfreq**2*partfunc)
                     !
                     ! split into a product of vib and rot parts 
                  else
  
positive.code:            !
             ! check of Jref is appropriate, integer/half integer
             !
             if (nint(Ji-QN%Jref)==0) then
               !
               ivib = ivib + 1 
               !           
  
negative.code:             !
              !if (info>0) cycle
              !
              if (nint(Ji-QN%Jref)<=0) then 
                Nvib_states = Nvib_states + 1
              endif 
              !  info = 1
  