)
wks_M01.from_list(i+1, M01_list[i])
wks_M02.from_list(i+1, M02_list[i])
wks_M03.from_list(i+1, M03_list[i])
wks_M10.from_list(i+1, M10_list[i])
wks_M11.from_list(i+1, M11_list[i])
wks_M12.from_list(i+1, M12_list[i])
wks_M13.from_list(i+1, M13_list[i])
wks_M20.from_list(i+1, M20_list[i])
wks_M21.from_list(i+1, M21_list[i])
wks_M22.from_list(i+1, M22_list[i])
wks_M23.from_list(i+1, M23_list[i])
wks_M30.from_list(i+1, M30_list[i])
wks_M31.from_list(i+1, M31_list[i])
wks_M32.from_list(i+1, M32_list[i])
wks_M33.from_list(i+1, M33_list[i])

path_template = os.getcwd()

new_path_template = path_template.replace('/', '\\')

gr = op.new_graph(template=(new_path_template + '\\Origin templates\\test_plot16.otpu'), lname='16MM plot')

poz = 0
neg = 0

for i in range(len(csv_file)):
if csv_file['M00'][i] > 0:
poz = poz + 1
else:
neg = neg + 1

gl_1 = gr[0]
p1 = gl_1.add_plot(f'{wks.lt_range()}!(?,1:end)')
p1.set_int('line.width', 3)
if (neg >= 1) and poz >=1:
gl_1.set_int("x.AtZero",1)
gl_1.group()
p1.colormap = 'Candy'
gl_1.rescale()
gl_1.xlim = (min(csv_file['Wavelength']), max(csv_file['Wavelength']))
lgnd_M00 = gr[0].label('Legend')
lgnd_M00.text = ' %(1)'

poz = 0
neg = 0

for i in range(len(csv_file)):
if csv_file['M01'][i] > 0:
poz = poz + 1
else:
neg = neg + 1

gl_2 = gr[1]
p2 = gl_2.add_plot(f'{wks_M01.lt_range()}!(?,1:end)')
p2.set_int('line.width', 3)
if (neg >= 1) and poz >=1:
gl_1.set_int("x.AtZero",1)
gl_2.group()
p2.colormap = 'Candy'
gl_2.rescale()
gl_2.xlim = (min(csv_file['Wavelength']), max(csv_file['Wavelength']))
lgnd_M01 = gr[1].label('Legend')
lgnd_M01.text = ' %(1)'

poz = 0
neg = 0

for i in range(len(csv_file)):
if csv_file['M02'][i] > 0:
poz = poz + 1
else:
neg = neg + 1

gl_3 = gr[2]
p3 = gl_3.add_plot(f'{wks_M02.lt_range()}!(?,1:end)')
p3.set_int('line.width', 3)
if (neg >= 1) and poz >=1:
gl_3.set_int("x.AtZero",1)
gl_3.group()
p3.colormap = 'Candy'
gl_3.rescale()
gl_3.xlim = (min(csv_file['Wavelength']), max(csv_file['Wavelength']))
lgnd_M02 = gr[2].label('Legend')
lgnd_M02.text = ' %(1)'

poz = 0
neg = 0

for i in range(len(csv_file)):
if csv_file['M03'][i] > 0:
poz = poz + 1
else:
neg = neg + 1

gl_4 = gr[3]
p4 = gl_4.add_plot(f'{wks_M03.lt_range()}!(?,1:end)')
p4.set_int('line.width', 3)
if (neg >= 1) and poz >=1:
gl_4.set_int("x.AtZero",1)
gl_4.group()
p4.colormap = 'Candy'
gl_4.rescale()
gl_4.xlim = (min(csv_file['Wavelength']), max(csv_file['Wavelength']))
lgnd_M03 = gr[3].label('Legend')
lgnd_M03.text = ' %(1)'

poz = 0
neg = 0

for i in range(len(csv_file)):
if csv_file['M10'][i] > 0:
poz = poz + 1
else:
neg = neg + 1

gl_5 = gr[4]
p5 = gl_5.add_plot(f'{wks_M10.lt_range()}!(?,1:end)')
p5.set_int('line.width', 3)
if (neg >= 1) and poz >=1:
gl_5.set_int("x.AtZero",1)
gl_5.group()
p5.colormap = 'Candy'
gl_5.rescale()
gl_5.xlim = (min(csv_file['Wavelength']), max(csv_file['Wavelength']))
lgnd_M10 = gr[4].label('Legend')
lgnd_M10.text = ' %(1)'

poz = 0
neg = 0

for i in range(len(csv_file)):
if csv_file['M11'][i] > 0:
poz = poz + 1
else:
neg = neg + 1

gl_6 = gr[5]
p6 = gl_6.add_plot(f'{wks_M11.lt_range()}!(?,1:end)')
p6.set_int('line.width', 3)
if (neg >= 1) and poz >=1:
gl_6.set_int("x.AtZero",1)
gl_6.group()
p6.colormap = 'Candy'
gl_6.rescale()
gl_6.xlim = (min(csv_file['Wavelength']), max(csv_file['Wavelength']))
lgnd_M11 = gr[5].label('Legend')
lgnd_M11.text = ' %(1)'

poz = 0
neg = 0

for i in range(len(csv_file)):
if csv_file['M12'][i] > 0:
poz = poz + 1
else:
neg = neg + 1

