dustem_compute_polext.pro 870 Bytes
FUNCTION dustem_compute_polext ,p_dim,sti,POLEXT_spec,_extra=extra,out_st=out_st,dustem_qext,dustem_uext,Q_ext,U_ext 

;CREATED AS A CONSEQUENCE OF THE INCLUSION OF THE PLUGINS


;written to homogonize DustEM

IF not keyword_set(sti) THEN BEGIN
  dustem_activate_plugins,p_dim/(*(*!dustem_fit).param_init_values)
  sti=dustem_run(p_dim)
ENDIF


POLEXT_spec = sti.polext.ext_tot * (*!dustem_HCD)/1.0e21  ; 



;ADDING PLUGIN TO SPECTRUM----------------
scopes=tag_names((*!dustem_plugin))
IF scopes[0] NE 'NONE' THEN BEGIN
;IF ptr_valid(!dustem_plugin) THEN BEGIN
  for i=0L,n_tags(*!dustem_plugin)-1 do begin
    if total(strsplit((*(*!dustem_plugin).(i).scope),'+',/extract) eq 'ADD_EXT') then POLEXT_spec+=(*(*!dustem_plugin).(i).spec)[*,0]
  endfor
ENDIF

dustem_polext = interpol(POLEXT_spec,sti.polext.wav,(*!dustem_data.polext).wav) 


return, dustem_polext



END