MAT1 = np.zeros([8, d.shape[0], d.shape[1]])
MAT1[0] = d.D1mag
MAT1[1] = d.D2mag
MAT1[2] = d.D1phase
MAT1[3] = d.D2phase
MAT1[4] = d.D1Pow
MAT1[5] = d.D2Pow
MAT1[6] = d.D1Pow2
MAT1[7] = d.D2Pow2
header1 = make_header(d.n1, d.n2, d.n2, meas_data='est. Photon #')
savemtx('output//S1_480_V.mtx', MAT1, header=header1)
y = (d.n2.lin*d.n2.lin/50.0) # position of the data
MAT2 = np.zeros([2, d.shape[0]*10, d.shape[1]])
y2, MAT2[0] = interp_y(y, d.D1Pow2, factor=10)
y2, MAT2[1] = interp_y(y, d.D2Pow2, factor=10)
d.dim_y2 = d.n2
d.dim_y2.start = y2[0]*1e3
d.dim_y2.stop = y2[-1]*1e3
d.dim_y2.pt = len(y2)
d.dim_y2.lin = y2*1e3
d.dim_y2.name = 'Pump power (mW)'
header2 = make_header(d.n1, d.dim_y2, d.n2, meas_data='est. Photon #')
savemtx('output//S1_480_mW.mtx', MAT2, header=header2)
d.shape = d.data.shape
xoff = 139.3e-3
x1flux = 479.6e-3
d.dim_1.lin = (d.dim_1.lin-xoff)/x1flux + 0.5
d.dim_1.start = d.dim_1.lin[0]
d.dim_1.stop = d.dim_1.lin[-1]
MAT1 = np.zeros([6, d.shape[0], d.shape[2]])
MAT1[0] = d.data[:, 3, :]
MAT1[1] = d.data[:, 4, :]
MAT1[2] = d.D1mag
MAT1[3] = d.D1ang
MAT1[4] = d.D1pow
MAT1[5] = d.D1lev
y = (d.dim_2.lin*d.dim_2.lin/50.0) # position of the data
y2, MAT2 = interp_y(y, d.D1pow)
d.dim_y2 = d.dim_2
d.dim_y2.start = y2[0]
d.dim_y2.stop = y2[-1]
d.dim_y2.pt = len(y2)
d.dim_y2.lin = y2
header1 = make_header(d.dim_1, d.dim_2, d.dim_3, meas_data='(a.u)')
savemtx('S1_471_DCE_MAP2_volt.mtx', MAT1, header=header1)
header2 = make_header(d.dim_1, d.dim_y2, d.dim_3, meas_data='(a.u)')
savemtx('S1_471_DCE_MAP2_wats.mtx', MAT2, header=header2)
