* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * Program: LaGriT V2.202 Linux m64 * * * date_compile: 2011/01/05 gfort 4.5 * * * Run Time: 2011/Feb 22 11:43:59 * * * Manual: http://lagrit.lanl.gov * * * * * * * * * * * * * * * * * * * * * * * * * * * * * -----oOo----- 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 hex points and write fehm files * identify node column along well path read avs brick_evclr.inp cmohex cmo/addatt/cmohex/in_obj/VINT/scalar/nelements/linear/permanent/gxaf/0.0 finish cmo/addatt/cmohex/itetpar/VINT/scalar/nelements/linear/permanent/gxaf/0.0 finish cmo/addatt/cmohex/itetkid/VINT/scalar/nelements/linear/permanent/gxaf/0.0 finish cmo/addatt/cmohex/itetlev/VINT/scalar/nelements/linear/permanent/gxaf/0.0 finish geniee finish cmo/status/brief The current-mesh-object(CMO) is: cmohex 1 Mesh Object name: cmohex number of nodes = 97156 number of elements = 87768 dimensions geometry = 3 element type = hex dimensions topology = 3 8 nodes 6 faces 12 edges boundary flag = 16000000 status = active finish cmo printatt cmohex -xyz- minmax ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH xic 1.448530000E+05 1.497170000E+05 4.864000000E+03 97156 yic 1.426160000E+05 1.518320000E+05 9.216000000E+03 97156 zic -3.840000000E+03 -4.480000000E+02 3.392000000E+03 97156 cmo create cmotet///tet copypts cmotet cmohex copypts added nodes 1 to 97156 filter 1,0,0 FILTER:Use internal epsilonl value = 0.243337806231E-08 FILTER:Dudding 0 duplicate points. rmpoint compress 0 points removed and 0 elements removed. The new point count is 97156. The new element count is 0. cmo setatt cmotet itetclr 1 Warning: nelements = 0 0 length attribute: itetclr cmo setatt cmotet imt 1 Warning: nelements = 0 97156 values reset for attribute imt1 connect noadd Coordinates of enclosing tetrahedron are: 0.14120D+06 0.13570D+06 -0.10624D+05 0.14120D+06 0.18178D+06 -0.21440D+04 0.14120D+06 0.13570D+06 0.63360D+04 0.16552D+06 0.14722D+06 -0.21440D+04 ntets= 29346 node index= 5000 nlstfail= 0 %done= 5.03 ntets= 59547 node index= 10000 nlstfail= 0 %done= 10.21 ntets= 90134 node index= 15000 nlstfail= 0 %done= 15.46 ntets= 124614 node index= 20000 nlstfail= 0 %done= 21.38 ntets= 157306 node index= 25000 nlstfail= 0 %done= 26.98 ntets= 186746 node index= 30000 nlstfail= 0 %done= 32.03 ntets= 216202 node index= 35000 nlstfail= 0 %done= 37.09 ntets= 245698 node index= 40000 nlstfail= 0 %done= 42.15 ntets= 274823 node index= 45000 nlstfail= 0 %done= 47.14 ntets= 304366 node index= 50000 nlstfail= 0 %done= 52.21 ntets= 333851 node index= 55000 nlstfail= 0 %done= 57.27 ntets= 363068 node index= 60000 nlstfail= 0 %done= 62.28 ntets= 392232 node index= 65000 nlstfail= 0 %done= 67.28 ntets= 421571 node index= 70000 nlstfail= 0 %done= 72.32 ntets= 450442 node index= 75000 nlstfail= 0 %done= 77.27 ntets= 479297 node index= 80000 nlstfail= 0 %done= 82.22 ntets= 508214 node index= 85000 nlstfail= 0 %done= 87.18 ntets= 536869 node index= 90000 nlstfail= 0 %done= 92.09 ntets= 558852 node index= 95000 nlstfail= 0 %done= 95.86 The mesh is now complete! rmpoint compress 0 points removed and 0 elements removed. The new point count is 97156. The new element count is 540932. quality epsilonl, epsilonaspect: 2.4333781E-09 1.4408832E-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 : 113834 elements with aspect ratio b/w .2 and .5 : 355578 elements with aspect ratio b/w .5 and 1. : 71520 min aspect ratio = 0.1310E+00 max aspect ratio = 0.1000E+01 epsilonvol: 3.3762306E-02 --------------------------------------- element volumes b/w 0.5461E+04 and 0.1255E+05: 50305 element volumes b/w 0.1255E+05 and 0.2883E+05: 25006 element volumes b/w 0.2883E+05 and 0.6622E+05: 10393 element volumes b/w 0.6622E+05 and 0.1521E+06: 25144 element volumes b/w 0.1521E+06 and 0.3495E+06: 430084 element volumes b/w 0.1521E+06 and 0.3495E+06: 430084 min volume = 5.4613333E+03 max volume = 3.4952533E+05 ----------------------------------------------------------- 540932 total elements evaluated. resetpts itp geniee finish ** add voronoi volumes as node attribute * use coordinates at zero for better precision cmo copy cmotmp cmotet cmo select cmotmp trans 1,0,0 zero cmo/addatt/cmotmp/transauto/VDOUBLE/scalar/vector/user/permanent/x/0.0 finish cmo addatt cmotmp vor_volume vorvol ADDATT/voronoi_volume: filling new attribute: vorvol AMatbld3d_stor: *****Zero Negative Coefficients ****** AMatbld3d_stor: Number of 'zero' (< 10e-8*max) coefs 0 AMatbld3d_stor: npoints = 97156 n connections = 681486 AMatbld3d_stor: Number of written coefs = 389321 AMatbld3d_stor: Maximum num. connections to a node = 9 AMatbld3d_stor: Volume min = 1.6384000E+04 AMatbld3d_stor: Volume max = 2.0971520E+06 AMatbld3d_stor: Total Volume: 1.5205191E+11 AMatbld3d_stor: abs(Aij/xij) min = 0.0000000E+00 AMatbld3d_stor: abs(Aij/xij) max = 2.0480000E+03 AMatbld3d_stor: (Aij/xij) max = 0.0000000E+00 AMatbld3d_stor: (Aij/xij) min = -2.0480000E+03 Compression used for graph, not coefficient values The area coefficient values were written as scalar area/distance cmotmp attribute with voronoi volumes created with name vorvol *** SPARSE COEFFICIENT MATRIX _gstor SUCCESSFUL *** cmo select cmotmp pset pneg attribute vorvol 1,0,0 lt 0. THE PSET pneg HAS 0 POINTS pset pout attribute itp 1,0,0 eq 10. THE PSET pout HAS 13366 POINTS pset pin not pout THE PSET pin HAS 83790 POINTS pset pBAD inter pin pneg THE PSET pBAD HAS 0 POINTS cmo printatt cmotmp vorvol minmax ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH vorvol 1.638399998E+04 2.097152000E+06 2.080768000E+06 97156 cmo/addatt/cmotet/vorvol/VDOUBLE/scalar/nnodes/linear/permanent/afgx/0. cmo copyatt/cmotet cmotmp/vorvol vorvol 97156 copied from cmotmp vorvol to -> cmotet vorvol cmo delete cmotmp Released Mesh Object: cmotmp * set materials imt (from EV) * set well locations * X 482800, Y 476300 ft (147157.44, 145176.24 meters) * The target interval Weber Sandstone at the well location for * CO2 injection ranges from -10600 to -11800 ft (elevation) * -3230.88 to -3596.64 meters cmo copyatt/cmotet cmohex/imt imt 97156 copied from cmohex imt1 to -> cmotet imt1 cmo delete cmohex Released Mesh Object: cmohex cmo/addatt/cmotet/iwell/vint/scalar/nnodes/linear/permanent/afgx/0. cmo select cmotet define MATNO 50 define/XMIN/147150.00 define/XMAX/147160.00 define/YMIN/145170.00 define/YMAX/145180.00 define/ZMIN/-5000. define/ZMAX/0. pset/pxy/geom/xyz/1,0,0/XMIN YMIN ZMIN/XMAX YMAX ZMAX/ THE PSET pxy HAS 107 POINTS cmo/setatt/cmotet/iwell pset,get,pxy MATNO 107 values reset for attribute iwell define MATNO 51 define/ZMIN/-3597. define/ZMAX/-3230. pset/pinj/geom/xyz/1,0,0/XMIN YMIN ZMIN/XMAX YMAX ZMAX/ THE PSET pinj HAS 12 POINTS cmo/setatt/cmotet/iwell pset,get,pinj MATNO 12 values reset for attribute iwell pset/pxy/well_nodes.zone/ascii/50 THE PSET pxy HAS 0 POINTS * write xyz and node information cmo copy cmotmp cmotet cmo/DELATT/cmotmp/itp1 cmo/DELATT/cmotmp/icr1 cmo/DELATT/cmotmp/isn1 cmo/addatt/cmotmp/xcoord/VDOUBLE/scalar/nnodes/linear/ cmo/addatt/cmotmp/ycoord/VDOUBLE/scalar/nnodes/linear/ cmo/addatt/cmotmp/zcoord/VDOUBLE/scalar/nnodes/linear/ cmo/copyatt/cmotmp cmotet/xcoord xic 97156 copied from cmotet xic to -> cmotmp xcoord cmo/copyatt/cmotmp cmotet/ycoord yic 97156 copied from cmotet yic to -> cmotmp ycoord cmo/copyatt/cmotmp cmotet/zcoord zic 97156 copied from cmotet zic to -> cmotmp zcoord dump/avs2/node_attributes.dat/cmotmp/0 0 1 0 WARNING: dump/avs iopt_points= 0 WARNING: dump/avs iopt_elements= 0 WARNING: dump/avs iopt_values_node= 1 WARNING: dump/avs iopt_values_elem= 0 WARNING: dump/avs will produce non-standard AVS output that WARNING: read/avs may not be able to read. cmo/modatt/-def-/-def-/ioflag/x finish cmo delete cmotmp Released Mesh Object: cmotmp cmo printatt cmotet -all- minmax ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH -def- 0.000000000E+00 0.000000000E+00 0.000000000E+00 97156 scalar 1 1 0 1 vector 3 3 0 1 nnodes 97156 97156 0 1 nedges 0 0 0 1 nfaces 0 0 0 1 nelements 540932 540932 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 2 2 97156 ialias 0 0 0 97156 imt1 1 13 12 97156 itp1 0 10 10 97156 icr1 0 0 0 97156 isn1 0 0 0 97156 xic 1.448530000E+05 1.497170000E+05 4.864000000E+03 97156 yic 1.426160000E+05 1.518320000E+05 9.216000000E+03 97156 zic -3.840000000E+03 -4.480000000E+02 3.392000000E+03 97156 xtetwd 0 0 0 540932 itetclr 1 1 0 540932 itettyp 5 5 0 540932 itetoff 0 2163724 2163724 540932 jtetoff 0 2163724 2163724 540932 itet 1 97156 97155 540932x4 jtet 1 16000000 15999999 540932x4 epsilon 1.000000004E-15 1.000000004E-15 0.000000000E+00 1 epsilonl 2.433378062E-09 2.433378062E-09 0.000000000E+00 1 epsilona 2.666729415E-05 2.666729415E-05 0.000000000E+00 1 epsilonv 3.376230597E-02 3.376230597E-02 0.000000000E+00 1 ipointi 1 1 0 1 ipointj 97156 97156 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 1.448530000E+05 1.448530000E+05 0.000000000E+00 1 ymin 1.426160000E+05 1.426160000E+05 0.000000000E+00 1 zmin -3.840000000E+03 -3.840000000E+03 0.000000000E+00 1 xmax 1.497170000E+05 1.497170000E+05 0.000000000E+00 1 ymax 1.518320000E+05 1.518320000E+05 0.000000000E+00 1 zmax -4.480000000E+02 -4.480000000E+02 0.000000000E+00 1 kdtree_level 0 0 0 1 max_number_sets 64 64 0 1 number_of_psets 2 2 0 1 number_of_eltsets 0 0 0 1 number_of_fsets 0 0 0 1 vorvol 1.638399998E+04 2.097152000E+06 2.080768000E+06 97156 iwell 0 51 51 97156 dump avs2 tet_connect.inp cmotet cmo/modatt/-def-/-def-/ioflag/x finish dump gmv tet_connect.gmv cmotet cmo/modatt/-def-/-def-/ioflag/x finish dump tecplot tet_grid cmotet ascii grid dump fehm tet cmotet ascii *** Write FEHMN GEOM AND ZONE FILES *** *********dump_material_lists******** Minimum material ID value = 1 Maximum material ID value = 13 Total possible materials = 13 Material 1 has 5254 nodes. #nodes/nnodes is 0.540779791772E-01 Material 2 has 4026 nodes. #nodes/nnodes is 0.414385125041E-01 Material 3 has 10764 nodes. #nodes/nnodes is 0.110790893435 Material 4 has 3051 nodes. #nodes/nnodes is 0.314031057060E-01 Material 5 has 3180 nodes. #nodes/nnodes is 0.327308662236E-01 Material 6 has 4695 nodes. #nodes/nnodes is 0.483243428171E-01 Material 7 has 6997 nodes. #nodes/nnodes is 0.720181986690E-01 Material 8 has 5329 nodes. #nodes/nnodes is 0.548499338329E-01 Material 9 has 12344 nodes. #nodes/nnodes is 0.127053394914 Material 10 has 9940 nodes. #nodes/nnodes is 0.102309688926 Material 11 has 3804 nodes. #nodes/nnodes is 0.391535274684E-01 Material 12 has 6347 nodes. #nodes/nnodes is 0.653279274702E-01 Material 13 has 21425 nodes. #nodes/nnodes is 0.220521628857 *********dump_interface_list******** *********dump_multi_material_list******** Material 1 has 4915 multi-material connections. Material 2 has 9742 multi-material connections. Material 3 has 9425 multi-material connections. Material 4 has 9170 multi-material connections. Material 5 has 9115 multi-material connections. Material 6 has 8983 multi-material connections. Material 7 has 8922 multi-material connections. Material 8 has 8927 multi-material connections. Material 9 has 8654 multi-material connections. Material 10 has 8248 multi-material connections. Material 11 has 8005 multi-material connections. Material 12 has 7947 multi-material connections. Material 13 has 3981 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 908 nodes. top Sum Voronoi Area_x Area_y Area_z 0.2949120E+06 0.1556480E+06 0.4482662E+08 cmo/addatt/cmotet/bottom/vint/scalar/nnodes/linear/permanent/afgx/-5.0/ finish Face bottom 2 has 908 nodes. bottom Sum Voronoi Area_x Area_y Area_z 0.2949120E+06 0.1556480E+06 0.4482662E+08 cmo/addatt/cmotet/left_w/vint/scalar/nnodes/linear/permanent/afgx/-5.0/ finish Face left_w 3 has 3959 nodes. left_w Sum Voronoi Area_x Area_y Area_z 0.3126067E+08 0.8683520E+06 0.2359296E+07 cmo/addatt/cmotet/right_e/vint/scalar/nnodes/linear/permanent/afgx/-5.0/ finish Face right_e 5 has 3959 nodes. right_e Sum Voronoi Area_x Area_y Area_z 0.3126067E+08 0.8683520E+06 0.2359296E+07 cmo/addatt/cmotet/back_n/vint/scalar/nnodes/linear/permanent/afgx/-5.0/ finish Face back_n 6 has 2140 nodes. back_n Sum Voronoi Area_x Area_y Area_z 0.8683520E+06 0.1649869E+08 0.1245184E+07 cmo/addatt/cmotet/front_s/vint/scalar/nnodes/linear/permanent/afgx/-5.0/ finish Face front_s 4 has 2140 nodes. front_s Sum Voronoi Area_x Area_y Area_z 0.8683520E+06 0.1649869E+08 0.1245184E+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 0 AMatbld3d_stor: npoints = 97156 n connections = 681486 AMatbld3d_stor: Number of written coefs = 28 AMatbld3d_stor: Maximum num. connections to a node = 9 AMatbld3d_stor: Volume min = 1.6383875E+04 AMatbld3d_stor: Volume max = 2.0971523E+06 AMatbld3d_stor: Total Volume: 1.5205191E+11 AMatbld3d_stor: abs(Aij/xij) min = 0.0000000E+00 AMatbld3d_stor: abs(Aij/xij) max = 2.0480000E+03 AMatbld3d_stor: (Aij/xij) max = 0.0000000E+00 AMatbld3d_stor: (Aij/xij) min = -2.0480000E+03 Compression used for graph and coefficient values The area coefficient values were written as scalar area/distance ascii STOR file written with name tet.stor *** SPARSE COEFFICIENT MATRIX _astor SUCCESSFUL *** 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: *****Zero Negative Coefficients ****** AMatbld3d_stor: Number of 'zero' (< 10e-8*max) coefs 0 AMatbld3d_stor: npoints = 97156 n connections = 681486 AMatbld3d_stor: Number of written coefs = 28 AMatbld3d_stor: Maximum num. connections to a node = 9 AMatbld3d_stor: Volume min = 1.6384000E+04 AMatbld3d_stor: Volume max = 2.0971520E+06 AMatbld3d_stor: Total Volume: 1.5205191E+11 AMatbld3d_stor: abs(Aij/xij) min = 0.0000000E+00 AMatbld3d_stor: abs(Aij/xij) max = 2.0480000E+03 AMatbld3d_stor: (Aij/xij) max = 0.0000000E+00 AMatbld3d_stor: (Aij/xij) min = -2.0480000E+03 Compression used for graph and coefficient values The area coefficient values were written as scalar area/distance ascii STOR file written with name tet.stor *** SPARSE COEFFICIENT MATRIX _astor SUCCESSFUL *** finish LaGriT successfully completed