Next step in the workflow is to obtain a Geometry_set which includes all the node boundaries for the initial formation in the forward simulation (final restored formation). This will be done by running another boundary extraction datafile. The data files for such operation are in: Geol_002b\Case01\05_Boundary_extract2\Data. It consists of:
1.Geol_002b_Case01_rest_formation01.geo geometry file for the initial formation. This was output at the end of the restoration simulation.
2.Geol_002b_Case01_boun_ext2.dat datafile defined to perform the identification and extraction of the boundary data. Note that this file contains different data and used different keywords than the data file used in workflow step 3 where the target was to generate data for couple freedoms.
Geol_002b_Case01_boun_ext.dat
This file is defined to identify and extract boundary nodes data. This requires:
1.Read .geo file with the geometry data. 2.Definition of the stratigraphy. 3.Usage of Util_create_hex_boundary data structure to generate the data.
Stratigraphy data
Data File
|
|
* Stratigraphy_definition
Unit_Names IDM=1
"formation01"
Basal_horizon "Basal_surface"
* Stratigraphy_horizon NUM=1
Name "Basal_surface"
Geometry_set "Basal_surface"
* Stratigraphy_horizon NUM=2
Name "formation01"
Geometry_set "formation01_1"
|
•The full stratigraphy is defined by the basal boundary and the top surface for formation01.
|
Util_create_hex_boundary
Data File
|
|
* Util_create_hex_boundary
File_name "Geol_002b_Case01_fsim_sides"
Boundary_type 1
Geometry_set_name "Sides_all"
|
•Util_create_hex_boundary data structure is used to identify and extract data for the boundary nodes.
•File_name keyword is used to define the output file name for the generated data.
•Boundary_type is set to 1 in order to use facet normals to find all nodes in the boundaries. If not set, the default would only output boundary nodes linked to stratigraphic surfaces (only valid for single element across the thickness/z direction). In this case it would make no difference as formation01 only has 1 element in the thickness direction.
•Geometry_set_name is used to define the name of the output geometry set which will contain all the nodes in the boundary. This has been set to "Sides_all" to avoid conflict with the geometry set named "Sides" carried from restoration.
|
|
Results
The results for the present workflow steps are provided in Geol_002b\Case01\05_Boundary_extract2\Results
The key generated datafile is Geol_002b_Case01_fsim_sides.dat which contains the required Geometry_set named "Sides_all" containing all boundary nodes.
Data File
|
|
(...)
* Geometry_set NUM=27
Name "Sides"
Set_type "Nodal"
Nodes IDM=126
1 2 ... 497
(...)
* Geometry_set NUM=154
Name "Sides_all"
Set_type "Nodal"
Nodes IDM=252
1 2 ... 994
|
•A Geometry_set named "Sides_all" is generated. Such geometry set contains all the nodes in the side boundaries
•Note that another Geometry_set named "Sides" is present as it was carried from restoration which includes half of the nodes (it does not include the side boundaries for the base) |
|
