;;; These variables have to be defined for each burst series. ;;; Cams(0) is used as 'master', and it is assumed that for each file ;;; in there, there is a pendant in the Cams(i) directory with the ;;; sequential number not further off than +/- 2 ;;; Id defines the wavelength of the camera ;;; Naming conventions: ;;; Target-ch-12345-yyyymmdd-hhmmss. ;;; e.g., AR8704-gb-0000-19991020-081730.fits ;;; AR8704-ca-2140-19991020-081730.fits ;;; AR8704-ha+0250-19991020-081730.fits base = 'AR0898' position = {lat: -6.6, lon: -16.4, mu: 0.946} fs = '-' Cams = ['C1', 'C2', 'C3', 'C4', 'C5'] Id = ['ca', 'gb', 'bc', 'rc', 'ha'] kvdl = [0, 0, 0, 0, 0] safe = [2, 2, 2, 3, 3] ;;; example for 2 Ca-positions with different extensions pext = [['', '', '', '', ''], $ ['1', '', '', '', ''], $ ['2', '', '', '', ''], $ ['3', '', '', '', ''], $ ['4', '', '', '', ''], $ ['5', '', '', '', ''], $ ['6', '', '', '', ''], $ ['7', '', '', '', ''], $ ['8', '', '', '', ''], $ ['9', '', '', '', '']] kvdl = replicate(0,5,10) lam = [['-0000', '-0000', '-0000', '-0000', '-0000'], $ ['-0500', '-0000', '-0000', '-0000', '-0000'], $ ['-1000', '-0000', '-0000', '-0000', '-0000'], $ ['-1500', '-0000', '-0000', '-0000', '-0000'], $ ['-2000', '-0000', '-0000', '-0000', '-0000'], $ ['-2500', '-0000', '-0000', '-0000', '-0000'], $ ['-3000', '-0000', '-0000', '-0000', '-0000'], $ ['-3500', '-0000', '-0000', '-0000', '-0000'], $ ['-4000', '-0000', '-0000', '-0000', '-0000'], $ ['-4500', '-0000', '-0000', '-0000', '-0000']] npmo = 10 pmod = 0 ;;; Fucking Cams are not linear! Try to compensate ;;; no measured values for this day - use those of July 13 gamm = [[.075, 313], [.0, 0], [.079, 212], [.0, 0], [.088, 641]] ;;; Decide if subimage tracking should be done notrack = 0 ;;; replace low-intensity parts (int lower than rli) with raw image part rli = 0.1 ;;; use drop-in-place for parts with rms lower than ldr ldr = 0.04 ;;; Setting this to 1 will cause the program to empty the original ;;; data file after successful reconstruction empty_file = 0 ;;; Filter the raw data against the electronic interference? filt = [1, 1, 1, 1, 1] ;filt = [0,0,0,0] conn_ord = 7 alwaysdrop = 0 ;;; lockfile for stopping computations stpf = 'STOP_COMPUTING' busy = ['|', '/', '-', '\'] nb = 0 sl = string(byte([27, 91, 68])) ;;; Old FITS format? oldfile = 0 ;;; Define and fill the control structure bname = 'Core-Test' dpos = strlen(bname)+8 filreg = 'RAW/'+Cams(0)+'/'+bname+'-*' files = findfile(filreg) fdat = long(strmid(files, dpos, 10)) num = n_elements(fdat) n_cams = n_elements(Cams) ntot = num*n_cams db = 0 ;;; bursts already done - for correct counting in logfile start = 0 wait_before_start = 1 restart_cond = 0