OV-GSO-DC / dustem-wrapper_idl

Blame view

src/idl/dustem_read_align.pro 3.82 KB
Edit Raw Normal View History
427f1205   Jean-Michel Glorian   version 4.2 merged 
1
2
 
FUNCTION dustem_read_align,dir_in,st_grains,silent=silent
66aff855   Jean-Philippe Bernard   modified to impro... 
3
4
5
6
7
8
 
;moved to dustem_init.pro
;defsysv, '!run_circ', 0.  
;defsysv, '!run_anis', 0.  
;defsysv, '!run_rrf', 0.    
;defsysv, '!run_univ', 0.    
;defsysv, '!run_lin', 0.
427f1205   Jean-Michel Glorian   version 4.2 merged 
9
10
11
12
13
14
15
16
17
18
19
 

file=dir_in+(*!dustem_inputs).align
openr,unit,file,/get_lun
Nlines=0L
str=''
WHILE not eof(unit) DO BEGIN
      readf,unit,str
      Nlines=Nlines+1
ENDWHILE
close,unit
free_lun,unit
b5ccb706   Jean-Philippe Bernard   improved to fit p... 
20
427f1205   Jean-Michel Glorian   version 4.2 merged 
21
22
23
24
25
26
27
28
29
30
31
 
;== now read the file
openr,unit,file,/get_lun
ncurrent=0L
REPEAT BEGIN
      readf,unit,str
      ncurrent=ncurrent+1
      first_char=strmid(str,0,1)
ENDREP UNTIL first_char NE '#'
key_str=str

IF stregex(key_str, 'lin', /bool) THEN !run_lin = 1
e69bf245   Jean-Philippe Bernard   improved 
32
33
 
IF stregex(key_str, 'univ', /bool) THEN BEGIN
  !run_univ = 1
8357b327   Jean-Philippe Bernard   removed a stop 
34
 
  ;stop
e69bf245   Jean-Philippe Bernard   improved 
35
 
ENDIF
427f1205   Jean-Michel Glorian   version 4.2 merged 
36
37
38
39
40
41
42
43
44
45
 
IF stregex(key_str, 'circ', /bool) THEN !run_circ = 1
IF stregex(key_str, 'anis', /bool) THEN !run_anis = 1
IF stregex(key_str, 'rrf', /bool) THEN !run_rrf = 1

readf,unit,str,format='(A100)'
xx=strsplit(str,' ',/regex,/extr)
anisG0 = float(xx(0))
ncurrent+=+1
nalig=nlines-ncurrent
b5ccb706   Jean-Philippe Bernard   improved to fit p... 
46
47
 
;stop
427f1205   Jean-Michel Glorian   version 4.2 merged 
48
 
one_st={aligned:0.,law:'',athresh:0.,plev:0.,pstiff:0.,TdsTg:0.,a0:0.,rvcut:0.}
b5ccb706   Jean-Philippe Bernard   improved to fit p... 
49
50
 
;Ngrains=st_grains.Ngrains
Ngrains=n_elements(st_grains)
427f1205   Jean-Michel Glorian   version 4.2 merged 
51
b5ccb706   Jean-Philippe Bernard   improved to fit p... 
52
53
54
 
st=replicate(one_st,Ngrains)

;st.aligned = stregex(st_grains.grains.type_keywords, 'pol', /bool)
36e8b879   Jean-Michel Glorian   update from deborah 
55
 
;count=0
b5ccb706   Jean-Philippe Bernard   improved to fit p... 
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
 

