Blame view

src/idl/dustem_write_mix.pro 1.11 KB
427f1205   Jean-Michel Glorian   version 4.2 merged
1
2
PRO dustem_write_mix,dir,st

68d2f391   Jean-Philippe Bernard   modified to cope ...
3
4
5
6
7
8
IF ptr_valid(st[0]) THEN BEGIN
  comments=(*st[0]).comments
  Ncomments=n_elements(comments)
ENDIF ELSE BEGIN
  Ncomments=0L
ENDELSE
427f1205   Jean-Michel Glorian   version 4.2 merged
9
  
68d2f391   Jean-Philippe Bernard   modified to cope ...
10
Nst=n_elements(st)
427f1205   Jean-Michel Glorian   version 4.2 merged
11
  
68d2f391   Jean-Philippe Bernard   modified to cope ...
12
13
14
15
16
17
18

frmts=[ $
         '(A)', $
         '(A)', $
         '(A)', $
         '(A)', $
         '(A)']
427f1205   Jean-Michel Glorian   version 4.2 merged
19
20

;Start loop over type of grain
68d2f391   Jean-Philippe Bernard   modified to cope ...
21
22
23
FOR i=0L,Nst-1 DO BEGIN
  ;Test if the pointer is pointing to something
  IF ptr_valid(st(i)) THEN BEGIN
e71a7ada   Jean-Philippe Bernard   modified to work ...
24
      frmts=replicate('(A)',n_elements((*st(i)).fmix))
68d2f391   Jean-Philippe Bernard   modified to cope ...
25
        ;Get filename
427f1205   Jean-Michel Glorian   version 4.2 merged
26
27
28
        ffile=(*st(i)).file
        fv=str_sep(ffile,'/')
        file=dir+fv(n_elements(fv)-1)
68d2f391   Jean-Philippe Bernard   modified to cope ...
29
30
        ;Open file
        message,'Will write MIX file '+file,/continue
427f1205   Jean-Michel Glorian   version 4.2 merged
31
        openw,unit,file,/get_lun
68d2f391   Jean-Philippe Bernard   modified to cope ...
32
33
34
35
36
37
38
        ;Print comments
        IF Ncomments NE 0 THEN BEGIN
          FOR ii=0,Ncomments-1 DO BEGIN
             printf,unit,comments[ii]
          ENDFOR
        ENDIF
        ;Print FMIX
427f1205   Jean-Michel Glorian   version 4.2 merged
39
        FOR k=0L,n_elements((*st(i)).fmix)-1 DO BEGIN
68d2f391   Jean-Philippe Bernard   modified to cope ...
40
            printf,unit,((*st(i)).fmix)(k),format=frmts[k]
427f1205   Jean-Michel Glorian   version 4.2 merged
41
42
43
        ENDFOR
        close,unit
        free_lun,unit
68d2f391   Jean-Philippe Bernard   modified to cope ...
44
45
46
47
   ENDIF ELSE BEGIN
      ;left blank intentionally
   ENDELSE
ENDFOR
427f1205   Jean-Michel Glorian   version 4.2 merged
48
49

END