diff --git a/src/idl/dustem_fit_intensity_example.pro b/src/idl/dustem_fit_intensity_example.pro
index d80991b..8ef7017 100644
--- a/src/idl/dustem_fit_intensity_example.pro
+++ b/src/idl/dustem_fit_intensity_example.pro
@@ -5,6 +5,7 @@ PRO dustem_fit_intensity_example,model=model $
                                 ,help=help $
                                 ,wait=wait $
                                 ,noobj=noobj $
+                                ,plot_all=plot_all $
                                 ,verbose=verbose
 
 ;+
@@ -57,8 +58,9 @@ PRO dustem_fit_intensity_example,model=model $
 ;    help      = if set, print this help
 ;    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
+;    verbose   = if set, subroutines will run in verbose mode
 ;    noobj     = if set, runs with no object graphics
+;    plot_all  = if set, will show all intermediate plotting outputs (default OFF)
 ;
 ; COMMON BLOCKS:
 ;    None
@@ -104,6 +106,8 @@ if exists ne 1 then $
 ;; ;===================================
 ;; ;=== PARSE OTHER INPUTS AND SET SOME BASIC RUN PARAMETERS
 ;; ;===================================
+do_all_plots=0
+if keyword_set(plot_all) then do_all_plots=1
 use_verbose=0
 if keyword_set(verbose) then use_verbose=1
 use_Nitermax=5        ; maximum number of iterations for the fit
@@ -268,7 +272,7 @@ fpd=[] & fiv=[]
 ;;===================================
 ;;=== INITIALISE DUSTEM
 ;;===================================
-dustem_init,model=use_model,polarization=use_polarization
+dustem_init,model=use_model,polarization=use_polarization,show=do_all_plots
 !dustem_nocatch=1
 !dustem_verbose=use_verbose
 IF keyword_set(noobj) THEN !dustem_noobj=1
--
libgit2 0.21.2