anchor.diff:             !
              nfiles = 1 
              call reada(intfilename(1))
  +           !
  +           ! for J-dependent histograms a J-value is expected in the last column 
  +           !
  +           if (histogramJ)   call readi(intJvalue(1))
              !
              cycle 
              !
  
positive.code:           intens(ipoint)=intens(ipoint)+abscoef
            !
          end select
          !
          cycle
       20  exit
  
negative.code:            if (i>nfiles_max) call report ("Too many files, increase nfiles_max"//trim(w),.true.)
             !
             intfilename(i) = trim(v)
             !
             ! for J-dependent histograms a J-value is expected in the last column 
             !
             if (histogramJ) then 
               call readi(intJvalue(i))
             endif 
             !
             call read_line(eof) ; if (eof) exit
             call reada(w)
  