This code snippet should allow the computation of the lightcurves using the simple, analytical model derived by Brian Metzger:

 import matplotlib
 import matplotlib.pyplot as plt
 matplotlib.use('TkAgg')
 %matplotlib inline
 from gwemlightcurves.sampler.model import Me2017_model
 t_Me, lbol_Me, mag_Me = Me2017_model(mej=1e-2,vej=0.1,beta=3,kappa_r=25)
 for i in range(1,8):
     plt.plot(t_Me, mag_Me[i])
 plt.xlim(1, 14)
 plt.ylim(-10, -15)
 plt.show();

Similarly you could also call other models, e.g., here the one derived by Dietrich and Ujevic in 2017:

 from gwemlightcurves.sampler.model import DiUj2017_model_ejecta
 import numpy as np
 t, lbol, mag = DiUj2017_model_ejecta(mej=1e-1,vej=0.1,th=0.3,ph=90)
 print(np.size(mag[1]))
 for i in range(1,8):
     plt.plot(t, mag[i])
 plt.xlim(4, 20)
 plt.ylim(-10, -18)
 plt.show();