87Sr/86Sr with Trace Element LASS

Ivan

Hello,

While attempting to process both trace elements and Sr-isotopes in the same experiment, the initial processing of the 87Sr86Sr works but after processing the TE DRS, it provides me with this error:

There was an error in your python runDRS function: : operands could not be broadcast together with shapes (28223,) (62723,) File "C:/Users/Steve French/Documents/iolite/Plugins/Data reduction schemes/Sr_isotopes_new.py", line 97, in runDRS PFract = (np.log(8.37520938 / (SrCaAr88/SrCaAr86))) / (np.log(87.9056/85.9093))*mask

1) I think what is happening is both DRS's are searching for a "Mass 88" and once the TE DRS uses it for Sr88 the "shape" becomes different sizes? I am unsure exactly. What modification could I make?

I thought maybe I could modify this portion of the 87Sr86Sr DRS?
##
SrCaAr88 = data.timeSeriesList(data.Intermediate, {'Mass': '88'})[0].data()
SrCaAr86 = data.timeSeriesList(data.Intermediate, {'Mass': '86'})[0].data()

Sr86channel = data.timeSeriesList(data.Input, {'Mass': '86'})[0]

SrRb87 = data.timeSeriesList(data.Intermediate, {'Mass': '87'})[0].data()
Rb85 = data.timeSeriesList(data.Intermediate, {'Mass': '85'})[0].data()
SrCaAr84 = data.timeSeriesList(data.Intermediate, {'Mass': '84'})[0].data()
CaAr83 = data.timeSeriesList(data.Intermediate, {'Mass': '83'})[0].data()
CaAr82 = data.timeSeriesList(data.Intermediate, {'Mass': '82'})[0].data()

Let me know if you have any work arounds you can think of?

Ivan Edgeworth

Bence

Hi @Ivan ,

I think you're right: it's finding the wrong Sr88, and that's why there's a different number of points for the different mass specs.

In the metadata for the 'right' Sr88, can you see if the machine name is recorded? You could perhaps then perhaps change
SrCaAr88 = data.timeSeriesList(data.Intermediate, {'Mass': '88'})[0].data()

SrCaAr88 = data.timeSeriesList(data.Intermediate, {'Mass': '88', 'Machine name': '{Your machine Name}' })[0].data()

Just make sure you get the name of the machine name property correct (i.e. 'Machine name' or 'Machine Name').

Please let us know how you go.