fixed final plot for wraptest mode

Annie Hughes
1 parent 2926480f
Showing 1 changed file with 28 additions and 21 deletions Show diff stats
src/idl/dustem_init.pro
@@ -11,20 +11,28 @@ PRO dustem_init,model=model $
 ;+
 ; NAME:
 ;   dustem_init
+;  
 ; PURPOSE:
 ;   Defines system variables to run the DUSTEM code
+;  
 ; CATEGORY:
 ;    DustEMWrap, Distributed, High-Level, Initialization
+;  
 ; CALLING SEQUENCE:
 ;   dustem_init,model=model,
+;
 ; INPUTS:
 ;    None
+;
 ; OPTIONAL INPUT PARAMETERS:
 ;    None
+;
 ; OUTPUTS:
 ;    None
+;
 ; OPTIONAL OUTPUT PARAMETERS:
-;    ;   st_model *** COMMENT AH -- what is this? END COMMENT AH ***
+;    ;   st_model
+;
 ; ACCEPTED KEY-WORDS:
 ;    model = specifies the interstellar dust mixture used by DustEM
 ;           'MC10' model from Compiegne et al 2010 (default)
@@ -46,8 +54,10 @@ PRO dustem_init,model=model $
 ;   polarization   = set this variables for runs taking polarization into account
 ;   noerase        : set this keyword to not erase the temporary memory used by dustemwrap
 ;   grain_keywords : Fortran grain keywords. Must be an array of Ngrains strings. Use '?' to keep default values.
+;
 ; COMMON BLOCKS:
 ;    None
+;
 ; SIDE EFFECTS:
 ;    1/ The following system variables are initialized to null values
 ;	or empty pointers:
@@ -63,13 +73,17 @@ PRO dustem_init,model=model $
 ;    2/ optionally runs the fortran code (if wraptest)
 ;       This step calls the Fortran code to run a standard model
 ;    3/ !dustem_filters is initialized calling dustem_filter_init.pro.
+;
 ; RESTRICTIONS:
-;    The dustem fortran code must be installed
-;    The dustem idl wrapper must be installed
+;    The DustEM fortran code must be installed
+;    The DustEMWrap idl code must be installed
+;
 ; PROCEDURES AND SUBROUTINES USED:
 ;    *** COMMENT AH --> is this really NONE? ****
+;
 ; EXAMPLES
 ;    dustem_init,model="MC10",/wrap,/plot_it
+;
 ; MODIFICATION HISTORY:
 ;    Written by JPB,NF,DP Jan-2007
 ;    Evolution details on the DustEMWrap gitlab.
@@ -321,21 +335,7 @@ defsysv,&#39;!dustem_filters&#39;,ptr_new()     ;filter information (filled by dustem_fi
 ;; ** Structure <1b95608>, 17 tags, length=8232, data length=8118, refs=1:
 ;;    IRAC            STRUCT    -> <Anonymous> Array[1]
 ;;    MIPS            STRUCT    -> <Anonymous> Array[1]
-;;    MSX             STRUCT    -> <Anonymous> Array[1]
-;;    IRAS            STRUCT    -> <Anonymous> Array[1]
-;;    DIRBE           STRUCT    -> <Anonymous> Array[1]
-;;    SPM             STRUCT    -> <Anonymous> Array[1]
-;;    ISOCAM          STRUCT    -> <Anonymous> Array[1]
-;;    ISOPHOTP        STRUCT    -> <Anonymous> Array[1]
-;;    ISOPHOTC        STRUCT    -> <Anonymous> Array[1]
-;;    EFIRAS          STRUCT    -> <Anonymous> Array[1]
-;;    ARCHEOPS        STRUCT    -> <Anonymous> Array[1]
-;;    HFI             STRUCT    -> <Anonymous> Array[1]
-;;    LFI             STRUCT    -> <Anonymous> Array[1]
-;;    WMAP            STRUCT    -> <Anonymous> Array[1]
-;;    SPIRE           STRUCT    -> <Anonymous> Array[1]
-;;    PACS            STRUCT    -> <Anonymous> Array[1]
-;;    AKARI           STRUCT    -> <Anonymous> Array[1]
+;; ... etc
  
 defsysv,'!dustem_verbose',0                                  ;1=verbose
 defsysv,'!dustem_show_plot',1           ;0=show no plot
@@ -534,9 +534,16 @@ IF keyword_set(wraptest) THEN BEGIN
   dustem_write_all,st_model,dir_out
   st=dustem_run()
   IF keyword_set(plot_it) THEN BEGIN
-    loadct,13
-    dustem_plot_nuinu_em,st,yr=[1e-14,1.e-6],/ysty,xr=[1,5e3],/xsty
-  ENDIF
+     window,0 & loadct,13
+     yr=[1e-14,1.e-6]
+     xr=[1,5e3]
+     tit='dustem_init wrapper test '+!dustem_model
+     xtit=textoidl('\lambda (\mum)')
+     ytit=textoidl('\nuI_\nu^{em} (W/m^2/sr for N_H=1.e20 H/cm^2)')
+     yr=[1e-11,6.e-7]
+     xr=[1,5e3]
+     dustem_plot_nuinu_em,st,yr=yr,/ysty,xr=xr,/xsty,/xlog,/ylog,title=tit,xtit=xtit,ytit=ytit
+ ENDIF
 ENDIF
  
 ;Initialize filters