anchor.diff:         do
            !
            ! Scan and find Jmax
  -         read(bunit,*,end=14) gamma_,n_,J_
  +         read(bunit,*,end=14) ch_broad(1:3),gamma_,n_,J_
            !
            Jmax = max(Jmax,J_)
            !
  
positive.code:    ! open and count number of lines (levels) in the Energy files
     !
     if (trim(enrfilename)/="none".and..not.hitran_do) then
        !
        write(ioname, '(a)') 'Energy file'
        call IOstart(trim(ioname),enunit)
  
negative.code:         !
        endif
        !
        if (verbose>=2) print('(1x,a/)'),'... reading energies done!'
        !
     endif  ! end of processing the states files 
     !   
     if (verbose>=2.or.(partfunc_do.and.verbose>0)) print('(1x,a,1x,es16.8/)'),'! partition function value is',partfunc
  