gl_7 = gr[6]
p7 = gl_7.add_plot(f'{wks_M12.lt_range()}!(?,1:end)')
p7.set_int('line.width', 3)
if (neg >= 1) and poz >=1:
gl_7.set_int("x.AtZero",1)
gl_7.group()
p7.colormap = 'Candy'
gl_7.rescale()
gl_7.xlim = (min(csv_file['Wavelength']), max(csv_file['Wavelength']))
lgnd_M12 = gr[6].label('Legend')
lgnd_M12.text = ' %(1)'

poz = 0
neg = 0

for i in range(len(csv_file)):
if csv_file['M13'][i] > 0:
poz = poz + 1
else:
neg = neg + 1

gl_8 = gr[7]
p8 = gl_8.add_plot(f'{wks_M13.lt_range()}!(?,1:end)')
p8.set_int('line.width', 3)
if (neg >= 1) and poz >=1:
gl_8.set_int("x.AtZero",1)
gl_8.group()
p8.colormap = 'Candy'
gl_8.rescale()
gl_8.xlim = (min(csv_file['Wavelength']), max(csv_file['Wavelength']))
lgnd_M13 = gr[7].label('Legend')
lgnd_M13.text = ' %(1)'

poz = 0
neg = 0

for i in range(len(csv_file)):
if csv_file['M20'][i] > 0:
poz = poz + 1
else:
neg = neg + 1

gl_9 = gr
p9 = gl_9.add_plot(f'{wks_M20.lt_range()}!(?,1:end)')
p9.set_int('line.width', 3)
if (neg >= 1) and poz >=1:
gl_9.set_int("x.AtZero",1)
gl_9.group()
p9.colormap = 'Candy'
gl_9.rescale()
gl_9.xlim = (min(csv_file['Wavelength']), max(csv_file['Wavelength']))
lgnd_M20 = gr.label('Legend')
lgnd_M20.text = ' %(1)'

poz = 0
neg = 0

for i in range(len(csv_file)):
if csv_file['M21'][i] > 0:
poz = poz + 1
else:
neg = neg + 1

gl_10 = gr[9]
p10 = gl_10.add_plot(f'{wks_M21.lt_range()}!(?,1:end)')
p10.set_int('line.width', 3)
if (neg >= 1) and poz >=1:
gl_10.set_int("x.AtZero",1)
gl_10.group()
p10.colormap = 'Candy'
gl_10.rescale()
gl_10.xlim = (min(csv_file['Wavelength']), max(csv_file['Wavelength']))
lgnd_M21 = gr[9].label('Legend')
lgnd_M21.text = ' %(1)'

poz = 0
neg = 0

for i in range(len(csv_file)):
if csv_file['M22'][i] > 0:
poz = poz + 1
else:
neg = neg + 1

gl_11 = gr[10]
p11 = gl_11.add_plot(f'{wks_M22.lt_range()}!(?,1:end)')
p11.set_int('line.width', 3)
if (neg >= 1) and poz >=1:
gl_11.set_int("x.AtZero",1)
gl_11.group()
p11.colormap = 'Candy'
gl_11.rescale()
gl_11.xlim = (min(csv_file['Wavelength']), max(csv_file['Wavelength']))
lgnd_M22 = gr[10].label('Legend')
lgnd_M22.text = ' %(1)'

poz = 0
neg = 0

for i in range(len(csv_file)):
if csv_file['M23'][i] > 0:
poz = poz + 1
else:
neg = neg + 1

gl_12 = gr[11]
p12 = gl_12.add_plot(f'{wks_M23.lt_range()}!(?,1:end)')
p12.set_int('line.width', 3)
if (neg >= 1) and poz >=1:
gl_12.set_int("x.AtZero",1)
gl_12.group()
p12.colormap = 'Candy'
gl_12.rescale()
gl_12.xlim = (min(csv_file['Wavelength']), max(csv_file['Wavelength']))
lgnd_M23 = gr[11].label('Legend')
lgnd_M23.text = ' %(1)'

poz = 0
neg = 0

for i in range(len(csv_file)):
if csv_file['M30'][i] > 0:
poz = poz + 1
else:
neg = neg + 1

gl_13 = gr[12]
p13 = gl_13.add_plot(f'{wks_M30.lt_range()}!(?,1:end)')
p13.set_int('line.width', 3)
if (neg >= 1) and poz >=1:
gl_13.set_int("x.AtZero",1)
gl_13.group()
p13.colormap = 'Candy'
gl_13.rescale()
gl_13.xlim = (min(csv_file['Wavelength']), max(csv_file['Wavelength']))
lgnd_M30 = gr[12].label('Legend')
lgnd_M30.text = ' %(1)'

poz = 0
neg = 0

for i in range(len(csv_file)):
if csv_file['M31'][i] > 0:
poz = poz + 1
else:
neg = neg + 1

gl_14 = gr[13]
p14 = gl_14.add_plot(f'{wks_M31.lt_range()}!(?,1:end)')
p14.set_int('line.width', 3)
if (neg >= 1) and poz >=1:
gl_14.set_int("x.AtZero",1)
gl_14.group()
p14.colormap = 'Candy'
gl_14.rescale()
gl_14.xlim = (min(csv_file['Wavelength']), max(csv_file['Wavelength']))
lgnd_M31 = gr[13].label('Legend')
lgnd_M31.text = ' %(1)'

poz = 0
neg = 0

for i in range(len(csv_file)):
if csv_file['M32'][i] > 0:
poz = poz + 1
else:
neg = neg + 1

gl_15 = gr[14]
p15 = gl_15.add_plot(f'{wks_M32.lt_range()}!(?,1:end)')
p15.set_int('line.width', 3)
if (neg >= 1) and poz >=1:
gl_15.set_int("x.AtZero",1)
gl_15.group()
p15.colormap = 'Candy'
gl_15.rescale()
gl_15.xlim = (min(csv_file['Wavelength']), max(csv_file['Wavelength']))
lgnd_M32 = gr[14].label('Legend')
lgnd_M32.text = ' %(1)'

