This example demonstrates how to set up Sedimentation_data to sub-divide a given depositional event into "n" smaller sub-events. This allows for thick sediment volumes to be deposited in several deposition sub-events which are later merged. This may be useful for example in an scenario in which a forward model is derived from a restoration analysis and the restoration model considered formation discretisation which may be too thick for forward modelling. The layer sub division algorithm works as follows:

 

1.The sedimentation volume for a given formation is divided into "n" user input sub-layers

2.The duration for each sub-layer is calculated as t / n where t is the total duration for formation deposition

3.At the end of the deposition of the sublayer n+1 it is merged with sublayer n into a single group. At the end of the deposition of the last sublayer there will be a single group for the whole formation.

 

For this example a simple 3D column geometry will be used. The simulation considers 3 stages as follows:

 

1.From 0 to 1 Ma: Initialisation of pre-existing sediment (Formation0)

2.From 1 to 2.5 Ma: Deposition of Formation1 (300m) sub-dividing it into 3 sub-layers

3.From 2.5 to 4.5 Ma: Deposition of Formation2 (200m) sub-dividing it into 2 sub-layers

 

 

The data file Ex_005_Case1.dat for the example is found in: ParaGeo Examples\General Examples\Ex_005\Data