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

 ;----What to do

 ;---------- Directories
  dirEuclid = '/observ/euclid/'
  dirEuclid = '/home/bcarry/work/data/euclid/'

  dirFilter = dirEuclid+'filters/'
  dirTaxo   = dirEuclid+'taxonomy/'
  dirStat = dirEuclid + 'stat/'



;--- survey limits
  limEcl = 15
  limGal = 20.
  vLimit = 24.  ;-- YJH -- 24.5 in V, but smearing...
  binH   = 0.25

 ;-geometry of obsevration
  elongation=90



;-----------------------------------------------------------------------------------------------;
;-----------------------------------------------------------------------------------------------;
;---  TAG ---   II -- SSO populations                                    -----------------------;
;-----------------------------------------------------------------------------------------------;
;-----------------------------------------------------------------------------------------------;
 ;--II.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)
  popPrint = where( popMain eq 1, nbPopPrint )
  print, 'Populations considered: ('+strtrim(string(nbPop),2)+'): ', popName


;-----------------------------------------------------------------------------------------------;
;-----------------------------------------------------------------------------------------------;
;---  TAG ---  III -- Compute LSST vs Euclid                             -----------------------;
;-----------------------------------------------------------------------------------------------;
;-----------------------------------------------------------------------------------------------;
  print, 'Population  H2V_L   H2V_E   diff  H_LSST  H_Euclid'
  for kPop=0, nbPop-1 do begin

    lsst1 = ssoHtoV( pop=popName[kPop], elong=180.) 
    lsst2 = ssoHtoV( pop=popName[kPop], elong=110.) 

    euclid = ssoHtoV( pop=popName[kPop], elong=90.) 

    diff = euclid - lsst1
 ;   print, popName[kPop], 'Euclid', 24.5 - dataBox(euclid), $
 ;          format='(A-10,1x,A-10,1x,5(F6.2,2x))'
 ;   print, popName[kPop], 'LSST', 24 - dataBox(lsst1), $
 ;          format='(A-10,1x,A-10,1x,5(F6.2,2x))'

   print, popName[kPop], $
          median(lsst1), median(euclid), median(lsst1)-median(euclid), $
          24 - median(lsst1), 24.5 - median(euclid), $
          format='(A-10,1x,5(F6.2,2x))'

;    print, popName[kPop], 24.5 - dataBox(euclid), $
;           format='(A-10,1x,5(F6.2,2x))'

  endfor


end