poz = 0
neg = 0

for i in range(len(csv_file)):
if csv_file['M33'][i] > 0:
poz = poz + 1
else:
neg = neg + 1

gl_16 = gr[15]
p16 = gl_16.add_plot(f'{wks_M33.lt_range()}!(?,1:end)')
p16.set_int('line.width', 3)
if (neg >= 1) and poz >=1:
gl_16.set_int("x.AtZero",1)
gl_16.group()
p16.colormap = 'Candy'
gl_16.rescale()
gl_16.xlim = (min(csv_file['Wavelength']), max(csv_file['Wavelength']))
lgnd_M33 = gr[15].label('Legend')
lgnd_M33.text = ' %(1)'


wks_ind = op.new_sheet(lname = 'Ind method') # for normal plots
wks_ind.name = 'CD-ind'

wks_ind.from_list(0, Wavelength_new, lname='Wavelength', units='nm')
wks_ind.from_list(1, CD_list_new[0], lname='CD-ind', units='mdeg')

wks_ind_CB = wks_ind.get_book().add_sheet('CB-ind')
wks_ind_CB.from_list(0, Wavelength_new, lname='Wavelength', units='nm')
wks_ind_CB.from_list(1, CB_list_new[0], lname='CB-ind', units='mdeg')

wks_ind_Abs = wks_ind.get_book().add_sheet('Abs')
wks_ind_Abs.from_list(0, Wavelength_new, lname='Wavelength', units='nm')
wks_ind_Abs.from_list(1, Abs_list[0], lname='Abs', units = 'au')

wks_ind_LD = wks_ind.get_book().add_sheet('LD-ind')
wks_ind_LD.from_list(0, Wavelength_new, lname='Wavelength', units='nm')
wks_ind_LD.from_list(1, LD_list_new[0], lname='LD-ind')

wks_ind_LB = wks_ind.get_book().add_sheet('LB-ind')
wks_ind_LB.from_list(0, Wavelength_new, lname='Wavelength', units='nm')
wks_ind_LB.from_list(1, LB_list_new[0], lname='LB-ind')

wks_ind_g_factor = wks_ind.get_book().add_sheet('G-factor-ind')
wks_ind_g_factor.from_list(0, Wavelength_new, lname='Wavelength', units='nm')
wks_ind_g_factor.from_list(1, g_factor_list_new[0], lname='G-factor-ind', units='CD/Abs')

wks_ind_LDp = wks_ind.get_book().add_sheet("LD'-ind")
wks_ind_LDp.from_list(0, Wavelength_new, lname='Wavelength', units='nm')
wks_ind_LDp.from_list(1, LDp_list_new[0], lname="LD'-ind")

wks_ind_LBp = wks_ind.get_book().add_sheet("LB'-ind")
wks_ind_LBp.from_list(0, Wavelength_new, lname='Wavelength', units='nm')
wks_ind_LBp.from_list(1, LBp_list_new[0], lname="LB'-ind")

wks_ind_T = wks_ind.get_book().add_sheet("Transmission")
wks_ind_T.from_list(0, Wavelength_new, lname='Wavelength', units='nm')
wks_ind_T.from_list(1, M00_avg, lname="Transmission_M00")
wks_ind_T.from_list(2, M11_avg, lname="Transmission_M11")
wks_ind_T.from_list(3, M22_avg, lname="Transmission_M22")
wks_ind_T.from_list(4, M33_avg, lname="Transmission_M33")

if len(M00_list) > 1:
for i in range(1, len(M00_list)):
wks_ind.from_list(i+1, CD_list_new[i])
wks_ind_CB.from_list(i+1, CB_list_new[i])
wks_ind_Abs.from_list(i+1, Abs_list[i])
wks_ind_LD.from_list(i+1, LD_list_new[i])
wks_ind_LB.from_list(i+1, LB_list_new[i])
wks_ind_g_factor.from_list(i+1, g_factor_list_new[i])
wks_ind_LDp.from_list(i+1, LDp_list_new[i])
wks_ind_LBp.from_list(i+1, LBp_list_new[i])

gr1 = op.new_graph(template=(new_path_template + '\\Origin templates\\test_plot9.otpu'), lname='Ind method plot')
# gr1 = op.new_graph(template=op.path('e') + 'test_plot9.otpu')
# gr1.name = 'Ind method plot'

poz = 0
neg = 0

for i in range(len(CD_list_new)):
for j in range(len(CD_list_new[i])):
if CD_list_new[i][j] > 0:
poz = poz + 1
else:
neg = neg + 1

gl_iCD = gr1[0]
p1_iCD = gl_iCD.add_plot(f'{wks_ind.lt_range()}!(?,1:end)')
p1_iCD.set_int('line.width', 3)
if (neg >= 1) and poz >=1:
gl_iCD.set_int("x.AtZero",1)
gl_iCD.group()
p1_iCD.colormap = 'Candy'
gl_iCD.rescale()
gl_iCD.xlim = (min(csv_file['Wavelength']), max(csv_file['Wavelength']))

poz = 0
neg = 0

for i in range(len(CD_list_new)):
for j in range(len(CD_list_new[i])):
if CB_list_new[i][j] > 0:
poz = poz + 1
else:
neg = neg + 1

