about OSL dating in OxCal

[Çiğdem Tepe]

Hi all;

I have 7 OSL samples taken from a trench which were excavated on a fault for a paloesismological study. How do I run these in the OxCal program. Can you show an example code? Because nothing I tried worked properly in the program. What should I do? 

I'll be happy if you can help me.

best wishes

çiğdem

[Çiğdem Tepe]

Hello All;

I have OSL ages taken from a trench. However, I could not figure out with which input code these ages should be entered into the Oxcal program.

It is mentioned in the manual that Oxcal is willing to treat OSL Dates as historical information. Is there a script I can use to input already calibrated OSL dates and not have Oxcal expand the uncertainty?

For example, how should input code that I need to enter into OxCal program for OSL dates in below be? I would be very appreciate if you can help?

best wishes

Sample1 (83,08 ± 4,16)   Sample 2 (60,02 ± 5,78)  

Sample 3(45,43 ± 2,31)   Sample 4(37,72 ± 6,25)

best wishes

Crocus Hill

[Erik Marsh]

Hi Crocus,

You can input these dates as normal distributions that are not calibrated with the radiocarbon curve. The code would be:

C_Date("Sample 1",83.08,4.16);

Also, you can use the AD format:

C_Date("Sample 1",AD(83.08),4.16);

Hope this helps

Erik

[Kate]

Hi Erik-

A variation on this question - is there a way to include non-normal distributions (i.e., multimodal pdfs) that are not calibrated? (Thinking a txt file with age and likelihood columns?).  

Thanks!
Kate

[Kate]

clarifying my question above -  the multimodal pdfs are not radiocarbon dates and do not need to be calibrated. 

kate

[Christopher Ramsey]

Yes - you can use either a Prior command or the P command - both at:

https://c14.arch.ox.ac.uk/oxcalhelp/hlp_analysis_inform.html#param

The Prior will probably be easier for what you want - and the format for these is given in:

https://c14.arch.ox.ac.uk/oxcalhelp/hlp_analysis_file.html#prior

It is very important to note that if being used for dates the first column should be in the fractional Gregorian years (effectively astronomical years for this purpose).

You might want to play around with using P and the inline array method - that is probably best if the number of points in the distribution is not too high. There you want the P command to be something like:

Date(P("Test",calBP(12000),calBP(10000),[0,0.1,0.3,0.6,0.8,1.0,0.8,0.7,0.3,0.2,0.1,0]));

You will get the idea of how this works from this example - and can modify it from there.

Best wishes

Christopher

> --
> You received this message because you are subscribed to the Google Groups "OxCal" group.
> To unsubscribe from this group and stop receiving emails from it, send an email to oxcal+un...@googlegroups.com.
> To view this discussion on the web visit https://groups.google.com/d/msgid/oxcal/16586d7c-2429-4451-b50c-50f4fe616c2fn%40googlegroups.com.

[Kate]

great- thanks! I remembered seeing it somewhere but then missed it on the re^n-read.

Appreciate it. 
Kate