phangs_make_hi_images.pro
2.38 KB
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
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
PRO phangs_make_hi_images,source_name=source_name,save=save,show_images=show_images,nostop=nostop,help=help
;+
; NAME:
; phangs_make_hi_images
; PURPOSE:
; makes HI images usable by the Phangs ISRF project
; CATEGORY:
; Dustem Phangs
; CALLING SEQUENCE:
; phangs_make_hi_images[,source_name=][,/save][,/show_images][,/nostop]
; INPUTS:
; None
; OPTIONAL INPUT PARAMETERS:
; source_name : source name (default = 'ngc0628')
; OUTPUTS:
; None
; OPTIONAL OUTPUT PARAMETERS:
; None
; ACCEPTED KEY-WORDS:
; help = If set, print this help
; save = If set, save result
; show_images= if set, show images
; nostop = if set, does not stop
; COMMON BLOCKS:
; None
; SIDE EFFECTS:
; A file is written
; RESTRICTIONS:
; None
; PROCEDURE:
; input HI images are asumed to be in units of K*km/s
; EXAMPLES
; phangs_make_hi_images,source_name='ngc0628',/save,/show,/nostop
; MODIFICATION HISTORY:
; Written by J.-Ph. Bernard (2023)
; Evolution details on the DustEMWrap gitlab.
; See http://dustemwrap.irap.omp.eu/ for FAQ and help.
;-
IF keyword_set(help) THEN BEGIN
doc_library,'phangs_make_hi_images'
goto,the_end
ENDIF
win=0L
dustem_define_la_common
use_source_name='ngc0628'
IF keyword_set(source_name) THEN use_source_name=source_name
save_data_dir=!phangs_data_dir+'/ISRF/WORK/'
file=save_data_dir+use_source_name+'_ref_header.sav'
restore,file,/verb
;% RESTORE: Restored variable: HREF.
Nx=sxpar(href,'NAXIS1')
Ny=sxpar(href,'NAXIS2')
data_dir=!phangs_data_dir+'/phangs_drive/PHANGS_HI/'
file=data_dir+use_source_name+'_21cm_strictmask_mom0.fits'
d=readfits(file,h)
;stop
;sxaddpar,h,'CTYPE1','RA---TAN'
;sxaddpar,h,'CTYPE2','RA---TAN'
sxaddpar,h,'EQUINOX',2000.
ind=where(finite(d) NE 1,count)
IF count NE 0 THEN d[ind]=la_undef()
IF sxpar(h,'NAXIS1') NE Nx OR sxpar(h,'NAXIS2') NE Ny THEN WHI=project2(h,d,href,/silent) ELSE WHI=d
fact=1.83e18/1.e21
NHI=la_mul(WHI,fact) ;NH from HI in 1e21 H/cm2
tit=source_name+' '+'NHI [1e21 H/cm2]'
;stop
IF keyword_set(show_images) THEN BEGIN
window,win & win=win+1
image_cont20,NHI,href,/square,imrange=[-0.1,5],image_color_table='jpbloadct',/silent,tit=tit,axis_color_table=1
IF not keyword_set(nostop) THEN stop
ENDIF
IF keyword_set(save) THEN BEGIN
save_file=save_data_dir+use_source_name+'_HI_images.sav'
save,NHI,href,file=save_file
message,'Saved '+save_file,/continue
ENDIF
the_end:
END