dir='/home/bcarry/data/binary/87_Sylvia/satellite/512x512/'
  nights=['20181015_1',$
          '20181015_2',$
          '20181019',$
          '20181112_1',$
          '20181112_2',$
          '20181125',$
          '20181126',$
          '20181129']
  nbNight=n_elements(nights)

  doZOOM =01
  doMERGE=01



  edgeSampling=200
  
  dirZOOM   = '/data/binary/87_Sylvia/satellite/512x512/zoom/'
  dirEDGES  = '/data/binary/87_Sylvia/satellite/512x512/edges/'
  file_mkdir, dirZOOM
  file_mkdir, dirEDGES

  
;##########################################################################################;
;##########################################################################################;
;######   EXTRACT ZOOMS AROUND THE SATELLITE     ##########################################;
;##########################################################################################;
;##########################################################################################;
  if doZOOM eq 1 then begin
    dimMZ=20
    dimWX=60
    edgeDefault={LoG: 1., $ ;tbd
                 bin: 1 ,$
                 sampling: 100 }
    if not keyword_set(edgeINFO) then edgeINFO=edgeDefault else edgeINFO=updateStructure(edgeDefault,edgeINFO)
  
    
    for kN=0, nbNight-1 do begin
  
       print, kN+1, nbNight, nights[kN], format='(I2,"/",I2,"  -  ",A-22)'
  
       spawn, 'ls -1 '+dir+nights[kN]+'/*cam1*fits', list
       nbFile=n_elements(list)
  ;     print, list
  
       x=fltarr(nbFile)
       y=fltarr(nbFile)
       
       for kF=0, nbFile-1 do begin
         print, list[kF]
  
         split = strSplit(list[kF],'/')
         fileName  = strMid(list[kF],split[-1],70)
         
         im=readfits(list[kF],h,/silent)
         xy = grabXY( im, subframe=dimMZ, /peak2D, title='Romulus Image', win=768 )
         x[kF]=xy[0]
         y[kF]=xy[1]
         print, xy

         iso=sxpar(h,'DATE-OBS')
         print, dir+nights[kN]+'/'+fileName+'-eph.dat'
         print, iso
         eph=sspEphemcc('87/1',iso,nbd=1,dump=dir+nights[kN]+'/'+fileName+'-eph.dat')

         stop

         
         zoom=cropFrame(im,xy,dimWX)
         writefits,  dirZOOM+fileName, zoom, h
       endfor
  
    endfor


  endif


;
;  commande du soft de laurent jorda:
;        for k in $(ls *fits); do ./img.calcon1.sh $k 0.5 ${k%%.fits};done
;
;
  

;##########################################################################################;
;##########################################################################################;
;######   MERGE CONTOURS                         ##########################################;
;##########################################################################################;
;##########################################################################################;
  if doMERGE eq 1 then begin

     epList=['d87Sylvia_2018-10-15', $
             'd87Sylvia_2018-10-19', $
             'd87Sylvia_2018-11-12', $
             'd87Sylvia_2018-11-26', $
             'd87Sylvia_2018-11-29' ]
     nbEp=n_elements(epList)
     for kE=0, nbEp-1 do begin

;       ;-1- read dat write into fits
;        spawn, ' ls -1 '+dirZoom+epList[kE]+'*dat', listF
;        nbF=n_elements(listF)
;        for kF=0, nbF-1 do begin
;           
;           split = strSplit(listF[kF],'/')
;           nameDAT  = strMid(listF[kF],split[-1],70)
;           nameFITS = strMid(nameDAT,0,strLen(nameDAT)-3)+'fits'
;
;           readcol, dirZOOM+nameDAT, x, y, z, /Silent
;           h=headfits(dirZoom+nameFITS)
;
;           writefits, dirEDGES+nameFITS, transpose([[x],[y]]), h
;;           cgPlot, x, y, /ISO
;        endfor

       ;-2- Merge profiles
        spawn, ' ls -1 '+dirEDGES+epList[kE]+'*fits', listE
        print, listE
        edge = edgesMerging( listE, sampling=edgeSampling ) ;, figure=figINFO )
        edgesDisplay, edge, dump=dirEDGES+epList[kE]+'.eps'
        print, dirEDGES+epList[kE]+'.eps'
        stop
     endfor

     
  endif



end