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], '.')