Commit 5c830aef577158b02194cd81a9d0895f804438e6
1 parent
63a74ca2
Exists in
master
removed modification by Ilyes which made the code to fail ...
Showing
1 changed file
with
4 additions
and
7 deletions
Show diff stats
src/idl/dustem_fit_sed_readme.pro
@@ -80,12 +80,10 @@ IF keyword_set(png) THEN BEGIN | @@ -80,12 +80,10 @@ IF keyword_set(png) THEN BEGIN | ||
80 | ENDIF | 80 | ENDIF |
81 | 81 | ||
82 | polarization=0. ;default is no polarization in models | 82 | polarization=0. ;default is no polarization in models |
83 | +;dustem_init,model=use_model,pol=polarization | ||
83 | 84 | ||
84 | ;COMMENT: for some reson polarization=1 doesn't work for ALL of the models except for the two last A and B models | 85 | ;COMMENT: for some reson polarization=1 doesn't work for ALL of the models except for the two last A and B models |
85 | 86 | ||
86 | -dustem_init,model=use_model,pol=polarization | ||
87 | -!dustem_nocatch=1 | ||
88 | - | ||
89 | ;=== Set which parameters you want to fit | 87 | ;=== Set which parameters you want to fit |
90 | CASE use_model OF | 88 | CASE use_model OF |
91 | 'DBP90':BEGIN | 89 | 'DBP90':BEGIN |
@@ -349,6 +347,9 @@ CASE use_model OF | @@ -349,6 +347,9 @@ CASE use_model OF | ||
349 | 347 | ||
350 | ENDCASE | 348 | ENDCASE |
351 | 349 | ||
350 | +dustem_init,model=use_model,pol=polarization | ||
351 | +!dustem_nocatch=1 | ||
352 | + | ||
352 | ;dustem_init,model=use_model,pol=polarization | 353 | ;dustem_init,model=use_model,pol=polarization |
353 | 354 | ||
354 | !dustem_verbose=1 | 355 | !dustem_verbose=1 |
@@ -362,10 +363,6 @@ file=dir+'Gal_composite_spectrum.xcat' | @@ -362,10 +363,6 @@ file=dir+'Gal_composite_spectrum.xcat' | ||
362 | ;file=dir+'Gal_composite_spectrum.xcat.old' | 363 | ;file=dir+'Gal_composite_spectrum.xcat.old' |
363 | 364 | ||
364 | spec=read_xcat(file,/silent) | 365 | spec=read_xcat(file,/silent) |
365 | -stop | ||
366 | -;spec.stokesQ=0. | ||
367 | -;spec.stokesU=0. | ||
368 | -;spec.stokesU=0. | ||
369 | 366 | ||
370 | ;ind=where(spec.error EQ 0.,count) | 367 | ;ind=where(spec.error EQ 0.,count) |
371 | ;IF count NE 0 THEN spec(ind).error=0.2*spec(ind).spec | 368 | ;IF count NE 0 THEN spec(ind).error=0.2*spec(ind).spec |