pro align_ba d_id = '1246871116' ;;; DARK f_id = '1246871071' ;;; FLAT r_id = '1246871035' ;;; Data par_ba = [-0.16, 0.997] p = assoc(1, intarr(1296, 1030), 2880) openr, 1, 'Al*'+d_id+'*john*' dr = float(p(16:*, 2:1025, 0)) FOR i=1, 49 DO dr += p(16:*, 2:1025, i) close, 1 dr /= 50. openr, 1, 'Al*'+d_id+'*gordon*' dh = float(p(16:*, 2:1025, 0)) FOR i=1, 49 DO dh += p(16:*, 2:1025, i) close, 1 dh /= 50. openr, 1, 'Al*'+f_id+'*john*' fr = float(p(16:*, 2:1025, 0)) FOR i=1, 49 DO fr += p(16:*, 2:1025, i) close, 1 fr /= 50. fr -= dr fr /= float(avg(fr)) openr, 1, 'Al*'+f_id+'*gordon*' fh = float(p(16:*, 2:1025, 0)) FOR i=1, 49 DO fh += p(16:*, 2:1025, i) close, 1 fh /= 50. fh -= dh fh /= float(avg(fh)) rms = fltarr(50) openr, 1, 'Al*'+r_id+'*john*' FOR i=0, 49 DO rms(i) = stdev(p(16:*, 2:1025, i)) ix = reverse(sort(rms)) pr = (float(p(16:*, 2:1025, ix(0)))-dr)/fr FOR i=1, 9 DO pr += (float(p(16:*, 2:1025, ix(i)))-dr)/fr close, 1 openr, 1, 'Al*'+r_id+'*gordon*' FOR i=0, 49 DO rms(i) = stdev(p(16:*, 2:1025, i)) ix = reverse(sort(rms)) ph = (float(p(16:*, 2:1025, ix(0)))-dh)/fh FOR i=1, 9 DO ph += (float(p(16:*, 2:1025, ix(i)))-dh)/fh close, 1 !X.Style = 1 !Y.Style = 1 pr = rotate(pr, 5) ph = rotate(ph, 5) redo = '' REPEAT BEGIN print, 'BaC-Ba:', (par_ba = rotscalenew(pic1=pr, pic2=ph, init=par_ba, /verb)) read, redo, Prompt = 'Redo with this value? (y/Enter)' ENDREP UNTIL (redo NE 'y') print print, 'BaC-Ba:', par_ba print, 'shift shc', shc(pr, rot(ph, par_ba(0), par_ba(1), cub=-.5), /filt) end