From 4f9a7092c88ffe480d0b89757e07662880ea6fed Mon Sep 17 00:00:00 2001 From: Jean-Philippe Bernard Date: Mon, 19 Sep 2022 11:43:09 +0200 Subject: [PATCH] polished the code --- src/idl/dustem_fit_polarisation_example.pro | 8 ++++---- 1 file changed, 4 insertions(+), 4 deletions(-) diff --git a/src/idl/dustem_fit_polarisation_example.pro b/src/idl/dustem_fit_polarisation_example.pro index 2ce5189..c156a17 100644 --- a/src/idl/dustem_fit_polarisation_example.pro +++ b/src/idl/dustem_fit_polarisation_example.pro @@ -145,13 +145,13 @@ pd = [ $ '(*!dustem_params).grains(0).mdust_o_mh',$ ;PAH0 mass fraction '(*!dustem_params).grains(1).mdust_o_mh',$ ;PAH0 mass fraction '(*!dustem_params).grains(2).mdust_o_mh', $ ;PAH1 mass fraction - ;'dustem_plugin_synchrotron_1', $ ;Synchrotron spectra index + 'dustem_plugin_modify_dust_pol_2', $ ;This will set the polarization angle 'dustem_plugin_synchrotron_2', $ ;Synchrotron amplitude at 10 mm 'dustem_plugin_synchrotron_3', $ ;Synchrotron polarization fraction 'dustem_plugin_synchrotron_4' $ ;Synchrotron polarization angle ] -p_truth=[1.,7.8000E-04,7.8000E-04,7.8000E-04, 1.e-2, 0.3, 45.] -iv=p_truth+[0.,1.e-4,1.e-4,1.e-4,0.,0.1,10] ;initial values shifted from solution +p_truth=[1.,7.8000E-04,7.8000E-04,7.8000E-04, 10.,1.e-2, 0.3, 45.] +iv=p_truth+[0.,1.e-4,1.e-4,1.e-4,-5.,0.,0.1,10] ;initial values shifted from solution ;Final parameters : 1.0000000 0.00075715537 0.00075280676 0.00092480417 0.010006239 0.29981309 45.000163 @@ -187,7 +187,7 @@ yrange=[1e-4,1.e3] loadct,13 ;!y.range=[1e-8,100] ;This is to ajust plot range from outside the routine t1=systime(0,/sec) -title='CACA' ;does not work if title is undefined .... +title='FIT' ;does not work if title is undefined .... 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) t2=systime(0,/sec) -- libgit2 0.21.2