if not keyword_set(ssos) then ssos=ssoDiscovery_read(/dump) 
  save=ssos


  doCSD   = 0
  doHisto = 01
  doSFD   = 0

;  t1 = '2007-01-01T12:00:00'
;  t2 = '2018-01-01T12:00:00'
;  t1 = '1990-01-01T12:00:00'
;  t2 = '2010-01-01T12:00:00'
;  
;  JD1 = date_conv(t1,'JULIAN')
;  JD2 = date_conv(t2,'JULIAN')
;  nbDay = JD2 - JD1
;  nbPage= 79.
;  step = nbDay / nbPage

  t1 = '1805-01-01T12:00:00'
  t2 = '1990-01-01T12:00:00'
  t3 = '2018-01-01T12:00:00'
  
  JD1 = date_conv(t1,'JULIAN')
  JD2 = date_conv(t2,'JULIAN')
  JD3 = date_conv(t3,'JULIAN')

  nbDay1  = JD2 - JD1
  nbPage1 = 15.
  step1 = nbDay1 / nbPage1

  nbDay2  = JD3 - JD2
  nbPage2 = 79. - nbPage1
  step2 = nbDay2 / nbPage2

  jdArr = [JD1 + step1*findgen(nbPage1), $
           JD2 + step2*findgen(nbPage2)  ]
           

  for k=0, nbPage-1 do begin
     print, k, jdArr[k],date_conv(jdArr[k],'FITS'), format='(I2,2x,D16.8,2x,A-20)'
  endfor