gl_iCB= gr1[1]
p1_iCB = gl_iCB.add_plot(f'{wks_ind_CB.lt_range()}!(?,1:end)')
p1_iCB.set_int('line.width', 3)
if (neg >= 1) and poz >=1:
gl_iCB.set_int("x.AtZero",1)
gl_iCB.group()
p1_iCB.colormap = 'Candy'
gl_iCB.rescale()
gl_iCB.xlim = (min(csv_file['Wavelength']), max(csv_file['Wavelength']))

poz = 0
neg = 0

for i in range(len(CD_list_new)):
for j in range(len(CD_list_new[i])):
if Abs_list[i][j] > 0:
poz = poz + 1
else:
neg = neg + 1

gl_iAbs= gr1[2]
p1_iAbs = gl_iAbs.add_plot(f'{wks_ind_Abs.lt_range()}!(?,1:end)')
p1_iAbs.set_int('line.width', 3)
if (neg >= 1) and poz >=1:
gl_iAbs.set_int("x.AtZero",1)
gl_iAbs.group()
p1_iAbs.colormap = 'Candy'
gl_iAbs.rescale()
gl_iAbs.xlim = (min(csv_file['Wavelength']), max(csv_file['Wavelength']))

poz = 0
neg = 0

for i in range(len(CD_list_new)):
for j in range(len(CD_list_new[i])):
if LD_list_new[i][j] > 0:
poz = poz + 1
else:
neg = neg + 1

gl_iLD= gr1[3]
p1_iLD = gl_iLD.add_plot(f'{wks_ind_LD.lt_range()}!(?,1:end)')
p1_iLD.set_int('line.width', 3)
if (neg >= 1) and poz >=1:
gl_iLD.set_int("x.AtZero",1)
gl_iLD.group()
p1_iLD.colormap = 'Candy'
gl_iLD.rescale()
gl_iLD.xlim = (min(csv_file['Wavelength']), max(csv_file['Wavelength']))

poz = 0
neg = 0

for i in range(len(CD_list_new)):
for j in range(len(CD_list_new[i])):
if LB_list_new[i][j] > 0:
poz = poz + 1
else:
neg = neg + 1

gl_iLB= gr1[4]
p1_iLB = gl_iLB.add_plot(f'{wks_ind_LB.lt_range()}!(?,1:end)')
p1_iLB.set_int('line.width', 3)
if (neg >= 1) and poz >=1:
gl_iLB.set_int("x.AtZero",1)
gl_iLB.group()
p1_iLB.colormap = 'Candy'
gl_iLB.rescale()
gl_iLB.xlim = (min(csv_file['Wavelength']), max(csv_file['Wavelength']))

poz = 0
neg = 0

for i in range(len(CD_list_new)):
for j in range(len(CD_list_new[i])):
if g_factor_list_new[i][j] > 0:
poz = poz + 1
else:
neg = neg + 1

gl_ig= gr1[5]
p1_ig = gl_ig.add_plot(f'{wks_ind_g_factor.lt_range()}!(?,1:end)')
p1_ig.set_int('line.width', 3)
if (neg >= 1) and poz >=1:
gl_ig.set_int("x.AtZero",1)
gl_ig.group()
p1_ig.colormap = 'Candy'
gl_ig.rescale()
gl_ig.xlim = (min(csv_file['Wavelength']), max(csv_file['Wavelength']))

poz = 0
neg = 0

for i in range(len(CD_list_new)):
for j in range(len(CD_list_new[i])):
if LDp_list_new[i][j] > 0:
poz = poz + 1
else:
neg = neg + 1

gl_iLDp= gr1[6]
p1_iLDp = gl_iLDp.add_plot(f'{wks_ind_LDp.lt_range()}!(?,1:end)')
p1_iLDp.set_int('line.width', 3)
if (neg >= 1) and poz >=1:
gl_iLDp.set_int("x.AtZero",1)
gl_iLDp.group()
p1_iLDp.colormap = 'Candy'
gl_iLDp.rescale()
gl_iLDp.xlim = (min(csv_file['Wavelength']), max(csv_file['Wavelength']))

poz = 0
neg = 0

for i in range(len(CD_list_new)):
for j in range(len(CD_list_new[i])):
if LBp_list_new[i][j] > 0:
poz = poz + 1
else:
neg = neg + 1

gl_iLBp= gr1[7]
p1_iLBp = gl_iLBp.add_plot(f'{wks_ind_LBp.lt_range()}!(?,1:end)')
p1_iLBp.set_int('line.width', 3)
if (neg >= 1) and poz >=1:
gl_iLBp.set_int("x.AtZero",1)
gl_iLBp.group()
p1_iLBp.colormap = 'Candy'
gl_iLBp.rescale()
gl_iLBp.xlim = (min(csv_file['Wavelength']), max(csv_file['Wavelength']))

poz = 0
neg = 0

for i in range(len(M00_avg)):
if M00_avg[i] > 0 or M11_avg[i] > 0 or M22_avg[i] > 0 or M33_avg[i] > 0:
poz = poz + 1
elif M00_avg[i] < 0 or M11_avg[i] < 0 or M22_avg[i] < 0 or M33_avg[i] < 0:
neg = neg + 1

gl_iT= gr1
p1_iT = gl_iT.add_plot(f'{wks_ind_T.lt_range()}!(?,1:end)')
p1_iT.set_int('line.width', 3)
if (neg >= 1) and poz >=1:
gl_iT.set_int("x.AtZero",1)
gl_iT.group()
p1_iT.colormap = 'Candy'
gl_iT.rescale()
gl_iT.xlim = (min(csv_file['Wavelength']), max(csv_file['Wavelength']))


wks_m = op.new_sheet(lname = 'Matrix ln method') # for normal plots
wks_m.name = 'mCD'

