* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * Program: LaGriT V3.101 Linux m64 * * * date_compile: 2013/04/25 RH gf4.5 * * * Run Time: 2013/Jun 19 08:57:39 * * * 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 25m # refine halves are 12.5 6.25 3.125 1.5625 # # Materials are Background 5, fault 4, repo 3 # Zones are Top left 1, Top right 2 # # read and add attribute for use with LEVEL values 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 = 39729 number of elements = 36160 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 = 39729 number of elements = 36160 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 3.973E+04 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 3.616E+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 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 3.973E+04 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 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 3.125m # 300 +- 12.5m define XMIN 100. define XMAX 300. define YMIN 400. define YMAX 600. define ZMIN 287. define ZMAX 313. # Refine Level 1 (12.5 m) # Select element sets # Refine octree based on element sets define LEVEL 1 infile get_eltsets_max.mlgi # macro infile for finding fault and repo # this version is for coarse mesh # make sure correct mo is current # make sure sets do not exist eltset/erefine/delete ERROR - THE NAME erefine 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/p1/region/rflt THE PSET p1 HAS 779 POINTS eltset/e1/inclusive/pset,get,p1 THE ELTSET e1 HAS 1760 ELEMENTS pset/p2/eltset e1 THE PSET p2 HAS 2993 POINTS eltset/eflt/inclusive/pset,get,p2 THE ELTSET eflt HAS 3840 ELEMENTS pset/p1/delete PSET p1 DELETED pset/p2/delete PSET p2 DELETED eltset/e1/delete ELTSET e1 DELETED pset/p1/geom/xyz/1,0,0/XMIN YMIN ZMIN/XMAX YMAX ZMAX THE PSET p1 HAS 81 POINTS eltset/e1/inclusive/pset,get,p1 THE ELTSET e1 HAS 200 ELEMENTS pset/p2/eltset e1 THE PSET p2 HAS 363 POINTS eltset/ebox/inclusive/pset,get,p2 THE ELTSET ebox HAS 576 ELEMENTS pset/p1/delete PSET p1 DELETED pset/p2/delete PSET p2 DELETED eltset/e1/delete ELTSET e1 DELETED eltset/erefine/union ebox eflt THE ELTSET erefine HAS 4416 ELEMENTS finish cmo/setatt/mo/itetclr/eltset,get,eflt 4 3840 values reset for attribute itetclr cmo/setatt/mo/itetclr/eltset,get,ebox 3 576 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 6052 POINTS finish refine/constant/imt1/linear/element/pset,get,refine_eltset_pset_tmp_abc_xyz/-1., 0.,0./exclusive/ refine elements: refine hex : 4416 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 4416 Octree refined elems: old= 36160 new= 71488 Dudding 48286 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 (6.25 m) # Select element sets # Refine octree based on element sets # move surfaces out just a tad to capture coarse elements 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 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 ERROR - THE NAME eboth DOES NOT EXIST ! eltset/elevel/delete ERROR - THE NAME elevel DOES NOT EXIST ! 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 14920 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 1496 ELEMENTS eltset/eboth/union ebox eflt THE ELTSET eboth HAS 16416 ELEMENTS # narrow refine area to current LEVEL eltset/elevel/itetlev/le/LEVEL THE ELTSET elevel HAS 71488 ELEMENTS eltset/erefine/inter eboth elevel THE ELTSET erefine HAS 16416 ELEMENTS finish cmo/setatt/mo/itetclr/eltset,get,eflt 4 14920 values reset for attribute itetclr cmo/setatt/mo/itetclr/eltset,get,ebox 3 1496 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 20690 POINTS finish refine/constant/imt1/linear/element/pset,get,refine_eltset_pset_tmp_abc_xyz/-1., 0.,0./exclusive/ refine elements: refine hex : 16416 Element iteration: 1 14096 Octree refined elems: old= 71488 new= 184256 Dudding 153411 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_lvl2.gmv mo cmo/modatt/-def-/-def-/ioflag/x finish # Refine Level 3 (3.125 m) # Select element sets # Refine octree based on element sets define LEVEL 3 infile get_eltsets_min.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/exclusive/pset,get,pflt THE ELTSET eflt HAS 54960 ELEMENTS 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 4608 ELEMENTS eltset/eboth/union ebox eflt THE ELTSET eboth HAS 59568 ELEMENTS # narrow refine area to current LEVEL eltset/elevel/itetlev/le/LEVEL THE ELTSET elevel HAS 184256 ELEMENTS eltset/erefine/inter eboth elevel THE ELTSET erefine HAS 59568 ELEMENTS finish cmo/setatt/mo/itetclr/eltset,get,eflt 4 54960 values reset for attribute itetclr cmo/setatt/mo/itetclr/eltset,get,ebox 3 4608 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 74353 POINTS finish refine/constant/imt1/linear/element/pset,get,refine_eltset_pset_tmp_abc_xyz/-1., 0.,0./exclusive/ refine elements: refine hex : 59568 Element iteration: 1 54816 Octree refined elems: old= 184256 new= 622784 Dudding 599863 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/motmp/mo cmo/copy/motmp/mo finish eltset/-parentset-/itetkid/gt/0 THE ELTSET -parentset- HAS 73328 ELEMENTS finish rmpoint/element/eltset,get,-parentset- Dudded 801560 points plus parents, and 73328 elements removed. finish rmpoint/compress 801560 points removed and 0 elements removed. The new point count is 631401. The new element count is 549456. 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. : 549456 min aspect ratio = 0.1000E+01 max aspect ratio = 0.1000E+01 epsilonvol: 1.3322676E-04 --------------------------------------- element volumes b/w 0.3052E+02 and 0.1063E+03: 438528 element volumes b/w 0.1063E+03 and 0.3700E+03: 57952 element volumes b/w 0.3700E+03 and 0.1289E+04: 0 element volumes b/w 0.1289E+04 and 0.4487E+04: 21232 element volumes b/w 0.4487E+04 and 0.1562E+05: 31744 min volume = 3.0517578E+01 max volume = 1.5625000E+04 ----------------------------------------------------------- 549456 total elements evaluated. finish cmo/select/motmp rmpoint/compress 0 points removed and 0 elements removed. The new point count is 631401. The new element count is 549456. 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 478130 POINTS eltset/eflt/inclusive/pset,get,pflt THE ELTSET eflt HAS 440000 ELEMENTS pset/pbox/geom/xyz/1,0,0/XMIN YMIN ZMIN/XMAX YMAX ZMAX THE PSET pbox HAS 38025 POINTS eltset/ebox/inclusive/pset,get,pbox THE ELTSET ebox HAS 35608 ELEMENTS eltset/eboth/union ebox eflt THE ELTSET eboth HAS 475608 ELEMENTS # narrow refine area to current LEVEL eltset/elevel/itetlev/le/LEVEL THE ELTSET elevel HAS 549456 ELEMENTS eltset/erefine/inter eboth elevel THE ELTSET erefine HAS 475608 ELEMENTS finish cmo/setatt/motmp/itetclr/eltset,get,eflt 4 440000 values reset for attribute itetclr cmo/setatt/motmp/itetclr/eltset,get,ebox 3 35608 values reset for attribute itetclr ## FINAL LEVEL 3 (3.125 m) BASE MESH cmo/setatt/motmp/ipolydat no 1 values reset for attribute ipolydat dump/avs/hex_lvl3.inp/motmp/1 1 0 0 cmo/modatt/-def-/-def-/ioflag/x finish dump/hex_lvl3.gmv/motmp dump/gmv/hex_lvl3.gmv/motmp cmo/modatt/-def-/-def-/ioflag/x finish finish cmo printatt motmp -all- minmax ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH -def- 0.000000000E+00 0.000000000E+00 0.000000000E+00 631401 scalar 1 1 0 1 vector 3 3 0 1 nnodes 631401 631401 0 1 nedges 0 0 0 1 nfaces 0 0 0 1 nelements 549456 549456 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 631401 ialias 0 1432961 1432961 631401 imt1 1 1 0 631401 itp1 0 10 10 631401 icr1 0 0 0 631401 isn1 0 0 0 631401 xic 0.000000000E+00 1.000000000E+03 1.000000000E+03 631401 yic 0.000000000E+00 1.000000000E+03 1.000000000E+03 631401 zic 0.000000000E+00 6.000000000E+02 6.000000000E+02 631401 xtetwd 8 30 22 549456 itetclr 1 4 3 549456 itettyp 8 8 0 549456 itetoff 0 4395640 4395640 549456 jtetoff 0 3296730 3296730 549456 itet 1 631401 631400 549456x8 jtet 2 16000000 15999998 549456x6 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 631401 631401 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 549456 itetpar 0 184224 184224 549456 itetkid 0 0 0 549456 cmo/delete/motmp Released Mesh Object: motmp # Refine once more in case higher resolution needed cmo select mo # Refine Level 4 (1.5625) # Select element sets # Refine octree based on element sets define LEVEL 4 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 478130 POINTS eltset/eflt/inclusive/pset,get,pflt THE ELTSET eflt HAS 505640 ELEMENTS pset/pbox/geom/xyz/1,0,0/XMIN YMIN ZMIN/XMAX YMAX ZMAX THE PSET pbox HAS 38025 POINTS eltset/ebox/inclusive/pset,get,pbox THE ELTSET ebox HAS 41200 ELEMENTS eltset/eboth/union ebox eflt THE ELTSET eboth HAS 546840 ELEMENTS # narrow refine area to current LEVEL eltset/elevel/itetlev/le/LEVEL THE ELTSET elevel HAS 622784 ELEMENTS eltset/erefine/inter eboth elevel THE ELTSET erefine HAS 546840 ELEMENTS finish cmo/setatt/mo/itetclr/eltset,get,eflt 4 505640 values reset for attribute itetclr cmo/setatt/mo/itetclr/eltset,get,ebox 3 41200 values reset for attribute itetclr dump gmv tmp_lvl3_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 560762 POINTS finish refine/constant/imt1/linear/element/pset,get,refine_eltset_pset_tmp_abc_xyz/-1., 0.,0./exclusive/ refine elements: refine hex : 546840 Element iteration: 1 475608 Octree refined elems: old= 622784 new= 4427648 Dudding5473313 duplicate points. filter/1,0,0 FILTER:Use internal epsilonl value = 0.341111394591E-09 FILTER:Dudding 0 duplicate points. geniee finish Element iteration: 1 3440 Octree refined elems: old= 4427648 new= 4455168 Dudding 40270 duplicate points. filter/1,0,0 FILTER:Use internal epsilonl value = 0.341111394591E-09 FILTER:Dudding 640 duplicate points. geniee finish finish pset/refine_eltset_pset_tmp_abc_xyz/delete/ PSET refine_e DELETED finish dump gmv hex_lvl4_octree.gmv mo cmo/modatt/-def-/-def-/ioflag/x finish # Clean up attributes for final hex mesh grid2grid/tree_to_fe/motmp/mo cmo/copy/motmp/mo finish eltset/-parentset-/itetkid/gt/0 THE ELTSET -parentset- HAS 552376 ELEMENTS finish rmpoint/element/eltset,get,-parentset- Dudded 6315783 points plus parents, and 552376 elements removed. finish rmpoint/compress 6315783 points removed and 0 elements removed. The new point count is 4219090. The new element count is 3902792. 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. : 3902792 min aspect ratio = 0.1000E+01 max aspect ratio = 0.1000E+01 epsilonvol: 1.3322676E-04 --------------------------------------- element volumes b/w 0.3815E+01 and 0.2013E+02: 3508224 element volumes b/w 0.2013E+02 and 0.1063E+03: 296640 element volumes b/w 0.1063E+03 and 0.5609E+03: 48392 element volumes b/w 0.5609E+03 and 0.2960E+04: 17792 element volumes b/w 0.2960E+04 and 0.1562E+05: 31744 min volume = 3.8146973E+00 max volume = 1.5625000E+04 ----------------------------------------------------------- 3902792 total elements evaluated. finish cmo/select/motmp rmpoint/compress 0 points removed and 0 elements removed. The new point count is 4219090. The new element count is 3902792. 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 3602527 POINTS eltset/eflt/inclusive/pset,get,pflt THE ELTSET eflt HAS 3450560 ELEMENTS pset/pbox/geom/xyz/1,0,0/XMIN YMIN ZMIN/XMAX YMAX ZMAX THE PSET pbox HAS 282897 POINTS eltset/ebox/inclusive/pset,get,pbox THE ELTSET ebox HAS 272936 ELEMENTS eltset/eboth/union ebox eflt THE ELTSET eboth HAS 3723496 ELEMENTS # narrow refine area to current LEVEL eltset/elevel/itetlev/le/LEVEL THE ELTSET elevel HAS 3902792 ELEMENTS eltset/erefine/inter eboth elevel THE ELTSET erefine HAS 3723496 ELEMENTS finish cmo/setatt/motmp/itetclr/eltset,get,eflt 4 3450560 values reset for attribute itetclr cmo/setatt/motmp/itetclr/eltset,get,ebox 3 272936 values reset for attribute itetclr ## FINAL LEVEL Level 4 (1.5625) BASE MESH cmo/setatt/motmp/ipolydat no 1 values reset for attribute ipolydat dump/avs/hex_lvl4.inp/motmp/1 1 0 0 cmo/modatt/-def-/-def-/ioflag/x finish dump/hex_lvl4.gmv/motmp dump/gmv/hex_lvl4.gmv/motmp cmo/modatt/-def-/-def-/ioflag/x finish finish cmo printatt motmp -all- minmax ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH -def- 0.000000000E+00 0.000000000E+00 0.000000000E+00 4219090 scalar 1 1 0 1 vector 3 3 0 1 nnodes 4219090 4219090 0 1 nedges 0 0 0 1 nfaces 0 0 0 1 nelements 3902792 3902792 0 1 mbndry 32000000 32000000 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 4219090 ialias 0 10534873 10534873 4219090 imt1 1 1 0 4219090 itp1 0 10 10 4219090 icr1 0 0 0 4219090 isn1 0 0 0 4219090 xic 0.000000000E+00 1.000000000E+03 1.000000000E+03 4219090 yic 0.000000000E+00 1.000000000E+03 1.000000000E+03 4219090 zic 0.000000000E+00 6.000000000E+02 6.000000000E+02 4219090 xtetwd 8 30 22 3902792 itetclr 1 4 3 3902792 itettyp 8 8 0 3902792 itetoff 0 31222328 31222328 3902792 jtetoff 0 23416746 23416746 3902792 itet 1 4219090 4219089 3902792x8 jtet 2 32000000 31999998 3902792x6 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 4219090 4219090 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 4 4 3902792 itetpar 0 622784 622784 3902792 itetkid 0 0 0 3902792 cmo/delete/motmp Released Mesh Object: motmp finish LaGriT successfully completed