In [ ]:

    

import numpy as np
import matplotlib.pyplot as plt

%matplotlib inline


    

In [ ]:

    

from newdust.graindist import composition


    

In [ ]:

    

ENERGY = np.logspace(-1, 1, 100)
LAMBDA = np.logspace(0.0, 2, 100)


    

In [ ]:

    

dtest = composition.CmDrude()


    

In [ ]:

    

ax = plt.subplot(111)
dtest.plot(ax, ENERGY)
plt.loglog()


    

In [ ]:

    

ax = plt.subplot(111)
dtest.plot(ax, LAMBDA, unit='angs')
plt.loglog()


    

In [ ]:

    

stest = composition.CmSilicate()


    

In [ ]:

    

ax = plt.subplot(111)
stest.plot(ax)
plt.loglog()


    

In [ ]:

    

ax = plt.subplot(111)
stest.plot(ax, lam=ENERGY)
plt.loglog()


    

In [ ]:

    

ax = plt.subplot(111)
stest.plot(ax, lam=LAMBDA, unit='angs')
plt.loglog()
plt.legend(loc='upper left', frameon=False)


    

In [ ]:

    

gtest = composition.CmGraphite()


    

In [ ]:

    

ax = plt.subplot(111)
gtest.plot(ax)
plt.loglog()
plt.xlim()


    

In [ ]:

    

ax = plt.subplot(111)
gtest.plot(ax, lam=LAMBDA, unit='angs')
plt.loglog()
plt.legend(loc='upper left', frameon=False)


    

In [ ]:

    

ax = plt.subplot(111)
gtest.plot(ax, lam=ENERGY, unit='kev')
plt.loglog()
plt.legend(loc='upper right', frameon=False)


    

In [ ]: