dirDATA='/data/colors/sdss-code/sdss-debias/data/'

  nbA = 187
  nbC = 14
  frac=fltarr(nbC,nbA,nbH-1) ;class, aArr, H

  minH=15.50
  binD = 1   &  hArr =findgen(14)+3   ;-LOG D
  nbH  = (hArr(n_elements(hArr)-1)-hArr(0))/binD + 1
  nbH = 13

 ;----MAKING THE LAST H BIN  15-15.5 instead of 15-16
  hArr(13) = 15.5

  outh= ['3','4','5','6','7','8','9','10','11','12','13','14','15'] 

  numHun  = fltarr(nbH)
  fracHun = fltarr(nbC)

  for z=0,nbH-2 do begin
     readfmt,dirDATA+'outhisto-num-'+outh(z)+'-ALL.txt','12F11.4',aArr,a,b,c,d,e,f,g,h,i,j,k,skipline=1, /silent


     frac( 0,*,z)=a             ;a
     frac( 1,*,z)=b             ;b
     frac( 2,*,z)=c             ;c
     frac( 3,*,z)=d             ;d
     frac( 4,*,z)=e             ;k
     frac( 5,*,z)=f             ;l
     frac( 6,*,z)=g             ;q
     frac( 7,*,z)=h             ;r
     frac( 8,*,z)=i             ;s
     frac( 9,*,z)=j             ;v
     frac(10,*,z)=k             ;x
     
  endfor


;--------------------- READ THE MPCORB
  if not keyword_set(H_mpc) then begin

   ;-Read MPCORB (need orbital elements for each classified object)
    readfmt,dirDATA+'MPCORB.DAT', $
            'A7,1x,F5,45x,'+$     	;-Absolute Mag
            '1x,F9,2x,F9,13x,F11,'+$ 	;-orbital elements
            '24x,F4,36x,'+$
            'A7,1x,A18',$       ;,1x,F8', $F4',$ ;      ;-Object Tag, Number, Name
            junk, H_mpc,  $
            i_mpc, e_mpc, a_mpc, $
            arc1,$
            num_mpc,name_mpc, $ ;,date_mpc, $
            /SILENT
;    nbMOC = n_elements( num_mpc )
  
    peri=a_mpc*(1-e_mpc)
    keep=where(peri gt 1.67 and a_mpc le 5.5, nbMPC)
    
    H_mpc   = H_mpc(keep)
    a_mpc   = a_mpc(keep)
    e_mpc   = e_mpc(keep)
    i_mpc   = i_mpc(keep)
    num_mpc = num_mpc(keep)
    name_mpc= name_mpc(keep)

  endif





 ;--- FOR EACH semi major axis BIN -- A in Hungaria
  hun = where( aArr lt 2, nbAinH)

  for kA=0,nbAinH-2 do begin ; loop through each semi-maj bin


   ;-current MPC pbject in that bin
    cA = where( a_mpc ge aArr(kA) and a_mpc lt aArr(kA+1), nbCurA ) ; find objs in each a bin
;    numMPC(kA) = nbCurA

    if cA(0) ne -1 then begin    

      for kH=0,nbH-2 do begin  ; kH=7 starts at H=10 loop only through H>10.

       ;-current objects with A bin and H bin
        cH = where( H_mpc(cA) ge hArr(kH) and H_mpc(cA) lt hArr(kH+1), nBin ) ; find objs in each H bin
        if cH(0) ne -1 then begin
;
;                massC(kH*10+kH2,kC,kA)=(classD(choice(kC))*vol) * nbin * frac(choice(kC),kA,kH) * frackeep(kA,kH)
           numHun[kH] += nBin


           for kC=0, 10 do begin
             fracHun[kC] += nBin * frac(kC,kA,kH)
           endfor

        endif
      endfor
    endif
  endfor


  cgPlot, hArr, numHun


end