;stop



  dirSFD='/home/bcarry/Documents/parutions/2017-HDR/misc/figures/mov-sfd/'
  dirHIS='/home/bcarry/Documents/parutions/2017-HDR/misc/figures/mov-histo/'
  rootCSD ='csd-'
  rootSFD ='sfd-'
  rootHIS ='hist-'
  sufPNG='.png'



 ;----------------------------------------------------------------------------------------------------
 ;----------------------------------------------------------------------------------------------------
 ;----------------------------------------------------------------------------------------------------
 ;----------------------------------------------------------------------------------------------------
 ;--I-- Cumulative distribution vs time  
  if doCSD eq 1 then begin
    kEp=0
    nameEp = rootCSD+string(kEp+1,format='(I02)')+sufPNG
    ssoDiscovery_plot, disco=ssos, /AbsoluteMagnitude, /Cumulative, $
                       ep=jdArr[kEp], dump=dirSFD+nameEp
  
    kEp=1
    nameEp = rootCSD+string(kEp+1,format='(I02)')+sufPNG
    ssoDiscovery_plot, disco=ssos, /AbsoluteMagnitude, /Cumulative, $
                       ep=jdArr[kEp-1:kEp], dump=dirSFD+nameEp, /Ancient
  
    kEp=2
    nameEp = rootCSD+string(kEp+1,format='(I02)')+sufPNG
    ssoDiscovery_plot, disco=ssos, /AbsoluteMagnitude, /Cumulative, $
                       ep=jdArr[kEp-2:kEp], dump=dirSFD+nameEp, /Ancient
  
    kEp=3
    nameEp = rootCSD+string(kEp+1,format='(I02)')+sufPNG
    ssoDiscovery_plot, disco=ssos, /AbsoluteMagnitude, /Cumulative, $
                       ep=jdArr[kEp-3:kEp], dump=dirSFD+nameEp, /Ancient
  
    kEp=4
    nameEp = rootCSD+string(kEp+1,format='(I02)')+sufPNG
    ssoDiscovery_plot, disco=ssos, /AbsoluteMagnitude, /Cumulative, $
                       ep=jdArr[kEp-4:kEp], dump=dirSFD+nameEp, /Ancient
  
    kEp=5
    nameEp = rootCSD+string(kEp+1,format='(I02)')+sufPNG
    ssoDiscovery_plot, disco=ssos, /AbsoluteMagnitude, /Cumulative, $
                       ep=jdArr[kEp-5:kEp], dump=dirSFD+nameEp, /Ancient
  
    for kEp=6, nbPage-1 do begin
       nameEp = rootCSD+string(kEp+1,format='(I02)')+sufPNG
       ssoDiscovery_plot, disco=ssos, /AbsoluteMagnitude, /Cumulative, $
                          ep=jdArr[kEp-6:kEp], dump=dirSFD+nameEp, /Ancient
    endfor

  endif


 ;----------------------------------------------------------------------------------------------------
 ;----------------------------------------------------------------------------------------------------
 ;----------------------------------------------------------------------------------------------------
 ;----------------------------------------------------------------------------------------------------
 ;--II-- 2d a vs e Histograms
  if doHisto eq 1 then begin

    for kEp=0, nbPage-1 do begin

      nameEp = rootHIS+string(kEp+1,format='(I02)')

      print, kEp, nameEp, date_conv(jdArr[kEp],'FITS'), $
              format='(I2,2x,A-20,2x,A-20)'

      timeFile = '/home/bcarry/Documents/parutions/2017-HDR/misc/figures/mov-histo/time.csv'
      forprint, 10, 0.55, strMid(date_conv(jdArr[kEp],'FITS'), 0, 10), $
                comment='a, e, year', textout=timeFile, $
                format='(F,",",F,",",A10)', /Silent


      spawn, "sed -e 's/ENTERJD/"+strTrim(string(jdArr[kEp],format='(D16.1)'),2)+"/' "+$
             dirHIS+"template.sh > "+dirHIS+"toto.sh"

      spawn, "sed -e 's/ENTERFILE/"+strTrim(nameEp,2)+"/' "+$
             dirHIS+"toto.sh > "+dirHIS+"run.sh"

      spawn, ".  run.sh"

    endfor

  endif


 ;----------------------------------------------------------------------------------------------------
 ;----------------------------------------------------------------------------------------------------
 ;----------------------------------------------------------------------------------------------------
 ;----------------------------------------------------------------------------------------------------
 ;--II--- Size distribution vs time  
  if doSFD eq 1 then begin

    kEp=0
    nameEp = rootSFD+string(kEp+1,format='(I02)')+sufPNG
    ssoDiscovery_plot, disco=ssos, /AbsoluteMagnitude, $
                       ep=jdArr[kEp], dump=dirSFD+nameEp
  
    kEp=1
    nameEp = rootSFD+string(kEp+1,format='(I02)')+sufPNG
    ssoDiscovery_plot, disco=ssos, /AbsoluteMagnitude, $
                       ep=jdArr[kEp-1:kEp], dump=dirSFD+nameEp, /Ancient
  
    kEp=2
    nameEp = rootSFD+string(kEp+1,format='(I02)')+sufPNG
    ssoDiscovery_plot, disco=ssos, /AbsoluteMagnitude, $
                       ep=jdArr[kEp-2:kEp], dump=dirSFD+nameEp, /Ancient
  
    kEp=3
    nameEp = rootSFD+string(kEp+1,format='(I02)')+sufPNG
    ssoDiscovery_plot, disco=ssos, /AbsoluteMagnitude, $
                       ep=jdArr[kEp-3:kEp], dump=dirSFD+nameEp, /Ancient
  
    kEp=4
    nameEp = rootSFD+string(kEp+1,format='(I02)')+sufPNG
    ssoDiscovery_plot, disco=ssos, /AbsoluteMagnitude, $
                       ep=jdArr[kEp-4:kEp], dump=dirSFD+nameEp, /Ancient
  
    kEp=5
    nameEp = rootSFD+string(kEp+1,format='(I02)')+sufPNG
    ssoDiscovery_plot, disco=ssos, /AbsoluteMagnitude, $
                       ep=jdArr[kEp-5:kEp], dump=dirSFD+nameEp, /Ancient
  
    for kEp=6, nbPage-1 do begin
       nameEp = rootSFD+string(kEp+1,format='(I02)')+sufPNG
       ssoDiscovery_plot, disco=ssos, /AbsoluteMagnitude, $
                          ep=jdArr[kEp-6:kEp], dump=dirSFD+nameEp, /Ancient
    endfor

  endif



end