* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * Program: LaGriT V3.101 Linux m64 * * * date_compile: 2013/04/25 RH gf4.5 * * * Run Time: 2013/Jun 17 10:38:47 * * * 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 # Create the top and bottom piece for mesh ################################# # Extents for full box # define X0 0. define X1 1000. define Y0 0. define Y1 1000. define Z0 0. define Z1 600. # Resolution for full box 25m define NX 41 define NY 41 define NZ 25 # Resolution for full box 50m define NX 21 define NY 21 define NZ 13 # # For top cutout for slope # Include a buffer of 25m elements # Do not include refined fault with top define Xcut 250. define Xcut2 700. define Zcut 500. define Zcut2 550. # # ####################################### # Make BASE PIECE # Create the brick mesh # We will use all but the line boundary # and top area below boundary # for the points of the mesh ####################################### # evenly space points along x,y and z axis cmo/create/mohex///hex createpts/brick/xyz/NX NY NZ/X0 Y0 Z0/X1 Y1 Z1/1,1,1/1,1,0/1. 1. 1. Number of nodes: 5733 Number of elements: 4800 Number of negative volume elements: 0 Total volume: 0.600000000000E+09 cmo/setatt/mohex/imt1/1,0,0/1 5733 values reset for attribute imt1 cmo/setatt/mohex/itetclr/1,0,0/1 4800 values reset for attribute itetclr filter 1,0,0 FILTER:Use internal epsilonl value = 0.341111394591E-09 FILTER:Dudding 0 duplicate points. rmpoint compress 0 points removed and 0 elements removed. The new point count is 5733. The new element count is 4800. resetpts itp geniee finish # cut out the place where the slope will go cmo select mohex pset/pxduds/attribute xic/1,0,0/gt Xcut THE PSET pxduds HAS 4095 POINTS pset/px2duds/attribute xic/1,0,0/lt Xcut2 THE PSET px2duds HAS 3822 POINTS pset/pzduds/attribute zic/1,0,0/gt Zcut THE PSET pzduds HAS 882 POINTS pset/pduds/inter pxduds px2duds pzduds THE PSET pduds HAS 336 POINTS eltset/eduds/inclusive pset,get,pduds THE ELTSET eduds HAS 360 ELEMENTS cmo setatt mohex itetclr eltset,get,eduds 2 360 values reset for attribute itetclr dump gmv tmp_box_full.gmv mohex cmo/modatt/-def-/-def-/ioflag/x finish rmpoint element eltset,get,eduds Dudded 336 points plus parents, and 360 elements removed. rmpoint compress 336 points removed and 0 elements removed. The new point count is 5397. The new element count is 4440. resetpts itp geniee finish # remove base top that will align with top slope cmo select mohex pset/pz2/attribute zic/1,0,0/gt Zcut2 THE PSET pz2 HAS 273 POINTS pset/px2/attribute xic/1,0,0/gt Xcut2 THE PSET px2 HAS 1638 POINTS pset/pduds2/inter pz2 px2 THE PSET pduds2 HAS 126 POINTS eltset/eduds2/inclusive pset,get,pduds2 THE ELTSET eduds2 HAS 120 ELEMENTS rmpoint element eltset,get,eduds2 Dudded 147 points plus parents, and 120 elements removed. cmo setatt mohex ipolydat no 1 values reset for attribute ipolydat dump gmv hex_base.gmv mohex cmo/modatt/-def-/-def-/ioflag/x finish dump avs hex_base.inp mohex cmo/modatt/-def-/-def-/ioflag/x finish ####################################### # Make TOP PIECE # This is from the tet mesh 25m version # make sure tip near x=700 is kept read avs tet_top_25m.inp motop geniee finish cmo/status/brief The current-mesh-object(CMO) is: motop 1 Mesh Object name: mohex 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 = inactive 2 Mesh Object name: motop number of nodes = 704630 number of elements = 4099282 dimensions geometry = 3 element type = tet dimensions topology = 3 4 nodes 4 faces 6 edges boundary flag = 25000000 status = active finish cmo printatt motop -xyz- minmax ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH xic 0.000000000E+00 1.000000000E+03 1.000000000E+03 704630 yic 0.000000000E+00 1.000000000E+03 1.000000000E+03 704630 zic 0.000000000E+00 6.000000000E+02 6.000000000E+02 704630 cmo select motop pset/px/attribute xic/1,0,0/ge Xcut THE PSET px HAS 672606 POINTS pset/px2/attribute xic/1,0,0/lt 701. THE PSET px2 HAS 130805 POINTS pset/pz/attribute zic/1,0,0/ge Zcut THE PSET pz HAS 68779 POINTS pset/pkeep/inter px px2 pz THE PSET pkeep HAS 5330 POINTS pset/pduds/not pkeep THE PSET pduds HAS 699300 POINTS eltset/eduds/inclusive pset,get,pduds THE ELTSET eduds HAS 4075328 ELEMENTS cmo setatt motop itetclr eltset,get,eduds 2 4075328 values reset for attribute itetclr dump gmv tmp_box_full.gmv motop cmo/modatt/-def-/-def-/ioflag/x finish cmo select motop rmpoint element eltset,get,eduds Dudded 699300 points plus parents, and 4075328 elements removed. rmpoint compress 699300 points removed and 0 elements removed. The new point count is 5330. The new element count is 23954. resetpts itp geniee finish cmo setatt motop itetclr 1 23954 values reset for attribute itetclr cmo setatt motop imt 1 5330 values reset for attribute imt1 cmo setatt motop ipolydat no 1 values reset for attribute ipolydat dump gmv tet_top.gmv motop cmo/modatt/-def-/-def-/ioflag/x finish dump avs tet_top.inp motop cmo/modatt/-def-/-def-/ioflag/x finish ### # merge for viewing addmesh/append/mov/mohex motop cmo setatt mov ipolydat no 1 values reset for attribute ipolydat dump gmv merge_base_top.gmv mov cmo/modatt/-def-/-def-/ioflag/x finish finish LaGriT successfully completed