from dgfit.dustmodel import DustModel, WD01DustModel
from dgfit.dustgrains import DustGrains
from dgfit.plotting.plot_effsize import plot

compnames = ["Silicates-DL84", "Carbonaceous-LD01"]

DustM = DustModel(componentnames=compnames, path="../../dgfit/data/indiv_grain/", dustmodel=None, obsdata=None, every_nth=2)
model = WD01DustModel(dustmodel=DustM)

plot(model)