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
  dustem_sed=dustem_compute_sed(p_min,st=st,cont=cont)
  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(p_min,waves=ww,spec=spec)
;       stop
       dustem_gb_plot_fit_sed,ww,spec,dustem_sed,_extra=extra
    END
  ENDCASE
ENDELSE

END