dat.analysis['d_in'] = 2.25
C.merge(C225)
C225_1 = lconfig.collection(afun=fsanalysis, asave=False)
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)
# look for the positive 0 crossing in voltage
if I[index]<0. and I[index+1]>=0.:
Vsample.append(\
V[index] - (V[index+1] - V[index])/(I[index+1] - I[index]) * I[index])
Is = I[index:index+100].tolist()
Is.sort(reverse=True)
tsample.append(t[index+50])
Isample.append(Is[4]) # Reject the 4 largest points
Tsample.append(np.mean(out['T'][index:index+100]))
index+=120
else:
index += 1
out['Tsample'] = Tsample
out['Isample'] = Isample
out['tsample'] = tsample
out['Vsample'] = Vsample
return out
datafile = 'data/a4_1.dat'
#datafile = 'data/b1_1.dat'
data = lconfig.dfile(datafile, afun=fsanalysis, afile=None)
plt.close('all')
lplot.set_defaults(screen_dpi=112)
ax = lplot.init_fig('Voltage (V)', 'Current ($\mu$A)', figure_size=(6., 4.5))
index = data.t()>15.
ax.plot(data.analysis['V'][index], data.analysis['I'][index], '.')