wks_m.from_list(0, Wavelength_new, lname='Wavelength', units='nm')
wks_m.from_list(1, mCD_list[0], lname='mCD', units='mdeg')

wks_m_CB = wks_m.get_book().add_sheet("mCB")
wks_m_CB.from_list(0, Wavelength_new, lname='Wavelength', units='nm')
wks_m_CB.from_list(1, mCB_list[0], lname='mCB', units='mdeg')

wks_m_Abs = wks_m.get_book().add_sheet("Abs")
wks_m_Abs.from_list(0, Wavelength_new, lname='Wavelength', units='nm')
wks_m_Abs.from_list(1, Abs_list[0], lname='Abs', units = 'au')

wks_m_LD = wks_m.get_book().add_sheet("mLD")
wks_m_LD.from_list(0, Wavelength_new, lname='Wavelength', units='nm')
wks_m_LD.from_list(1, mLD_list[0], lname='mLD')

wks_m_LB = wks_m.get_book().add_sheet("mLB")
wks_m_LB.from_list(0, Wavelength_new, lname='Wavelength', units='nm')
wks_m_LB.from_list(1, mLB_list[0], lname='mLB')

wks_m_g_factor = wks_m.get_book().add_sheet("mG-factor")
wks_m_g_factor.from_list(0, Wavelength_new, lname='Wavelength', units='nm')
wks_m_g_factor.from_list(1, mg_factor_list[0], lname='mG-factor', units='CD/Abs')

wks_m_LDp = wks_m.get_book().add_sheet("mLD'")
wks_m_LDp.from_list(0, Wavelength_new, lname='Wavelength', units='nm')
wks_m_LDp.from_list(1, mLDp_list[0], lname="mLD'")

wks_m_LBp = wks_m.get_book().add_sheet("mLB'")
wks_m_LBp.from_list(0, Wavelength_new, lname='Wavelength', units='nm')
wks_m_LBp.from_list(1, mLBp_list[0], lname="mLB'")

wks_m_T = wks_m.get_book().add_sheet("Ts")
wks_m_T.from_list(0, Wavelength_new, lname='Wavelength', units='nm')
wks_m_T.from_list(1, M00_avg, lname="Transmission_M00")
wks_m_T.from_list(2, M11_avg, lname="Transmission_M11")
wks_m_T.from_list(3, M22_avg, lname="Transmission_M22")
wks_m_T.from_list(4, M33_avg, lname="Transmission_M33")

if len(M00_list) > 1:
for i in range(1, len(M00_list)):
wks_m.from_list(i+1, mCD_list[i])
wks_m_CB.from_list(i+1, mCB_list[i])
wks_m_Abs.from_list(i+1, Abs_list[i])
wks_m_LD.from_list(i+1, mLD_list[i])
wks_m_LB.from_list(i+1, mLB_list[i])
wks_m_g_factor.from_list(i+1, mg_factor_list[i])
wks_m_LDp.from_list(i+1, mLDp_list[i])
wks_m_LBp.from_list(i+1, mLBp_list[i])

gr2 = op.new_graph(template=(new_path_template + '\\Origin templates\\test_plot9.otpu'), lname='Matrix method plot')

# gr2 = op.new_graph(template=op.path('e') + 'test_plot9.otpu')
# gr2.name = 'Matrix ln method plot'

poz = 0
neg = 0

for i in range(len(CD_list_new)):
for j in range(len(CD_list_new[i])):
if mCD_list[i][j] > 0:
poz = poz + 1
else:
neg = neg + 1

gl_mCD = gr2[0]
p1_mCD = gl_mCD.add_plot(f'{wks_m.lt_range()}!(?,1:end)')
p1_mCD.set_int('line.width', 3)
if (neg >= 1) and poz >=1:
gl_mCD.set_int("x.AtZero",1)
gl_mCD.group()
p1_mCD.colormap = 'Candy'
gl_mCD.rescale()
gl_mCD.xlim = (min(csv_file['Wavelength']), max(csv_file['Wavelength']))

poz = 0
neg = 0

for i in range(len(CD_list_new)):
for j in range(len(CD_list_new[i])):
if mCB_list[i][j] > 0:
poz = poz + 1
else:
neg = neg + 1

gl_mCB= gr2[1]
p1_mCB = gl_mCB.add_plot(f'{wks_m_CB.lt_range()}!(?,1:end)')
p1_mCB.set_int('line.width', 3)
if (neg >= 1) and poz >=1:
gl_mCB.set_int("x.AtZero",1)
gl_mCB.group()
p1_mCB.colormap = 'Candy'
gl_mCB.rescale()
gl_mCB.xlim = (min(csv_file['Wavelength']), max(csv_file['Wavelength']))

poz = 0
neg = 0

for i in range(len(CD_list_new)):
for j in range(len(CD_list_new[i])):
if Abs_list[i][j] > 0:
poz = poz + 1
else:
neg = neg + 1

gl_mAbs= gr2[2]
p1_mAbs = gl_mAbs.add_plot(f'{wks_m_Abs.lt_range()}!(?,1:end)')
p1_mAbs.set_int('line.width', 3)
if (neg >= 1) and poz >=1:
gl_mAbs.set_int("x.AtZero",1)
gl_mAbs.group()
p1_mAbs.colormap = 'Candy'
gl_mAbs.rescale()
gl_mAbs.xlim = (min(csv_file['Wavelength']), max(csv_file['Wavelength']))

poz = 0
neg = 0

for i in range(len(CD_list_new)):
for j in range(len(CD_list_new[i])):
if mLD_list[i][j] > 0:
poz = poz + 1
else:
neg = neg + 1

