PRO dustem_sed_plot,p_min,_EXTRA=extra,function_name=function_name

IF not keyword_set(function_name) THEN BEGIN
  ;Run dustem with as an interface to mpfitfun

;stop
  ;== JPB: st used to come out from here
  ;even when not set. Had to replavce by out_st ...
  dustem_sed=dustem_compute_sed(p_min,st=st,cont=cont,out_st=out_st)
  st=out_st

  dustem_plot_fit_sed,st,dustem_sed,cont,_extra=extra
ENDIF ELSE BEGIN
  CASE function_name OF
    'dustem_greybody_mpfit':BEGIN

       dustem_sed=dustem_compute_gb_sed_fast(p_min,waves=ww,spec=spec)

;       stop
       dustem_gb_plot_fit_sed,ww,spec,dustem_sed,_extra=extra
    END
  ENDCASE
ENDELSE

END