OV-GSO-DC / dustem-wrapper_idl

Commit fdf4c603391e26c5b476e633f96cce638b4d8cbd

Authored by Annie Hughes
1 parent f9089dc2
Exists in master

cosmetic changes to help

Showing 1 changed file with 22 additions and 8 deletions   Show diff stats
Inline Side-by-side
src/idl/dustem_write_fits_table.pro
  Diff comments   View file @fdf4c60
@@ -3,30 +3,43 @@ PRO dustem_write_fits_table,filename=filename,help=help @@ -3,30 +3,43 @@ PRO dustem_write_fits_table,filename=filename,help=help
3 ;+ 3 ;+
4 ; NAME: 4 ; NAME:
5 ; dustem_write_fits_table 5 ; dustem_write_fits_table
  6 +;
6 ; PURPOSE: 7 ; PURPOSE:
7 ; Write a FITS table with content of current DustEMWrap system variables 8 ; Write a FITS table with content of current DustEMWrap system variables
  9 +;
8 ; CATEGORY: 10 ; CATEGORY:
9 ; DustEMWrap, Distributed, High-Level, User Convenience 11 ; DustEMWrap, Distributed, High-Level, User Convenience
  12 +;
10 ; CALLING SEQUENCE: 13 ; CALLING SEQUENCE:
11 ; dustem_write_fits_table,file[,/help] 14 ; dustem_write_fits_table,file[,/help]
  15 +;
12 ; INPUTS: 16 ; INPUTS:
13 ; file = File name to be written (default='./dustemwrap_results.fits') 17 ; file = File name to be written (default='./dustemwrap_results.fits')
  18 +;
14 ; OPTIONAL INPUT PARAMETERS: 19 ; OPTIONAL INPUT PARAMETERS:
15 ; None 20 ; None
  21 +;
16 ; OUTPUTS: 22 ; OUTPUTS:
17 ; None 23 ; None
  24 +;
18 ; OPTIONAL OUTPUT PARAMETERS: 25 ; OPTIONAL OUTPUT PARAMETERS:
19 ; None 26 ; None
  27 +;
20 ; ACCEPTED KEY-WORDS: 28 ; ACCEPTED KEY-WORDS:
21 ; help = If set, print this help 29 ; help = If set, print this help
  30 +;
22 ; COMMON BLOCKS: 31 ; COMMON BLOCKS:
23 ; None 32 ; None
  33 +;
24 ; SIDE EFFECTS: 34 ; SIDE EFFECTS:
25 ; Output FITS file is written to disk. 35 ; Output FITS file is written to disk.
  36 +;
26 ; RESTRICTIONS: 37 ; RESTRICTIONS:
27 ; The DustEMWrap IDL code must be installed 38 ; The DustEMWrap IDL code must be installed
  39 +;
28 ; PROCEDURE: 40 ; PROCEDURE:
29 ; None 41 ; None
  42 +;
30 ; EXAMPLES 43 ; EXAMPLES
31 ; 44 ;
32 ; MODIFICATION HISTORY: 45 ; MODIFICATION HISTORY:
@@ -43,14 +56,15 @@ ENDIF @@ -43,14 +56,15 @@ ENDIF
43 56
44 ;===== run dustem to get predicted spectrum and sed 57 ;===== run dustem to get predicted spectrum and sed
45 dustem_predicted_sed=dustem_compute_sed(*(*!dustem_fit).current_param_values,out_st=dustem_spectra_st) 58 dustem_predicted_sed=dustem_compute_sed(*(*!dustem_fit).current_param_values,out_st=dustem_spectra_st)
46 -;stop 59 +
47 IF !run_pol THEN BEGIN 60 IF !run_pol THEN BEGIN
48 ;This is needed to compute Q,U of both SED and spectra 61 ;This is needed to compute Q,U of both SED and spectra
49 - dustem_predicted_polsed=dustem_compute_stokes(*(*!dustem_fit).current_param_values,bidon, $  
50 - dustem_predicted_Qsed,dustem_predicted_Used, $  
51 - dustem_predicted_Q_spec,dustem_predicted_U_spec,dustem_predicted_PSI_spec, $  
52 - dustem_psi_em,out_st=out_st)  
53 - ;Same in extinction, but does not work in polarization extinction was not fitted. 62 + dustem_predicted_polsed=dustem_compute_stokes(*(*!dustem_fit).current_param_values,bidon, $
  63 + dustem_predicted_Qsed,dustem_predicted_Used, $
  64 + dustem_predicted_Q_spec,dustem_predicted_U_spec,dustem_predicted_PSI_spec, $
  65 + dustem_psi_em,out_st=out_st)
  66 +
  67 +;Same in extinction, but does not work in polarization extinction was not fitted.
