Blame view

src/idl/dustem_write_mix.pro 1.06 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
24
FOR i=0L,Nst-1 DO BEGIN
  ;Test if the pointer is pointing to something
  IF ptr_valid(st(i)) THEN BEGIN
        ;Get filename
427f1205   Jean-Michel Glorian   version 4.2 merged
25
26
27
        ffile=(*st(i)).file
        fv=str_sep(ffile,'/')
        file=dir+fv(n_elements(fv)-1)
68d2f391   Jean-Philippe Bernard   modified to cope ...
28
29
        ;Open file
        message,'Will write MIX file '+file,/continue
427f1205   Jean-Michel Glorian   version 4.2 merged
30
        openw,unit,file,/get_lun
68d2f391   Jean-Philippe Bernard   modified to cope ...
31
32
33
34
35
36
37
        ;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
38
        FOR k=0L,n_elements((*st(i)).fmix)-1 DO BEGIN
68d2f391   Jean-Philippe Bernard   modified to cope ...
39
            printf,unit,((*st(i)).fmix)(k),format=frmts[k]
427f1205   Jean-Michel Glorian   version 4.2 merged
40
41
42
        ENDFOR
        close,unit
        free_lun,unit
68d2f391   Jean-Philippe Bernard   modified to cope ...
43
44
45
46
   ENDIF ELSE BEGIN
      ;left blank intentionally
   ENDELSE
ENDFOR
427f1205   Jean-Michel Glorian   version 4.2 merged
47
48

END