finished for tutorial

Annie Hughes
1 parent bb69363f
Showing 1 changed file with 25 additions and 39 deletions Show diff stats
src/idl/dustem_myisrf_example.pro
@@ -5,6 +5,7 @@ PRO dustem_myisrf_example,model=model $
                          ,postscript=postscript $
                          ,fits_save_and_restore=fits_save_and_restore $
                          ,wait=wait $
+                         ,noobj=noobj $
                          ,verbose=verbose $
                          ,help=help
  
@@ -64,6 +65,7 @@ PRO dustem_myisrf_example,model=model $
 ;    wait      = if set, wait this many seconds between each step of
 ;                the code (for illustration purposes)
 ;    verbose      = if set, subroutines will run in verbose mode
+;    noobj     = if set, runs with no object graphics
 ;
 ; COMMON BLOCKS:
 ;    None
@@ -125,53 +127,35 @@ 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]
+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 
+; To suppress the Mathis field altogether
 ;fpd=['(*!dustem_params).G0']
-;fiv=[1.e-10]
+;fiv=[1.e-12]
  
  
-;; ;=== AN EXAMPLE FOR DL07 (or DL01 or WD01_RV5P5B)
+;; ;=== AN EXAMPLE FOR DL07 
 ;; ;=== 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='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(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]
-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.,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.]
+;; pd = [ $
+;;      '(*!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]
+;; Npar=n_elements(pd)
+;; ulimed=replicate(0,Npar)
+;; llimed=replicate(1,Npar)
+;; llims=replicate(1.e-15,Npar)
+;; ; To suppress the Mathis field altogether
+;; ;fpd=['(*!dustem_params).gas.G0']
+;; ;fiv=[1.e-12]
  
  
 if keyword_set(wait) then begin
@@ -179,10 +163,10 @@ if keyword_set(wait) then begin
    wait,wait
 end
  
-
 dustem_init,model=use_model,polarization=use_polarization
 !dustem_nocatch=1
 !dustem_verbose=use_verbose
+IF keyword_set(noobj) THEN !dustem_noobj=1
  
 ;== Fill a structure with default inputs to the model. This
 ;== includes the default ISRF, which we want here just for plotting a
@@ -194,7 +178,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=1,fname=my_isrf_file,x=mywaves)
+   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([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)
@@ -221,8 +205,10 @@ al_legend,/top,/right,clear=0,box=0 $
           ,thick=2,charsize=1.3
  
 ;=== READ THE ISRF
-if keyword_set(isrf_file) then !dustem_isrf_file=ptr_new(isrf_file)
-
+if keyword_set(isrf_file) then begin
+   message,'Setting ISRF component from file: '+isrf_file,/info
+   !dustem_isrf_file=ptr_new(isrf_file)
+end
  
 ;=== READ EXAMPLE SED DATA: 
 dir=!dustem_wrap_soft_dir+'/Data/EXAMPLE_OBSDATA/'