54 ; dustem_predicted_polext=dustem_compute_stokext(*(*!dustem_fit).current_param_values,bidon, $ 68 ; dustem_predicted_polext=dustem_compute_stokext(*(*!dustem_fit).current_param_values,bidon, $
55 ; dustem_predicted_Qext,dustem_predicted_Uext, $ 69 ; dustem_predicted_Qext,dustem_predicted_Uext, $
56 ; dustem_predicted_Qext_spec,dustem_predicted_Uext_spec,dustem_predicted_PSI_ext, $ 70 ; dustem_predicted_Qext_spec,dustem_predicted_Uext_spec,dustem_predicted_PSI_ext, $
@@ -80,7 +94,7 @@ str_predicted_EXT=replicate(one_str_predicted_SED,N_predicted_sed) @@ -80,7 +94,7 @@ str_predicted_EXT=replicate(one_str_predicted_SED,N_predicted_sed)
80 fully_undefined_sed=str_predicted_SED 94 fully_undefined_sed=str_predicted_SED
81 str_predicted_SED.FILTER=(*(*!dustem_data).sed).filt_names ;taken from input SED 95 str_predicted_SED.FILTER=(*(*!dustem_data).sed).filt_names ;taken from input SED
82 str_predicted_SED.Wavelength=(*(*!dustem_data).sed).wav ;taken from input SED 96 str_predicted_SED.Wavelength=(*(*!dustem_data).sed).wav ;taken from input SED
83 -str_predicted_SED.I=dustem_predicted_sed ;but where is Q,U ?? 97 +str_predicted_SED.I=dustem_predicted_sed
84 ;stop 98 ;stop
85 IF !run_pol THEN BEGIN 99 IF !run_pol THEN BEGIN
86 str_predicted_SED.Q=dustem_predicted_Qsed 100 str_predicted_SED.Q=dustem_predicted_Qsed
@@ -215,7 +229,7 @@ ENDFOR @@ -215,7 +229,7 @@ ENDFOR
215 sxaddpar,header_predicted_SED,'COMMENT','===================================================',after='PARF'+strtrim(i,2) 229 sxaddpar,header_predicted_SED,'COMMENT','===================================================',after='PARF'+strtrim(i,2)
216 230
217 ;===== complement extension header for predicted spectrum total 231 ;===== complement extension header for predicted spectrum total
218 -sxaddpar,header_predicted_spectrum_tot,'TITLE','OUTPUT BEST FIT TOTAL SPECTRUM BY DUSTEMWRAP' 232 +;sxaddpar,header_predicted_spectrum_tot,'TITLE','OUTPUT BEST FIT TOTAL SPECTRUM BY DUSTEMWRAP'
219 sxaddpar,header_predicted_spectra,'TITLE','OUTPUT BEST FIT PER-GRAIN AND TOTAL EMISSION SPECTRA BY DUSTEMWRAP' 233 sxaddpar,header_predicted_spectra,'TITLE','OUTPUT BEST FIT PER-GRAIN AND TOTAL EMISSION SPECTRA BY DUSTEMWRAP'
220 sxaddpar,header_predicted_extinctions,'TITLE','OUTPUT BEST FIT PER-GRAIN AND TOTAL EXTINCTION BY DUSTEMWRAP' 234 sxaddpar,header_predicted_extinctions,'TITLE','OUTPUT BEST FIT PER-GRAIN AND TOTAL EXTINCTION BY DUSTEMWRAP'
221 235