* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * Program: LaGriT V3.101 Linux m64 * * * date_compile: 2013/04/25 RH gf4.5 * * * Run Time: 2013/Jun 19 16:31:43 * * * 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 # refine base box shape mesh in area of repo and fault # Background spacing is 50m # refine halves are 25.0 12.5 6.25 # # Materials are Background 5, fault 4, repo 3 # Zones are Top left 1, Top right 2 # read avs hex_base.inp mo geniee finish cmo/status/brief The current-mesh-object(CMO) is: mo 1 Mesh Object name: mo number of nodes = 5397 number of elements = 4320 dimensions geometry = 3 element type = hex dimensions topology = 3 8 nodes 6 faces 12 edges boundary flag = 16000000 status = active finish cmo/addatt/mo/itetlev/vint/scalar/nelements cmo/setatt/mo/ipolydat no 1 values reset for attribute ipolydat cmo status mo The current-mesh-object(CMO) is: mo 1 Mesh Object name: mo number of nodes = 5397 number of elements = 4320 dimensions geometry = 3 element type = hex dimensions topology = 3 8 nodes 6 faces 12 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 5.397E+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 4.320E+03 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 8.000E+00 12 edges_per_eleme INT scalar scalar const perma l 1.200E+01 13 faces_per_eleme INT scalar scalar const perma l 6.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 no 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.411E-10 35 epsilona REAL scalar scalar const perma l 5.240E-07 36 epsilonv REAL scalar scalar const perma l 1.332E-04 37 ipointi INT scalar scalar const perma l 1.000E+00 38 ipointj INT scalar scalar const perma l 5.397E+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 0.000E+00 52 ymin REAL scalar scalar const perma l 0.000E+00 53 zmin REAL scalar scalar const perma l 0.000E+00 54 xmax REAL scalar scalar const perma l 1.000E+03 55 ymax REAL scalar scalar const perma l 1.000E+03 56 zmax REAL scalar scalar const perma l 6.000E+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 itetlev VINT scalar nelements linea tempo agl 0.000E+00 cmo select mo # PLANES for fault 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 # Expand fault width to capture elements of 50m spacing surface/sright/intrface/plane/720.0 0. 0./840.0 1000.0 550.0/840.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 6.25m # this sets refinement slightly larger than final set define XMIN 100. define XMAX 300. define YMIN 400. define YMAX 600. define ZMIN 287. define ZMAX 313. # Refine Level 1 (25 m) # Select element sets # Refine octree based on element sets define LEVEL 1 infile get_eltsets.mlgi # macro infile for finding fault and repo # make sure correct mo is current # make sure sets do not exist # Ignore errors on first level of refinement eltset/erefine/delete ERROR - THE NAME erefine DOES NOT EXIST ! eltset/eboth/delete ERROR - THE NAME eboth DOES NOT EXIST ! eltset/elevel/delete ERROR - THE NAME elevel DOES NOT EXIST ! eltset/ebox/delete ERROR - THE NAME ebox DOES NOT EXIST ! eltset/eflt/delete ERROR - THE NAME eflt DOES NOT EXIST ! region/rflt/delete REGION ERROR: delete IS AN INVALID OPERATOR, OR A NON-DEFINED SURFACE OR REGI ON error in command : region/rflt/delete pset/pbox/delete ERROR - THE NAME pbox DOES NOT EXIST ! pset/pflt/delete ERROR - THE NAME pflt DOES NOT EXIST ! region/rflt/ge sleft and le sright pset/pflt/region/rflt THE PSET pflt HAS 210 POINTS eltset/eflt/inclusive/pset,get,pflt THE ELTSET eflt HAS 440 ELEMENTS pset/pbox/geom/xyz/1,0,0/XMIN YMIN ZMIN/XMAX YMAX ZMAX THE PSET pbox HAS 25 POINTS eltset/ebox/inclusive/pset,get,pbox THE ELTSET ebox HAS 72 ELEMENTS eltset/eboth/union ebox eflt THE ELTSET eboth HAS 512 ELEMENTS # narrow refine area to current LEVEL eltset/elevel/itetlev/le/LEVEL THE ELTSET elevel HAS 4320 ELEMENTS eltset/erefine/inter eboth elevel THE ELTSET erefine HAS 512 ELEMENTS finish cmo/setatt/mo/itetclr/eltset,get,eflt 4 440 values reset for attribute itetclr cmo/setatt/mo/itetclr/eltset,get,ebox 3 72 values reset for attribute itetclr dump gmv tmp_lvl0_clrs.gmv mo cmo/modatt/-def-/-def-/ioflag/x finish refine/eltset/eltset,get,erefine pset/refine_eltset_pset_tmp_abc_xyz/eltset/erefine/ THE PSET refine_eltset_pset_tmp_abc_xyz HAS 945 POINTS finish refine/constant/imt1/linear/element/pset,get,refine_eltset_pset_tmp_abc_xyz/-1., 0.,0./exclusive/ refine elements: refine hex : 512 cmo/addatt/mo/itetpar/VINT/scalar/nelements/linear/permanent//0.0 finish cmo/addatt/mo/itetkid/VINT/scalar/nelements/linear/permanent//0.0 finish Element iteration: 1 512 Octree refined elems: old= 4320 new= 8416 Dudding 4949 duplicate points. filter/1,0,0 FILTER:Use internal epsilonl value = 0.341111394591E-09 FILTER:Dudding 0 duplicate points. geniee finish finish pset/refine_eltset_pset_tmp_abc_xyz/delete/ PSET refine_e DELETED finish dump gmv tmp_lvl1.gmv mo cmo/modatt/-def-/-def-/ioflag/x finish # Refine Level 2 (12.5 m) # Select element sets # Refine octree based on element sets define LEVEL 2 infile get_eltsets.mlgi # macro infile for finding fault and repo # make sure correct mo is current # make sure sets do not exist # Ignore errors on first level of refinement eltset/erefine/delete ELTSET erefine DELETED eltset/eboth/delete ELTSET eboth DELETED eltset/elevel/delete ELTSET elevel DELETED eltset/ebox/delete ELTSET ebox DELETED eltset/eflt/delete ELTSET eflt DELETED region/rflt/delete pset/pbox/delete PSET pbox DELETED pset/pflt/delete PSET pflt DELETED region/rflt/ge sleft and le sright pset/pflt/region/rflt THE PSET pflt HAS 1517 POINTS eltset/eflt/inclusive/pset,get,pflt THE ELTSET eflt HAS 2920 ELEMENTS pset/pbox/geom/xyz/1,0,0/XMIN YMIN ZMIN/XMAX YMAX ZMAX THE PSET pbox HAS 81 POINTS eltset/ebox/inclusive/pset,get,pbox THE ELTSET ebox HAS 272 ELEMENTS eltset/eboth/union ebox eflt THE ELTSET eboth HAS 3192 ELEMENTS # narrow refine area to current LEVEL eltset/elevel/itetlev/le/LEVEL THE ELTSET elevel HAS 8416 ELEMENTS eltset/erefine/inter eboth elevel THE ELTSET erefine HAS 3192 ELEMENTS finish cmo/setatt/mo/itetclr/eltset,get,eflt 4 2920 values reset for attribute itetclr cmo/setatt/mo/itetclr/eltset,get,ebox 3 272 values reset for attribute itetclr dump gmv tmp_lvl1_clrs.gmv mo cmo/modatt/-def-/-def-/ioflag/x finish refine/eltset/eltset,get,erefine pset/refine_eltset_pset_tmp_abc_xyz/eltset/erefine/ THE PSET refine_eltset_pset_tmp_abc_xyz HAS 4551 POINTS finish refine/constant/imt1/linear/element/pset,get,refine_eltset_pset_tmp_abc_xyz/-1., 0.,0./exclusive/ refine elements: refine hex : 3192 Element iteration: 1 2680 Octree refined elems: old= 8416 new= 29856 Dudding 27983 duplicate points. filter/1,0,0 FILTER:Use internal epsilonl value = 0.341111394591E-09 FILTER:Dudding 0 duplicate points. geniee finish Element iteration: 1 220 Octree refined elems: old= 29856 new= 31616 Dudding 2506 duplicate points. filter/1,0,0 FILTER:Use internal epsilonl value = 0.341111394591E-09 FILTER:Dudding 40 duplicate points. geniee finish finish pset/refine_eltset_pset_tmp_abc_xyz/delete/ PSET refine_e DELETED finish dump gmv tmp_lvl2.gmv mo cmo/modatt/-def-/-def-/ioflag/x finish # Refine Level 3 (6.25 m) # Select element sets # Refine octree based on element sets # reset surface back to thinner width cmo select mo surface/sright/remove 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 define LEVEL 3 infile get_eltsets.mlgi # macro infile for finding fault and repo # make sure correct mo is current # make sure sets do not exist # Ignore errors on first level of refinement eltset/erefine/delete ELTSET erefine DELETED eltset/eboth/delete ELTSET eboth DELETED eltset/elevel/delete ELTSET elevel DELETED eltset/ebox/delete ELTSET ebox DELETED eltset/eflt/delete ELTSET eflt DELETED region/rflt/delete pset/pbox/delete PSET pbox DELETED pset/pflt/delete PSET pflt DELETED region/rflt/ge sleft and le sright pset/pflt/region/rflt THE PSET pflt HAS 8829 POINTS eltset/eflt/inclusive/pset,get,pflt THE ELTSET eflt HAS 15280 ELEMENTS pset/pbox/geom/xyz/1,0,0/XMIN YMIN ZMIN/XMAX YMAX ZMAX THE PSET pbox HAS 867 POINTS eltset/ebox/inclusive/pset,get,pbox THE ELTSET ebox HAS 1568 ELEMENTS eltset/eboth/union ebox eflt THE ELTSET eboth HAS 16848 ELEMENTS # narrow refine area to current LEVEL eltset/elevel/itetlev/le/LEVEL THE ELTSET elevel HAS 31616 ELEMENTS eltset/erefine/inter eboth elevel THE ELTSET erefine HAS 16848 ELEMENTS finish cmo/setatt/mo/itetclr/eltset,get,eflt 4 15280 values reset for attribute itetclr cmo/setatt/mo/itetclr/eltset,get,ebox 3 1568 values reset for attribute itetclr dump gmv tmp_lvl2_clrs.gmv mo cmo/modatt/-def-/-def-/ioflag/x finish refine/eltset/eltset,get,erefine pset/refine_eltset_pset_tmp_abc_xyz/eltset/erefine/ THE PSET refine_eltset_pset_tmp_abc_xyz HAS 21127 POINTS finish refine/constant/imt1/linear/element/pset,get,refine_eltset_pset_tmp_abc_xyz/-1., 0.,0./exclusive/ refine elements: refine hex : 16888 Element iteration: 1 14096 Octree refined elems: old= 31616 new= 144384 Dudding 153411 duplicate points. filter/1,0,0 FILTER:Use internal epsilonl value = 0.341111394591E-09 FILTER:Dudding 0 duplicate points. geniee finish Element iteration: 1 720 Octree refined elems: old= 144384 new= 150144 Dudding 8278 duplicate points. filter/1,0,0 FILTER:Use internal epsilonl value = 0.341111394591E-09 FILTER:Dudding 80 duplicate points. geniee finish Element iteration: 1 152 Octree refined elems: old= 150144 new= 151360 Dudding 1677 duplicate points. filter/1,0,0 FILTER:Use internal epsilonl value = 0.341111394591E-09 FILTER:Dudding 0 duplicate points. geniee finish finish pset/refine_eltset_pset_tmp_abc_xyz/delete/ PSET refine_e DELETED finish dump gmv hex_lvl3_octree.gmv mo cmo/modatt/-def-/-def-/ioflag/x finish # Clean up attributes for final hex mesh grid2grid/tree_to_fe/mo_clean/mo cmo/copy/mo_clean/mo finish eltset/-parentset-/itetkid/gt/0 THE ELTSET -parentset- HAS 18380 ELEMENTS finish rmpoint/element/eltset,get,-parentset- Dudded 199071 points plus parents, and 18380 elements removed. finish rmpoint/compress 199071 points removed and 0 elements removed. The new point count is 155546. The new element count is 132980. finish resetpts/itp geniee finish finish eltset/-parentset-/delete ELTSET -parentset- DELETED 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: 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. : 132980 min aspect ratio = 0.1000E+01 max aspect ratio = 0.1000E+01 epsilonvol: 1.3322676E-04 --------------------------------------- element volumes b/w 0.2441E+03 and 0.8501E+03: 112768 element volumes b/w 0.8501E+03 and 0.2960E+04: 13104 element volumes b/w 0.2960E+04 and 0.1031E+05: 0 element volumes b/w 0.1031E+05 and 0.3590E+05: 3672 element volumes b/w 0.3590E+05 and 0.1250E+06: 3436 min volume = 2.4414062E+02 max volume = 1.2500000E+05 ----------------------------------------------------------- 132980 total elements evaluated. finish cmo/select/mo_clean rmpoint/compress 0 points removed and 0 elements removed. The new point count is 155546. The new element count is 132980. resetpts/itp geniee finish infile get_eltsets.mlgi # macro infile for finding fault and repo # make sure correct mo is current # make sure sets do not exist # Ignore errors on first level of refinement eltset/erefine/delete ELTSET erefine DELETED eltset/eboth/delete ELTSET eboth DELETED eltset/elevel/delete ELTSET elevel DELETED eltset/ebox/delete ELTSET ebox DELETED eltset/eflt/delete ELTSET eflt DELETED region/rflt/delete pset/pbox/delete PSET pbox DELETED pset/pflt/delete PSET pflt DELETED region/rflt/ge sleft and le sright pset/pflt/region/rflt THE PSET pflt HAS 69069 POINTS eltset/eflt/inclusive/pset,get,pflt THE ELTSET eflt HAS 84960 ELEMENTS pset/pbox/geom/xyz/1,0,0/XMIN YMIN ZMIN/XMAX YMAX ZMAX THE PSET pbox HAS 5445 POINTS eltset/ebox/inclusive/pset,get,pbox THE ELTSET ebox HAS 6936 ELEMENTS eltset/eboth/union ebox eflt THE ELTSET eboth HAS 91896 ELEMENTS # narrow refine area to current LEVEL eltset/elevel/itetlev/le/LEVEL THE ELTSET elevel HAS 132980 ELEMENTS eltset/erefine/inter eboth elevel THE ELTSET erefine HAS 91896 ELEMENTS finish cmo/setatt/mo_clean/itetclr/eltset,get,eflt 4 84960 values reset for attribute itetclr cmo/setatt/mo_clean/itetclr/eltset,get,ebox 3 6936 values reset for attribute itetclr ## WRITE FINAL LEVEL 3 (6.25 m) BASE MESH cmo/setatt/mo_clean/ipolydat no 1 values reset for attribute ipolydat dump/avs/hex_lvl3.inp/mo_clean/1 1 0 0 cmo/modatt/-def-/-def-/ioflag/x finish dump/hex_lvl3.gmv/mo_clean dump/gmv/hex_lvl3.gmv/mo_clean cmo/modatt/-def-/-def-/ioflag/x finish finish cmo printatt mo_clean -all- minmax ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH -def- 0.000000000E+00 0.000000000E+00 0.000000000E+00 155546 scalar 1 1 0 1 vector 3 3 0 1 nnodes 155546 155546 0 1 nedges 0 0 0 1 nfaces 0 0 0 1 nelements 132980 132980 0 1 mbndry 16000000 16000000 0 1 ndimensions_topo 3 3 0 1 ndimensions_geom 3 3 0 1 nodes_per_element 8 8 0 1 edges_per_element 12 12 0 1 faces_per_element 6 6 0 1 isetwd 0 2 2 155546 ialias 0 354617 354617 155546 imt1 1 1 0 155546 itp1 0 12 12 155546 icr1 0 0 0 155546 isn1 0 0 0 155546 xic 0.000000000E+00 1.000000000E+03 1.000000000E+03 155546 yic 0.000000000E+00 1.000000000E+03 1.000000000E+03 155546 zic 0.000000000E+00 6.000000000E+02 6.000000000E+02 155546 xtetwd 8 30 22 132980 itetclr 1 4 3 132980 itettyp 8 8 0 132980 itetoff 0 1063832 1063832 132980 jtetoff 0 797874 797874 132980 itet 1 155546 155545 132980x8 jtet 2 16796891 16796889 132980x6 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 155546 155546 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 2 2 0 1 number_of_eltsets 0 0 0 1 number_of_fsets 0 0 0 1 itetlev 0 3 3 132980 itetpar 0 31136 31136 132980 itetkid 0 0 0 132980 # check point distribution of top to fault refine area cmo select mo_clean pset/psub/attribute zic/1,0,0 lt 400. THE PSET psub HAS 115606 POINTS rmpoint pset,get,psub Dudded 113205 points plus parents, plus 98332 elements removed. rmpoint compress 113205 points removed and 0 elements removed. The new point count is 42341. The new element count is 34648. read avs tet_top.inp motop geniee finish cmo/status/brief The current-mesh-object(CMO) is: motop 1 Mesh Object name: mo number of nodes = 354617 number of elements = 151360 dimensions geometry = 3 element type = hex dimensions topology = 3 8 nodes 6 faces 12 edges boundary flag = 16000000 status = inactive 2 Mesh Object name: mo_clean number of nodes = 42341 number of elements = 34648 dimensions geometry = 3 element type = hex dimensions topology = 3 8 nodes 6 faces 12 edges boundary flag = 16000000 status = inactive 3 Mesh Object name: motop number of nodes = 5330 number of elements = 23954 dimensions geometry = 3 element type = tet dimensions topology = 3 4 nodes 4 faces 6 edges boundary flag = 16000000 status = active finish addmesh append mov mo_clean motop searching for mov attributes in motop attribute not found: 66 itetlev attribute not found: 67 itetpar attribute not found: 68 itetkid ADDMESH: 3 attributes not copied from motop to mov dump gmv merge_octree_top.gmv mov cmo/modatt/-def-/-def-/ioflag/x finish finish LaGriT successfully completed