Hi Samantha,
The function `prepsources()` job is to take single measurements through time and to aggregate them for each location.
We developed that to handle weather station data.
In your case, you should skip this function altogether and directly create a data frame or tibble similar to what this function usually produces.
Here is an example of output for terrestrial data:
source_ID mean_source_value var_source_value n_source_value lat long elev
1 BAD SALZUFLEN -49.57500 123.61929 8 52.10 8.75 135
2 BERLIN -50.63125 132.25067 8 52.46 13.40 48
3 BRAUNSCHWEIG -36.10000 85.51143 8 52.29 10.44 81
4 CUXHAVEN -45.90000 156.42000 8 53.87 8.70 5
5 EMMERICH -42.67500 208.13071 8 51.83 6.25 43
6 GARMISCH-PARTENKIRCHEN -61.80000 304.84000 8 47.48 11.06 719
So as you can see, it should not be difficult to generate yourself.
The first column needs some kind of ID that would be unique for each location.
Then comes mean_source_value, var_source_value and n_source_value, which indicate the mean, variance and number of point for a given isotope at the corresponding location.
Finally, lat, long, elev are needed for fitting the model, for you elev should be the depth.
Just make sure the column names are correct.
For your other issues, I already replied directly from my private email box where you sent more questions.
Best,
Alex