gl_mLD= gr2[3]
p1_mLD = gl_mLD.add_plot(f'{wks_m_LD.lt_range()}!(?,1:end)')
p1_mLD.set_int('line.width', 3)
if (neg >= 1) and poz >=1:
gl_mLD.set_int("x.AtZero",1)
gl_mLD.group()
p1_mLD.colormap = 'Candy'
gl_mLD.rescale()
gl_mLD.xlim = (min(csv_file['Wavelength']), max(csv_file['Wavelength']))

poz = 0
neg = 0

for i in range(len(CD_list_new)):
for j in range(len(CD_list_new[i])):
if mLB_list[i][j] > 0:
poz = poz + 1
else:
neg = neg + 1

gl_mLB= gr2[4]
p1_mLB = gl_mLB.add_plot(f'{wks_m_LB.lt_range()}!(?,1:end)')
p1_mLB.set_int('line.width', 3)
if (neg >= 1) and poz >=1:
gl_mLB.set_int("x.AtZero",1)
gl_mLB.group()
p1_mLB.colormap = 'Candy'
gl_mLB.rescale()
gl_mLB.xlim = (min(csv_file['Wavelength']), max(csv_file['Wavelength']))

poz = 0
neg = 0

for i in range(len(CD_list_new)):
for j in range(len(CD_list_new[i])):
if mg_factor_list[i][j] > 0:
poz = poz + 1
else:
neg = neg + 1

gl_mg= gr2[5]
p1_mg = gl_mg.add_plot(f'{wks_m_g_factor.lt_range()}!(?,1:end)')
p1_mg.set_int('line.width', 3)
if (neg >= 1) and poz >=1:
gl_mg.set_int("x.AtZero",1)
gl_mg.group()
p1_mg.colormap = 'Candy'
gl_mg.rescale()
gl_mg.xlim = (min(csv_file['Wavelength']), max(csv_file['Wavelength']))

poz = 0
neg = 0

for i in range(len(CD_list_new)):
for j in range(len(CD_list_new[i])):
if mLDp_list[i][j] > 0:
poz = poz + 1
else:
neg = neg + 1

gl_mLDp= gr2[6]
p1_mLDp = gl_mLDp.add_plot(f'{wks_m_LDp.lt_range()}!(?,1:end)')
p1_mLDp.set_int('line.width', 3)
if (neg >= 1) and poz >=1:
gl_mLDp.set_int("x.AtZero",1)
gl_mLDp.group()
p1_mLDp.colormap = 'Candy'
gl_mLDp.rescale()
gl_mLDp.xlim = (min(csv_file['Wavelength']), max(csv_file['Wavelength']))

poz = 0
neg = 0

for i in range(len(CD_list_new)):
for j in range(len(CD_list_new[i])):
if mLBp_list[i][j] > 0:
poz = poz + 1
else:
neg = neg + 1

gl_mLBp= gr2[7]
p1_mLBp = gl_mLBp.add_plot(f'{wks_m_LBp.lt_range()}!(?,1:end)')
p1_mLBp.set_int('line.width', 3)
if (neg >= 1) and poz >=1:
gl_mLBp.set_int("x.AtZero",1)
gl_mLBp.group()
p1_mLBp.colormap = 'Candy'
gl_mLBp.rescale()
gl_mLBp.xlim = (min(csv_file['Wavelength']), max(csv_file['Wavelength']))

poz = 0
neg = 0

for i in range(len(M00_avg)):
if M00_avg[i] > 0 or M11_avg[i] > 0 or M22_avg[i] > 0 or M33_avg[i] > 0:
poz = poz + 1
elif M00_avg[i] < 0 or M11_avg[i] < 0 or M22_avg[i] < 0 or M33_avg[i] < 0:
neg = neg + 1

gl_mT= gr2
p1_mT = gl_mT.add_plot(f'{wks_m_T.lt_range()}!(?,1:end)')
p1_mT.set_int('line.width', 3)
if (neg >= 1) and poz >=1:
gl_mT.set_int("x.AtZero",1)
gl_mT.group()
p1_mT.colormap = 'Candy'
gl_mT.rescale()
gl_mT.xlim = (min(csv_file['Wavelength']), max(csv_file['Wavelength']))


wks_a = op.new_sheet(lname='Analytical method') # for normal plots
wks_a.name = 'aCD'

wks_a.from_list(0, Wavelength_new, lname='Wavelength', units='nm')
wks_a.from_list(1, aCD_list[0], lname='aCD', units='mdeg')

wks_a_CB = wks_a.get_book().add_sheet("aCB")
wks_a_CB.from_list(0, Wavelength_new, lname='Wavelength', units='nm')
wks_a_CB.from_list(1, aCB_list[0], lname='aCB', units='mdeg')

wks_a_Abs = wks_a.get_book().add_sheet("Abs")
wks_a_Abs.from_list(0, Wavelength_new, lname='Wavelength', units='nm')
wks_a_Abs.from_list(1, Abs_list[0], lname='Abs', units = 'au')

wks_a_LD = wks_a.get_book().add_sheet("aLD")
wks_a_LD.from_list(0, Wavelength_new, lname='Wavelength', units='nm')
wks_a_LD.from_list(1, aLD_list[0], lname='aLD')

wks_a_LB = wks_a.get_book().add_sheet("aLB")
wks_a_LB.from_list(0, Wavelength_new, lname='Wavelength', units='nm')
wks_a_LB.from_list(1, aLB_list[0], lname='aLB')