FOR i=0L,Ngrains-1 DO BEGIN
  ;stop
  is_aligned = stregex(st_grains(i).type_keywords, 'pol', /bool) 

  ;IF st(i).aligned then begin
  IF is_aligned THEN BEGIN
    ;stop
    readf,unit,str,format='(A100)'
    strv=str_sep(strcompress(strtrim(str,2)),' ')
    Nstrv=n_elements(strv)
  
    ii=0L
    st(i).law=strv(ii) & ii=ii+1

    ; IDG
    IF stregex(st(i).law, 'idg', /bool) THEN BEGIN
      st[i].TdsTg = strv[ii] & ii=ii+1
      st[i].a0    = strv[ii] & ii=ii+1
      st[i].rvcut = strv[ii] & ii=ii+1
    ENDIF

    ; RAT
    IF stregex(st(i).law, 'rat', /bool) THEN BEGIN
    ENDIF

    ; PARAMETRIC
    IF stregex(st(i).law, 'par', /bool) THEN BEGIN
      ; Grain radius Threshold for alignment, given in microns and converted in cm
      st[i].athresh = strv(ii) & ii=ii+1
      st[i].pstiff  = strv(ii) & ii=ii+1
      st[i].plev    = strv(ii) & ii=ii+1
    ENDIF
    st[i].aligned = stregex(st_grains(i).type_keywords, 'pol', /bool)
4ae78eb8   Jean-Philippe Bernard   removed stop 
90
 
    ;stop
b5ccb706   Jean-Philippe Bernard   improved to fit p... 
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
 

    IF !run_univ eq 1 THEN BEGIN
      st = replicate(st(i),Ngrains)
      ;st.aligned = stregex(st_grains.grains.type_keywords, 'pol', /bool) 
      st.aligned = stregex(st_grains(i).type_keywords, 'pol', /bool) 
      break
    ENDIF
  ENDIF
ENDFOR

; FOR i=0L,st_grains.Ngrains-1 DO BEGIN
;    IF st(i).aligned then begin
;         ; if count EQ 0. then begin   
;          readf,unit,str,format='(A100)'   
;       ;endif
36e8b879   Jean-Michel Glorian   update from deborah 
106
b5ccb706   Jean-Philippe Bernard   improved to fit p... 
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
 
;       	strv=str_sep(strcompress(strtrim(str,2)),' ')
;       	Nstrv=n_elements(strv)

; 	ii=0L
;       	st(i).law=strv(ii) & ii=ii+1 

; 	; IDG
;         IF stregex(st(i).law, 'idg', /bool) THEN begin
;            ii=ii+1 
;       		st(i).TdsTg = strv(ii) & ii=ii+1
;       		st(i).a0    = strv(ii) & ii=ii+1
;       		st(i).rvcut = strv(ii) & ii=ii+1
; 	endif

; 	; RAT
;       	IF stregex(st(i).law, 'rat', /bool) THEN begin
; 	endif

; 	; PARAMETRIC
;         IF stregex(st(i).law, 'par', /bool) THEN begin
;                                 ; Grain radius Threshold for
;                                 ; alignment, given in microns and
;                                 ; converted in cm
;       		st(i).athresh = strv(ii) & ii=ii+1
;       		st(i).pstiff  = strv(ii) & ii=ii+1
;                 st(i).plev    = strv(ii) & ii=ii+1
;                 ; count=1
; 	endif

;         IF !run_univ eq 1 then begin
; 		st = replicate(st(i),st_grains.Ngrains)
;                 st.aligned = stregex(st_grains.grains.type_keywords, 'pol', /bool)
;                 ;stop
; 		break
; 	endif
36e8b879   Jean-Michel Glorian   update from deborah 
142
143
144
b5ccb706   Jean-Philippe Bernard   improved to fit p... 
145
 
;       ENDIF
427f1205   Jean-Michel Glorian   version 4.2 merged 
146
b5ccb706   Jean-Philippe Bernard   improved to fit p... 
147
 
; ENDFOR
427f1205   Jean-Michel Glorian   version 4.2 merged 
148
149
150
151
152
 

close,unit
free_lun,unit

full_st={keywords:key_str,anisG0:anisG0,gamma:0,grains:st}
36e8b879   Jean-Michel Glorian   update from deborah 
153
427f1205   Jean-Michel Glorian   version 4.2 merged 
154
155
156
 
RETURN,full_st

END