dustem_write_all_web3p8.pro
1.34 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
PRO dustem_write_all_web3p8,st,dir_out
dir_out_dat=dir_out+'/data/'
dir_out_qabs=dir_out+'/oprop/'
dir_out_capa=dir_out+'/hcap/'
;== ISRF
file_out=dir_out_dat+'ISRF.DAT'
dustem_write_isrf_lv,file_out,st.isrf
;== LAMBDA
file_out=dir_out_qabs+'LAMBDA.DAT'
dustem_write_lambda,file_out,st.lambda
;== GRAIN
file_out=dir_out_dat+'GRAIN.DAT'
dustem_write_grain_web3p8,file_out,st.grain
;== MIX
file_out=dir_out_dat
dustem_write_mix,file_out,st.mix
;== SIZE
file_out=dir_out_dat
dustem_write_size_lv,file_out,st.size
;== QABS
file_out=dir_out_qabs
dustem_write_qabs_lv,file_out,st.qabs
;== CALOR
file_out=dir_out_capa
dustem_write_calor_lv,file_out,st.calor
;== POL
if !run_pol then begin
file_out=dir_out_dat+'ALIGN.DAT'
; Pol file
dustem_write_align,file_out,st.align
; Linear polarization
file_out=dir_out_qabs
IF not stregex(st.align.keywords, 'rrf', /bool) THEN for i_axis = 1, 2 do dustem_write_qpol,file_out,st,i_axis
; Circular polarization
if !run_circ then dustem_write_qcirc,file_out,st
; Anisotropic heating
IF (st.align.anisG0 > 0) THEN for i_axis = 1, 2 do dustem_write_qh,file_out,st,i_axis
; RRF Files
IF stregex(st.align.keywords, 'rrf', /bool) THEN for i_axis = 1, 2 do dustem_write_qpol_rrf,file_out,st,i_axis
endif
;DP : add TLS
file_out=dir_out_dat
;stop
dustem_write_tls,file_out,st.tls
file_out=dir_out_qabs
sortie:
END