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.