Full extents - 4 layer 3 material mesh for Insalah Site
May 4 2016

Source layers are from corner point geometry in AVS file mesh.inp
Middle layer is 4m thick, top and bottom are 8m with aprox 500m lateral width. 

This hex mesh does not have any zero or negative volumes.
It is very thin with elements aprox 530 wide and  4 high


Here are the numbers for this mesh:

number of nodes =     1608010        
number of elements =  8852802

ATTRIBUTE        MIN               MAX         DIFFERENCE    LENGTH  
 xic     0.000000000E+00  2.000000000E+04 2.000000000E+04   1608010  
 yic     0.000000000E+00  2.650000000E+04 2.650000000E+04   1608010  
 zic    -1.540000000E+03 -1.240000000E+03 3.000000000E+02   1608010  

---------------------------------------                                         
element volumes b/w  0.6739E-01 and  0.1426E+01:     16451                      
element volumes b/w  0.1426E+01 and  0.3018E+02:    160195                      
element volumes b/w  0.3018E+02 and  0.6386E+03:     37336                      
element volumes b/w  0.6386E+03 and  0.1351E+05:   6814244                      
element volumes b/w  0.1351E+05 and  0.2860E+06:   1824576                      
min volume =   6.7392260E-02  max volume =   2.8597917E+05                      
-----------------------------------------------------------         

AMatbld3d_stor: npoints =  1608010  ncoefs =   21110542                         
AMatbld3d_stor: Number of unique coefs =   1551877                              
AMatbld3d_stor: Maximum num. connections to a node =  21  
AMatbld3d_stor: Volume min =    1.6562500E+03                                    
AMatbld3d_stor: Volume max =    4.3577781E+05                                    
AMatbld3d_stor: Total Volume:   1.5900002E+11 
AMatbld3d_stor: abs(Aij/xij) min =   0.0000000E+00                              
AMatbld3d_stor: abs(Aij/xij) max =   1.6652937E+03                              
AMatbld3d_stor: (Aij/xij) max =   0.0000000E+00                                 
AMatbld3d_stor: (Aij/xij) min =  -1.6652937E+03             

=====================================================================
Tet Mesh with 1 meter buffers
Full mesh 20000x26500, Subset is  14400x12190 meters
Refine horizontal 3 times (near 66m width)

Connect into tet mesh
  1) refine laterally until connect works without neg ccoefs
  2) create buffer version and connect
     move each layer 1m up and down, add flat top and bottom
     the voronoi interfaces will be between these 2m spaced nodes

Project directory /scratch/fwo/tam/CO2/insalah
mesh directory tet_buff1m_full/
FEHM directory fehm_tet_buff1m_full/

Mesh file: tet_buff1m_full/pyramid/tet_bufflayers.inp 

Material  1 has Lower    321602 nodes. #nodes/nnodes is   0.20
Material  2 has Middle   321602 nodes. #nodes/nnodes is   0.20
Material  3 has Upper    321602 nodes. #nodes/nnodes is   0.20
Material  4 has Top      321602 nodes. #nodes/nnodes is   0.20
Material  5 has Bottom   321602 nodes. #nodes/nnodes is   0.20

Mesh Files:

-rw-r--r-- 1 tamiller sft      477133651 May  4 08:26 tet.fehmn
-rw-r--r-- 1 tamiller sft      543690083 May  4 12:06 tet.stor
-rw-r--r-- 1 tamiller sft       17688287 May  4 08:26 tet_material.zone
-rw-r--r-- 1 tamiller sft       84788342 May  4 12:04 tet_multi_mat.zone
-rw-r--r-- 1 tamiller sft        3714319 May  4 12:04 tet_outside.zone
-rw-r--r-- 1 tamiller sft       21440921 May  4 12:04 tet_outside_vor.area