OV-GSO-DC / dustem-wrapper_idl

Blame view

src/idl/dustem_write_all_release.pro 1.62 KB
Edit Raw Normal View History
662437c0   Annie Hughes   Added 
1
2
 
PRO dustem_write_all_release,st,dir_out
b05b5475   Ilyes Choubani   adding test so th... 
3
4
5
6
7
 
if strmid(!dustem_dat,0,/reverse) ne '/' then !dustem_dat=!dustem_dat+'/' 

dir_out_dat=dir_out+'data/'
dir_out_qabs=dir_out+'oprop/'
dir_out_capa=dir_out+'hcap/'
662437c0   Annie Hughes   Added 
8
9
10
 

;== ISRF
file_out=dir_out_dat+'ISRF.DAT'
29421054   Ilyes Choubani   write_isrf_lv is ... 
11
 
dustem_write_isrf_release,file_out,st.isrf
662437c0   Annie Hughes   Added 
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
 

;== 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_release,file_out,st.grains

;== GAS
file_out=dir_out_dat+'GAS.DAT'

dustem_write_gas,file_out,st.GAS
29421054   Ilyes Choubani   write_isrf_lv is ... 
27
28
29
662437c0   Annie Hughes   Added 
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
 

;== SPIN
dustem_write_spin,dir_out_dat,st.spin

;stop
;== CHRG
dustem_write_chrg,dir_out_dat,st.chrg

;== MIX
dustem_write_mix,dir_out_dat,st.mix

;== SIZE
dustem_write_size_lv,dir_out_dat,st.size

;== QABS
dustem_write_qabs_lv,dir_out_qabs,st.qabs

;== CALOR
dustem_write_calor_lv,dir_out_capa,st.calor

;== POL
;dustem_write_pol,dir_out_dat,st.pol
;FOR i_axis = 1, 3 DO dustem_write_qabspol,dir_out_qabs,st,i_axis

;== POL
if !run_pol then begin
    ;stop
	file_out=dir_out_dat+'ALIGN.DAT'
    ;stop
	; 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

;== TLS
file_out=dir_out_dat
dustem_write_tls,dir_out_dat,st.tls

sortie:

END