# IPython log file

fn = '/mbari/LRAUV/tethys/missionlogs/2013/20130920_20131007/20130920T204656/201309202047_201309210606.mat'
import okeanidanalysis as oa
s = oa.logs.OkeanidLog(fn)
s.print_tree()
s.plot_timeseries('volume_scattering_470_nm')
s.plot_timeseries('volume_scattering_650_nm')
v,t = s.timeseries('volume_scattering_650_nm')
dv = np.diff(v)
plot(dv)
np.any(dv = 0)
(dv == 0).any()
np.find(dv == 0)
idx = np.where(dv == 0)[0]
plot(idx, dv[idx], 'r.')
figure)(
figure()
plot(idx)
clf()
hist(diff(idx))
hist(diff(idx), np.arange(20))
clf()
hist(diff(idx), np.arange(20))
clf()
plot(diff(idx))
fn
s.print_tree('CTD_NeilBrown')
s.plot_timeseries('CTD_NeilBrown/sea_water_salinity')
s.plot_timeseries('CTD_NeilBrown/sea_water_conductivity')
s.plot_timeseries('CTD_NeilBrown/sea_water_electrical_ conductivity')
s.plot_timeseries('CTD_NeilBrown/sea_water_electrical_conductivity')
clf()
s.plot_timeseries('CTD_NeilBrown/sea_water_electrical_conductivity')
figure()
s.plot_timeseries('CTD_NeilBrown/sea_water_salinity')
fig = plt.figure()
ax = fig.add_subplot(2,1,1)
s.plot_timeseries('CTD_NeilBrown/sea_water_electrical_conductivity', ax = ax)
s.plot_timeseries('CTD_NeilBrown/sea_water_electrical_conductivity', axes = ax)
ax2 = fig.add_subplot(2,1,2, sharex = ax)
s.plot_timeseries('CTD_NeilBrown/sea_water_salinity', axes=ax2)
draw()
ax1 = fig.add_subplot(5,1,1)
ax2 = fig.add_subplot(5,1,2, sharex = ax1)
ax3 = fig.add_subplot(5,1,3, sharex = ax1)
ax4 = fig.add_subplot(5,1,4, sharex = ax1)
ax5 = fig.add_subplot(5,1,5, sharex = ax1)
s.plot_timeseries('CTD_NeilBrown/sea_water_electrical_conductivity', axes = ax1)
s.print_tree('CTD_NeilBrown')
s.plot_timeseries('CTD_NeilBrown/sea_water_temperature', axes = ax2)
s.plot_timeseries('CTD_NeilBrown/sea_water_pressure', axes = ax3)
s.plot_timeseries('CTD_NeilBrown/sea_water_salinity', axes = ax4)
s.plot_timeseries('CTD_NeilBrown/depth', axes = ax5)
for ax in fig.get_axes(): ax.legend()
for ax in fig.get_axes(): ax.grid(True)
draw()
ax1 = fig.add_subplot(5,1,1)
s.plot_timeseries('CTD_NeilBrown/sea_water_electrical_conductivity', axes = ax1)
draw()
clf()
fig = plt.figure()
ax1 = fig.add_subplot(5,1,1)
s.plot_timeseries('CTD_NeilBrown/sea_water_electrical_conductivity', axes = ax1)
draw()
ax2 = fig.add_subplot(5,1,2, sharex = ax1)
ax3 = fig.add_subplot(5,1,3, sharex = ax1)
ax4 = fig.add_subplot(5,1,4, sharex = ax1)
ax5 = fig.add_subplot(5,1,5, sharex = ax1)
s.plot_timeseries('CTD_NeilBrown/depth', axes = ax5)
s.plot_timeseries('CTD_NeilBrown/sea_water_temperature', axes = ax2)
s.plot_timeseries('CTD_NeilBrown/sea_water_pressure', axes = ax3)
s.plot_timeseries('CTD_NeilBrown/sea_water_salinity', axes = ax4)
draw()
for ax in fig.get_axes(): ax.grid(True)
for ax in fig.get_axes(): ax.legend()
draw()
fig.suptitle(fn.rsplit('/')[-1])
draw()
get_ipython().magic(u'logstart /home/squall/Desktop/ctd.py')
exit()