forked from akleroy/phangs_imaging_scripts
-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathprocess_singledish.pro
47 lines (36 loc) · 1.76 KB
/
process_singledish.pro
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
tol = 1d-6
; in cycle 3
cube = readfits('ALMA_TP.M74.v0p2.fits', hdr)
blank = where(abs(cube - sxpar(hdr,'BLANK')) lt tol)
cube[blank] = !values.f_nan
writefits, 'ALMA_TP.M74.v0p2.blanked.fits', cube, hdr
cube = readfits('ALMA_TP.NGC_1087.CO21.v0p2.gildas.fits', hdr)
blank = where(abs(cube - sxpar(hdr,'BLANK')) lt tol)
cube[blank] = !values.f_nan
writefits, 'ALMA_TP.NGC_1087.CO21.v0p2.gildas.blanked.fits', cube, hdr
cube = readfits('ALMA_TP.NGC_1300.CO21.v0p2.image.VLSRK.gildas.fits', hdr)
blank = where(abs(cube - sxpar(hdr,'BLANK')) lt tol)
cube[blank] = !values.f_nan
writefits, 'ALMA_TP.NGC_1300.CO21.v0p2.image.VLSRK.gildas.blanked.fits', cube, hdr
cube = readfits('ALMA_TP.NGC_1385.CO21.v0p2.image.VLSRK.gildas.fits', hdr)
blank = where(abs(cube - sxpar(hdr,'BLANK')) lt tol)
cube[blank] = !values.f_nan
writefits, 'ALMA_TP.NGC_1385.CO21.v0p2.image.VLSRK.gildas.blanked.fits', cube, hdr
cube = readfits('ALMA_TP.NGC_1433.CO21.v0p2.gildas.NuclearFix.fits', hdr)
blank = where(abs(cube - sxpar(hdr,'BLANK')) lt tol)
cube[blank] = !values.f_nan
writefits, 'ALMA_TP.NGC_1433.CO21.v0p2.gildas.NuclearFix.blanked.fits', cube, hdr
cube = readfits('ALMA_TP.NGC_1566.CO21.v0p2.gildas.fits', hdr)
blank = where(abs(cube - sxpar(hdr,'BLANK')) lt tol)
cube[blank] = !values.f_nan
writefits, 'ALMA_TP.NGC_1566.CO21.v0p2.gildas.blanked.fits', cube, hdr
; in cycle 5
tol = 1d-6
cube = readfits('ALMA_TP.NGC4293.CO21.blank.image.VLSRK.fits', hdr)
blank = where(abs(cube - sxpar(hdr,'BLANK')) lt tol)
cube[blank] = !values.f_nan
writefits, 'ALMA_TP.NGC4293.CO21.processed.fits', cube, hdr
cube = readfits('ALMA_TP.IC1954.CO21.blank.image.VLSRK.fits', hdr)
blank = where(abs(cube - sxpar(hdr,'BLANK')) lt tol)
cube[blank] = !values.f_nan
writefits, 'ALMA_TP.IC1954.CO21.processed.fits', cube, hdr