Mesh 3D octree with unstructured top conforming to slope.
This method uses hex with octree refine joined with triangulated slope piece.
  /scratch/sft/tam/grid_gen/generic_repository/3Dgrid_octree_coarse
  based on ../3Dgrid_octree
  directory images/gallery.html all images

Gridding Notes:
  The triangulated slope piece must be Delaunay (recon 1).
  Start refinement of fault wider than the fault, 
  otherwise the fault will stair step outside the smallest refine area.
  Repository box is centered at Z=300, select height based on resolution.
  Refined resolution is 3.25m so repository size is multiples of 3.25

The 25m octree mesh:
number of nodes    =     634148
number of elements =    3686223

This 50m octree mesh:
number of nodes    =  160583
number of elements =  921085

              MIN               MAX         DIFFERENCE    
 xic      0.000000000E+00  1.000000000E+03 1.000000000E+03 
 yic      0.000000000E+00  1.000000000E+03 1.000000000E+03
 zic      0.000000000E+00  6.000000000E+02 6.000000000E+02


FEHM Files:

-rw-r--r-- 1 tam sftus 49165932 Jun 19 16:47 tet.fehmn
-rw-r--r-- 1 tam sftus 34312186 Jun 19 16:47 tet.stor
-rw-r--r-- 1 tam sftus  1766590 Jun 19 16:47 tet_material.zone
-rw-r--r-- 1 tam sftus  4705462 Jun 19 16:47 tet_multi_mat.zone
-rw-r--r-- 1 tam sftus    94923 Jun 19 16:47 tet_outside.zone
-rw-r--r-- 1 tam sftus   547136 Jun 19 16:47 tet_outside_vor.area
-rw-r--r-- 1 tam sftus 38467720 Jun 19 16:47 tet_final_zones.gmv
-rw-r--r-- 1 tam sftus 91938827 Jun 19 16:47 tet_final_zones.inp

From file write_fehm.out

Material Top left     1 has     1417 nodes. #nodes/nnodes is   0.882409699261E-02
Material Top Right    2 has     1173 nodes. #nodes/nnodes is   0.730463396758E-02
Material Repository   3 has     5445 nodes. #nodes/nnodes is   0.339077003300E-01
Material Fault        4 has    46207 nodes. #nodes/nnodes is   0.287745267153
Material Background   5 has   106341 nodes. #nodes/nnodes is   0.662218272686


*** Construct and Compress Sparse Matrix:3D ***                                 
   *** Compress Area Coefficient Values ***                                     
AMatbld3d_stor: *****Zero Negative Coefficients ******                          
AMatbld3d_stor: Number of 'zero' (< 10e-8 *max) coefs         0                 
AMatbld3d_stor: npoints =   160583  ncoefs =    1281451                         
AMatbld3d_stor: Number of unique coefs =       632                              
AMatbld3d_stor: Maximum num. connections to a node =         23                 
AMatbld3d_stor: Volume min =   6.1035156E+01                                    
AMatbld3d_stor: Volume max =   1.2500000E+05                                    
AMatbld3d_stor: Total Volume:   5.7500000E+08                                   
AMatbld3d_stor: abs(Aij/xij) min =   0.0000000E+00                              
AMatbld3d_stor: abs(Aij/xij) max =   1.3333333E+02                              
AMatbld3d_stor: (Aij/xij) max =   0.0000000E+00                                 
AMatbld3d_stor: (Aij/xij) min =  -1.3333333E+02                                 
AMatbld3d_stor Matrix coefficient values stored as scalar area/distance         
AMatbld3d_stor Matrix compression used for graph and coefficient values         
ascii STOR file written with name tet.stor                                      
 
*** SPARSE COEFFICIENT MATRIX _astor SUCCESSFUL ***                             
3D Matrix Coefficient file written with name tet.stor        


Exodus Files:

-rw-r--r-- 1 tam sftus    36243 Jun 19 16:49 exo_1_top_left.faceset
-rw-r--r-- 1 tam sftus    49663 Jun 19 16:49 exo_2_top_right.faceset
-rw-r--r-- 1 tam sftus 26013176 Jun 19 16:49 tet_exodus.exo
-rw-r--r-- 1 tam sftus 36906997 Jun 19 16:49 tet_exodus_ascii.exo


number of dimension: 3 
number of nodes:     160583 
number of elements:  921085 

Material 1 background num_el_in_blk1 = 531503 ;
Material 2 repository num_el_in_blk2 = 365103 ;
Material 3 fault      num_el_in_blk3 = 24479 ;

Faceset  1 Top left   num_side_ss1 = 2780 ;
Faceset  2 Top right  num_side_ss2 = 3540 ;


Method ------------------------------------------------------------

See sequence of calls in file RUNS.

LaGriT input files *.lgi are commented with explanations.

Project Files -----------------------------------------------------

This is the 50m version of the 25m mesh in ../3Dgrid_octree
/scratch/sft/tam/grid_gen/generic_repository/3Dgrid_octree_coarse

This replaces 3Dgrid using stacked 2D triangulations.
/scratch/sft/spchu/reference_case/grid/gen_grid_v1_y25
/scratch/sft/spchu/reference_case/grid/gen_grid_v1_y25_new

Quan's 2D slices 
/home/quanb/src/grid_gen/generic_repository2D

Mesh 2D octree with unstructured top conforming to slope.
The mesh uses a combination of octree for repo and fault, and unstructured 
along top slope. This is similar to the 2D unstructured slope in 
/scratch/fwo/tam/scott/2D_proto/2D_vers02_unstructured
and /scratch/fwo/tam/scott/2D_proto/scripts/work_2d_unstructured