wks_a_g_factor = wks_a.get_book().add_sheet("aG-factor")
wks_a_g_factor.from_list(0, Wavelength_new, lname='Wavelength', units='nm')
wks_a_g_factor.from_list(1, ag_factor_list[0], lname='aG-factor', units='CD/Abs')

wks_a_LDp = wks_a.get_book().add_sheet("aLD'")
wks_a_LDp.from_list(0, Wavelength_new, lname='Wavelength', units='nm')
wks_a_LDp.from_list(1, aLDp_list[0], lname="aLD'")

wks_a_LBp = wks_a.get_book().add_sheet("aLB'")
wks_a_LBp.from_list(0, Wavelength_new, lname='Wavelength', units='nm')
wks_a_LBp.from_list(1, aLBp_list[0], lname="aLB'")

wks_a_T = wks_a.get_book().add_sheet("Ts")
wks_a_T.from_list(0, Wavelength_new, lname='Wavelength', units='nm')
wks_a_T.from_list(1, M00_avg, lname="Transmission_M00")
wks_a_T.from_list(2, M11_avg, lname="Transmission_M11")
wks_a_T.from_list(3, M22_avg, lname="Transmission_M22")
wks_a_T.from_list(4, M33_avg, lname="Transmission_M33")

if len(M00_list) > 1:
for i in range(1, len(M00_list)):
wks_a.from_list(i+1, aCD_list[i])
wks_a_CB.from_list(i+1, aCB_list[i])
wks_a_Abs.from_list(i+1, Abs_list[i])
wks_a_LD.from_list(i+1, aLD_list[i])
wks_a_LB.from_list(i+1, aLB_list[i])
wks_a_g_factor.from_list(i+1, ag_factor_list[i])
wks_a_LDp.from_list(i+1, aLDp_list[i])
wks_a_LBp.from_list(i+1, aLBp_list[i])


gr3 = op.new_graph(template=(new_path_template + '\\Origin templates\\test_plot9.otpu'), lname='Analytical method plot')

# gr3 = op.new_graph(template=op.path('e') + 'test_plot9.otpu')

poz = 0
neg = 0

for i in range(len(CD_list_new)):
for j in range(len(CD_list_new[i])):
if aCD_list[i][j] > 0:
poz = poz + 1
else:
neg = neg + 1

gl_aCD = gr3[0]
p1_aCD = gl_aCD.add_plot(f'{wks_a.lt_range()}!(?,1:end)')
p1_aCD.set_int('line.width', 3)
if (neg >= 1) and poz >=1:
gl_aCD.set_int("x.AtZero",1)
gl_aCD.group()
p1_aCD.colormap = 'Candy'
gl_aCD.rescale()
gl_aCD.xlim = (min(csv_file['Wavelength']), max(csv_file['Wavelength']))

poz = 0
neg = 0

for i in range(len(CD_list_new)):
for j in range(len(CD_list_new[i])):
if aCB_list[i][j] > 0:
poz = poz + 1
else:
neg = neg + 1

gl_aCB= gr3[1]
p1_aCB = gl_aCB.add_plot(f'{wks_a_CB.lt_range()}!(?,1:end)')
p1_aCB.set_int('line.width', 3)
if (neg >= 1) and poz >=1:
gl_aCB.set_int("x.AtZero",1)
gl_aCB.group()
p1_aCB.colormap = 'Candy'
gl_aCB.rescale()
gl_aCB.xlim = (min(csv_file['Wavelength']), max(csv_file['Wavelength']))

poz = 0
neg = 0

for i in range(len(CD_list_new)):
for j in range(len(CD_list_new[i])):
if Abs_list[i][j] > 0:
poz = poz + 1
else:
neg = neg + 1

gl_aAbs= gr3[2]
p1_aAbs = gl_aAbs.add_plot(f'{wks_a_Abs.lt_range()}!(?,1:end)')
p1_mAbs.set_int('line.width', 3)
if (neg >= 1) and poz >=1:
gl_aAbs.set_int("x.AtZero",1)
gl_aAbs.group()
p1_aAbs.colormap = 'Candy'
gl_aAbs.rescale()
gl_aAbs.xlim = (min(csv_file['Wavelength']), max(csv_file['Wavelength']))

poz = 0
neg = 0

for i in range(len(CD_list_new)):
for j in range(len(CD_list_new[i])):
if aLD_list[i][j] > 0:
poz = poz + 1
else:
neg = neg + 1

gl_aLD= gr3[3]
p1_aLD = gl_aLD.add_plot(f'{wks_a_LD.lt_range()}!(?,1:end)')
p1_aLD.set_int('line.width', 3)
if (neg >= 1) and poz >=1:
gl_aLD.set_int("x.AtZero",1)
gl_aLD.group()
p1_aLD.colormap = 'Candy'
gl_aLD.rescale()
gl_aLD.xlim = (min(csv_file['Wavelength']), max(csv_file['Wavelength']))

poz = 0
neg = 0

for i in range(len(CD_list_new)):
for j in range(len(CD_list_new[i])):
if aLB_list[i][j] > 0:
poz = poz + 1
else:
neg = neg + 1

gl_aLB= gr3[4]
p1_aLB = gl_aLB.add_plot(f'{wks_a_LB.lt_range()}!(?,1:end)')
p1_aLB.set_int('line.width', 3)
if (neg >= 1) and poz >=1:
gl_aLB.set_int("x.AtZero",1)
gl_aLB.group()
p1_aLB.colormap = 'Candy'
gl_aLB.rescale()
gl_aLB.xlim = (min(csv_file['Wavelength']), max(csv_file['Wavelength']))

poz = 0
neg = 0

