OV-GSO-DC / dustem-wrapper_idl

Blame view

src/idl/dustem_str_inst.pro 1.43 KB
Edit Raw Normal View History
52a3cc37   Annie Hughes   First commit of I... 
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
 
FUNCTION DUSTEM_STR_INST, n1, n2=n2, n3=n3
; returns structure containing flux and CC in band 
; N1   (I): number of data points
; N2   (I): nr of grain types or different models
; N3   (I): nr of points in transmission (band flux data)

  if n_elements(n2) EQ 0 then n2 = 1

  if n_elements(n3) NE 0 then begin
     strct = { NAME :  strarr(n1),       $
               X :     dblarr(n1),       $  
               YD :    dblarr(n1),       $
               ERR :   dblarr(n1),       $
               YM  :   dblarr(n1,n2),    $
               FLX :   dblarr(n1,n2),    $
               CC :    dblarr(n1,n2),    $                   
               RR :    dblarr(n1),       $                   
               ISEL :  intarr(n1)+1,     $
               UNIT :  '',               $ 
               NPAR :  0,                $
               CHI2 :  0d                 }
     s1 = CREATE_STRUCT('TRANS', {X : dblarr(n1,n3), Y : dblarr(n1,n3) })
     strct = CREATE_STRUCT( strct, s1 )
  endif else if n_elements(BAND) EQ 0 then begin
        strct = { X :     dblarr(n1),           $  
                  YD :    dblarr(n1),           $
                  ERR :   dblarr(n1) ,          $
                  YM :    dblarr(n1,n2),        $
                  ISEL :  intarr(n1)+1,         $
                  UNIT :  '',                   $ 
                  NPAR :  0,                    $
                  CHI2 :  0d                      }
  endif

  RETURN, strct
END