dustem_make_grid_description_table.pro
1.19 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
PRO dustem_make_grid_description_table,table_name,nolatex=nolatex
;dir=!phangs_data_dir+'/ISRF/GRIDS/'
;isrf_class_str='_isrfclass'+strtrim(isrf_class,2)
;model='DBP90'
;table_name=dir+model+'_MuseISRF_JWST_G0_YPAH_YVSG_4Phangs_noionis'+isrf_class_str+'.fits'
;dustem_make_grid_description_table,
;table_name='/Volumes/PILOT_FLIGHT1/PHANGS/ISRF/GRIDS/DL07_logn_MuseISRF_JWST_adaptedG0_Ypah1added_Yvsgadded_4Phangs_noionis_F0upext_isrfclass8.fits'
;dustem_make_grid_description_table,table_name
dustem_define_grid,table_name,filters=filters,status=status
Npar=(*!dustem_grid).NPARAMS
one_st={parname:'',Pmin:0.,Pmax:0.,pNvals:0L,pLog:0L}
st=replicate(one_st,Npar)
st.parname=(*!dustem_grid).PARAMS
st.pmin=(*!dustem_grid).PMIN_VALUES
st.pmax=(*!dustem_grid).PMAX_VALUES
st.pNvals=(*!dustem_grid).PNVALS
st.plog=(*!dustem_grid).PLOG
write_xcat,st,'/tmp/bidon.xcat'
IF not keyword_set(nolatex) THEN BEGIN
struct2latex_table,st,'/tmp/bidon.tex'
ENDIF
;stop
;===== check parameter values
;stop
;count parameters uniq values
n_diffs=lonarr(Npar)
FOR i=0L,Npar-1 DO BEGIN
vv=(*!dustem_grid).st_params.(3)
order=sort(vv)
vv=vv[order]
un=uniq(vv)
vv=vv[un]
n_diffs[i]=n_elements(vv)
ENDFOR
print,n_diffs
stop
END