#pragma rtGlobals=1		// Use modern global access method.

//E.g., ExtractOrgMatrices(HighSO4_conV, 1.4, 0.5)
Function ExtractOrgMatrices(wTodo, vOrg_RIE, vCE)
	wave wTodo
	variable vOrg_RIE, vCE
	svar ms0_sv_speclist=root:panel:ms0_sv_speclist
	ms0_sv_speclist = "Org"
	nvar ms0_ck_plotNew=root:panel:ms0_ck_plotNew
	ms0_ck_plotnew = 0
	nvar trunc2zero=root:panel:ms0_ck_trunc2zero
	trunc2zero = 0
	nvar ion2mass=root:panel:ms0_ck_convert2ug
	ion2mass = 1
	nvar ms0_ck_addControls=root:panel:ms0_ck_addControls
	ms0_ck_addControls = 0
	make/o/d/n=(0,0) organics_MS
	newdatafolder root:tempFolder
	Make/o/d/n=0 :tempFolder:t_series
	wave Myt_series = :tempFolder:t_series
	Make/o/d/n=0 :tempFolder:rn_series
	wave Myrn_series = :tempFolder:rn_series
	Make/o/d/n=0 :tempFolder:SO4
	wave MySO4 = :tempFolder:SO4
	Make/o/d/n=0 :tempFolder:NO3
	wave MyNO3 = :tempFolder:NO3
	Make/o/d/n=0 :tempFolder:Org
	wave MyOrg = :tempFolder:Org
	Make/o/d/n=0 :tempFolder:NH4
	wave MyNH4 = :tempFolder:NH4
	Make/o/d/n=0 :tempFolder:Chl
	wave MyChl = :tempFolder:Chl
	svar toDo_makeFromName_sv=root:panel:toDo_makeFromName_sv
	wave t_series = :index:t_series
	wave rn_series = :index:rn_series
	wave org, SO4, NO3, NH4, Chl
	variable i, j
	For (i=0; i<numpnts(wTodo); i+=1)
		svar toDo_makeFromRuns_sv=root:panel:toDo_makeFromRuns_sv
		toDo_makeFromRuns_sv = num2str(wTodo[i])
		toDo_makeFromName_sv = "ForOrgMat"
		sq_butt_ToDoFromRuns("ams_toDoNewToDo")
		sq_butt_MS_avg("ms0_avcalc") 
		variable pntnum = BinarySearch(:index:rn_series, wTodo[i])
		if (waveexists(MSSDiff_p_ForOrgMat_org))
			redimension/n=(i+1) :tempFolder:t_series
			Myt_series[i] = t_series[pntnum]
			redimension/n=(i+1) :tempFolder:rn_series
			Myrn_series[i] = rn_series[pntnum]
			redimension/n=(i+1) :tempFolder:SO4
			MySO4[i] = SO4[pntnum]
			redimension/n=(i+1) :tempFolder:NO3
			MyNO3[i] = NO3[pntnum]
			redimension/n=(i+1) :tempFolder:Org
			MyOrg[i] = Org[pntnum]
			redimension/n=(i+1) :tempFolder:NH4
			MyNH4[i] = NH4[pntnum]
			redimension/n=(i+1) :tempFolder:Chl
			MyChl[i] = Chl[pntnum]
			wave MSSDiff_p_ForOrgMat_org
			redimension/n=(i+1, numpnts(MSSDiff_p_ForOrgMat_org)) organics_MS 
			For (j=0; j<numpnts(MSSDiff_p_ForOrgMat_org); j+=1)
				organics_MS[i][j] = MSSDiff_p_ForOrgMat_org[j]/(vOrg_RIE*vCE)
			Endfor
		elseif (waveexists(MSSD_ForOrgMat_org))
			redimension/n=(i+1) :tempFolder:t_series
			Myt_series[i] = t_series[pntnum]
			redimension/n=(i+1) :tempFolder:rn_series
			Myrn_series[i] = rn_series[pntnum]
			redimension/n=(i+1) :tempFolder:SO4
			MySO4[i] = SO4[pntnum]
			redimension/n=(i+1) :tempFolder:NO3
			MyNO3[i] = NO3[pntnum]
			redimension/n=(i+1) :tempFolder:Org
			MyOrg[i] = Org[pntnum]
			redimension/n=(i+1) :tempFolder:NH4
			MyNH4[i] = NH4[pntnum]
			redimension/n=(i+1) :tempFolder:Chl
			MyChl[i] = Chl[pntnum]
			wave MSSD_ForOrgMat_org
			redimension/n=(i+1, numpnts(MSSD_ForOrgMat_org)) organics_MS 
			For (j=0; j<numpnts(MSSD_ForOrgMat_org); j+=1)
				organics_MS[i][j] = MSSD_ForOrgMat_org[j]/(vOrg_RIE*vCE)
			Endfor
		endif
	endfor
	make/o/d/n=(dimsize(organics_MS, 0), dimsize(organics_MS, 1)) organics_MS_err = 1
	Save/T/M="\r\n" root:TempFolder:t_series,root:TempFolder:rn_series, root:TempFolder:org, root:TempFolder:SO4, root:TempFolder:NO3, root:TempFolder:NH4, root:TempFolder:Chl, Organics_MS, Organics_MS_err as "WavesMatricesForOrganicAnalysis.itx"
	KillDatafolder/Z "root:TempFolder"
End