* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * Program: LaGriT V3.101 Linux m64 * * * date_compile: 2013/04/25 RH gf4.5 * * * Run Time: 2013/Dec 4 13:46:07 * * * 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 # read subset truncated at top and bottom read avs tet_trunc.inp cmotet cmo/addatt/cmotet/xval/VDOUBLE/scalar/nnodes/linear/permanent/gxaf/0.0 finish cmo/addatt/cmotet/nlayer/VDOUBLE/scalar/nnodes/linear/permanent/gxaf/0.0 finish cmo/addatt/cmotet/imat/VDOUBLE/scalar/nnodes/linear/permanent/gxaf/0.0 finish geniee finish cmo/status/brief The current-mesh-object(CMO) is: cmotet 1 Mesh Object name: cmotet number of nodes = 117450 number of elements = 662382 dimensions geometry = 3 element type = tet dimensions topology = 3 4 nodes 4 faces 6 edges boundary flag = 16000000 status = active finish resetpts itp geniee finish cmo printatt cmotet -all- minmax ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH -def- 0.000000000E+00 0.000000000E+00 0.000000000E+00 117450 scalar 1 1 0 1 vector 3 3 0 1 nnodes 117450 117450 0 1 nedges 0 0 0 1 nfaces 0 0 0 1 nelements 662382 662382 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 0 0 117450 ialias 0 0 0 117450 imt1 1 1 0 117450 itp1 0 12 12 117450 icr1 0 0 0 117450 isn1 0 0 0 117450 xic 4.201000000E+05 4.299000000E+05 9.800000000E+03 117450 yic 4.307100000E+06 4.316900000E+06 9.800000000E+03 117450 zic -2.208199951E+03 1.081999969E+02 2.316399948E+03 117450 xtetwd 0 0 0 662382 itetclr 1 9 8 662382 itettyp 5 5 0 662382 itetoff 0 2649524 2649524 662382 jtetoff 0 2649524 2649524 662382 itet 1 117450 117449 662382x4 jtet 1 18649416 18649415 662382x4 epsilon 1.000000004E-15 1.000000004E-15 0.000000000E+00 1 epsilonl 3.120068113E-09 3.120068113E-09 0.000000000E+00 1 epsilona 4.384175439E-05 4.384175439E-05 0.000000000E+00 1 epsilonv 4.939760845E-02 4.939760845E-02 0.000000000E+00 1 ipointi 1 1 0 1 ipointj 117450 117450 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 4.201000000E+05 4.201000000E+05 0.000000000E+00 1 ymin 4.307100000E+06 4.307100000E+06 0.000000000E+00 1 zmin -2.208199951E+03 -2.208199951E+03 0.000000000E+00 1 xmax 4.299000000E+05 4.299000000E+05 0.000000000E+00 1 ymax 4.316900000E+06 4.316900000E+06 0.000000000E+00 1 zmax 1.081999969E+02 1.081999969E+02 0.000000000E+00 1 kdtree_level 0 0 0 1 max_number_sets 64 64 0 1 number_of_psets 0 0 0 1 number_of_eltsets 0 0 0 1 number_of_fsets 0 0 0 1 xval 1.000000000E+00 9.000000000E+00 8.000000000E+00 117450 nlayer 1.000000000E+00 4.700000000E+01 4.600000000E+01 117450 imat 1.000000000E+00 9.000000000E+00 8.000000000E+00 117450 ********* Write FEHM setup files ******** cmo/addatt/cmotet/matsav/vint/scalar/nelements/linear/permanent/afgx/0. cmo/addatt/cmotet/cc_neg/vint/scalar/nnodes/linear/permanent/afgx/0. cmo copyatt/cmotet cmotet/matsav itetclr 662382 copied from cmotet itetclr to -> cmotet matsav cmo setatt cmotet itetclr 1 662382 values reset for attribute itetclr cmo copyatt/cmotet cmotet/imt imat Warning: Attribute type VDOUBLE written to VINT 117450 copied from cmotet imat to -> cmotet imt1 filter 1,0,0 FILTER:Use internal epsilonl value = 0.312006811336E-08 FILTER:Dudding 0 duplicate points. rmpoint compress 0 points removed and 0 elements removed. The new point count is 117450. The new element count is 662382. resetpts itp geniee finish quality epsilonl, epsilonaspect: 3.1200681E-09 3.0373317E-26 -------------------------------------------- elements with aspect ratio < .01: 0 elements with aspect ratio b/w .01 and .02: 0 elements with aspect ratio b/w .02 and .05: 0 elements with aspect ratio b/w .05 and .1 : 0 elements with aspect ratio b/w .1 and .2 : 131064 elements with aspect ratio b/w .2 and .5 : 321436 elements with aspect ratio b/w .5 and 1. : 209882 min aspect ratio = 0.1002E+00 max aspect ratio = 0.7849E+00 epsilonvol: 4.9397608E-02 --------------------------------------- element volumes b/w 0.1267E+06 and 0.1952E+06: 345744 element volumes b/w 0.1952E+06 and 0.3007E+06: 86436 element volumes b/w 0.3007E+06 and 0.4633E+06: 57330 element volumes b/w 0.4633E+06 and 0.7139E+06: 172578 element volumes b/w 0.7139E+06 and 0.1100E+07: 294 min volume = 1.2666667E+05 max volume = 1.1000000E+06 ----------------------------------------------------------- 662382 total elements evaluated. * Must be single defined mesh cmo select cmotet resetpts/parent rmpoint compress 0 points removed and 0 elements removed. The new point count is 117450. The new element count is 662382. cmo status cmotet The current-mesh-object(CMO) is: cmotet 1 Mesh Object name: cmotet number of nodes = 117450 number of elements = 662382 dimensions geometry = 3 element type = tet dimensions topology = 3 4 nodes 4 faces 6 edges boundary flag = 16000000 status = active NAME TYPE RANK LENGTH INTER PERSI IO VALUE 1 -def- VDOU scalar nnodes linea tempo agl 0.000E+00 2 scalar INT scalar scalar const perma l 1.000E+00 3 vector INT scalar scalar const perma l 3.000E+00 4 nnodes INT scalar scalar const perma l 1.175E+05 5 nedges INT scalar scalar const perma l 0.000E+00 6 nfaces INT scalar scalar const perma l 0.000E+00 7 nelements INT scalar scalar const perma l 6.624E+05 8 mbndry INT scalar scalar const perma l 1.600E+07 9 ndimensions_top INT scalar scalar const perma l 3.000E+00 10 ndimensions_geo INT scalar scalar const perma l 3.000E+00 11 nodes_per_eleme INT scalar scalar const perma l 4.000E+00 12 edges_per_eleme INT scalar scalar const perma l 6.000E+00 13 faces_per_eleme INT scalar scalar const perma l 4.000E+00 14 isetwd VDOU scalar nnodes or perma l 0.000E+00 15 ialias VINT scalar nnodes seque perma l 0.000E+00 16 imt1 VINT scalar nnodes max perma gal 0.000E+00 17 itp1 VINT scalar nnodes min perma gal 0.000E+00 18 icr1 VINT scalar nnodes min perma gal 0.000E+00 19 isn1 VINT scalar nnodes user perma gal 0.000E+00 20 xic VDOU scalar nnodes linea perma l 0.000E+00 21 yic VDOU scalar nnodes linea perma l 0.000E+00 22 zic VDOU scalar nnodes linea perma l 0.000E+00 23 xtetwd VDOU scalar nelements or perma l 0.000E+00 24 itetclr VINT scalar nelements user perma l 0.000E+00 25 itettyp VINT scalar nelements user perma l 0.000E+00 26 itetoff VINT scalar nelements user perma l 0.000E+00 27 jtetoff VINT scalar nelements user perma l 0.000E+00 28 itet VINT nodes_per_ nelements user perma l 0.000E+00 29 jtet VINT faces_per_ nelements user perma l 0.000E+00 30 ipolydat CHAR scalar scalar const perma l yes 31 vor2d CHAR scalar scalar const perma l yes 32 vor3d CHAR scalar scalar const perma l no 33 epsilon REAL scalar scalar const perma l 1.000E-15 34 epsilonl REAL scalar scalar const perma l 3.120E-09 35 epsilona REAL scalar scalar const perma l 4.384E-05 36 epsilonv REAL scalar scalar const perma l 4.940E-02 37 ipointi INT scalar scalar const perma l 1.000E+00 38 ipointj INT scalar scalar const perma l 1.175E+05 39 idebug INT scalar scalar const perma l 0.000E+00 40 itypconv_sm INT scalar scalar const perma l 1.000E+00 41 maxiter_sm INT scalar scalar const perma l 2.500E+01 42 tolconv_sm REAL scalar scalar const perma l 1.000E+00 43 nnfreq INT scalar scalar const perma l 1.000E+00 44 ivoronoi INT scalar scalar const perma l 1.000E+00 45 iopt2to2 INT scalar scalar const perma l 2.000E+00 46 dumptype CHAR scalar scalar const perma l binary 47 velname CHAR scalar scalar const perma l vels 48 densname CHAR scalar scalar const perma l ric 49 presname CHAR scalar scalar const perma l pic 50 enername CHAR scalar scalar const perma l eic 51 xmin REAL scalar scalar const perma l 4.201E+05 52 ymin REAL scalar scalar const perma l 4.307E+06 53 zmin REAL scalar scalar const perma l -2.208E+03 54 xmax REAL scalar scalar const perma l 4.299E+05 55 ymax REAL scalar scalar const perma l 4.317E+06 56 zmax REAL scalar scalar const perma l 1.082E+02 57 kdtree_level INT scalar scalar const perma l 0.000E+00 58 max_number_sets INT scalar scalar const perma l 6.400E+01 59 number_of_psets INT scalar scalar const perma l 0.000E+00 60 number_of_eltse INT scalar scalar const perma l 0.000E+00 61 psetnames VCHA scalar max_number_sets const perma l 62 eltsetnames VCHA scalar max_number_sets const perma l 63 geom_name CHAR scalar max_number_sets const perma l -defaultge 64 fsetnames VCHA scalar max_number_sets const perma l 65 number_of_fsets INT scalar scalar const perma l 0.000E+00 66 xval VDOU scalar nnodes linea perma gxa 0.000E+00 67 nlayer VDOU scalar nnodes linea perma gxa 0.000E+00 68 imat VDOU scalar nnodes linea perma gxa 0.000E+00 69 matsav VINT scalar nelements linea perma afg 0.000E+00 70 cc_neg VINT scalar nnodes linea perma afg 0.000E+00 cmo setatt cmotet ipolydat no 1 values reset for attribute ipolydat dump fehm tet cmotet ascii *** Write FEHMN GEOM AND ZONE FILES *** *********dump_material_lists******** Minimum material ID value = 1 Maximum material ID value = 9 Total possible materials = 9 Material 1 has 12500 nodes. #nodes/nnodes is 0.106428265572 Material 2 has 25000 nodes. #nodes/nnodes is 0.212856531143 Material 3 has 12500 nodes. #nodes/nnodes is 0.106428265572 Material 4 has 5000 nodes. #nodes/nnodes is 0.425713062286E-01 Material 5 has 7500 nodes. #nodes/nnodes is 0.638569593430E-01 Material 6 has 7500 nodes. #nodes/nnodes is 0.638569593430E-01 Material 7 has 25000 nodes. #nodes/nnodes is 0.212856531143 Material 8 has 15000 nodes. #nodes/nnodes is 0.127713918686 Material 9 has 7450 nodes. #nodes/nnodes is 0.634312480688E-01 *********dump_interface_list******** *********dump_multi_material_list******** Material 1 has 9801 multi-material connections. Material 2 has 19602 multi-material connections. Material 3 has 19602 multi-material connections. Material 4 has 19602 multi-material connections. Material 5 has 19602 multi-material connections. Material 6 has 19602 multi-material connections. Material 7 has 19602 multi-material connections. Material 8 has 19602 multi-material connections. Material 9 has 9801 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 cmo/addatt/cmotet/top/vint/scalar/nnodes/linear/permanent/afgx/-5.0/ finish Face top 1 has 2500 nodes. top Sum Voronoi Area_x Area_y Area_z 0.2673440E+07 0.9800000E+06 0.9604000E+08 cmo/addatt/cmotet/bottom/vint/scalar/nnodes/linear/permanent/afgx/-5.0/ finish Face bottom 2 has 2500 nodes. bottom Sum Voronoi Area_x Area_y Area_z 0.2673440E+07 0.9800000E+06 0.9604000E+08 cmo/addatt/cmotet/left_w/vint/scalar/nnodes/linear/permanent/afgx/-5.0/ finish Face left_w 3 has 2350 nodes. left_w Sum Voronoi Area_x Area_y Area_z 0.2097200E+08 0.4280000E+06 0.1960000E+07 cmo/addatt/cmotet/right_e/vint/scalar/nnodes/linear/permanent/afgx/-5.0/ finish Face right_e 5 has 2350 nodes. right_e Sum Voronoi Area_x Area_y Area_z 0.2097200E+08 0.4280000E+06 0.1960000E+07 cmo/addatt/cmotet/back_n/vint/scalar/nnodes/linear/permanent/afgx/-5.0/ finish Face back_n 6 has 2349 nodes. back_n Sum Voronoi Area_x Area_y Area_z 0.4625600E+06 0.2097200E+08 0.1960000E+07 cmo/addatt/cmotet/front_s/vint/scalar/nnodes/linear/permanent/afgx/-5.0/ finish Face front_s 4 has 2349 nodes. front_s Sum Voronoi Area_x Area_y Area_z 0.4625600E+06 0.2097200E+08 0.1960000E+07 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 1 AMatbld3d_stor: npoints = 117450 ncoefs = 1033054 AMatbld3d_stor: Number of unique coefs = 384 AMatbld3d_stor: Maximum num. connections to a node = 10 AMatbld3d_stor: Volume min = 1.8999575E+05 AMatbld3d_stor: Volume max = 5.1941486E+06 AMatbld3d_stor: Total Volume: 2.0552558E+11 AMatbld3d_stor: abs(Aij/xij) min = 1.0941969E-12 AMatbld3d_stor: abs(Aij/xij) max = 2.1023453E+03 AMatbld3d_stor: (Aij/xij) max = -1.0941969E-12 AMatbld3d_stor: (Aij/xij) min = -2.1023453E+03 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 cmo copyatt cmotet cmotet itetclr matsav 662382 copied from cmotet matsav to -> cmotet itetclr cmo/printatt/cmotet/ccoef/minmax cannot attribute ccoef in cmotet # IF ccoef exists # check to see if neg ccoefs are on outside of mesh # pset ncoef attribute ccoef 1,0,0 lt -0.001 # pset pout attribute itp 1,0,0 eq 10. # pset pin not pout # pset pBAD inter pin ncoef # cmo setatt cmotet cc_neg pset,get,ncoef 1 quality epsilonl, epsilonaspect: 3.1200681E-09 3.0373317E-26 -------------------------------------------- elements with aspect ratio < .01: 0 elements with aspect ratio b/w .01 and .02: 0 elements with aspect ratio b/w .02 and .05: 0 elements with aspect ratio b/w .05 and .1 : 0 elements with aspect ratio b/w .1 and .2 : 131064 elements with aspect ratio b/w .2 and .5 : 321436 elements with aspect ratio b/w .5 and 1. : 209882 min aspect ratio = 0.1002E+00 max aspect ratio = 0.7849E+00 epsilonvol: 4.9397608E-02 --------------------------------------- element volumes b/w 0.1267E+06 and 0.1952E+06: 345744 element volumes b/w 0.1952E+06 and 0.3007E+06: 86436 element volumes b/w 0.3007E+06 and 0.4633E+06: 57330 element volumes b/w 0.4633E+06 and 0.7139E+06: 172578 element volumes b/w 0.7139E+06 and 0.1100E+07: 294 min volume = 1.2666667E+05 max volume = 1.1000000E+06 ----------------------------------------------------------- 662382 total elements evaluated. dump gmv tet_fehm.gmv cmotet cmo/modatt/-def-/-def-/ioflag/x finish * translate to 0,0 for better voronoi volume precision trans 1,0,0 zero cmo/addatt/cmotet/transauto/VDOUBLE/scalar/vector/user/permanent/x/0.0 finish dump stor tet cmotet ascii *** Construct and Compress Sparse Matrix:3D *** *** Compress Area Coefficient Values *** AMatbld3d_stor: *****Negative Coefficients ****** AMatbld3d_stor: Total Number of Negative Coefficients 1 AMatbld3d_stor: Number of Significant Negative Coefs 0 AMatbld3d_stor: Number of 'zero' (< 10e-8 *max) coefs 1 cmo/addatt/cmotet/ccoef/VDOUBLE/scalar/nnodes/linear/permanent/afgx/0.0/ finish cmo/addatt/cmotet/ij_ccoef/VINT/scalar/nnodes/linear/permanent/afgx/0.0/ finish AMatbld3d_stor: npoints = 117450 ncoefs = 1033054 AMatbld3d_stor: Number of unique coefs = 384 AMatbld3d_stor: Maximum num. connections to a node = 10 AMatbld3d_stor: Volume min = 1.9000000E+05 AMatbld3d_stor: Volume max = 5.1941482E+06 AMatbld3d_stor: Total Volume: 2.0552560E+11 AMatbld3d_stor: abs(Aij/xij) min = 1.7116691E-15 AMatbld3d_stor: abs(Aij/xij) max = 2.1023453E+03 AMatbld3d_stor: (Aij/xij) max = 1.7116691E-15 AMatbld3d_stor: (Aij/xij) min = -2.1023453E+03 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 # Extract surface mesh for viewing voronoi polygons cmo select cmotet extract/surfmesh/1 0 0/cmos/cmotet/external cmo/addatt/cmos/itetclr0/vint/scalar/nelements////0 finish cmo/addatt/cmos/itetclr1/vint/scalar/nelements////0 finish cmo/addatt/cmos/idnode0/vint/scalar/nnodes////0 finish cmo/addatt/cmos/idelem0/vint/scalar/nelements////0 finish cmo/addatt/cmos/idelem1/vint/scalar/nelements////0 finish cmo/addatt/cmos/facecol/vint/scalar/nelements////0 finish cmo/addatt/cmos/idface0/vint/scalar/nelements////0 finish cmo/addatt/cmos/idface1/vint/scalar/nelements////0 finish cmo/addatt/cmos/xval/VDOUBLE/scalar/nnodes/linear/permanent/gxaf finish cmo/addatt/cmos/nlayer/VDOUBLE/scalar/nnodes/linear/permanent/gxaf finish cmo/addatt/cmos/imat/VDOUBLE/scalar/nnodes/linear/permanent/gxaf finish cmo/addatt/cmos/cc_neg/VINT/scalar/nnodes/linear/permanent/afgx finish cmo/addatt/cmos/ccoef/VDOUBLE/scalar/nnodes/linear/permanent/afgx finish cmo/addatt/cmos/ij_ccoef/VINT/scalar/nnodes/linear/permanent/afgx finish cmo/addatt/cmos/ncon50/int/scalar/scalar/constant/permanent/x/0 finish cmo/addatt/cmos/nconbnd/int/scalar/scalar/constant/permanent/x/0 finish cmo/addatt/cmos/icontab/vint/scalar/ncon50/constant/permanent/x/0 finish resetpts/itp geniee finish finish cmo copyatt/cmos cmos/itetclr itetclr1 27632 copied from cmos itetclr1 to -> cmos itetclr cmo setatt cmos ipolydat yes 1 values reset for attribute ipolydat dump gmv tet_surfmesh.gmv cmos cmo/modatt/-def-/-def-/ioflag/x finish finish LaGriT successfully completed