* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * Program: LaGriT V3.101 Linux m64 * * * date_compile: 2013/04/25 RH gf4.5 * * * Run Time: 2013/Jun 19 16:46:52 * * * 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 = 160583 number of elements = 921085 dimensions geometry = 3 element type = tet dimensions topology = 3 4 nodes 4 faces 6 edges boundary flag = 16000000 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: 2 elements with aspect ratio b/w .02 and .05: 6 elements with aspect ratio b/w .05 and .1 : 0 elements with aspect ratio b/w .1 and .2 : 165 elements with aspect ratio b/w .2 and .5 : 4485 elements with aspect ratio b/w .5 and 1. : 916427 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.2106E+02: 8 element volumes b/w 0.2106E+02 and 0.1181E+03: 743938 element volumes b/w 0.1181E+03 and 0.6624E+03: 106458 element volumes b/w 0.6624E+03 and 0.3715E+04: 42964 element volumes b/w 0.3715E+04 and 0.2083E+05: 27717 min volume = 3.7560096E+00 max volume = 2.0833333E+04 ----------------------------------------------------------- 921085 total elements evaluated. cmo printatt mo -xyz- minmax ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH xic 0.000000000E+00 1.000000000E+03 1.000000000E+03 160583 yic 0.000000000E+00 1.000000000E+03 1.000000000E+03 160583 zic 0.000000000E+00 6.000000000E+02 6.000000000E+02 160583 # PLANES for fault cmo/select/mo 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 surface/sright/intrface/plane/700.0 0. 0./820.0 1000.0 550.0/820.0 0.0 550.0 The unit vector for sright is 0.9770158i -0.0000000j -0.2131671k # BOX for repository middle at 300. (200x200x20) # select based on spacing of 6.25m # repo needs 12.5 on either side of 300m define XMIN 100. define XMAX 300. define YMIN 400. define YMAX 600. define ZMIN 287. define ZMAX 313. # select set of repo zones # height is 12.50 with Z min 300. max 312.25 # 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 3234 nodes. top Sum Voronoi Area_x Area_y Area_z 0.1000000E+06 0.2574584E+05 0.1000000E+07 cmo/addatt/mo/bottom/vint/scalar/nnodes/linear/permanent/afgx/-5.0/ finish Face bottom 2 has 2273 nodes. bottom Sum Voronoi Area_x Area_y Area_z 0.5000000E+05 0.4159180E+05 0.1000000E+07 cmo/addatt/mo/left_w/vint/scalar/nnodes/linear/permanent/afgx/-5.0/ finish Face left_w 3 has 273 nodes. left_w Sum Voronoi Area_x Area_y Area_z 0.6000000E+06 0.3000000E+05 0.5000000E+05 cmo/addatt/mo/right_e/vint/scalar/nnodes/linear/permanent/afgx/-5.0/ finish Face right_e 5 has 252 nodes. right_e Sum Voronoi Area_x Area_y Area_z 0.5500000E+06 0.2750000E+05 0.5000000E+05 cmo/addatt/mo/back_n/vint/scalar/nnodes/linear/permanent/afgx/-5.0/ finish Face back_n 6 has 1287 nodes. back_n Sum Voronoi Area_x Area_y Area_z 0.2937500E+05 0.5750000E+06 0.3612305E+05 cmo/addatt/mo/front_s/vint/scalar/nnodes/linear/permanent/afgx/-5.0/ finish Face front_s 4 has 1287 nodes. front_s Sum Voronoi Area_x Area_y Area_z 0.2937500E+05 0.5750000E+06 0.3612305E+05 ############################################# # Set the zones for FEHM files cmo select mo # set all node imt and element itetclr to 5 cmo setatt mo imt 5 160583 values reset for attribute imt1 cmo setatt mo itetclr 5 921085 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 3234 POINTS eltset ezone1 inclusive pset,get,pzone1 THE ELTSET ezone1 HAS 18782 ELEMENTS cmo/setatt/mo/imt/pset get pzone1/1 3234 values reset for attribute imt1 cmo setatt mo itetclr eltset,get,ezone1 1 18782 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 1817 POINTS eltset ezone2 inclusive pset,get,pzone2 THE ELTSET ezone2 HAS 10674 ELEMENTS cmo setatt mo imt pset,get,pzone2 2 1817 values reset for attribute imt1 cmo setatt mo itetclr eltset,get,ezone2 2 10674 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 5445 POINTS eltset/ebox/exclusive/pset,get,pbox THE ELTSET ebox HAS 24479 ELEMENTS cmo setatt mo imt pset,get,pbox 3 5445 values reset for attribute imt1 cmo setatt mo itetclr eltset,get,ebox 3 24479 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 46207 POINTS eltset/eflt/inclusive/pset get pflt THE ELTSET eflt HAS 365103 ELEMENTS cmo/setatt/mo/imt/pset get pflt/4 46207 values reset for attribute imt1 cmo/setatt/mo/itetclr/eltset get eflt/4 365103 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 921085 copied from mo itetclr to -> mo matclr cmo setatt mo itetclr 1 921085 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 160583. The new element count is 921085. 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: 2 elements with aspect ratio b/w .02 and .05: 6 elements with aspect ratio b/w .05 and .1 : 0 elements with aspect ratio b/w .1 and .2 : 165 elements with aspect ratio b/w .2 and .5 : 4485 elements with aspect ratio b/w .5 and 1. : 916427 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.2106E+02: 8 element volumes b/w 0.2106E+02 and 0.1181E+03: 743938 element volumes b/w 0.1181E+03 and 0.6624E+03: 106458 element volumes b/w 0.6624E+03 and 0.3715E+04: 42964 element volumes b/w 0.3715E+04 and 0.2083E+05: 27717 min volume = 3.7560096E+00 max volume = 2.0833333E+04 ----------------------------------------------------------- 921085 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 160583. The new element count is 921085. cmo status mo brief The current-mesh-object(CMO) is: mo 1 Mesh Object name: mo number of nodes = 160583 number of elements = 921085 dimensions geometry = 3 element type = tet dimensions topology = 3 4 nodes 4 faces 6 edges boundary flag = 16000000 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 1417 nodes. #nodes/nnodes is 0.882409699261E-02 Material 2 has 1173 nodes. #nodes/nnodes is 0.730463396758E-02 Material 3 has 5445 nodes. #nodes/nnodes is 0.339077003300E-01 Material 4 has 46207 nodes. #nodes/nnodes is 0.287745267153 Material 5 has 106341 nodes. #nodes/nnodes is 0.662218272686 *********dump_interface_list******** *********dump_multi_material_list******** Material 1 has 5622 multi-material connections. Material 2 has 5289 multi-material connections. Material 3 has 10944 multi-material connections. Material 4 has 126058 multi-material connections. Material 5 has 146159 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 3234 nodes. top Sum Voronoi Area_x Area_y Area_z 0.1000000E+06 0.2574584E+05 0.1000000E+07 Face bottom 2 has 2273 nodes. bottom Sum Voronoi Area_x Area_y Area_z 0.5000000E+05 0.4159180E+05 0.1000000E+07 Face left_w 3 has 273 nodes. left_w Sum Voronoi Area_x Area_y Area_z 0.6000000E+06 0.3000000E+05 0.5000000E+05 Face right_e 5 has 252 nodes. right_e Sum Voronoi Area_x Area_y Area_z 0.5500000E+06 0.2750000E+05 0.5000000E+05 Face back_n 6 has 1287 nodes. back_n Sum Voronoi Area_x Area_y Area_z 0.2937500E+05 0.5750000E+06 0.3612305E+05 Face front_s 4 has 1287 nodes. front_s Sum Voronoi Area_x Area_y Area_z 0.2937500E+05 0.5750000E+06 0.3612305E+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 = 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 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: 2 elements with aspect ratio b/w .02 and .05: 6 elements with aspect ratio b/w .05 and .1 : 0 elements with aspect ratio b/w .1 and .2 : 165 elements with aspect ratio b/w .2 and .5 : 4485 elements with aspect ratio b/w .5 and 1. : 916427 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.2106E+02: 8 element volumes b/w 0.2106E+02 and 0.1181E+03: 743938 element volumes b/w 0.1181E+03 and 0.6624E+03: 106458 element volumes b/w 0.6624E+03 and 0.3715E+04: 42964 element volumes b/w 0.3715E+04 and 0.2083E+05: 27717 min volume = 3.7560096E+00 max volume = 2.0833333E+04 ----------------------------------------------------------- 921085 total elements evaluated. cmo printatt mo -all- minmax ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH -def- 0.000000000E+00 0.000000000E+00 0.000000000E+00 160583 scalar 1 1 0 1 vector 3 3 0 1 nnodes 160583 160583 0 1 nedges 0 0 0 1 nfaces 0 0 0 1 nelements 921085 921085 0 1 mbndry 16000000 16000000 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 160583 ialias 0 0 0 160583 imt1 1 5 4 160583 itp1 0 10 10 160583 icr1 0 0 0 160583 isn1 0 0 0 160583 xic 0.000000000E+00 1.000000000E+03 1.000000000E+03 160583 yic 0.000000000E+00 1.000000000E+03 1.000000000E+03 160583 zic 0.000000000E+00 6.000000000E+02 6.000000000E+02 160583 xtetwd 0 11 11 921085 itetclr 1 1 0 921085 itettyp 5 5 0 921085 itetoff 0 3684336 3684336 921085 jtetoff 0 3684336 3684336 921085 itet 1 160583 160582 921085x4 jtet 1 16000000 15999999 921085x4 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 160583 160583 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 921085 finish LaGriT successfully completed