code now conforms to user tutorials

Annie Hughes
1 parent 37926d83
Showing 1 changed file with 48 additions and 25 deletions Show diff stats
src/idl/dustem_myisrf_example.pro
@@ -95,7 +95,7 @@ END
 IF keyword_set(model) THEN BEGIN
   use_model=strupcase(model)
 ENDIF ELSE BEGIN
-  use_model='DL07'    ;Default is last dustem model
+  use_model='DBP90'    ;default dust model for this example
 ENDELSE
  
  
@@ -112,6 +112,29 @@ dustem_define_la_common
 ;=== Refer to the DustEM and DustEMWrap userguides for an explanation
 ;    of the different grain types
  
+;=== AN EXAMPLE FOR DBP90
+;=== Here we fit the PAH abundance of the model, and the amplitude
+;=== of a user-supplied ISRF. A synchrotron component is included
+;=== The other grain parameters are fixed to their default values.
+;=== The free parameters are all lower-bounded at zero.
+;=== Note that if '(*!dustem_params).G0' is not set explicitly as a fixed parameter, then
+;=== the ISRF used by DustEM will include a standard Mathis field with
+;=== amplitude, G0=1.
+use_model='DBP90' ; you should specify this above, or in the command line
+pd = [ $
+     '(*!dustem_params).grains(0).mdust_o_mh',$ ;PAH0 mass fraction
+     'dustem_plugin_modify_isrf_1', $           ; Amplitude of user-defined ISRF
+     'dustem_plugin_synchrotron_2']             ;Synchrotron amplitude at 10 mm
+iv =   [5.4e-4,0.1, 0.01]
+Npar=n_elements(pd)
+ulimed=replicate(0,Npar)
+llimed=replicate(1,Npar)
+llims=replicate(1.e-15,Npar)
+; To suppress the 
+;fpd=['(*!dustem_params).G0']
+;fiv=[1.e-10]
+
+
 ;; ;=== AN EXAMPLE FOR DL07 (or DL01 or WD01_RV5P5B)
 ;; ;=== Here we fit the PAH abundances of the model, and the amplitude
 ;; ;=== of a user-supplied ISRF. A synchrotron component is included
@@ -119,36 +142,36 @@ dustem_define_la_common
 ;; ;=== The free parameters are all lower-bounded at zero.
 ;; use_model='DL07' ; you should specify this above, or in the command line
 pd = [ $
-         '(*!dustem_params).grains(0).mdust_o_mh',$     ;PAH0 mass fraction
+     '(*!dustem_params).grains(0).mdust_o_mh',$     ;PAH0 mass fraction
      '(*!dustem_params).grains(1).mdust_o_mh',$         ;PAH1 mass fraction
      'dustem_plugin_modify_isrf_1', $                   ; Ampitude of user-defined ISRF
      'dustem_plugin_synchrotron_2']                     ;Synchrotron amplitude at 10 mm
-iv =   [5.4e-4, 5.4e-4,0.1, 0.01]
+iv =   [5.4e-4, 5.4e-4, 0.1, 0.01]
 Npar=n_elements(pd)
 ulimed=replicate(0,Npar)
 llimed=replicate(1,Npar)
 llims=replicate(1.e-15,Npar)
-;   fpd=['(*!dustem_params).G0', $
-;       '(*!dustem_params).gas.G0']
-;   fiv=[1.e-10,1.e-10]
-
-;; ;=== AN EXAMPLE FOR DBP90
-;; ;=== Here we fit the PAH abundances of the model, and the amplitude
-;; ;=== of a user-supplied ISRF. A synchrotron component is included
-;; ;=== The other grain parameters are fixed to their default values.
-;; ;=== The free parameters are all lower-bounded at zero.
-;; use_model='DBP90' ; you should specify this above, or in the command line
-;; pd = [ $
-;;      '(*!dustem_params).grains(0).mdust_o_mh',$ ;PAH0 mass fraction
-;;      'dustem_plugin_modify_isrf_1', $           ; Ampitude of user-defined ISRF
-;;      'dustem_plugin_synchrotron_2']             ;Synchrotron amplitude at 10 mm
-;; iv =   [5.4e-4,0.1, 0.01]
-;; Npar=n_elements(pd)
-;; ulimed=replicate(0,Npar)
-;; llimed=replicate(1,Npar)
-;; llims=replicate(1.e-15,Npar)
-;; ;   fpd=['(*!dustem_params).G0']
-;; ;   fiv=[1.e-10]
+fpd=['(*!dustem_params).G0']
+     '(*!dustem_params).gas.G0']
+   fiv=[1.,1.]
+
+;=== AN EXAMPLE FOR DBP90
+;=== Here we fit the PAH abundances of the model, and the amplitude
+;=== of a user-supplied ISRF. A synchrotron component is included
+;=== The other grain parameters are fixed to their default values.
+;=== The free parameters are all lower-bounded at zero.
+use_model='DBP90' ; you should specify this above, or in the command line
+pd = [ $
+     '(*!dustem_params).grains(0).mdust_o_mh',$ ;PAH0 mass fraction
+     'dustem_plugin_modify_isrf_1', $           ; Ampitude of user-defined ISRF
+     'dustem_plugin_synchrotron_2']             ;Synchrotron amplitude at 10 mm
+iv =   [5.4e-4,0.1, 0.01]
+Npar=n_elements(pd)
+ulimed=replicate(0,Npar)
+llimed=replicate(1,Npar)
+llims=replicate(1.e-15,Npar)
+fpd=['(*!dustem_params).G0']
+fiv=[1.]
  
  
 if keyword_set(wait) then begin
@@ -171,7 +194,7 @@ st_model=dustem_read_all(dir_in)
 if not keyword_set(isrf_file) then begin
    mywaves=[]
    my_isrf_file='./myisrf_3bb_8000K.dat'
-   myisrf=dustem_create_rfield([8000,8000,8000],wdil=[1.d-11,1.d-11,1.d-11],isrf=0,fname=my_isrf_file,x=mywaves)
+   myisrf=dustem_create_rfield([8000,8000,8000],wdil=[1.d-11,1.d-11,1.d-11],isrf=1,fname=my_isrf_file,x=mywaves)
 ;myisrf=dustem_create_rfield([3500,4000,5500,6500,8000],wdil=[1.d-13,1.d-13,1.d-13,1.d-13,1.d-13],isrf=1,fname='./myisrf_mathis_5bb.dat',x=mywaves)
 ;myisrf=dustem_create_rfield([3500,4000,5500,6500,8000],wdil=[1.d-13,1.d-13,1.d-13,1.d-13,1.d-13],isrf=2,fname='./myisrf_habing_5bb.dat',x=mywaves)
 ;myisrf=dustem_create_rfield([3500,4000,5500,6500,8000],wdil=[1.d-13,1.d-13,1.d-13,1.d-13,1.d-13],isrf=0,fname='./myisrf_5bb.dat',x=mywaves)