* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * Program: LaGriT V3.101 Linux m64 * * * date_compile: 2013/04/25 RH gf4.5 * * * Run Time: 2013/Oct 16 09:09:13 * * * 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 # stack the 2D cylinder surfaces into a 3D tet mesh # read the 2D Delaunay triangulation # and use as top and bottom of stacked mesh define INPUT_FILE cyl_2D_delaunay.inp define Z_BOTTOM 0.0 define Z_TOP 15.0 define NLAYERS 15 read/avs/INPUT_FILE/mo_tri cmo/addatt/mo_tri/idnode0/VINT/scalar/nnodes/linear/permanent/gxaf/0.0 finish cmo/addatt/mo_tri/wdist/VDOUBLE/scalar/nnodes/linear/permanent/gxaf/0.0 finish geniee finish cmo/status/brief The current-mesh-object(CMO) is: mo_tri 1 Mesh Object name: mo_tri number of nodes = 494 number of elements = 894 dimensions geometry = 3 element type = tri dimensions topology = 2 3 nodes 3 faces 3 edges boundary flag = 16000000 status = active finish # Write the top and bottom surfaces for stack command cmo/setatt/mo_tri/zic/Z_BOTTOM 494 values reset for attribute zic dump/avs/tmp_bottom.inp/mo_tri cmo/modatt/-def-/-def-/ioflag/x finish cmo/setatt/mo_tri/zic/Z_TOP 494 values reset for attribute zic dump/avs/tmp_top.inp/mo_tri cmo/modatt/-def-/-def-/ioflag/x finish ############# STACK LAYERS ####################### # assume itetclr and imt are set with materials # all triangle faces are upward normal cmo/create/mo_stack stack/layers/avs/tmp_bottom.inp tmp_top.inp/1 NLAYERS syntax warning: file1 integer / file2 integer integer refine starts at second surface, number ignored for file1. Layers to create: 17 Max material number: 1 Reading 2 surface files... cmo/addatt/mo_stack/layertyp/VINT/scalar/nnodes//permanent/agfx finish Read first surface: tmp_bottom.inp read/avs/tmp_bottom.inp/def1 cmo/addatt/def1/idnode0/VINT/scalar/nnodes/linear/permanent/gxaf/0.0 finish cmo/addatt/def1/wdist/VDOUBLE/scalar/nnodes/linear/permanent/gxaf/0.0 finish geniee finish cmo/status/brief The current-mesh-object(CMO) is: def1 1 Mesh Object name: mo_tri number of nodes = 494 number of elements = 894 dimensions geometry = 3 element type = tri dimensions topology = 2 3 nodes 3 faces 3 edges boundary flag = 16000000 status = inactive 2 Mesh Object name: mo_stack number of nodes = 0 number of elements = 0 dimensions geometry = 3 element type = tet dimensions topology = 3 4 nodes 4 faces 6 edges boundary flag = 16000000 status = inactive 3 Mesh Object name: def1 number of nodes = 494 number of elements = 894 dimensions geometry = 3 element type = tri dimensions topology = 2 3 nodes 3 faces 3 edges boundary flag = 16000000 status = active finish finish Read surface: tmp_top.inp read/avs/tmp_top.inp/def2 cmo/addatt/def2/idnode0/VINT/scalar/nnodes/linear/permanent/gxaf/0.0 finish cmo/addatt/def2/wdist/VDOUBLE/scalar/nnodes/linear/permanent/gxaf/0.0 finish geniee finish cmo/status/brief The current-mesh-object(CMO) is: def2 1 Mesh Object name: mo_tri number of nodes = 494 number of elements = 894 dimensions geometry = 3 element type = tri dimensions topology = 2 3 nodes 3 faces 3 edges boundary flag = 16000000 status = inactive 2 Mesh Object name: mo_stack number of nodes = 0 number of elements = 0 dimensions geometry = 3 element type = tri dimensions topology = 2 3 nodes 3 faces 3 edges boundary flag = 16000000 status = inactive 3 Mesh Object name: def1 number of nodes = 494 number of elements = 894 dimensions geometry = 3 element type = tri dimensions topology = 2 3 nodes 3 faces 3 edges boundary flag = 16000000 status = inactive 4 Mesh Object name: cmoprev number of nodes = 494 number of elements = 894 dimensions geometry = 3 element type = tri dimensions topology = 2 3 nodes 3 faces 3 edges boundary flag = 16000000 status = inactive 5 Mesh Object name: def2 number of nodes = 494 number of elements = 894 dimensions geometry = 3 element type = tri dimensions topology = 2 3 nodes 3 faces 3 edges boundary flag = 16000000 status = active finish finish trilayer/derive/proportional/cmonxt/def1/def2/6.250000E-02/ cmo/copy/cmonxt/def2/ finish Created new trilayer with average distance of 9.3750000E-01 New height: min/max 9.3750000E-01 9.3750000E-01 Old height: min/max 1.5000000E+01 1.5000000E+01 finish cmo/copyatt/cmonxt/def1/itetclr/itetclr/ 894 copied from def1 itetclr to -> cmonxt itetclr finish Layers compared: tmp_bottom.inp and refine cmo/delete/cmoprev/ Released Mesh Object: cmoprev finish cmo/copy/cmoprev/cmonxt/ finish addmesh/merge/mo_stack/def1/cmonxt/ searching for mo_stack attributes in cmonxt attribute not found: 66 layertyp ADDMESH: 1 attributes not copied from cmonxt to mo_stack finish trilayer/derive/proportional/cmonxt/def1/def2/1.250000E-01/ cmo/copy/cmonxt/def2/ Released Mesh Object: cmonxt finish Created new trilayer with average distance of 1.8750000E+00 New height: min/max 1.8750000E+00 1.8750000E+00 Old height: min/max 1.5000000E+01 1.5000000E+01 finish cmo/copyatt/cmonxt/def1/itetclr/itetclr/ 894 copied from def1 itetclr to -> cmonxt itetclr finish Layers compared: refine and refine cmo/delete/cmoprev/ Released Mesh Object: cmoprev finish cmo/copy/cmoprev/cmonxt/ finish addmesh/merge/mo_stack/mo_stack/cmonxt/ cmo/move/-tmp_source_internal-/mo_stack finish searching for mo_stack attributes in cmonxt attribute not found: 66 layertyp ADDMESH: 1 attributes not copied from cmonxt to mo_stack cmo/delete/-tmp_source_internal-/ Released Mesh Object: -tmp_source_internal- finish finish trilayer/derive/proportional/cmonxt/def1/def2/1.875000E-01/ cmo/copy/cmonxt/def2/ Released Mesh Object: cmonxt finish Created new trilayer with average distance of 2.8125000E+00 New height: min/max 2.8125000E+00 2.8125000E+00 Old height: min/max 1.5000000E+01 1.5000000E+01 finish cmo/copyatt/cmonxt/def1/itetclr/itetclr/ 894 copied from def1 itetclr to -> cmonxt itetclr finish Layers compared: refine and refine cmo/delete/cmoprev/ Released Mesh Object: cmoprev finish cmo/copy/cmoprev/cmonxt/ finish addmesh/merge/mo_stack/mo_stack/cmonxt/ cmo/move/-tmp_source_internal-/mo_stack finish searching for mo_stack attributes in cmonxt attribute not found: 66 layertyp ADDMESH: 1 attributes not copied from cmonxt to mo_stack cmo/delete/-tmp_source_internal-/ Released Mesh Object: -tmp_source_internal- finish finish trilayer/derive/proportional/cmonxt/def1/def2/2.500000E-01/ cmo/copy/cmonxt/def2/ Released Mesh Object: cmonxt finish Created new trilayer with average distance of 3.7500000E+00 New height: min/max 3.7500000E+00 3.7500000E+00 Old height: min/max 1.5000000E+01 1.5000000E+01 finish cmo/copyatt/cmonxt/def1/itetclr/itetclr/ 894 copied from def1 itetclr to -> cmonxt itetclr finish Layers compared: refine and refine cmo/delete/cmoprev/ Released Mesh Object: cmoprev finish cmo/copy/cmoprev/cmonxt/ finish addmesh/merge/mo_stack/mo_stack/cmonxt/ cmo/move/-tmp_source_internal-/mo_stack finish searching for mo_stack attributes in cmonxt attribute not found: 66 layertyp ADDMESH: 1 attributes not copied from cmonxt to mo_stack cmo/delete/-tmp_source_internal-/ Released Mesh Object: -tmp_source_internal- finish finish trilayer/derive/proportional/cmonxt/def1/def2/3.125000E-01/ cmo/copy/cmonxt/def2/ Released Mesh Object: cmonxt finish Created new trilayer with average distance of 4.6875000E+00 New height: min/max 4.6875000E+00 4.6875000E+00 Old height: min/max 1.5000000E+01 1.5000000E+01 finish cmo/copyatt/cmonxt/def1/itetclr/itetclr/ 894 copied from def1 itetclr to -> cmonxt itetclr finish Layers compared: refine and refine cmo/delete/cmoprev/ Released Mesh Object: cmoprev finish cmo/copy/cmoprev/cmonxt/ finish addmesh/merge/mo_stack/mo_stack/cmonxt/ cmo/move/-tmp_source_internal-/mo_stack finish searching for mo_stack attributes in cmonxt attribute not found: 66 layertyp ADDMESH: 1 attributes not copied from cmonxt to mo_stack cmo/delete/-tmp_source_internal-/ Released Mesh Object: -tmp_source_internal- finish finish trilayer/derive/proportional/cmonxt/def1/def2/3.750000E-01/ cmo/copy/cmonxt/def2/ Released Mesh Object: cmonxt finish Created new trilayer with average distance of 5.6250000E+00 New height: min/max 5.6250000E+00 5.6250000E+00 Old height: min/max 1.5000000E+01 1.5000000E+01 finish cmo/copyatt/cmonxt/def1/itetclr/itetclr/ 894 copied from def1 itetclr to -> cmonxt itetclr finish Layers compared: refine and refine cmo/delete/cmoprev/ Released Mesh Object: cmoprev finish cmo/copy/cmoprev/cmonxt/ finish addmesh/merge/mo_stack/mo_stack/cmonxt/ cmo/move/-tmp_source_internal-/mo_stack finish searching for mo_stack attributes in cmonxt attribute not found: 66 layertyp ADDMESH: 1 attributes not copied from cmonxt to mo_stack cmo/delete/-tmp_source_internal-/ Released Mesh Object: -tmp_source_internal- finish finish trilayer/derive/proportional/cmonxt/def1/def2/4.375000E-01/ cmo/copy/cmonxt/def2/ Released Mesh Object: cmonxt finish Created new trilayer with average distance of 6.5625000E+00 New height: min/max 6.5625000E+00 6.5625000E+00 Old height: min/max 1.5000000E+01 1.5000000E+01 finish cmo/copyatt/cmonxt/def1/itetclr/itetclr/ 894 copied from def1 itetclr to -> cmonxt itetclr finish Layers compared: refine and refine cmo/delete/cmoprev/ Released Mesh Object: cmoprev finish cmo/copy/cmoprev/cmonxt/ finish addmesh/merge/mo_stack/mo_stack/cmonxt/ cmo/move/-tmp_source_internal-/mo_stack finish searching for mo_stack attributes in cmonxt attribute not found: 66 layertyp ADDMESH: 1 attributes not copied from cmonxt to mo_stack cmo/delete/-tmp_source_internal-/ Released Mesh Object: -tmp_source_internal- finish finish trilayer/derive/proportional/cmonxt/def1/def2/5.000000E-01/ cmo/copy/cmonxt/def2/ Released Mesh Object: cmonxt finish Created new trilayer with average distance of 7.5000000E+00 New height: min/max 7.5000000E+00 7.5000000E+00 Old height: min/max 1.5000000E+01 1.5000000E+01 finish cmo/copyatt/cmonxt/def1/itetclr/itetclr/ 894 copied from def1 itetclr to -> cmonxt itetclr finish Layers compared: refine and refine cmo/delete/cmoprev/ Released Mesh Object: cmoprev finish cmo/copy/cmoprev/cmonxt/ finish addmesh/merge/mo_stack/mo_stack/cmonxt/ cmo/move/-tmp_source_internal-/mo_stack finish searching for mo_stack attributes in cmonxt attribute not found: 66 layertyp ADDMESH: 1 attributes not copied from cmonxt to mo_stack cmo/delete/-tmp_source_internal-/ Released Mesh Object: -tmp_source_internal- finish finish trilayer/derive/proportional/cmonxt/def1/def2/5.625000E-01/ cmo/copy/cmonxt/def2/ Released Mesh Object: cmonxt finish Created new trilayer with average distance of 8.4375000E+00 New height: min/max 8.4375000E+00 8.4375000E+00 Old height: min/max 1.5000000E+01 1.5000000E+01 finish cmo/copyatt/cmonxt/def1/itetclr/itetclr/ 894 copied from def1 itetclr to -> cmonxt itetclr finish Layers compared: refine and refine cmo/delete/cmoprev/ Released Mesh Object: cmoprev finish cmo/copy/cmoprev/cmonxt/ finish addmesh/merge/mo_stack/mo_stack/cmonxt/ cmo/move/-tmp_source_internal-/mo_stack finish searching for mo_stack attributes in cmonxt attribute not found: 66 layertyp ADDMESH: 1 attributes not copied from cmonxt to mo_stack cmo/delete/-tmp_source_internal-/ Released Mesh Object: -tmp_source_internal- finish finish trilayer/derive/proportional/cmonxt/def1/def2/6.250000E-01/ cmo/copy/cmonxt/def2/ Released Mesh Object: cmonxt finish Created new trilayer with average distance of 9.3750000E+00 New height: min/max 9.3750000E+00 9.3750000E+00 Old height: min/max 1.5000000E+01 1.5000000E+01 finish cmo/copyatt/cmonxt/def1/itetclr/itetclr/ 894 copied from def1 itetclr to -> cmonxt itetclr finish Layers compared: refine and refine cmo/delete/cmoprev/ Released Mesh Object: cmoprev finish cmo/copy/cmoprev/cmonxt/ finish addmesh/merge/mo_stack/mo_stack/cmonxt/ cmo/move/-tmp_source_internal-/mo_stack finish searching for mo_stack attributes in cmonxt attribute not found: 66 layertyp ADDMESH: 1 attributes not copied from cmonxt to mo_stack cmo/delete/-tmp_source_internal-/ Released Mesh Object: -tmp_source_internal- finish finish trilayer/derive/proportional/cmonxt/def1/def2/6.875000E-01/ cmo/copy/cmonxt/def2/ Released Mesh Object: cmonxt finish Created new trilayer with average distance of 1.0312500E+01 New height: min/max 1.0312500E+01 1.0312500E+01 Old height: min/max 1.5000000E+01 1.5000000E+01 finish cmo/copyatt/cmonxt/def1/itetclr/itetclr/ 894 copied from def1 itetclr to -> cmonxt itetclr finish Layers compared: refine and refine cmo/delete/cmoprev/ Released Mesh Object: cmoprev finish cmo/copy/cmoprev/cmonxt/ finish addmesh/merge/mo_stack/mo_stack/cmonxt/ cmo/move/-tmp_source_internal-/mo_stack finish searching for mo_stack attributes in cmonxt attribute not found: 66 layertyp ADDMESH: 1 attributes not copied from cmonxt to mo_stack cmo/delete/-tmp_source_internal-/ Released Mesh Object: -tmp_source_internal- finish finish trilayer/derive/proportional/cmonxt/def1/def2/7.500000E-01/ cmo/copy/cmonxt/def2/ Released Mesh Object: cmonxt finish Created new trilayer with average distance of 1.1250000E+01 New height: min/max 1.1250000E+01 1.1250000E+01 Old height: min/max 1.5000000E+01 1.5000000E+01 finish cmo/copyatt/cmonxt/def1/itetclr/itetclr/ 894 copied from def1 itetclr to -> cmonxt itetclr finish Layers compared: refine and refine cmo/delete/cmoprev/ Released Mesh Object: cmoprev finish cmo/copy/cmoprev/cmonxt/ finish addmesh/merge/mo_stack/mo_stack/cmonxt/ cmo/move/-tmp_source_internal-/mo_stack finish searching for mo_stack attributes in cmonxt attribute not found: 66 layertyp ADDMESH: 1 attributes not copied from cmonxt to mo_stack cmo/delete/-tmp_source_internal-/ Released Mesh Object: -tmp_source_internal- finish finish trilayer/derive/proportional/cmonxt/def1/def2/8.125000E-01/ cmo/copy/cmonxt/def2/ Released Mesh Object: cmonxt finish Created new trilayer with average distance of 1.2187500E+01 New height: min/max 1.2187500E+01 1.2187500E+01 Old height: min/max 1.5000000E+01 1.5000000E+01 finish cmo/copyatt/cmonxt/def1/itetclr/itetclr/ 894 copied from def1 itetclr to -> cmonxt itetclr finish Layers compared: refine and refine cmo/delete/cmoprev/ Released Mesh Object: cmoprev finish cmo/copy/cmoprev/cmonxt/ finish addmesh/merge/mo_stack/mo_stack/cmonxt/ cmo/move/-tmp_source_internal-/mo_stack finish searching for mo_stack attributes in cmonxt attribute not found: 66 layertyp ADDMESH: 1 attributes not copied from cmonxt to mo_stack cmo/delete/-tmp_source_internal-/ Released Mesh Object: -tmp_source_internal- finish finish trilayer/derive/proportional/cmonxt/def1/def2/8.750000E-01/ cmo/copy/cmonxt/def2/ Released Mesh Object: cmonxt finish Created new trilayer with average distance of 1.3125000E+01 New height: min/max 1.3125000E+01 1.3125000E+01 Old height: min/max 1.5000000E+01 1.5000000E+01 finish cmo/copyatt/cmonxt/def1/itetclr/itetclr/ 894 copied from def1 itetclr to -> cmonxt itetclr finish Layers compared: refine and refine cmo/delete/cmoprev/ Released Mesh Object: cmoprev finish cmo/copy/cmoprev/cmonxt/ finish addmesh/merge/mo_stack/mo_stack/cmonxt/ cmo/move/-tmp_source_internal-/mo_stack finish searching for mo_stack attributes in cmonxt attribute not found: 66 layertyp ADDMESH: 1 attributes not copied from cmonxt to mo_stack cmo/delete/-tmp_source_internal-/ Released Mesh Object: -tmp_source_internal- finish finish trilayer/derive/proportional/cmonxt/def1/def2/9.375000E-01/ cmo/copy/cmonxt/def2/ Released Mesh Object: cmonxt finish Created new trilayer with average distance of 1.4062500E+01 New height: min/max 1.4062500E+01 1.4062500E+01 Old height: min/max 1.5000000E+01 1.5000000E+01 finish cmo/copyatt/cmonxt/def1/itetclr/itetclr/ 894 copied from def1 itetclr to -> cmonxt itetclr finish Layers compared: refine and refine cmo/delete/cmoprev/ Released Mesh Object: cmoprev finish cmo/copy/cmoprev/cmonxt/ finish addmesh/merge/mo_stack/mo_stack/cmonxt/ cmo/move/-tmp_source_internal-/mo_stack finish searching for mo_stack attributes in cmonxt attribute not found: 66 layertyp ADDMESH: 1 attributes not copied from cmonxt to mo_stack cmo/delete/-tmp_source_internal-/ Released Mesh Object: -tmp_source_internal- finish finish Released Mesh Object: cmonxt cmo/delete/def2/ Released Mesh Object: def2 finish Layers compared: refine and tmp_top.inp cmo/delete/cmoprev/ Released Mesh Object: cmoprev finish cmo/copy/cmoprev/cmonxt/ finish addmesh/merge/mo_stack/mo_stack/cmonxt/ cmo/move/-tmp_source_internal-/mo_stack finish searching for mo_stack attributes in cmonxt attribute not found: 66 layertyp ADDMESH: 1 attributes not copied from cmonxt to mo_stack cmo/delete/-tmp_source_internal-/ Released Mesh Object: -tmp_source_internal- finish finish ................................................................ surface name layer color type llcorner zic tmp_bottom.inp 1 1 -1 1) 0.000000E+00 refine 2 1 2 495) 9.375000E-01 refine 3 1 2 989) 1.875000E+00 refine 4 1 2 1483) 2.812500E+00 refine 5 1 2 1977) 3.750000E+00 refine 6 1 2 2471) 4.687500E+00 refine 7 1 2 2965) 5.625000E+00 refine 8 1 2 3459) 6.562500E+00 refine 9 1 2 3953) 7.500000E+00 refine 10 1 2 4447) 8.437500E+00 refine 11 1 2 4941) 9.375000E+00 refine 12 1 2 5435) 1.031250E+01 refine 13 1 2 5929) 1.125000E+01 refine 14 1 2 6423) 1.218750E+01 refine 15 1 2 6917) 1.312500E+01 refine 16 1 2 7411) 1.406250E+01 tmp_top.inp 17 1 -2 7905) 1.500000E+01 Elements per layer: 894 stacked total: 15198 Nodes per layer: 494 stacked total: 8398 No Truncating layer specified. files read: 2 from total: 2 ................................................................ cmo/addatt/mo_stack/nlayers/INT/scalar/scalar/constant// finish cmo/addatt/mo_stack/nnperlayer/INT/scalar/scalar/constant// finish cmo/addatt/mo_stack/neperlayer/INT/scalar/scalar/constant// finish cmo/setatt/mo_stack/nlayers 17 1 values reset for attribute nlayers finish cmo/setatt/mo_stack/nnperlayer 494 1 values reset for attribute nnperlayer finish cmo/setatt/mo_stack/neperlayer 894 1 values reset for attribute neperlayer finish Released Mesh Object: def1 Released Mesh Object: cmonxt Released Mesh Object: cmoprev stack done. dump gmv tmp_layers_stack.gmv mo_stack cmo/modatt/-def-/-def-/ioflag/x finish dump avs tmp_layers_stack.inp mo_stack cmo/modatt/-def-/-def-/ioflag/x finish cmo printatt mo_stack -xyz- minmax ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH xic -1.500000000E+01 1.500000000E+01 3.000000000E+01 8398 yic -1.500000000E+01 1.500000000E+01 3.000000000E+01 8398 zic 0.000000000E+00 1.500000000E+01 1.500000000E+01 8398 # fill layers with prism elements # this can be used with hextotet if connect does # not preserve the central cylinder smoothly stack/fill/mo_prism/mo_stack cmo/copy/mo_prism/mo_stack/ finish geniee finish cmo select mo_prism filter/1,0,0 FILTER:Use internal epsilonl value = 0.999999993923E-08 FILTER:Dudding 0 duplicate points. rmpoint/compress 0 points removed and 0 elements removed. The new point count is 8398. The new element count is 14304. resetpts itp geniee finish define WELL_ID 1 define MAT_ID 2 define CMO mo_prism infile color_well.mlgi # Select center cylinder set of nodes and cells for well # MUST DEFINE MAT_ID, WELL_ID, CMO cmo setatt CMO imt MAT_ID 8398 values reset for attribute imt1 cmo setatt CMO itetclr MAT_ID 14304 values reset for attribute itetclr # pset select the center node pset/px/attribute xic/1,0,0/eq 0. THE PSET px HAS 17 POINTS pset/pcenter/attribute yic/pset,get,px/eq 0. THE PSET pcenter HAS 17 POINTS # eltset select cells connected to center node eltset/ewell/inclusive pset,get,pcenter THE ELTSET ewell HAS 112 ELEMENTS # pset select all node vertices of well cells pset/pwell/eltset/ewell THE PSET pwell HAS 136 POINTS # overwrite material MAT_ID with WELL_ID for well cmo setatt CMO imt pset,get,pwell WELL_ID 136 values reset for attribute imt1 cmo setatt CMO itetclr eltset,get,ewell WELL_ID 112 values reset for attribute itetclr finish dump/gmv/prism_stack.gmv/mo_prism cmo/modatt/-def-/-def-/ioflag/x finish dump/avs/prism_stack.inp/mo_prism cmo/modatt/-def-/-def-/ioflag/x finish cmo status mo_prism The current-mesh-object(CMO) is: mo_prism 3 Mesh Object name: mo_prism number of nodes = 8398 number of elements = 14304 dimensions geometry = 3 element type = pri dimensions topology = 3 6 nodes 5 faces 9 edges boundary flag = 16000000 status = active NAME TYPE RANK LENGTH INTER PERSI IO VALUE 1 -def- VDOU scalar nnodes linea tempo x 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 8.398E+03 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 1.430E+04 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 6.000E+00 12 edges_per_eleme INT scalar scalar const perma l 9.000E+00 13 faces_per_eleme INT scalar scalar const perma l 5.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 1.000E-08 35 epsilona REAL scalar scalar const perma l 1.000E-08 36 epsilonv REAL scalar scalar const perma l 1.000E-08 37 ipointi INT scalar scalar const perma l 1.000E+00 38 ipointj INT scalar scalar const perma l 8.398E+03 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 -1.000E+00 52 ymin REAL scalar scalar const perma l -1.000E+00 53 zmin REAL scalar scalar const perma l -1.000E+00 54 xmax REAL scalar scalar const perma l 1.000E+00 55 ymax REAL scalar scalar const perma l 1.000E+00 56 zmax REAL scalar scalar const perma l 1.000E+00 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 3.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 layertyp VINT scalar nnodes linea perma agf 0.000E+00 67 nlayers INT scalar scalar const tempo agl 1.700E+01 68 nnperlayer INT scalar scalar const tempo agl 4.940E+02 69 neperlayer INT scalar scalar const tempo agl 8.940E+02 cmo printatt mo_prism -all- minmax ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH -def- 0.000000000E+00 0.000000000E+00 0.000000000E+00 8398 scalar 1 1 0 1 vector 3 3 0 1 nnodes 8398 8398 0 1 nedges 0 0 0 1 nfaces 0 0 0 1 nelements 14304 14304 0 1 mbndry 16000000 16000000 0 1 ndimensions_topo 3 3 0 1 ndimensions_geom 3 3 0 1 nodes_per_element 6 6 0 1 edges_per_element 9 9 0 1 faces_per_element 5 5 0 1 isetwd 0 7 7 8398 ialias 0 7904 7904 8398 imt1 1 2 1 8398 itp1 0 10 10 8398 icr1 0 0 0 8398 isn1 0 0 0 8398 xic -1.500000000E+01 1.500000000E+01 3.000000000E+01 8398 yic -1.500000000E+01 1.500000000E+01 3.000000000E+01 8398 zic 0.000000000E+00 1.500000000E+01 1.500000000E+01 8398 xtetwd 0 1 1 14304 itetclr 1 2 1 14304 itettyp 7 7 0 14304 itetoff 0 85818 85818 14304 jtetoff 0 71515 71515 14304 itet 1 8398 8397 14304x6 jtet 2 16000000 15999998 14304x5 epsilon 1.000000004E-15 1.000000004E-15 0.000000000E+00 1 epsilonl 9.999999939E-09 9.999999939E-09 0.000000000E+00 1 epsilona 9.999999939E-09 9.999999939E-09 0.000000000E+00 1 epsilonv 9.999999939E-09 9.999999939E-09 0.000000000E+00 1 ipointi 1 1 0 1 ipointj 8398 8398 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.000000000E+00 -1.000000000E+00 0.000000000E+00 1 ymin -1.000000000E+00 -1.000000000E+00 0.000000000E+00 1 zmin -1.000000000E+00 -1.000000000E+00 0.000000000E+00 1 xmax 1.000000000E+00 1.000000000E+00 0.000000000E+00 1 ymax 1.000000000E+00 1.000000000E+00 0.000000000E+00 1 zmax 1.000000000E+00 1.000000000E+00 0.000000000E+00 1 kdtree_level 0 0 0 1 max_number_sets 64 64 0 1 number_of_psets 3 3 0 1 number_of_eltsets 0 0 0 1 number_of_fsets 0 0 0 1 layertyp -2 2 4 8398 nlayers 17 17 0 1 nnperlayer 494 494 0 1 neperlayer 894 894 0 1 cmo printatt mo_prism itetclr minmax ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH itetclr 1 2 1 14304 ############# STACK LAYERS ####################### # Connect points into 3D Delaunay mesh cmo create cmotet///tet copypts cmotet mo_prism cmo/addatt/cmotet/layertyp/VINT/scalar/nnodes/linear/permanent/agfx/0.0 finish copypts added nodes 1 to 8398 cmo select cmotet cmo printatt cmotet -xyz- minmax Warning: nelements = 0 ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH xic -1.500000000E+01 1.500000000E+01 3.000000000E+01 8398 yic -1.500000000E+01 1.500000000E+01 3.000000000E+01 8398 zic 0.000000000E+00 1.500000000E+01 1.500000000E+01 8398 cmo/setatt/cmotet/imt 1 Warning: nelements = 0 8398 values reset for attribute imt1 filter 1,0,0 FILTER:Use internal epsilonl value = 0.999200722163E-11 FILTER:Dudding 0 duplicate points. rmpoint compress 0 points removed and 0 elements removed. The new point count is 8398. The new element count is 0. dump gmv points_for_tet.gmv cmotet cmo/modatt/-def-/-def-/ioflag/x finish dump avs2 points_for_tet.inp cmotet cmo/modatt/-def-/-def-/ioflag/x finish connect noadd Coordinates of enclosing tetrahedron are: -0.37500D+02 -0.37500D+02 -0.30000D+02 -0.37500D+02 0.11250D+03 0.75000D+01 -0.37500D+02 -0.37500D+02 0.45000D+02 0.11250D+03 0.00000D+00 0.75000D+01 ntets= 30870 node index= 5000 nlstfail= 0 %done= 61.24 The mesh is now complete! LaGriT FINISH: connect rmpoint compress 0 points removed and 0 elements removed. The new point count is 8398. The new element count is 42912. quality epsilonl, epsilonaspect: 9.9920072E-12 9.9760408E-34 -------------------------------------------- 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 : 0 elements with aspect ratio b/w .2 and .5 : 0 elements with aspect ratio b/w .5 and 1. : 42912 min aspect ratio = 0.4114E+00 max aspect ratio = 0.8150E+00 epsilonvol: 2.9976022E-09 --------------------------------------- element volumes b/w 0.1118E+00 and 0.1473E+00: 816 element volumes b/w 0.1473E+00 and 0.1940E+00: 4944 element volumes b/w 0.1940E+00 and 0.2556E+00: 18720 element volumes b/w 0.2556E+00 and 0.3368E+00: 17136 element volumes b/w 0.3368E+00 and 0.4438E+00: 1296 min volume = 1.1176980E-01 max volume = 4.4378813E-01 ----------------------------------------------------------- 42912 total elements evaluated. resetpts itp geniee finish cmo addatt cmotet vor_volume vorvol ADDATT/voronoi_volume: creating new attribute: vorvol AMatbld3d_stor: *****Zero Negative Coefficients ****** AMatbld3d_stor: Number of 'zero' (< 10e-8 *max) coefs 0 AMatbld3d_stor: npoints = 8398 ncoefs = 71364 AMatbld3d_stor: Number of unique coefs = 39881 AMatbld3d_stor: Maximum num. connections to a node = 13 AMatbld3d_stor: Volume min = 1.6537376E-01 AMatbld3d_stor: Volume max = 2.2975411E+00 AMatbld3d_stor: Total Volume: 1.0594635E+04 AMatbld3d_stor: abs(Aij/xij) min = 0.0000000E+00 AMatbld3d_stor: abs(Aij/xij) max = 2.6140912E+00 AMatbld3d_stor: (Aij/xij) max = 0.0000000E+00 AMatbld3d_stor: (Aij/xij) min = -2.6140912E+00 AMatbld3d_stor Matrix coefficient values stored as scalar area/distance AMatbld3d_stor Matrix compression used for graph, not coefficient values cmotet attribute with voronoi volumes created with name vorvol *** SPARSE COEFFICIENT MATRIX _gstor SUCCESSFUL *** 3D Matrix Coefficient file written with name -notset- cmo select cmotet cmo printatt cmotet ccoef minmax cannot attribute ccoef in cmotet define WELL_ID 1 define MAT_ID 2 define CMO cmotet infile color_well.mlgi # Select center cylinder set of nodes and cells for well # MUST DEFINE MAT_ID, WELL_ID, CMO cmo setatt CMO imt MAT_ID 8398 values reset for attribute imt1 cmo setatt CMO itetclr MAT_ID 42912 values reset for attribute itetclr # pset select the center node pset/px/attribute xic/1,0,0/eq 0. THE PSET px HAS 17 POINTS pset/pcenter/attribute yic/pset,get,px/eq 0. THE PSET pcenter HAS 17 POINTS # eltset select cells connected to center node eltset/ewell/inclusive pset,get,pcenter THE ELTSET ewell HAS 336 ELEMENTS # pset select all node vertices of well cells pset/pwell/eltset/ewell THE PSET pwell HAS 136 POINTS # overwrite material MAT_ID with WELL_ID for well cmo setatt CMO imt pset,get,pwell WELL_ID 136 values reset for attribute imt1 cmo setatt CMO itetclr eltset,get,ewell WELL_ID 336 values reset for attribute itetclr finish dump/fehm/cyl_tet/cmotet/ascii *** Write FEHMN GEOM AND ZONE FILES *** *********dump_material_lists******** Minimum material ID value = 1 Maximum material ID value = 2 Total possible materials = 2 Material 1 has 136 nodes. #nodes/nnodes is 0.161943323910E-01 Material 2 has 8262 nodes. #nodes/nnodes is 0.983805656433 *********dump_interface_list******** *********dump_multi_material_list******** Material 1 has 660 multi-material connections. Material 2 has 660 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 494 nodes. top Sum Voronoi Area_x Area_y Area_z 0.2809222E+02 0.2809222E+02 0.7063090E+03 cmo/addatt/cmotet/bottom/vint/scalar/nnodes/linear/permanent/afgx/-5.0/ finish Face bottom 2 has 494 nodes. bottom Sum Voronoi Area_x Area_y Area_z 0.2809222E+02 0.2809222E+02 0.7063090E+03 cmo/addatt/cmotet/left_w/vint/scalar/nnodes/linear/permanent/afgx/-5.0/ finish Face left_w 3 has 425 nodes. left_w Sum Voronoi Area_x Area_y Area_z 0.3389731E+03 0.1539062E+03 0.3240273E+02 cmo/addatt/cmotet/right_e/vint/scalar/nnodes/linear/permanent/afgx/-5.0/ finish Face right_e 5 has 425 nodes. right_e Sum Voronoi Area_x Area_y Area_z 0.3389731E+03 0.1539062E+03 0.3531194E+02 cmo/addatt/cmotet/back_n/vint/scalar/nnodes/linear/permanent/afgx/-5.0/ finish Face back_n 6 has 425 nodes. back_n Sum Voronoi Area_x Area_y Area_z 0.1539062E+03 0.3389731E+03 0.3277769E+02 cmo/addatt/cmotet/front_s/vint/scalar/nnodes/linear/permanent/afgx/-5.0/ finish Face front_s 4 has 425 nodes. front_s Sum Voronoi Area_x Area_y Area_z 0.1539062E+03 0.3389731E+03 0.3351008E+02 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 = 8398 ncoefs = 71364 AMatbld3d_stor: Number of unique coefs = 3255 AMatbld3d_stor: Maximum num. connections to a node = 13 AMatbld3d_stor: Volume min = 1.6537376E-01 AMatbld3d_stor: Volume max = 2.2975411E+00 AMatbld3d_stor: Total Volume: 1.0594635E+04 AMatbld3d_stor: abs(Aij/xij) min = 0.0000000E+00 AMatbld3d_stor: abs(Aij/xij) max = 2.6140912E+00 AMatbld3d_stor: (Aij/xij) max = 0.0000000E+00 AMatbld3d_stor: (Aij/xij) min = -2.6140912E+00 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 cyl_tet.stor *** SPARSE COEFFICIENT MATRIX _astor SUCCESSFUL *** 3D Matrix Coefficient file written with name cyl_tet.stor dump/zone_outside/cyl_tet/cmotet/keepatt *********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 494 nodes. top Sum Voronoi Area_x Area_y Area_z 0.2809222E+02 0.2809222E+02 0.7063090E+03 cmo/addatt/cmotet/bottom/vint/scalar/nnodes/linear/permanent/afgx/-5.0/ finish Face bottom 2 has 494 nodes. bottom Sum Voronoi Area_x Area_y Area_z 0.2809222E+02 0.2809222E+02 0.7063090E+03 cmo/addatt/cmotet/left_w/vint/scalar/nnodes/linear/permanent/afgx/-5.0/ finish Face left_w 3 has 425 nodes. left_w Sum Voronoi Area_x Area_y Area_z 0.3389731E+03 0.1539062E+03 0.3240273E+02 cmo/addatt/cmotet/right_e/vint/scalar/nnodes/linear/permanent/afgx/-5.0/ finish Face right_e 5 has 425 nodes. right_e Sum Voronoi Area_x Area_y Area_z 0.3389731E+03 0.1539062E+03 0.3531194E+02 cmo/addatt/cmotet/back_n/vint/scalar/nnodes/linear/permanent/afgx/-5.0/ finish Face back_n 6 has 425 nodes. back_n Sum Voronoi Area_x Area_y Area_z 0.1539062E+03 0.3389731E+03 0.3277769E+02 cmo/addatt/cmotet/front_s/vint/scalar/nnodes/linear/permanent/afgx/-5.0/ finish Face front_s 4 has 425 nodes. front_s Sum Voronoi Area_x Area_y Area_z 0.1539062E+03 0.3389731E+03 0.3351008E+02 dump gmv cyl_tet_delaunay.gmv cmotet cmo/modatt/-def-/-def-/ioflag/x finish dump avs cyl_tet_delaunay.inp cmotet cmo/modatt/-def-/-def-/ioflag/x finish finish LaGriT successfully completed