PRO dustem_run_plugins $
	;,p_dim $

                      ,param_descs $
                      ,param_values $
                      ,run_order $
                      ,ask_for_st=ask_for_st $
                      ,st=st $
                      ,N2be_run=N2be_run $
                      ,help=help

;+
; NAME:
;    dustem_run_plugins

;  

; PURPOSE:

;    runs the plugins according to their scopes

;    a N_ELEMENTS(PLUGINS)X N_ELEMENTS(SCOPES) matrix is created to achieve this 

;

; CATEGORY:

;    DustEMWrap, Distributed, Mid-Level, Plugins

;

; CALLING SEQUENCE:

;    dustem_run_plugins,p_dim,param_descs,param_values,run_order,[,/ask_for_st][,st][N2be_run=][,/help] 

;

; INPUTS:

;    p_dim  = parameter values  (actually not used)

;    param_descs = parameter description vector
;    param_values = current parameter values

;    run_order    = run order of the plugins to be run

;

; OPTIONAL INPUT PARAMETERS:

;

; OUTPUTS:
;    None

;

; OPTIONAL OUTPUT PARAMETERS:
;    st         = Dustem output structure

;    N2be_run   = Number of plugins which have been run

; ACCEPTED KEY-WORDS:
;    help                  = if set, print this help

;

; COMMON BLOCKS:
;    None

;

; SIDE EFFECTS:
;    None

;

; RESTRICTIONS:
;    The DustEM fortran code must be installed

;    The DustEMWrap IDL code must be installed

;    CAUTION: p_dim is actually not used (!!)

; PROCEDURES AND SUBROUTINES USED:

;

; EXAMPLES

;    dustem_run_plugins, p_min, param_descs, param_values, scopes, st=st

;

; MODIFICATION HISTORY:
;    Written by IC 
;    Evolution details on the DustEMWrap gitlab.
;    See http://dustemwrap.irap.omp.eu/ for FAQ and help.  
;-

dew_prefix='dustem_plugin_' ; we assume that this prefix is at the start of all plugin routines

;=== if no plugins, do nothing

IF (*!dustem_plugin)[0].name EQ 'NONE' THEN BEGIN
	N2be_run=0L
	goto, the_end 
ENDIF

;LOCATING THE PLUGINS THAT SHOULD BE RUN

Nparams=n_elements(param_descs)
Nplugins=n_elements(*!dustem_plugin)
must_be_run=intarr(Nplugins)

FOR i=0L,Nplugins-1 DO BEGIN

		IF (*!dustem_plugin)[i].run_order EQ run_order THEN BEGIN
			must_be_run[i]=1
		ENDIF

ENDFOR

ind2be_run=where(must_be_run EQ 1,count)
N2be_run=count

;=== run plugins that need be run
FOR i=0L,Nplugins-1 DO BEGIN
	IF must_be_run[i] THEN BEGIN
		keys=[0L]
		vals=[0.d0]
		FOR j=0L,Nparams-1 DO BEGIN
			parameter_type=dustem_parameter_description2type(param_descs[j],string_name=string_name,plugin_name=plugin_name,key=key)
			IF parameter_type EQ 'PLUGIN' THEN BEGIN
				IF plugin_name EQ (*!dustem_plugin)[i].name THEN BEGIN
					ind=where(param_descs EQ 'dustem_'+string_name,count)
					IF count NE 0 THEN BEGIN
						keys=[keys,key]
						vals=[vals,param_values[ind[0]]]
					ENDIF
				ENDIF
			ENDIF
		ENDFOR
		keys=keys[1:*]
		vals=vals[1:*]

		IF keyword_set(ask_for_st) THEN BEGIN
			str='spectrum='+(*!dustem_plugin)[i].name+'(key=keys,val=vals,st=st)'
			message,'executing '+str,/info
			toto=execute(str)

		ENDIF ELSE BEGIN
			str='spectrum='+(*!dustem_plugin)[i].name+'(key=keys,val=vals)'
			message,'executing '+str,/info
			toto=execute(str)
		ENDELSE

		(*!dustem_plugin)[i].spec=ptr_new(spectrum)

	ENDIF
ENDFOR

;== JPB: This is now done in dustem_activate_plugins

;IF keyword_set(force_dustem_run) and N2be_run EQ 0 THEN BEGIN

;   st = dustem_run(p_dim,use_previous_fortran=use_previous_fortran)
;   !dustem_current = ptr_new(st)
;ENDIF

;=== JPB: This was initially above the previous test. I moved it below.
the_end:

END