* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * Program: LaGriT V3.101 Linux m64 * * * date_compile: 2013/04/25 RH gf4.5 * * * Run Time: 2013/Jun 19 14:29:02 * * * Manual: http://lagrit.lanl.gov * * * * * * * * * * * * * * * * * * * * * * * * * * * * * -----oOo----- LaGriT V3 LACC-2012-084 LaGriT Copyright: This program was prepared by Los Alamos National Security, LLC at Los Alamos National Laboratory (LANL) under contract No. DE-AC52-06NA25396 with the U.S. Department of Energy (DOE). All rights in the program are reserved by the DOE and Los Alamos National Security, LLC. Permission is granted to the public to copy and use this software without charge, provided that this Notice and any statement of authorship are reproduced on all copies. Neither the U.S. Government nor LANS makes any warranty, express or implied, or assumes any liability or responsibility for the use of this software. -----oOo----- Output log file: outx3dgen Command log file: logx3dgen # write FEHM files with imt set to zones read avs tet_final.inp mo geniee finish cmo/status/brief The current-mesh-object(CMO) is: mo 1 Mesh Object name: mo number of nodes = 634148 number of elements = 3686223 dimensions geometry = 3 element type = tet dimensions topology = 3 4 nodes 4 faces 6 edges boundary flag = 22000000 status = active finish quality epsilonl, epsilonaspect: 3.4111139E-10 3.9690693E-29 -------------------------------------------- elements with aspect ratio < .01: 0 elements with aspect ratio b/w .01 and .02: 10 elements with aspect ratio b/w .02 and .05: 8 elements with aspect ratio b/w .05 and .1 : 14 elements with aspect ratio b/w .1 and .2 : 124 elements with aspect ratio b/w .2 and .5 : 10923 elements with aspect ratio b/w .5 and 1. : 3675144 min aspect ratio = 0.1866E-01 max aspect ratio = 0.1000E+01 epsilonvol: 1.3322676E-04 --------------------------------------- element volumes b/w 0.3756E+01 and 0.1596E+02: 2922656 element volumes b/w 0.1596E+02 and 0.6784E+02: 349929 element volumes b/w 0.6784E+02 and 0.2883E+03: 107848 element volumes b/w 0.2883E+03 and 0.1225E+04: 122024 element volumes b/w 0.1225E+04 and 0.5208E+04: 183766 min volume = 3.7560096E+00 max volume = 5.2083333E+03 ----------------------------------------------------------- 3686223 total elements evaluated. cmo printatt mo -xyz- minmax ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH xic 0.000000000E+00 1.000000000E+03 1.000000000E+03 634148 yic 0.000000000E+00 1.000000000E+03 1.000000000E+03 634148 zic 0.000000000E+00 6.000000000E+02 6.000000000E+02 634148 cmo/select/mo # PLANES for fault # 800 to min top x 818.75 and 825 surface/sleft/intrface/plane/680.0 0. 0./800.0 1000.0 550.0/800.0 0.0 550.0 The unit vector for sleft is 0.9770158i -0.0000000j -0.2131671k # 700.0 0. 0. / 820.0 1000.0 550.0 / 820.0 0.0 550.0 surface/sright/intrface/plane/710.0 0. 0./830.0 1000.0 550.0/830.0 0.0 550.0 The unit vector for sright is 0.9770158i -0.0000000j -0.2131671k # BOX for repository middle at 300. (200x200x20) # choose height to include 20m at resolution 3.125m # 300 middle, z min 290.625 max 309.375 define XMIN 100. define XMAX 300. define YMIN 400. define YMAX 600. define ZMIN 287. define ZMAX 313. define ZMIN 290. define ZMAX 310. # create outside node attributes dump zone_outside tmp mo keepatt *********dump_outside_list******** Voronoi Areas used for outside faces. assign_quadrant_id epsilon for zero: 1.0000000E-15 cmo/addatt/mo/top/vint/scalar/nnodes/linear/permanent/afgx/-5.0/ finish Face top 1 has 5907 nodes. top Sum Voronoi Area_x Area_y Area_z 0.7500000E+05 0.1798949E+05 0.1000000E+07 cmo/addatt/mo/bottom/vint/scalar/nnodes/linear/permanent/afgx/-5.0/ finish Face bottom 2 has 4852 nodes. bottom Sum Voronoi Area_x Area_y Area_z 0.2500000E+05 0.2291748E+05 0.1000000E+07 cmo/addatt/mo/left_w/vint/scalar/nnodes/linear/permanent/afgx/-5.0/ finish Face left_w 3 has 1025 nodes. left_w Sum Voronoi Area_x Area_y Area_z 0.6000000E+06 0.1500000E+05 0.2500000E+05 cmo/addatt/mo/right_e/vint/scalar/nnodes/linear/permanent/afgx/-5.0/ finish Face right_e 5 has 943 nodes. right_e Sum Voronoi Area_x Area_y Area_z 0.5500000E+06 0.1375000E+05 0.2500000E+05 cmo/addatt/mo/back_n/vint/scalar/nnodes/linear/permanent/afgx/-5.0/ finish Face back_n 6 has 2952 nodes. back_n Sum Voronoi Area_x Area_y Area_z 0.1500000E+05 0.5750000E+06 0.2290771E+05 cmo/addatt/mo/front_s/vint/scalar/nnodes/linear/permanent/afgx/-5.0/ finish Face front_s 4 has 2952 nodes. front_s Sum Voronoi Area_x Area_y Area_z 0.1500000E+05 0.5750000E+06 0.2290771E+05 ############################################# # Set the zones for FEHM files cmo select mo # set all node imt and element itetclr to 5 cmo setatt mo imt 5 634148 values reset for attribute imt1 cmo setatt mo itetclr 5 3686223 values reset for attribute itetclr # set all top nodes to 1 and this will be overwritten by top zone 2 pset/pzone1/attribute top/1,0,0/eq 1 THE PSET pzone1 HAS 5907 POINTS eltset ezone1 inclusive pset,get,pzone1 THE ELTSET ezone1 HAS 34609 ELEMENTS cmo/setatt/mo/imt/pset get pzone1/1 5907 values reset for attribute imt1 cmo setatt mo itetclr eltset,get,ezone1 1 34609 values reset for attribute itetclr # set top section right to 2 pset/pzone2/attribute/xic/pset get pzone1/ge/700.0 THE PSET pzone2 HAS 4185 POINTS eltset ezone2 inclusive pset,get,pzone2 THE ELTSET ezone2 HAS 24778 ELEMENTS cmo setatt mo imt pset,get,pzone2 2 4185 values reset for attribute imt1 cmo setatt mo itetclr eltset,get,ezone2 2 24778 values reset for attribute itetclr # set repository box shape to 3 pset/pbox/geom/xyz/1,0,0/XMIN YMIN ZMIN/XMAX YMAX ZMAX THE PSET pbox HAS 29575 POINTS eltset/ebox/exclusive/pset,get,pbox THE ELTSET ebox HAS 146894 ELEMENTS cmo setatt mo imt pset,get,pbox 3 29575 values reset for attribute imt1 cmo setatt mo itetclr eltset,get,ebox 3 146894 values reset for attribute itetclr dump gmv tmp.gmv mo cmo/modatt/-def-/-def-/ioflag/x finish # set fault imt to 4 # set last so overwrite the nodes set to top region/rflt/ge sleft and le sright pset/pflt/region/rflt THE PSET pflt HAS 478130 POINTS eltset/eflt/exclusive/pset get pflt THE ELTSET eflt HAS 2458156 ELEMENTS cmo/setatt/mo/imt/pset get pflt/4 478130 values reset for attribute imt1 cmo/setatt/mo/itetclr/eltset get eflt/4 2458156 values reset for attribute itetclr resetpts itp geniee finish cmo setatt mo ipolydat no 1 values reset for attribute ipolydat dump gmv tet_final_zones.gmv mo cmo/modatt/-def-/-def-/ioflag/x finish dump avs tet_final_zones.inp mo cmo/modatt/-def-/-def-/ioflag/x finish # do some FEHM preparation and checking cmo/addatt/mo/matclr/vint/scalar/nelements////0 cmo/copyatt/mo mo/matclr itetclr 3686223 copied from mo itetclr to -> mo matclr cmo setatt mo itetclr 1 3686223 values reset for attribute itetclr filter 1,0,0 FILTER:Use internal epsilonl value = 0.341111394591E-09 FILTER:Dudding 0 duplicate points. rmpoint compress 0 points removed and 0 elements removed. The new point count is 634148. The new element count is 3686223. resetpts itp geniee finish quality epsilonl, epsilonaspect: 3.4111139E-10 3.9690693E-29 -------------------------------------------- elements with aspect ratio < .01: 0 elements with aspect ratio b/w .01 and .02: 10 elements with aspect ratio b/w .02 and .05: 8 elements with aspect ratio b/w .05 and .1 : 14 elements with aspect ratio b/w .1 and .2 : 124 elements with aspect ratio b/w .2 and .5 : 10923 elements with aspect ratio b/w .5 and 1. : 3675144 min aspect ratio = 0.1866E-01 max aspect ratio = 0.1000E+01 epsilonvol: 1.3322676E-04 --------------------------------------- element volumes b/w 0.3756E+01 and 0.1596E+02: 2922656 element volumes b/w 0.1596E+02 and 0.6784E+02: 349929 element volumes b/w 0.6784E+02 and 0.2883E+03: 107848 element volumes b/w 0.2883E+03 and 0.1225E+04: 122024 element volumes b/w 0.1225E+04 and 0.5208E+04: 183766 min volume = 3.7560096E+00 max volume = 5.2083333E+03 ----------------------------------------------------------- 3686223 total elements evaluated. * Must be single defined mesh cmo select mo resetpts/parent rmpoint compress 0 points removed and 0 elements removed. The new point count is 634148. The new element count is 3686223. cmo status mo brief The current-mesh-object(CMO) is: mo 1 Mesh Object name: mo number of nodes = 634148 number of elements = 3686223 dimensions geometry = 3 element type = tet dimensions topology = 3 4 nodes 4 faces 6 edges boundary flag = 22000000 status = active dump fehm tet mo *** Write FEHMN GEOM AND ZONE FILES *** *********dump_material_lists******** Minimum material ID value = 1 Maximum material ID value = 5 Total possible materials = 5 Material 1 has 1722 nodes. #nodes/nnodes is 0.271545443684E-02 Material 2 has 1296 nodes. #nodes/nnodes is 0.204368703999E-02 Material 3 has 29575 nodes. #nodes/nnodes is 0.466373786330E-01 Material 4 has 478130 nodes. #nodes/nnodes is 0.753972232342 Material 5 has 123425 nodes. #nodes/nnodes is 0.194631218910 *********dump_interface_list******** *********dump_multi_material_list******** Material 1 has 6795 multi-material connections. Material 2 has 6744 multi-material connections. Material 3 has 37394 multi-material connections. Material 4 has 287342 multi-material connections. Material 5 has 334265 multi-material connections. log/tty/off finish *********dump_outside_list******** Voronoi Areas used for outside faces. assign_quadrant_id epsilon for zero: 1.0000000E-15 Face top 1 has 5907 nodes. top Sum Voronoi Area_x Area_y Area_z 0.7500000E+05 0.1798949E+05 0.1000000E+07 Face bottom 2 has 4852 nodes. bottom Sum Voronoi Area_x Area_y Area_z 0.2500000E+05 0.2291748E+05 0.1000000E+07 Face left_w 3 has 1025 nodes. left_w Sum Voronoi Area_x Area_y Area_z 0.6000000E+06 0.1500000E+05 0.2500000E+05 Face right_e 5 has 943 nodes. right_e Sum Voronoi Area_x Area_y Area_z 0.5500000E+06 0.1375000E+05 0.2500000E+05 Face back_n 6 has 2952 nodes. back_n Sum Voronoi Area_x Area_y Area_z 0.1500000E+05 0.5750000E+06 0.2290771E+05 Face front_s 4 has 2952 nodes. front_s Sum Voronoi Area_x Area_y Area_z 0.1500000E+05 0.5750000E+06 0.2290771E+05 cmo/DELATT//top finish cmo/DELATT//bottom finish cmo/DELATT//left_w finish cmo/DELATT//right_e finish cmo/DELATT//back_n finish cmo/DELATT//front_s finish log/tty/on finish *********dump_parent_list******** *** 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 = 634148 ncoefs = 5033908 AMatbld3d_stor: Number of unique coefs = 610 AMatbld3d_stor: Maximum num. connections to a node = 20 AMatbld3d_stor: Volume min = 7.6293943E+00 AMatbld3d_stor: Volume max = 1.5625000E+04 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 dump gmv tet_final_ccoef.gmv mo cmo/modatt/-def-/-def-/ioflag/x finish quality epsilonl, epsilonaspect: 3.4111139E-10 3.9690693E-29 -------------------------------------------- elements with aspect ratio < .01: 0 elements with aspect ratio b/w .01 and .02: 10 elements with aspect ratio b/w .02 and .05: 8 elements with aspect ratio b/w .05 and .1 : 14 elements with aspect ratio b/w .1 and .2 : 124 elements with aspect ratio b/w .2 and .5 : 10923 elements with aspect ratio b/w .5 and 1. : 3675144 min aspect ratio = 0.1866E-01 max aspect ratio = 0.1000E+01 epsilonvol: 1.3322676E-04 --------------------------------------- element volumes b/w 0.3756E+01 and 0.1596E+02: 2922656 element volumes b/w 0.1596E+02 and 0.6784E+02: 349929 element volumes b/w 0.6784E+02 and 0.2883E+03: 107848 element volumes b/w 0.2883E+03 and 0.1225E+04: 122024 element volumes b/w 0.1225E+04 and 0.5208E+04: 183766 min volume = 3.7560096E+00 max volume = 5.2083333E+03 ----------------------------------------------------------- 3686223 total elements evaluated. cmo printatt mo -all- minmax ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH -def- 0.000000000E+00 0.000000000E+00 0.000000000E+00 634148 scalar 1 1 0 1 vector 3 3 0 1 nnodes 634148 634148 0 1 nedges 0 0 0 1 nfaces 0 0 0 1 nelements 3686223 3686223 0 1 mbndry 22000000 22000000 0 1 ndimensions_topo 3 3 0 1 ndimensions_geom 3 3 0 1 nodes_per_element 4 4 0 1 edges_per_element 6 6 0 1 faces_per_element 4 4 0 1 isetwd 0 11 11 634148 ialias 0 0 0 634148 imt1 1 5 4 634148 itp1 0 10 10 634148 icr1 0 0 0 634148 isn1 0 0 0 634148 xic 0.000000000E+00 1.000000000E+03 1.000000000E+03 634148 yic 0.000000000E+00 1.000000000E+03 1.000000000E+03 634148 zic 0.000000000E+00 6.000000000E+02 6.000000000E+02 634148 xtetwd 0 11 11 3686223 itetclr 1 1 0 3686223 itettyp 5 5 0 3686223 itetoff 0 14744888 14744888 3686223 jtetoff 0 14744888 14744888 3686223 itet 1 634148 634147 3686223x4 jtet 2 22000000 21999998 3686223x4 epsilon 1.000000004E-15 1.000000004E-15 0.000000000E+00 1 epsilonl 3.411113946E-10 3.411113946E-10 0.000000000E+00 1 epsilona 5.240252676E-07 5.240252676E-07 0.000000000E+00 1 epsilonv 1.332267630E-04 1.332267630E-04 0.000000000E+00 1 ipointi 1 1 0 1 ipointj 634148 634148 0 1 idebug 0 0 0 1 itypconv_sm 1 1 0 1 maxiter_sm 25 25 0 1 tolconv_sm 1.000000000E+00 1.000000000E+00 0.000000000E+00 1 nnfreq 1 1 0 1 ivoronoi 1 1 0 1 iopt2to2 2 2 0 1 xmin 0.000000000E+00 0.000000000E+00 0.000000000E+00 1 ymin 0.000000000E+00 0.000000000E+00 0.000000000E+00 1 zmin 0.000000000E+00 0.000000000E+00 0.000000000E+00 1 xmax 1.000000000E+03 1.000000000E+03 0.000000000E+00 1 ymax 1.000000000E+03 1.000000000E+03 0.000000000E+00 1 zmax 6.000000000E+02 6.000000000E+02 0.000000000E+00 1 kdtree_level 0 0 0 1 max_number_sets 64 64 0 1 number_of_psets 4 4 0 1 number_of_eltsets 0 0 0 1 number_of_fsets 0 0 0 1 matclr 1 5 4 3686223 finish LaGriT successfully completed