Worked on reducing the Arecibo continuum data. There were some serious problems with the off positions previously, in particular generating weird slopes.
I suspect the problem was with these lines
sorted1 = sort(totals[0:-1:2])
sorted2 = sort(totals[1:-1:2])
since the indexing is wrong, but I'm not sure yet. I re-wrote the whole off-position measurement to average over all scans. Still in progress, but the fitting works now. Here's an example:
I've edited my idlang.vim to auto-identify files that start with a semicolon. Add Line 19: `` syn match idlangStatement "^s*;s"`` Line 61/2 (allow spaces before ;): syn match idlangComment "\s*[\;].*$" contains=idlangTodo AG
`` map = readfits('MOSAIC.fits',hdr) crpix1 = sxpar(hdr,'CRPIX1') crpix2 = sxpar(hdr,'CRPIX2') crval1 = sxpar(hdr,'CRVAL1') crval2 = sxpar(hdr,'CRVAL2') cd1_1 = sxpar(hdr,'CD1_1') cd2_2 = sxpar(hdr,'CD2_2') x = lindgen(n_e(map[,0])) y = lindgen(n_e(map[0,])) l = (x-crpix1)*cd1_1+crval1-360 b = (y-crpix2)*cd2_2+crval2 imdisp,map,/axis,xrange=[max(l),min(l)],yrange=[min(b),max(b)],range=[-1,8]`` This code makes use of imdisp.pro.
Page 1 / 1