;pro demo_for_desaturation
;;; Demo for the desaturation procedure DESAT for AIA images

;;; Cutout data request line.
;ssw_cutout_service,'2011/09/06 22:00:00','2011/09/06 22:30:00',ref_helio='N15W18',fovx=350,fovy=350,email='torre@dima.unige.it',waves=[94,131,171,193,211,304,335],max_frames=1000,instrument='aia',aec=1
;ssw_cutout_service,'2017/09/10 15:45:00','2017/09/10 16:15:00',ref_helio='S09W91',fovx=350,fovy=350,email='guastavino@dima.unige.it',waves=[94,131,171,193],max_frames=1000,instrument='aia',aec=1

;
; NOTE: in the code images will be resized to 499 * 499 pixels FOV by ctrl_data function.
; So, a bigger FOV for the input data is needed.

;cd,'personal_path' ;

add_path, './', /EXPAND
RESOLVE_ALL, /CONTINUE_ON_ERROR 


tstart =   '25-sep-2011 03:28:55';'10-sep-2017 16:06:20'; ; start time selection
tend   =   '25-sep-2011 03:28:58';'10-sep-2017 16:06:22'; ; end time selection

;cd, current=current_dir
;path = current_dir+'\25sep2011'                    ; data storage folder path
path = './25sep2011'
;wav  = ['94','131','171','193','211','304','335']  ; wavelength to process
wav = ['193']

obj = obj_new('desat')

result = obj ->desaturation( wav , tstart , tend  , path , save_fts = 1 , aec = 1, loud = 1, lev=1.0)

end