Case1b Base Case Grid defined at Cell centres

 

The data files for the project are in directory: SpatVar_001\Exercises\Case 1b\Data. Look for ! TODO  text in the datafile and follow the instructions.

 

Case1b example shows how to define the variables at grid cell centres as opposed to grid points. Same grid as Case1 will be used (note that there are slight differences in porosity distribution as in Case 1 grid points with Z=1000m were assigned a porosity of 0.2 whereas in this case porosity of 0.2 has been assigned at cells with Z=1250m) Grid data can be copied from excel file "SpatVar_001_case1b.xlsx".

 

Include

Data File

 

 

* Include

! ---------------------------------

 Filename "porosity_SpatVar_001_case1b.spat"          

 

1The included file name should point to the edited .spat file (ending with suffix _case1b).

 

 

 

Spatial_grid

Data File

 

 

* Spatial_grid    NUM=1

! ----------------------------------------

 Name        "Porosity_Spatial"

 Type                   "Grid1"

 Num_cells_x                 10

 Num_cells_y                  4

 Num_cells_z                  8

 Cell_division_X            800              

 Cell_division_Y           1000

 Cell_division_Z            500

 Grid_origin    IDM=3

   0.0   0.0  0.0  

 Element_variables   IDM=2

   "Porosity"  

  "P_strnv"

 Cell_values   IDM=2   JDM=320

  0.052528889   -0.552797207

  0.051586667   -0.553105649

     .

     .

     .

  0.353779067   -0.21182561

 

1Grid topology definition is identical as in case 1

2Either Element_variables or Cell_variables keywords can be used to define the variables that will be defined at grid cell centres.

3Either Element_values or Cell_values keywords can be used to define values for Element_variables.

 

 

 

 

Results

The result files for the project are in directory: SpatVar_001\Case1b\Results. Note that element values for a grid inputting values at cell centre are not interpolated (Element variables are assigned values of the closest cell).

 

SpatVar_001_Case1b_Results

 

Element porosity values (left). Nodal porosity values (right).