# <codecell>
map(lambda x:x*x,[-1,1,24])
# <codecell>
reload(Sp)
if 'meas' in locals():
del meas
import gc; gc.collect()
# <codecell>
# first convert measurements to Sp class, with inherent parameters defined
meas = Sp.Sp(m)
meas.params()
# <markdowncell>
# Plot the parameters for the specified time
# <codecell>
fig2,ax2 = plt.subplots(5,3,sharex=True,figsize=(15,8))
ax2 = ax2.ravel()
for i in range(meas.npar-1):
ax2[i].plot(meas.utc,Sp.smooth(meas.par[:,i],3))
ax2[i].set_title('Parameter '+str(i))
ax2[i].grid()
ax2[i].set_xlim([17,19])
print 'sp (wp, wvl, z, re, ta)'
# create custom key for sorting via wavelength
iwvls = np.argsort(s.zenlambda)
s.wv = np.sort(s.zenlambda)
# <codecell>
if 'Sp' in locals():
reload(Sp)
if 'lut' in locals():
del lut
import gc; gc.collect()
# <codecell>
lut = Sp.Sp(s,irrad=True)
lut.params()
lut.param_hires()
# <codecell>
lut.sp_hires()
# <codecell>
print lut.tau.shape
print lut.ref.shape
print lut.sp.ndim
print lut.par.size
print lut.par.shape
print lut.ref
