OV-GSO-DC / dustem-wrapper_idl

Commit 4f9a7092c88ffe480d0b89757e07662880ea6fed

Authored by Jean-Philippe Bernard
1 parent 5cb00808
Exists in master

polished the code

Showing 1 changed file with 4 additions and 4 deletions   Show diff stats
Inline Side-by-side
src/idl/dustem_fit_polarisation_example.pro
  Diff comments   View file @4f9a709
@@ -145,13 +145,13 @@ pd = [ $ @@ -145,13 +145,13 @@ pd = [ $
145 '(*!dustem_params).grains(0).mdust_o_mh',$ ;PAH0 mass fraction 145 '(*!dustem_params).grains(0).mdust_o_mh',$ ;PAH0 mass fraction
146 '(*!dustem_params).grains(1).mdust_o_mh',$ ;PAH0 mass fraction 146 '(*!dustem_params).grains(1).mdust_o_mh',$ ;PAH0 mass fraction
147 '(*!dustem_params).grains(2).mdust_o_mh', $ ;PAH1 mass fraction 147 '(*!dustem_params).grains(2).mdust_o_mh', $ ;PAH1 mass fraction
148 - ;'dustem_plugin_synchrotron_1', $ ;Synchrotron spectra index 148 + 'dustem_plugin_modify_dust_pol_2', $ ;This will set the polarization angle
149 'dustem_plugin_synchrotron_2', $ ;Synchrotron amplitude at 10 mm 149 'dustem_plugin_synchrotron_2', $ ;Synchrotron amplitude at 10 mm
150 'dustem_plugin_synchrotron_3', $ ;Synchrotron polarization fraction 150 'dustem_plugin_synchrotron_3', $ ;Synchrotron polarization fraction
151 'dustem_plugin_synchrotron_4' $ ;Synchrotron polarization angle 151 'dustem_plugin_synchrotron_4' $ ;Synchrotron polarization angle
152 ] 152 ]
153 -p_truth=[1.,7.8000E-04,7.8000E-04,7.8000E-04, 1.e-2, 0.3, 45.]  
154 -iv=p_truth+[0.,1.e-4,1.e-4,1.e-4,0.,0.1,10] ;initial values shifted from solution 153 +p_truth=[1.,7.8000E-04,7.8000E-04,7.8000E-04, 10.,1.e-2, 0.3, 45.]
  154 +iv=p_truth+[0.,1.e-4,1.e-4,1.e-4,-5.,0.,0.1,10] ;initial values shifted from solution
155 ;Final parameters : 1.0000000 0.00075715537 0.00075280676 0.00092480417 0.010006239 0.29981309 45.000163 155 ;Final parameters : 1.0000000 0.00075715537 0.00075280676 0.00092480417 0.010006239 0.29981309 45.000163
156 156
157 157
@@ -187,7 +187,7 @@ yrange=[1e-4,1.e3] @@ -187,7 +187,7 @@ yrange=[1e-4,1.e3]
187 loadct,13 187 loadct,13
188 ;!y.range=[1e-8,100] ;This is to ajust plot range from outside the routine 188 ;!y.range=[1e-8,100] ;This is to ajust plot range from outside the routine
189 t1=systime(0,/sec) 189 t1=systime(0,/sec)
190 -title='CACA' ;does not work if title is undefined .... 190 +title='FIT' ;does not work if title is undefined ....
191 res=dustem_mpfit_data(tol=tol,xtol=xtol,Nitermax=Nitermax,xrange=xrange,/xstyle,yrange=yrange,/ysty,/ylog,/xlog,xtit='wavelength [mic]',ytit='Brightness []',title=title) 191 res=dustem_mpfit_data(tol=tol,xtol=xtol,Nitermax=Nitermax,xrange=xrange,/xstyle,yrange=yrange,/ysty,/ylog,/xlog,xtit='wavelength [mic]',ytit='Brightness []',title=title)
192 t2=systime(0,/sec) 192 t2=systime(0,/sec)
193 193