for i in range(len(CD_list_new)):
for j in range(len(CD_list_new[i])):
if ag_factor_list[i][j] > 0:
poz = poz + 1
else:
neg = neg + 1

gl_ag= gr3[5]
p1_ag = gl_ag.add_plot(f'{wks_a_g_factor.lt_range()}!(?,1:end)')
p1_ag.set_int('line.width', 3)
if (neg >= 1) and poz >=1:
gl_ag.set_int("x.AtZero",1)
gl_ag.group()
p1_ag.colormap = 'Candy'
gl_ag.rescale()
gl_ag.xlim = (min(csv_file['Wavelength']), max(csv_file['Wavelength']))

poz = 0
neg = 0

for i in range(len(CD_list_new)):
for j in range(len(CD_list_new[i])):
if mLDp_list[i][j] > 0:
poz = poz + 1
else:
neg = neg + 1

gl_aLDp= gr3[6]
p1_aLDp = gl_aLDp.add_plot(f'{wks_a_LDp.lt_range()}!(?,1:end)')
p1_aLDp.set_int('line.width', 3)
if (neg >= 1) and poz >=1:
gl_aLDp.set_int("x.AtZero",1)
gl_aLDp.group()
p1_aLDp.colormap = 'Candy'
gl_aLDp.rescale()
gl_aLDp.xlim = (min(csv_file['Wavelength']), max(csv_file['Wavelength']))

poz = 0
neg = 0

for i in range(len(CD_list_new)):
for j in range(len(CD_list_new[i])):
if aLBp_list[i][j] > 0:
poz = poz + 1
else:
neg = neg + 1

gl_aLBp= gr3[7]
p1_aLBp = gl_aLBp.add_plot(f'{wks_a_LBp.lt_range()}!(?,1:end)')
p1_aLBp.set_int('line.width', 3)
if (neg >= 1) and poz >=1:
gl_aLBp.set_int("x.AtZero",1)
gl_aLBp.group()
p1_aLBp.colormap = 'Candy'
gl_aLBp.rescale()
gl_aLBp.xlim = (min(csv_file['Wavelength']), max(csv_file['Wavelength']))

poz = 0
neg = 0

for i in range(len(M00_avg)):
if M00_avg[i] > 0 or M11_avg[i] > 0 or M22_avg[i] > 0 or M33_avg[i] > 0:
poz = poz + 1
elif M00_avg[i] < 0 or M11_avg[i] < 0 or M22_avg[i] < 0 or M33_avg[i] < 0:
neg = neg + 1

gl_aT= gr3
p1_aT = gl_aT.add_plot(f'{wks_a_T.lt_range()}!(?,1:end)')
p1_aT.set_int('line.width', 3)
if (neg >= 1) and poz >=1:
gl_aT.set_int("x.AtZero",1)
gl_aT.group()
p1_aT.colormap = 'Candy'
gl_aT.rescale()
gl_aT.xlim = (min(csv_file['Wavelength']), max(csv_file['Wavelength']))


path, name = os.path.split(file)

new_path = path.replace('/', '\\')

if str(reply5) == 'Yes':
if os.path.isfile('' + path + "/" + str(session_name2) + str('_all_data_baseline_substracted') + '_' + str(name.split('.')[0]) + '.opju') == False:
op.save('' + new_path + "\\" + str(session_name2) + str('_all_data_baseline_substracted') + '_' + str(name.split('.')[0]) + '.opju')
else:
for i in range(1, 99):
if os.path.isfile('' + path + "/" + str(session_name2) + str('_all_data_baseline_substracted') + '_' + str(name.split('.')[0]) + 'Copy' + str(i) + '.opju') == False:
op.save('' + new_path + "\\" + str(session_name2) + str('_all_data_baseline_substracted') + '_' + str(name.split('.')[0]) + 'Copy' + str(i) + '.opju')
break
else:
if os.path.isfile('' + path + "/" + str(session_name2) + str('_all_data') + '_' + str(name.split('.')[0]) + '.opju') == False:
op.save('' + new_path + "\\" + str(session_name2) + str('_all_data') + '_' + str(name.split('.')[0]) + '.opju')
else:
for i in range(1, 99):
if os.path.isfile('' + path + "/" + str(session_name2) + str('_all_data') + '_' + str(name.split('.')[0]) + 'Copy' + str(i) + '.opju') == False:
op.save('' + new_path + "\\" + str(session_name2) + str('_all_data') + '_' + str(name.split('.')[0]) + 'Copy' + str(i) + '.opju')
break


if op.oext:
op.exit()


Here is the code. Sorry for the long code. Every sheet takes 289 spectra, so it makes in total about that number. I think the error comes after this line: gl_iLDp= gr1[6]

quote:

---

[i]Originally posted by cpyang

Origin worksheet has limitation of the number of columns, so you need to put 11560 or more columns into a worksheet, there might be memory issues.

Can you show the code snippet, like are you using wks.from_df, or wks.from_list?

CP

---

quote:

---

Originally posted by cpyang

Can you add some print to see where exactly it had trouble?
So each spectra is two column? How many sheets in the book?

If 289x2 is the total number of columns, then that should be no problem at all.

CP

---

quote:

---

Originally posted by cpyang

Can you find out exactly which line failed? Can you try inside Origin, does that also fail?

Can you also check Task Manager to see if computer memory is still plenty available?

Is there a way that you can prepare the package so we can run it to see the problem?

CP

---

quote:

---

Originally posted by cpyang

Can you find out exactly which line failed? Can you try inside Origin, does that also fail?

Can you also check Task Manager to see if computer memory is still plenty available?

Is there a way that you can prepare the package so we can run it to see the problem?

CP

---