included an help section

Jean-Philippe Bernard
1 parent 9eb53891
Showing 2 changed files with 154 additions and 6 deletions Show diff stats
src/idl/dustem_init_plugins.pro
src/idl/dustemwrap_plot.pro
-PRO dustem_init_plugins, pd, fpd
+PRO dustem_init_plugins,pd,fpd,help=help
+
+;+
+; NAME:
+;    dustem_init_plugins
+; PURPOSE:
+;    initializes variable !dustem_plugin
+; CATEGORY:
+;    DustEMWrap
+; CALLING SEQUENCE:
+;    dustem_init_plugins,pd[,fpd]
+; INPUTS:
+;    pd : variable parameters description
+; OPTIONAL INPUT PARAMETERS:
+;    fpd : fixed parameters description
+; OUTPUTS:
+;    None
+; OPTIONAL OUTPUT PARAMETERS:
+;    None
+; ACCEPTED KEY-WORDS:
+;    help
+; COMMON BLOCKS:
+;    None
+; SIDE EFFECTS:
+;    !dustem_plugin is set
+; RESTRICTIONS:
+;    The DustEM fortran code must be installed
+;    The DustEMWrap IDL code must be installed
+; PROCEDURES AND SUBROUTINES USED:
+;
+; EXAMPLES
+;    
+; MODIFICATION HISTORY:
+;    Written by Ilyes Choubani
+;    Evolution details on the DustEMWrap gitlab.
+;-
+
+IF keyword_set(help) THEN BEGIN
+  doc_library,'dustem_init_plugins'
+  goto,the_end
+END
 plugin_names=['']
@@ -63,5 +103,6 @@ ENDIF
 defsysv, '!dustem_plugin', ptr_new(strct)  ; Need the use of the pointer because dustem_init_plugins is called in dustem_init (avoid conflicting structures error). 
+the_end:
-end
+END
-PRO dustemwrap_plot,p_dim,st,dustem_sed,SED_spec,dustem_qsed,Q_spec,dustem_used,U_spec,dustem_polsed,P_spec,dustem_polfrac,SP_spec,dustem_psi_em,PSI_spec,dustem_ext,EXT_spec,dustem_qext,QEXT_spec,dustem_uext,UEXT_spec,dustem_polext,POLEXT_spec,dustem_fpolext,SPEXT_spec,dustem_psi_ext,PSIEXT_spec,_extra=_extra 
-
+PRO dustemwrap_plot,p_dim,st, $
+                    dustem_sed,SED_spec, $
+                    dustem_qsed,Q_spec, $
+                    dustem_used,U_spec, $
+                    dustem_polsed,P_spec, $
+                    dustem_polfrac,SP_spec, $
+                    dustem_psi_em,PSI_spec, $
+                    dustem_ext,EXT_spec, $
+                    dustem_qext,QEXT_spec, $
+                    dustem_uext,UEXT_spec, $
+                    dustem_polext,POLEXT_spec, $
+                    dustem_fpolext,SPEXT_spec, $
+                    dustem_psi_ext,PSIEXT_spec, $
+                    _extra=_extra,help=help
+
+;+
+; NAME:
+;    dustemwrap_plot
+; PURPOSE:
+;    plots dustemwrap observed sed, predicted sed, predicted spectra for total Intensity, polarization in emission and extinction
+; CATEGORY:
+;    DustEMWrap
+; CALLING SEQUENCE:
+;    dustemwrap_plot,p_dim,st,dustem_sed,SED_spec,dustem_qsed,Q_spec,dustem_used,U_spec,dustem_polsed,P_spec,dustem_polfrac,SP_spec,dustem_psi_em,PSI_spec,dustem_ext,EXT_spec,dustem_qext,QEXT_spec,dustem_uext,UEXT_spec,dustem_polext,POLEXT_spec,dustem_fpolext,SPEXT_spec,dustem_psi_ext,PSIEXT_spec,_extra=_extra 
+; INPUTS:
+;    p_dim
+; OPTIONAL INPUT PARAMETERS:
+;    st
+;    dustem_sed
+;    dustem_qsed
+;    dustem_used
+;    dustem_polsed
+;    dustem_polfrac
+;    dustem_psi_em
+;    dustem_ext
+;    dustem_qext
+;    dustem_uext
+;    dustem_polext
+;    dustem_fpolext
+;    dustem_psi_ext
+; OUTPUTS:
+;    dustem_sed
+;    dustem_qsed
+;    dustem_used
+;    dustem_polsed
+;    dustem_polfrac
+;    dustem_psi_em
+;    dustem_ext
+;    dustem_qext
+;    dustem_uext
+;    dustem_polext
+;    dustem_fpolext
+;    dustem_psi_ext
+;    SED_spec
+;    Q_spec
+;    U_spec
+;    P_spec
+;    SP_spec
+;    PSI_spec
+;    EXT_spec
+;    QEXT_spec
+;    UEXT_spec
+;    POLEXT_spec
+;    SPEXT_spec
+;    PSIEXT_spec
+; OPTIONAL OUTPUT PARAMETERS:
+;    
+
+; ACCEPTED KEY-WORDS:
+;    help
+; COMMON BLOCKS:
+;    None
+; SIDE EFFECTS:
+;    Does a plot
+; RESTRICTIONS:
+;    The DustEM fortran code must be installed
+;    The DustEMWrap IDL code must be installed
+; PROCEDURES AND SUBROUTINES USED:
+;
+; EXAMPLES
+;    
+; MODIFICATION HISTORY:
+;    Written by Ilyes Choubani
+;    Evolution details on the DustEMWrap gitlab.
+;-
+
+IF keyword_set(help) THEN BEGIN
+  doc_library,'dustemwrap_plot'
+  goto,the_end
+END
 ;This is also necessary for the plotting of the results of the fit (Last iteration) 
 IF not keyword_set(st) THEN BEGIN
@@ -83,8 +171,25 @@ new_mouchard=mouchard[0]
 new_mouchard.iteration=!dustem_iter.act
 new_mouchard.chi2=(*!dustem_fit).chi2
 new_mouchard.rchi2=(*!dustem_fit).rchi2
-new_mouchard.recomputed_chi2=recomputed_chi2_sed
-new_mouchard.recomputed_rchi2=recomputed_rchi2_sed
+;new_mouchard.recomputed_chi2=recomputed_chi2_sed
+;new_mouchard.recomputed_rchi2=recomputed_rchi2_sed
+
+;JPB: need here to store values of current parameter values (including plugins) into !iteration_mouchard
+help,(*!dustem_fit).CURRENT_PARAM_VALUES   ;This is not currently filled up ...
+;Below will work only for non plugin parameters
+goto,skip1
+param_descs=*(*!dustem_fit).param_descs
+params_current=[0.D0]
+FOR i=0L,n_elements(param_descs)-1 DO BEGIN
+    str='params_current=[params_current,'+param_descs[i]+']'
+    print,str
+    toto=execute(str)
+ENDFOR
+params_current=params_current[1:*]
+print,params_current
+
+stop
+skip1:
 mouchard=[mouchard,new_mouchard]
 !iteration_mouchard=ptr_new(mouchard)
@@ -1066,4 +1171,6 @@ IF (tag_names(*!dustem_plugin))(0) NE &#39;NONE&#39; THEN BEGIN
 ENDIF
+the_end:
+
 END