added dust polarization angle

Jean-Philippe Bernard
1 parent e0d00de0
Showing 1 changed file with 16 additions and 2 deletions Show diff stats
src/idl/dustem_make_polarization_sed_example.pro
@@ -126,6 +126,19 @@ pd = [ $
              ]                   
 p_truth=[1.,7.8000E-04,7.8000E-04,7.8000E-04,3., 1.e-2, 0.3, 45.]
  
+pd = [ $
+             '(*!dustem_params).G0', $                   ;G0
+             '(*!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_modify_dust_pol_2',  $     ;This will set the polarization angle
+             'dustem_plugin_synchrotron_1', $   ;Synchrotron spectra index
+             '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,10.,3., 1.e-2, 0.3, 45.]
+
  
 ;iv=p_truth+[0.,8.e-4,8.e-4,8.e-4]      ;shifted from solution
 iv=p_truth
@@ -193,9 +206,10 @@ message,&#39;Wrote &#39;+filename,/continue
  
 filename_final=!dustem_wrap_soft_dir+'/Data/EXAMPLE_OBSDATA/example_sed3.xcat'
 str='cp '+filename+' '+filename_final
-spawn,str
+message,'Do '+str+' to make change permanent',/continue
+;spawn,str
  
-message,'Wrote '+filename_final,/continue
+;message,'Wrote '+filename_final,/continue
  
 the_end: