Some small updates

Mostly updated the default sample/standard/tail numbering.

Author: tiger-b

Front End/MIT_MasterUThDataReduction.m

@@ -34,7 +34,7 @@
 % ----------------------------------------------------------------------- %
 % AUTHORS: Christine Y Chen [CYC], David McGee [DM], Ben Hardt [BH]
 % Created:      04/30/15
-% Last updated: 09/09/15
+% Last updated: 02/15/24
 %
 %% ---------------------- REVISION HISTORY ------------------------------ %
 % 02/15/24   [BHT]     Removed signifiant figures output. Updated units
@@ -53,7 +53,7 @@
 %
 %% --------------------------------------------------------------------- %%
 function MIT_MasterUThDataReduction(filename_sample_info)
-versionNum = '1.1.090915';
+versionNum = '1.1.021524';
 
 fprintf('----------------- BEGIN U-TH DATA REDUCTION ----------------- \n')
 

Front End/MIT_readRawSummaryData.asv

@@ -122,8 +122,8 @@
         std1_U_tail_num(i) = samp_U_num(i) - 2;
         std2_U_num(i)      = samp_U_num(i) + 3;
         std2_U_tail_num(i) = samp_U_num(i) + 4;
-        std1_Th_num(i)     = samp_Th_num(i) - 1;
-        std2_Th_num(i)     = samp_Th_num(i) + 1;
+        std1_Th_num(i)     = samp_Th_num(i) - 2;
+        std2_Th_num(i)     = samp_Th_num(i) + 2;
 
         % If the user has defined values for analysis numbers for the
         % U tail/stds or Th stds, then set those values
@@ -426,7 +426,7 @@
     % U standard 1
     index1 = find(allUstd_analyNum == std1_U_num(i));
     if isempty(index1)
-        error('Error. Cannot find data for Th standard #1 in text file, based on analysis number.')
+        error('Error. Cannot find data for U standard #1 in text file, based on analysis number.')
     end
     Ustd1_data(i, :)  = allUstd_data(index1, :);  % Find U standard 1 data
     Ustd1_datetime(i) = allUstd_date(index1);     % Find the date of analysis
@@ -436,7 +436,7 @@
     % U standard 2
     index2 = find(allUstd_analyNum == std2_U_num(i));
     if isempty(index2)
-        error('Error. Cannot find data for Th standard #2 in text file, based on analysis number.')
+        error('Error. Cannot find data for U standard #2 in text file, based on analysis number.')
     end
     Ustd2_data(i, :) = allUstd_data(index2, :);   % Find U standard 2 data 
     Ustd2_datetime(i) = allUstd_date(index2);     % Find the date of analysis
@@ -557,7 +557,7 @@
     % U sample tail
     index = find(allUtail_analyNum == samp_U_tail_num(i));
     if isempty(index)
-        error('Error. Cannot find data for Th standard #1 in text file, based on analysis number.')
+        error('Error. Cannot find data for U tail in text file, based on analysis number.')
     end
     Usamp_tail_data(i, :)  = allUtail_data(index, :);
     Usamp_tail_datetime(i) = allUtail_date(index);   % Find analysis date
@@ -567,7 +567,7 @@
     % U standard 1 tail
     index1 = find(allUtail_analyNum == std1_U_tail_num(i));
     if isempty(index1)
-        error('Error. Cannot find data for Th standard #1 in text file, based on analysis number.')
+        error('Error. Cannot find data for U standard tail #1 in text file, based on analysis number.')
     end
     Ustd1_tail_data(i, :)  = allUtail_data(index1, :);
     Ustd1_tail_datetime(i) = allUtail_date(index1);   % Find analysis date
@@ -577,7 +577,7 @@
     % U standard 2 tail 
     index2 = find(allUtail_analyNum == std2_U_tail_num(i));
     if isempty(index2)
-        error('Error. Cannot find data for Th standard #2 in text file, based on analysis number.')
+        error('Error. Cannot find data for U standard tail #2 in text file, based on analysis number.')
     end
     Ustd2_tail_data(i, :) = allUtail_data(index2, :);
     Ustd2_tail_datetime(i) = allUtail_date(index2);   % Find analysis date

Front End/MIT_readRawSummaryData.m

@@ -0,0 +1,72 @@
+Date of reduction: 12-Mar-2024 16:29:47
+Version of script: 1.1.021524 
+
+FILENAMES OF RAW DATA FILES
+User input of sample data	mad_input_3.csv
+U samp file	U samp_no tail_meas.csv
+Th samp file	Th samp_meas.csv
+U std file	U std_no tail_meas.csv
+U tail file	U tail-237-236.5-234.5-233.5_meas.csv
+Th std file	Th std_meas.csv
+
+USER INPUT OF VALUES OF CONSTANTS
+230Th/232Th Tail Correction	6.3900e-07
+U239/U238 ratio	2.2900e-06
+Detrital Th230/Th232 ratio	4.400000e+00
+Detrital Th230/Th232 uncertainty 	2.2000e+00
+
+USER INPUT OF SAMPLE DATA
+Sample Name	Spike Type	Powder Weight (g)	Spike Weight (g)	U frac vol (mL)	Th frac vol (mL)	Date of Chem	U std conc	Th std conc
+Procedural blank	Aragonite		0.022790000000	1.00	0.30	1/1/2021	12.50	0.50
+AB11-2A (138mm)	Aragonite	0.052520000000	0.089840000000	1.00	0.30	1/1/2021	12.50	0.50
+AB11-2B (270mm)	Aragonite	0.050150000000	0.097750000000	1.00	0.30	1/1/2021	12.50	0.50
+AB13-2A (34mm)	Aragonite	0.050150000000	0.081240000000	2.00	0.30	1/1/2021	12.50	0.50
+AB13-2B (85.5mm)	Aragonite	0.058380000000	0.076400000000	2.00	0.30	1/1/2021	12.50	0.50
+AB13-2C (130mm)	Aragonite	0.050430000000	0.081990000000	2.00	0.30	1/1/2021	12.50	0.50
+AB13-2D (177mm)	Aragonite	0.057130000000	0.104260000000	2.00	0.30	1/1/2021	12.50	0.50
+AB13-2E (229mm)	Aragonite	0.048520000000	0.083110000000	2.00	0.30	1/1/2021	12.50	0.50
+AB13-2F (265mm)	Aragonite	0.054670000000	0.080880000000	2.00	0.30	1/1/2021	12.50	0.50
+CC5A	Aragonite	0.002460000000	0.099890000000	1.00	0.30	1/1/2021	12.50	0.50
+
+ANALYSIS NUMBERS AND DATE OF ANALYSIS
+Sample Name	U samp #	Date	Th samp #	Date	U samp tail #	Date	U std1 #	Date	U std1 tail #	Date	U std2 #	Date	U std2 tail #	Date	Th std1 #	Date	Th std2 #	Date
+Procedural blank	526	Analysed on Monday	592	Analysed on Wednesday	530	Analysed on Monday	523	Analysed on Monday	524	Analysed on Monday	529	Analysed on Monday	530	Analysed on Monday	590	Analysed on Wednesday	594	Analysed on Wednesday
+AB11-2A (138mm)	532	Analysed on Monday	596	Analysed on Wednesday	536	Analysed on Monday	529	Analysed on Monday	530	Analysed on Monday	535	Analysed on Monday	536	Analysed on Monday	594	Analysed on Wednesday	598	Analysed on Wednesday
+AB11-2B (270mm)	538	Analysed on Monday	600	Analysed on Wednesday	542	Analysed on Monday	535	Analysed on Monday	536	Analysed on Monday	541	Analysed on Monday	542	Analysed on Monday	598	Analysed on Wednesday	602	Analysed on Wednesday
+AB13-2A (34mm)	544	Analysed on Monday	604	Analysed on Wednesday	548	Analysed on Monday	541	Analysed on Monday	542	Analysed on Monday	547	Analysed on Monday	548	Analysed on Monday	602	Analysed on Wednesday	606	Analysed on Wednesday
+AB13-2B (85.5mm)	550	Analysed on Monday	608	Analysed on Wednesday	554	Analysed on Monday	547	Analysed on Monday	548	Analysed on Monday	553	Analysed on Monday	554	Analysed on Monday	606	Analysed on Wednesday	610	Analysed on Wednesday
+AB13-2C (130mm)	556	Analysed on Monday	612	Analysed on Wednesday	560	Analysed on Monday	553	Analysed on Monday	554	Analysed on Monday	559	Analysed on Monday	560	Analysed on Monday	610	Analysed on Wednesday	614	Analysed on Wednesday
+AB13-2D (177mm)	562	Analysed on Monday	616	Analysed on Wednesday	566	Analysed on Monday	559	Analysed on Monday	560	Analysed on Monday	565	Analysed on Monday	566	Analysed on Monday	614	Analysed on Wednesday	618	Analysed on Wednesday
+AB13-2E (229mm)	568	Analysed on Monday	620	Analysed on Wednesday	572	Analysed on Monday	565	Analysed on Monday	566	Analysed on Monday	571	Analysed on Monday	572	Analysed on Monday	618	Analysed on Wednesday	622	Analysed on Wednesday
+AB13-2F (265mm)	574	Analysed on Monday	624	Analysed on Wednesday	578	Analysed on Monday	571	Analysed on Monday	572	Analysed on Monday	577	Analysed on Monday	578	Analysed on Monday	622	Analysed on Wednesday	626	Analysed on Wednesday
+
+SAMPLE RAW INTENSITIES, IMPACT OF BLANK CORRECTION, STANDARD FLAGS, AND PERCENT YIELD
+	Raw Intensity	Raw Intensity	Raw Intensity	Raw Intensity	Raw Intensity	Raw Intensity	Blk Corr (%)	Blk Corr (%)	Blk Corr (%)	Blk Corr (%)	Flag	Flag	Flag	Flag	Flag	Flag	Yield	Yield	
+Sample Name	U238 (V)	U236 (mV)	Th230 (cps)	Th229 (mV)	Th230bg (cps)	234tail/234 (per mil)	238/236	238/234	230/229	232/229	Th mass bias	U mass bias	Th IC yield	U IC yield	Th low spike	U low spike	U236 (%)	Th229 (%)
+Procedural blank	0.003490000000	0.032500000000	7.290000000000	0.001240000000	1.630000000000	0.939745751944								4.4248	1.2400		91.91	45.13
+AB11-2A (138mm)	13.800000000000	93.900000000000	1290.000000000000	13.500000000000	1.270000000000	1.833818126161	0.02	1.09	0.78	0.75			2.3202				67.92	122.34
+AB11-2B (270mm)	44.800000000000	92.600000000000	5310.000000000000	15.000000000000	1.170000000000	0.594926884811	0.00	0.31	0.19	1.38				4.4053			59.57	121.83
+AB13-2A (34mm)	24.200000000000	44.800000000000	2320.000000000000	8.950000000000	2.050000000000	0.947405518786	0.01	0.33	0.31	0.87			2.3256		8.9500		68.53	86.30
+AB13-2B (85.5mm)	38.400000000000	40.400000000000	5040.000000000000	9.920000000000	1.030000000000	0.626788958837	0.00	0.20	0.17	1.09			2.3202		9.9200		65.19	101.48
+AB13-2C (130mm)	28.300000000000	48.000000000000	4660.000000000000	14.200000000000	1.280000000000	0.818757838571	0.00	0.30	0.24	9.78							71.61	134.47
+AB13-2D (177mm)	32.500000000000	62.600000000000	3280.000000000000	10.800000000000	1.580000000000	0.723317067944	0.00	0.27	0.20	16.67			2.3256				74.03	78.35
+AB13-2E (229mm)	39.600000000000	55.600000000000	6030.000000000000	13.500000000000	1.680000000000	0.657977910974	0.00	0.24	0.17	12.71			2.3202				82.48	121.56
+AB13-2F (265mm)	34.000000000000	53.300000000000	3800.000000000000	8.980000000000	0.950000000000	0.690203166576	0.00	0.28	0.19	2.36			2.3256		8.9800		80.61	82.56
+CC5A	27.200000000000	131.000000000000	13900.000000000000	11.700000000000	1.380000000000	0.551917022761	0.01	0.37	0.05	1.24							80.53	86.18
+
+PROCEDURAL BLANK MASSES
+230Th (fg)	232Th (pg)	234U (fg)	238U (pg)
+0.102785782285	1.154928566203	-12.849382564355	3.872415518239
+
+FINAL SAMPLE CONCENTRATIONS AND AGES
+Sample Name	238U	±(2sig)	232Th	±(2sig)	d234U	±(2sig)	(230Th/238U)	±(2sig)	230Th/232Th	±(2sig)	Age (uncorr)	±(2sig)	Age (corr)	±(2sig)	d234init	±(2sig)	Age (corr)	±(2sig)
+	ng/g		pg/g		per mille		activity		ppm atomic		yr		yr		per mille		yr BP	
+AB11-2A (138mm)	415.815590252216	8.316345437700	2947.265016028193	60.032205291056	14.744371077077	5.764751762493	0.038256749087	0.000395203738	85.697376343710	0.928732327451	4190.228458754335	50.425754248846	3978.981722376706	108.240323810428	14.910867587862	5.829850320203	3907.981722376706	108.240323810428
+AB11-2B (270mm)	1560.215632358844	31.207055155869	1668.211101741723	35.381198838656	8.148656583258	2.048585635387	0.044650497839	0.000213818694	663.037065906486	5.485480576451	4939.494266389836	26.312428813714	4907.443581747300	19.036948775464	8.262293233618	2.077154123172	4836.443581747300	26.312428813714
+AB13-2A (34mm)	1446.771902539385	28.937580033713	2662.154153466994	54.609918502921	11.943643091294	2.255686203595	0.028976962394	0.000235052088	250.033158592556	2.269538810002	3167.909875824807	27.047515974276	3112.953079576941	28.380761746096	12.049029533201	2.275589777849	3041.953079576941	28.380761746096
+AB13-2B (85.5mm)	2056.240050764869	41.138622880899	1811.661436302398	37.633796186283	11.525573346476	1.654636328815	0.032563093319	0.000203365131	586.813791236596	4.775137782224	3567.923647079982	23.421399098620	3541.602205118457	14.431757591214	11.641344605631	1.671256769343	3470.602205118457	23.421399098620
+AB13-2C (130mm)	1586.473200663519	31.734273398731	214.077214510091	12.369800871351	12.785390879016	2.068114738545	0.033654962981	0.000180795376	3959.946297519251	215.562774951635	3684.907024128820	21.551734969783	3680.881227688973	7.942905933703	12.918893486522	2.089709617369	3609.881227688973	21.551734969783
+AB13-2D (177mm)	1566.481507135785	31.331878517427	102.416706158462	10.447107981641	10.952235733609	1.987891688051	0.035619072909	0.000214779006	8649.977216846106	866.611291503199	3911.083992778754	25.267629813824	3909.129890110423	7.928908153488	11.073727525960	2.009943139219	3838.129890110423	25.267629813824
+AB13-2E (229mm)	2017.076531120247	40.346838729769	165.587625518018	12.507546717109	9.622056178380	1.874872263148	0.038167872188	0.000205586016	7381.936945085129	538.943612259973	4202.053166764154	24.422259722981	4199.596292102374	8.090887817203	9.736770548547	1.897224545462	4128.596292102374	24.422259722981
+AB13-2F (265mm)	1558.615266274400	31.177371636032	885.284593054883	20.726452922484	8.844301442828	2.005110751332	0.040580464651	0.000248654785	1134.360039982433	15.275836661865	4476.711437797218	29.453827582442	4459.698230380256	12.460769236902	8.956315136405	2.030505643034	4388.698230380256	29.453827582442
+CC5A	13967.303946012087	279.350179059111	37939.220128165129	795.863263914517	934.907028560216	2.270485133737	0.348351112809	0.001190701969	2036.194551695861	14.026369855286	21306.915922617529	83.884811294486	21265.717813496216	34.407358215472	992.731498238354	2.412842127493	21194.717813496216	83.884811294486

UTh-reduced-data-12-Mar-2024.xls

@@ -0,0 +1,11 @@
+Sample Name,Spike used,Sample Wt (g),Spike Wt (g),U fraction vol (mL),Th fraction vol (mL),Date of chemistry,U samp #,Th samp #,U samp tail #,U std1 #,U std1 tail #,U std2 #,U std2 tail #,Th std 1 #,Th std 2 #,U std conc (pg/g),Th std conc (ng/g), ,Enter Names or Values Here,
+Procedural blank,Aragonite,0,0.02279,1,0.3,1/1/2021,526,592,,,,,,,,12.5,0.5,Filename for U samp,U samp_no tail_meas.csv,
+AB11-2A (138mm),Aragonite,0.05252,0.08984,1,0.3,1/1/2021,532,596,,,,,,,,12.5,0.5,Filename for Th samp,Th samp_meas.csv,
+AB11-2B (270mm),Aragonite,0.05015,0.09775,1,0.3,1/1/2021,538,600,,,,,,,,12.5,0.5,Filename for U standards,U std_no tail_meas.csv,
+AB13-2A (34mm),Aragonite,0.05015,0.08124,2,0.3,1/1/2021,544,604,,,,,,,,12.5,0.5,Filename for U tail,U tail-237-236.5-234.5-233.5_meas.csv,
+AB13-2B (85.5mm),Aragonite,0.05838,0.0764,2,0.3,1/1/2021,550,608,,,,,,,,12.5,0.5,Filename for Th standards,Th std_meas.csv,
+AB13-2C (130mm),Aragonite,0.05043,0.08199,2,0.3,1/1/2021,556,612,,,,,,,,12.5,0.5,230Th/232Th Tail Correction,6.39E-07,
+AB13-2D (177mm),Aragonite,0.05713,0.10426,2,0.3,1/1/2021,562,616,,,,,,,,12.5,0.5,U239/U238 ratio,2.29E-06,
+AB13-2E (229mm),Aragonite,0.04852,0.08311,2,0.3,1/1/2021,568,620,,,,,,,,12.5,0.5,Detrital Th230/Th232 ratio,4.4,
+AB13-2F (265mm),Aragonite,0.05467,0.08088,2,0.3,1/1/2021,574,624,,,,,,,,12.5,0.5,Detrital Th230/Th232 uncertainty ,2.2,
+CC5A,Aragonite,0.00246,0.09989,1,0.3,1/1/2021,580,628,,,,,,,,12.5,0.5,Counts Per Second Per Volt,62415093,