out['tsample'] = tsample
out['Vsample'] = Vsample
return out
files = ['data/b1_1.dat', 'data/b2_3.dat']
tstart = [24., 24.]
tstop = [344., 600.]
styles = ['o', 's']
LAB = ['20scfh', '25scfh']
plt.close('all')
lplot.set_defaults(screen_dpi=112)
ax0 = lplot.init_fig('Floating Potential (V)',
'Saturation Current ($\mu$A)',
figure_size=(6., 4.5))
for index in range(len(files)):
data = lconfig.dfile(files[index], afun=fsanalysis, afile=None)
t = np.array(data.analysis['tsample'])
V = np.array(data.analysis['Vsample'])
I = np.array(data.analysis['Isample'])
ii = (t > tstart[index]) * (t < tstop[index])
ax0.plot(V[ii], I[ii], styles[index], mec='k', mfc='w', label=LAB[index])
ax0.legend(loc=0)
ax0.get_figure().savefig('float_sat.png')
C225_1.add_dir('data/225_1')
for dat in C225_1.data:
dat.analysis['d_in'] = 2.25
C.merge(C225_1)
C.table(('d_in', 'flow_scfh', 'ratio_fto', 'plate_t_c', 'tip_t_c',
'plate_q_kw', 'standoff_in', 'b1', 'R1', 'b2', 'R2', 'b3', 'R3'),
fileout='table.wsv')
# Construct comparison plots
plt.close('all')
lplot.set_defaults()
# Figure 1
# Overplot the IV curves for tests
ax = lplot.init_fig('Voltage (V)', 'Current ($\mu$A)')
for dat in C150_1:
ax.plot(dat.analysis['V'], dat.analysis['I'], 'y.')
for dat in C105:
ax.plot(dat.analysis['V'], dat.analysis['I'], 'b.')
for dat in C105_1:
ax.plot(dat.analysis['V'], dat.analysis['I'], 'r.')
ax.legend(loc=0)
ax.set_xlim([-10, 10])
# Figure 2
# Show tests with 20scfh
out['Tsample'] = Tsample
out['Isample'] = Isample
out['tsample'] = tsample
out['Vsample'] = Vsample
return out
datafile = 'data/b2_1.dat'
#datafile = 'data/b1_1.dat'
data = lconfig.dfile(datafile, afun=fsanalysis, afile=None)
plt.close('all')
lplot.set_defaults(screen_dpi=112)
ax0 = lplot.init_fig('time (s)', 'Current ($\mu$A)', figure_size=(6., 4.5))
I = median_filter(data.analysis['I'], 20)
ax0.plot(data.t(), data.analysis['I'], color=[.45, .45, .45])
ax0.plot(data.t(), I, 'k')
ax0.set_ylim([0, 120])
zax = lplot.zoom_ax(ax0, [.27, .65, .4, .27],
xlim=[64.2, 65.6],
ylim=[20., 65.],
box=False)
zax.set_xticks([64.5, 65., 65.5])
zax.set_yticks([20, 40, 60])
ax0.get_figure().savefig('ivst.png')
ax0.get_figure().savefig('ivst.')
ax0.get_figure().savefig('ivst.png')
index = data.analysis['T'] > 100.