;-----------------------------------------------------------------------------------------------;
;-----------------------------------------------------------------------------------------------;
;---  TAG ---    I -- Initialization And Input Verification              -----------------------;
;-----------------------------------------------------------------------------------------------;
;-----------------------------------------------------------------------------------------------;

 ;----What to do
  kPopStart = 0
  readRS    = 0

  debug     = 0
      ;- 0 
      ;- 1 
      ;- 2 
      ;- 3 
;  if debug ne 0 then set_plot,'X'


 ;--I.1-- Define Directories  ------------------------------------------------------------------;
  spawn, 'hostname', host
  case strTrim(host) of
     'hyperion': dirEuclid = '/data/euclid/'
;     'endymion': dirEuclid = '/home/bcarry/work/data/euclid/'
     'endymion': dirEuclid = '/data/euclid/'
     else: stop
  endcase

  dirFilter = dirEuclid+'filters/'
  dirTaxo   = dirEuclid+'taxonomy/'
  dirStat   = dirEuclid + 'stat/'
  dirSFD    = dirStat + 'SFD/'
  dirFrac   = dirEuclid + 'fraction/'
  dirPop    = dirEuclid + 'synthetic/'
  dirRS     = dirEuclid + 'rs/'




 ;--I.3-- Solar Elongation of the Survey  ------------------------------------------------------;
  if readRS eq 1 or not keyword_set(rs) then $
    rs=euclidReadRS(dirRS+'wide.csv')
  cgHistoPlot, rs.SAA, binSize=binElong, loc=elong, histData=elongPDF
  euclidElong = [ [elong],[elongPDF] ]
  wdelete

;-----------------------------------------------------------------------------------------------;
;-----------------------------------------------------------------------------------------------;
;---  TAG ---  III -- SSO populations                                    -----------------------;
;-----------------------------------------------------------------------------------------------;
;-----------------------------------------------------------------------------------------------;
 ;--III.1-- Read SSO populations  --------------------------------------------------------------;
  readcol, dirEuclid+'pop', delimiter=',', $
           format='(A,I1,F4.1,F4.1,F4.1,F4.1,F3.1,F3,A,A)', $
           popName, popMain, popHmin, popHmax, xMin, xMax, yMin, yMax, popC, poS, /Silent
  popName = strTrim(popName,2)
  nbPop=n_elements(popName)
;  nbPop=1



  for kPop=kPopStart, nbPop-1 do $
    if popMain[kPop] eq 1 then begin

     print, ''
     print, popName[kPop]
     print, '   '+dirSFD+ popName[kPop]+'.csv'

    ;-- H2V -----------------------------------------------------------
     H2V = ssoHtoV(pop=popName[kpop], elong=91.5, stepDist=0.25)

    ;-- Cumulative SFD
     readcol, dirSFD+popName[kPop]+'.csv', delimiter=',', H, $
              cumL, cumM, cumH, $
              format='(F,F,F,F)', /Silent

    ;-- SFD ------------------------------------------------------------
     sfdL = cumL
     sfdM = cumM
     sfdH = cumH

     nbH=n_elements( sfdM )
     sfdL[0] = 0
     sfdM[0] = 0
     sfdH[0] = 0
     for kH=1, nbH-1 do begin
       sfdL[kH] = cumL[kH]-cumL[kH-1]
       sfdM[kH] = cumM[kH]-cumM[kH-1]
       sfdH[kH] = cumH[kH]-cumH[kH-1]
     endfor

     V=H
     PDF = V*0.
     nBin = n_elements(V)

     for kH=0, nbH-1 do begin

       y = sfdL[kH]*H2V.pdf
       x = H[kH] + H2V.mag

       for kB=0, nBin-2 do begin
         pp = where( x ge V[kB] and x lt V[kB+1], nbPP )
         if pp[0] ne -1 then begin
;           print, kH, kB, nbPP
           PDF[kB] += total(y[pp])
         endif
       endfor
       
     endfor
     PDF/=total(PDF)

     forprint, V, PDF, format='(F5.2,",",E8.2)', textout=dirSFD+popName[kPop]+'-V.csv', $
               comment='H, Probability', /Silent



  endif
        


end