Plot axes scales and page breaks

50 views
Skip to first unread message

Peter Almond

unread,
Feb 28, 2022, 6:04:35 PM2/28/22
to OxCal
Kia ora koutou

I am trying to get a consistent axis amongst multiple plots generated by quite a big model. I initially defined the x axis range for the plot
  Axis( BC(500), AD(2000));

but I found some plots had different axis ranges. I also wanted to put in page breaks, so I inserted these commands at the appropriate places,
 Page( );
   Axis( BC(500), AD(2000));  .

However, I still end up with different axis ranges and one page break doesn't work. I inserted the Page and axis commands at the end of one sequence and before another as so:
  ....
Before("Big matai germination on Blackburn south",C_date(2001)-U(251,271));
  };
  Page( );
  Axis( BC(500), AD(2000));

  Sequence("Littleman")
  {
   Boundary ("LME_VI")
.....

but the page break doesn't happen. I'd be grateful for some advice on the axis  and the page break problems. Thanks in anticipation.

Peter Almond


The full model specification is below.

 Plot("Littleman" )
 {
  Label("Both fans");
  Curve("SHCal20","shcal20.14c");
  Axis( BC(500), AD(2000));
  Sequence( "Blackburn")
  {
   Phase( "SSU S3")
   {
    First( "oldest OM in SSU S3");
    R_Date( "E8582/101-105T", 1562, 45)+U( 200, 240);
    R_Date( "Blackburn 3", 1291, 47);
    R_Date( "PT (a)", 1142, 45);
    R_Date( "PT3", 1119, 47);
    R_Date( "PT (b)", 1184, 55);
    R_Date("MLG_S2_4", 1429, 39);
    R_Date("MLG_S2_5", 1580, 37);
    R_Date("MLG_S2_3", 1208, 39);
    R_Date("Bb/Wh61/7Bb", 1241, 49);
    R_Date("Wh 61, Bb/Wh61/7Bc", 1758, 47);
    R_Date("GG_S1_1", 1002, 40);
    R_Date("N1_GG_S1_1", 1223, 40);
    Last( "Burial of S3");
   };
   Boundary( "BBE_III")
   {
    color="red";
   };
   Page( );
   Axis( BC(500), AD(2000));
   Phase( "Sediment package L3")
   {
    R_Date( "Wh14/27a", 996, 55);
    R_Date( "Wh 65- PT4 ", 967, 50);
    R_Date("S_SECT18_1", 948, 35);
   };
   //SSU S4 soil residence time based on Wh61 PS2
   Before("BBE_II minus S4 soil residence time",BBE_II-N(240,24));
   Phase( "SSU S4")
   {
    First( "begin_S4");
    R_Date( "Wh64, OT (b)", 774, 45);
    R_Date( "Wh64, OT (a)", 686, 50);
    R_Date( "OT (f)", 648, 50);
    R_Date( "OT (d)", 673, 45);
    R_Date( "OT(e)", 881, 50);
    Last( "end_S4");
   };
   Boundary( "BBE_II")
   {
    color="red";
   };
   Page( );
   Axis( BC(500), AD(2000));
   Before("CE2001 - Blackburn north soil residence time (Wh61 and Wh 67)",C_date(2001)-N(566,17));
   After("Difference in soil residence time of surface soils on Blackburn N and S",C_date(2001)-(N(566,17)-N(199,20)));
   Boundary( "BBE_I")
   {
    color="red";
   };
   Before("Surface soil residence time at Blackburn south based on Wh73",C_date(2001)-N(199,20));
   Before("Big matai germination on Blackburn south",C_date(2001)-U(251,271));
  };
  Page( );
  Axis( BC(500), AD(2000));
  Sequence("Littleman")
  {
   Boundary ("LME_VI")
   {
    color="red";
   };
   //Residence time of SSU S1 from mean of 3 estimates +/- s.e.
   Before(LME_V-N(746,57));
   Combine( "SSU S1 burial")
   {
    R_Date( "Wh57_Kelly soil", 1633, 58);
    R_Date( "Wh41, Kelly soil", 1689, 36);
    R_Date( "Wh55", 1696, 54);
   };
   Boundary( "LME_V")
   {
    color="red";
   };
   //SSU S2 soil residence time based on Wh57b
   Before(LME_IV-N(691,70));
   Phase( "SSU S2 at Wh50")
   {
    R_Date( "S_Wh50a_2", 967, 35);
    R_Date( "Wh50a.2", 1129, 55);
    R_Date( "Wh50a-Sb", 1057, 50);
   };
   Boundary( "LME_IV")
   {
    color="red";
   };
   Page( );
  Axis( BC(500), AD(2000));
   Phase( "Germination of Wh40 and Wh45")
   {
    First( "minimum age of start of SSU S3 soil formation");
    // tree in growth position rooted in Sm
    Combine( "Germination Wh40")
    {
     R_Date( "Wh40-40T", 797, 40)+U( -185, -184);
     //offset includes 28 y establishment lag
     R_Date( "Wh40-174T", 955, 40)+U( -51, -50);
     //offset includes 28 y establishment lag
    };
   };
   Combine( "Wh40_tree death & SSU S3 burial")
   {
    R_Date( "Wh40-174T", 955, 40)+U( 174, 175);
    R_Date( "Wh40-40T", 797, 40)+U( 40, 41);
    R_Date( "Wh45", 643, 43);
    R_Date( "Wh45-Sm1", 677, 50);
   };
   Boundary( "LME_III")
   {
    color="red";
   };
   Page( );
  Axis( BC(500), AD(2000));
   Phase( "Later trees growth period - germination to death. SSU S4 formation period.")
   {
    First( "establishment of later trees");
    //Germination of later trees - 28 years establishment time has been included
    R_Date( "Wh 50b-0T", 685, 60)+U( -33,-28);
    R_Date( "Wh50d-0T", 433, 65)+U( -33,-28);
    R_Date( "Wh50c - 43T", 328, 60)+U( -76,-71);
    R_Date( "Stuart's folly", 422, 42)+U( -198,-178);
    //Death of trees
    R_Date( "Wh 50b-0T", 685, 60)+U( 270, 275);
    R_Date( "Wh50d-0T", 433, 65)+U( 140, 141);
    R_Date( "Stuart's folly", 422, 42);
    R_Date( "Wh46", 427, 47);
    //death of prone tree in St
    R_Date( "Wh50c - 43T", 328, 60)+U( 133, 134);
    Last( "Death of later trees");
   };
   Boundary( "LME_II")
   {
    color="red";
   };
   Page( );
  Axis( BC(500), AD(2000));
   R_Date( "Wh59", 316, 46);
   Phase( "SSU S5")
   {
    R_Date( "St19", 267, 49);
    R_Date( "Wh47", 244, 45);
    R_Date("Wh50a.1",203,46 );
    R_Date("Wh50a-St",260,55 );
    R_Date( "Wh70", 226, 34);
   };
   Boundary("LME_I")
   {
    color="red";
   };
   //CE2001-Wh57c residence time
   Before(C_Date("Wh57c soil formation start",1706,30));
   Page( );
  };
  Difference("V-VI interval", "LME_V", "LME_VI");
  Difference("IV-V interval", "LME_IV", "LME_V");
  Difference("III-IV interval", "LME_III", "LME_IV");
  Difference("II-III interval", "LME_II", "LME_III");
  Difference("I-II interval", "LME_I", "LME_II");
  Difference("I-II interval", "BBE_I", "BBE_II");
  Difference("test", "BBE_I", "LME_II");
  Difference("LME_IV-BBE_III", "LME_IV", "BBE_III");
  Difference("LME_III-BBE_II", "LME_III", "BBE_II");
 };

Christopher Ramsey

unread,
Mar 1, 2022, 5:43:56 AM3/1/22
to OxCal group
The Page and Axis commands are interpreted by the graphing interface in a way which can cause issues.

In general the Page command works after a normal plottable item (like a date or a boundary). The Axis command is best used when the interface already has digested something from the page and knows what it is plotting.

It can be useful to shift the Page statements after plottable items - this can be a label if you wish and put the Axis statement after something already there in the page.

Your model makes this difficult because of all the 'Before' statements.

In such cases you may have to rely on the the output interface which allows you to select items after which you want a page break manually. But you will have to do that after the run. The following works by adding in a dummy label:

Plot("Littleman" )
{
Label("Both fans");
Curve("SHCal20","shcal20.14c");
Axis( BC(500), AD(2000));
Sequence( "Blackburn")
{
Phase( "SSU S3")
{
First( "oldest OM in SSU S3");
R_Date( "E8582/101-105T", 1562, 45)+U( 200, 240);
R_Date( "Blackburn 3", 1291, 47);
R_Date( "PT (a)", 1142, 45);
R_Date( "PT3", 1119, 47);
R_Date( "PT (b)", 1184, 55);
R_Date("MLG_S2_4", 1429, 39);
R_Date("MLG_S2_5", 1580, 37);
R_Date("MLG_S2_3", 1208, 39);
R_Date("Bb/Wh61/7Bb", 1241, 49);
R_Date("Wh 61, Bb/Wh61/7Bc", 1758, 47);
R_Date("GG_S1_1", 1002, 40);
R_Date("N1_GG_S1_1", 1223, 40);
Last( "Burial of S3");
};
Boundary( "BBE_III")
{
color="red";
};
Page( );
Phase( "Sediment package L3")
{
Axis( BC(500), AD(2000));
R_Date( "Wh14/27a", 996, 55);
R_Date( "Wh 65- PT4 ", 967, 50);
R_Date("S_SECT18_1", 948, 35);
};
//SSU S4 soil residence time based on Wh61 PS2
Before("BBE_II minus S4 soil residence time",BBE_II-N(240,24));
Phase( "SSU S4")
{
First( "begin_S4");
R_Date( "Wh64, OT (b)", 774, 45);
R_Date( "Wh64, OT (a)", 686, 50);
R_Date( "OT (f)", 648, 50);
R_Date( "OT (d)", 673, 45);
R_Date( "OT(e)", 881, 50);
Last( "end_S4");
};
Boundary( "BBE_II")
{
color="red";
};
Page( );
Before("CE2001 - Blackburn north soil residence time (Wh61 and Wh 67)",C_date(2001)-N(566,17));
Axis( BC(500), AD(2000));
After("Difference in soil residence time of surface soils on Blackburn N and S",C_date(2001)-(N(566,17)-N(199,20)));
Boundary( "BBE_I")
{
color="red";
};
Before("Surface soil residence time at Blackburn south based on Wh73",C_date(2001)-N(199,20));
Before("Big matai germination on Blackburn south",C_date(2001)-U(251,271));
Label(" ");
Page( );
};
Sequence("Littleman")
{
Axis( BC(500), AD(2000));
Boundary ("LME_VI")
{
color="red";
};
//Residence time of SSU S1 from mean of 3 estimates +/- s.e.
Before(LME_V-N(746,57));
Combine( "SSU S1 burial")
{
R_Date( "Wh57_Kelly soil", 1633, 58);
R_Date( "Wh41, Kelly soil", 1689, 36);
R_Date( "Wh55", 1696, 54);
};
Boundary( "LME_V")
{
color="red";
};
//SSU S2 soil residence time based on Wh57b
Before(LME_IV-N(691,70));
Phase( "SSU S2 at Wh50")
{
R_Date( "S_Wh50a_2", 967, 35);
R_Date( "Wh50a.2", 1129, 55);
R_Date( "Wh50a-Sb", 1057, 50);
};
Boundary( "LME_IV")
{
color="red";
};
Page( );
Phase( "Germination of Wh40 and Wh45")
{
First( "minimum age of start of SSU S3 soil formation");
Axis( BC(500), AD(2000));
// tree in growth position rooted in Sm
Combine( "Germination Wh40")
{
R_Date( "Wh40-40T", 797, 40)+U( -185, -184);
//offset includes 28 y establishment lag
R_Date( "Wh40-174T", 955, 40)+U( -51, -50);
//offset includes 28 y establishment lag
};
};
Combine( "Wh40_tree death & SSU S3 burial")
{
R_Date( "Wh40-174T", 955, 40)+U( 174, 175);
R_Date( "Wh40-40T", 797, 40)+U( 40, 41);
R_Date( "Wh45", 643, 43);
R_Date( "Wh45-Sm1", 677, 50);
};
Boundary( "LME_III")
{
color="red";
};
Page( );
Phase( "Later trees growth period - germination to death. SSU S4 formation period.")
{
Axis( BC(500), AD(2000));
First( "establishment of later trees");
//Germination of later trees - 28 years establishment time has been included
R_Date( "Wh 50b-0T", 685, 60)+U( -33,-28);
R_Date( "Wh50d-0T", 433, 65)+U( -33,-28);
R_Date( "Wh50c - 43T", 328, 60)+U( -76,-71);
R_Date( "Stuart's folly", 422, 42)+U( -198,-178);
//Death of trees
R_Date( "Wh 50b-0T", 685, 60)+U( 270, 275);
R_Date( "Wh50d-0T", 433, 65)+U( 140, 141);
R_Date( "Stuart's folly", 422, 42);
R_Date( "Wh46", 427, 47);
//death of prone tree in St
R_Date( "Wh50c - 43T", 328, 60)+U( 133, 134);
Last( "Death of later trees");
};
Boundary( "LME_II")
{
color="red";
};
Page( );
R_Date( "Wh59", 316, 46);
Axis( BC(500), AD(2000));
Phase( "SSU S5")
{
R_Date( "St19", 267, 49);
R_Date( "Wh47", 244, 45);
R_Date("Wh50a.1",203,46 );
R_Date("Wh50a-St",260,55 );
R_Date( "Wh70", 226, 34);
};
Boundary("LME_I")
{
color="red";
};
//CE2001-Wh57c residence time
Before(C_Date("Wh57c soil formation start",1706,30));
Page( );
};
Difference("V-VI interval", "LME_V", "LME_VI");
Difference("IV-V interval", "LME_IV", "LME_V");
Difference("III-IV interval", "LME_III", "LME_IV");
Difference("II-III interval", "LME_II", "LME_III");
Difference("I-II interval", "LME_I", "LME_II");
Difference("I-II interval", "BBE_I", "BBE_II");
Difference("test", "BBE_I", "LME_II");
Difference("LME_IV-BBE_III", "LME_IV", "BBE_III");
Difference("LME_III-BBE_II", "LME_III", "BBE_II");
};


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/4a2957da-7fb6-4000-b736-3f5d86f2dd65n%40googlegroups.com.

Reply all
Reply to author
Forward
0 new messages