pro go_corr, t0, t1, hours_grid=hours_grid, xrt=xrt, level0=level0, $ dir_top=dir_top, search_array=search_array, percentd_min=percentd_min, $ cat_sot=cat_sot, files_sot=files_sot, refresh_cat=refresh_cat, $ every_nth=every_nth, n_max_files=n_max_files, $ force_invert=force_invert, nsig_in=nsig_in, interactive=interactive, $ do_raw_cube=do_raw_cube, do_cubes=do_cubes, do_reg=do_reg, wcs_ref=wcs_ref, $ do_xstep=do_xstep, do_movie_map=do_movie_map, $ do_write=do_write, do_test=do_test, $ verbose=verbose, no_log_time=no_log_time, $ xdisp_ref=xdisp_ref, ydisp_ref=ydisp_ref, $ lil_reg_cube=lil_reg_cube, do_debug=do_debug, do_composite=do_composite, $ add_waves=add_waves, add_map_waves=add_map_waves, $ _extra=_extra ;+ ; name: go_corr ; call: go_corr, t0, t1 ; - if not exist(do_write) then do_write = 1 if not exist(hrs_forward) then hrs_forward = 6 if not exist(minutes_backoff) then minutes_backoff = 0 if not exist(level0) and not keyword_set(quicklook) then level0 = 1 if not exist(mode_select) then mode_select = 'fg' if not exist(hours_grid) then hours_grid = 1.0 ; 0.5 if not exist(dir_top) then dir_top = get_logenv('SOT_CORR_DB') file_log_exec_times = concat_dir(dir_top, 'sot_l1_corr_exec_times.log') if ( (not exist(search_array)) and (not keyword_set(interactive)) ) then $ ; search_array = ['naxis1=1024 && naxis2= 512', 'wave=Ca*II*H*line'] ; search_array = ['naxis1=1024 && naxis2=1024', 'wave=Ca*II*H*line'] search_array = ['naxis1=1024 && naxis2= 512 || naxis1=512 && naxis2=512 || naxis1=864 && naxis2=512', $ 'wave=Ca*II*H*line'] if not exist(dir_top) then begin if keyword_set(do_test) then $ dir_top = '$SSWDB/hinode/gen/hinode_sdo_alignment/sot_aia_align_test' else $ dir_top = '$SSWDB/hinode/gen/hinode_sdo_alignment/sot_aia_align' endif if not exist(percentd_min) then percentd_min = 98 if not exist(peak_str_min) then peak_str_min = -1 use_shimizu = 0 ; TODO - Fix the cat_sot refresh logic: delvarx, cat_sot delvarx, xdisp_rot, ydisp_rot, xdisp_ref_arr, ydisp_ref_arr, $ indexp_arr, lil_compare_cube delvarx, lil_reg_cube, cutoutp_add_wave_cube, $ n_max_files, index_str_arr, wcs_ref if not exist(cat_sot) or keyword_set(refresh_cat) then $ sot_cat, t0, t1, cat_sot, files_sot, xrt=keyword_set(xrt), search_array=search_array, $ level0=level0, quicklook=quicklook, count=count, tcount=tcount, refresh=refresh_cat help, cat_sot, files_sot prstr, all_vals(cat_sot.wave) if keyword_set(do_debug) then $ STOP, 'Stopping after call to sot_cat.' if count gt 0 then begin if keyword_set(interactive) then begin ss_mode = sot_umodes(incat=cat_sot, /int) files_sot_select = files_sot[ss_mode] ; [last_nelem(ss_mode,200)] endif else begin files_sot_select = files_sot endelse n_files_select = n_elements(files_sot_select) ; Optionally only process every 'nth' file: if keyword_set(do_select_sub_sampling) then begin every_nth = '' print, "Enter sub_sample frequency ('every_nth' - default is 1)", every_nth read, every_nth if every_nth eq '' then every_nth = 1 else every_nth = fix(every_nth) endif if keyword_set(every_nth) then begin ss_nth = indgen(n_files_select/every_nth) * every_nth files_sot_select = files_sot_select[ss_nth] n_files_select = n_elements(files_sot_select) endif t_files = file2time(files_sot_select) t_sec_files = anytim(t_files) ss_sort = sort(t_sec_files) files_sot_select = files_sot_select[ss_sort] if not exist(n_max_files_in) then $ n_max_files = n_files_select else n_max_files = n_max_files_in if keyword_set(debug) then STOP mreadfits_header, files_sot_select, index_raw error_corr = 0 for i=0, n_files_select-1 do begin t0_exec_full_loop_sec = anytim(!stime) delvarx, index_sot, data_sot file_sot = reform(files_sot_select[i]) read_sot, file_sot, index_sot, data_sot ; If read_sot failed, try once more: if not exist(index_sot) then begin print, 'read_sot failed. Waiting 100 seconds and trying again.' wait, 100 read_sot, file_sot, index_sot, data_sot endif ; Define percentd for images: if tag_exist(index_sot, 'percentd') then percentd = index_sot.percentd else $ percentd = (n_elements(where(data_sot ne 0))/n_elements(data_sot))*100 if percentd ge percentd_min then begin ; TOCHECK - NOTA BENA: ASSUMING SIGN IS NEGATIVE HERE! if not exist(xdisp_ref_arr) then begin sot_offsets = get_shimizu(index_sot.date_obs) xdisp_ref_arr = sot_offsets[0] ydisp_ref_arr = sot_offsets[1] xdisp_rot = sot_offsets[0] ydisp_rot = sot_offsets[1] if keyword_set(do_debug) then STOP, 'GO_CORR: Stopping on request prior to updating XCEN/YCEN.' endif else begin if ~ finite(peak_str) then begin print, 'PEAK_STR is not finite. Excluding.' endif else begin if (peak_str ge peak_str_min) then begin xdisp_ref_arr = [xdisp_ref_arr, xdisp_rot] ydisp_ref_arr = [ydisp_ref_arr, ydisp_rot] endif else begin print, 'PEAK_STR = ' + strtrim(peak_str,2) print, 'PEAK_STR exceeds minimum value. Excluding.' endelse endelse endelse print, 'Updating level0 CRVAL1/2 and X/YCEN (initial guesses for sdo cutout location).' xcen_orig = index_sot.xcen ycen_orig = index_sot.ycen index_sot.xcen = index_sot.xcen + total(xdisp_ref_arr) index_sot.ycen = index_sot.ycen + total(ydisp_ref_arr) index_sot.crval1 = $ comp_fits_crval(index_sot.xcen, index_sot.cdelt1, index_sot.naxis1, index_sot.crpix1) index_sot.crval2 = $ comp_fits_crval(index_sot.ycen, index_sot.cdelt2, index_sot.naxis2, index_sot.crpix2) index_sot.crpix1 = $ comp_fits_crpix(index_sot.xcen, index_sot.cdelt1, index_sot.naxis1, index_sot.crval1) index_sot.crpix2 = $ comp_fits_crpix(index_sot.ycen, index_sot.cdelt2, index_sot.naxis2, index_sot.crval2) new_corr, index_sot, data_sot, indexp_sot, datap_sot, $ xcen_orig=xcen_orig, ycen_orig=ycen_orig, $ file_sot=file_sot, file_sdo=file_sdo, $ dt_max=dt_max, use_maps=use_map, $ xquery=xquery, $ xdisp_ref_arr=xdisp_ref_arr, ydisp_ref_arr=ydisp_ref_arr, $ use_shimizu=use_shimizu, $ do_prep_sot=do_prep_sot, $ norollcorrect=norollcorrect, nocrop2disk=nocrop2disk, $ xdisp_rot=xdisp_rot, ydisp_rot=ydisp_rot, $ peak_str=peak_str, $ last_good_xcen=last_good_xcen, last_good_ycen=last_good_ycen, $ do_img_plot=do_img_plot, do_hist_plot=do_hist_plot, $ do_plot_map=do_plot_map, do_blink=do_blink, $ error_corr=error_corr, verbose=verbose, _extra=_extra if ( (keyword_set(do_write)) and (error_corr ne 1) ) then begin if not exist(indexp_sot_arr) then indexp_sot_arr = indexp_sot else $ indexp_sot_arr = concat_struct(indexp_sot_arr, indexp_sot) endif endif print, ' Finished ' + strtrim(i,2) + ' of ' + strtrim(n_files_select,2) + $ ' for time interval from ' + anytim(t0, /yoh) + ' to ' + anytim(t1, /yoh) + ' .' t1_exec_full_loop_sec = anytim(!stime) full_loop_time_sec = t1_exec_full_loop_sec - t0_exec_full_loop_sec sbuff = 'Full loop time for this image = ' + strtrim(full_loop_time_sec,2) + ' seconds.' print, sbuff file_append, file_log_exec_times, sbuff ; Error handling and clean up: catch, Error_status catch, /cancel error_corr = 0 ; Reset ERROR_CORR delvarx, indexp_sot endfor if ( exist(indexp_sot_arr) and keyword_set(do_write) ) then begin write_genxcat, indexp_sot_arr, topdir=dir_top, /day_round, $ prefix='sot_corr_offset_', catname=catname catname = catname + '.genx' endif ; Optionally create a composite co-aligned movie of SOT and other waves (for now SDO): if ( exist(indexp_sot_arr) and keyword_set(do_composite) ) then begin nn_test_movie, files_genx=catname, add_map_waves=add_map_waves, $ do_sot_map_movie=do_sot_map_movie, do_sdo_map_movie=do_sdo_map_movie, $ do_interleave_movie=do_interleave_movie, _extra=_extra endif endif end