* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * Program: LaGriT V3.101 Linux m64 * * * date_compile: 2013/04/25 RH gf4.5 * * * Run Time: 2013/Jun 19 16:48:55 * * * 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 # Read final tet mesh created for FEHM modeling # Define mesh cell materials and outside zones for exodus file # Conventions were setup by ealier meshing so there is no material 1 or 2 # # FEHM zone id numbers are changed to fit Exodus conventions # Exodus is based on element (cell) colors in itetclr # and has boundaries defined as outside surfaces (tri face of tet) # Note the node imt zone values are unchanged since Exodus does not use them. # Only element itetclr and facesets are used. # # For Exodus, we have 3 mesh materials plus top facesets (side sets) # 1 = background (fehm zone 5) # 2 = fault (fehm zone 4) # 3 = Repo, same as fehm zone # # Faceset integer id is assigned by order they are defined: # 1 = Top Left and middle (fehm zone 1) # 2 = Top Right (fehm zone 2) # 3 = Top Fault (fehm zone overwrites top fault to 4) # read avs tet_final_zones.inp mo cmo/addatt/mo/top/VINT/scalar/nnodes/linear/permanent/gxaf/0.0 finish cmo/addatt/mo/bottom/VINT/scalar/nnodes/linear/permanent/gxaf/0.0 finish cmo/addatt/mo/left_w/VINT/scalar/nnodes/linear/permanent/gxaf/0.0 finish cmo/addatt/mo/right_e/VINT/scalar/nnodes/linear/permanent/gxaf/0.0 finish cmo/addatt/mo/back_n/VINT/scalar/nnodes/linear/permanent/gxaf/0.0 finish cmo/addatt/mo/front_s/VINT/scalar/nnodes/linear/permanent/gxaf/0.0 finish geniee finish cmo/status/brief The current-mesh-object(CMO) is: mo 1 Mesh Object name: mo number of nodes = 160583 number of elements = 921085 dimensions geometry = 3 element type = tet dimensions topology = 3 4 nodes 4 faces 6 edges boundary flag = 16000000 status = active finish cmo printatt mo -all- minmax ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH -def- 0.000000000E+00 0.000000000E+00 0.000000000E+00 160583 scalar 1 1 0 1 vector 3 3 0 1 nnodes 160583 160583 0 1 nedges 0 0 0 1 nfaces 0 0 0 1 nelements 921085 921085 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 0 0 160583 ialias 0 0 0 160583 imt1 1 5 4 160583 itp1 0 12 12 160583 icr1 0 0 0 160583 isn1 0 0 0 160583 xic 0.000000000E+00 1.000000000E+03 1.000000000E+03 160583 yic 0.000000000E+00 1.000000000E+03 1.000000000E+03 160583 zic 0.000000000E+00 6.000000000E+02 6.000000000E+02 160583 xtetwd 0 0 0 921085 itetclr 1 5 4 921085 itettyp 5 5 0 921085 itetoff 0 3684336 3684336 921085 jtetoff 0 3684336 3684336 921085 itet 1 160583 160582 921085x4 jtet 1 19684194 19684193 921085x4 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 160583 160583 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 0 0 0 1 number_of_eltsets 0 0 0 1 number_of_fsets 0 0 0 1 top -5 1 6 160583 bottom -5 2 7 160583 left_w -5 3 8 160583 right_e -5 5 10 160583 back_n -5 6 11 160583 front_s -5 4 9 160583 ### CELL MATERIAL COLORS ################## # for Exodus set cells (tets) using node imt colors # Node colors will not be used. # Repository should be about 20 vertical # set everything to 1, then overwrite with fault and repo cmo select mo cmo setatt mo itetclr 1 921085 values reset for attribute itetclr cmo printatt mo imt minmax ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH imt1 1 5 4 160583 # repository stays 3 # use exclusive # inclusive selects too many elements making # the repository more than 20 m high pset/p3/attribute imt/1,0,0/eq 3 THE PSET p3 HAS 5445 POINTS eltset/e3/exclusive pset,get,p3 THE ELTSET e3 HAS 24479 ELEMENTS cmo setatt mo itetclr eltset,get,e3 3 24479 values reset for attribute itetclr # fault find imt 4, set to 2 pset/p4/attribute imt/1,0,0/eq 4 THE PSET p4 HAS 46207 POINTS eltset/e4/inclusive pset,get,p4 THE ELTSET e4 HAS 365103 ELEMENTS cmo setatt mo itetclr eltset,get,e4 2 365103 values reset for attribute itetclr # write element sets for each material eltset/e_mat3/itetclr/eq/3 THE ELTSET e_mat3 HAS 24479 ELEMENTS eltset/e_mat2/itetclr/eq/2 THE ELTSET e_mat2 HAS 365103 ELEMENTS eltset/e_mat1/itetclr/eq/1 THE ELTSET e_mat1 HAS 531503 ELEMENTS eltset/e_mat3/write/output_mat_3.cellset/ascii THERE ARE 5 ELTSETS DEFINED e3 e4 e_mat3 e_mat2 e_mat1 ELTSET: OUTPUT e_mat3 1 ELTSETS TO FILE eltset/e_mat2/write/output_mat_2.cellset/ascii THERE ARE 5 ELTSETS DEFINED e3 e4 e_mat3 e_mat2 e_mat1 ELTSET: OUTPUT e_mat2 1 ELTSETS TO FILE eltset/e_mat1/write/output_mat_1.cellset/ascii THERE ARE 5 ELTSETS DEFINED e3 e4 e_mat3 e_mat2 e_mat1 ELTSET: OUTPUT e_mat1 1 ELTSETS TO FILE dump gmv tet_colors.gmv mo cmo/modatt/-def-/-def-/ioflag/x finish ###### SORT CELLS for Exodus ################# # use sort to order element blocks as exodus will order # if this is not done, lagrit faceset numbers will not # correlate to exodus faceset numbers # itetclr must be ordered correctly # additional sort keys can be based on user preference # sort based on element itetclr number and median location createpts/median cmo/addatt//xmed/VDOUBLE/scalar/nelements///gal/0. finish cmo/addatt//ymed/VDOUBLE/scalar/nelements///gal/0. finish cmo/addatt//zmed/VDOUBLE/scalar/nelements///gal/0. finish sort/mo/index/ascending/ikey/itetclr xmed ymed zmed cmo/addatt/mo/ikey/vint/scalar/nelements///gax/0 finish SORT: order key written to attribute: ikey reorder/mo/ikey geniee finish dump/gmv/out_tmp_sort1.gmv/mo cmo/modatt/-def-/-def-/ioflag/x finish cmo/DELATT/mo/ikey # sort nodes based on mesh convention z, y, then x sort/mo/index/ascending/ikey/zic yic xic cmo/addatt/mo/ikey/vint/scalar/nnodes///gax/0 finish SORT: order key written to attribute: ikey reorder/mo/ikey dump/gmv/out_tmp_sort2.gmv/mo cmo/modatt/-def-/-def-/ioflag/x finish cmo/DELATT/mo/xmed cmo/DELATT/mo/ymed cmo/DELATT/mo/zmed cmo/DELATT/mo/ikey cmo/DELATT/mo/matclr DELATT no action, attribute does not exist: matclr ###### TOP FACESETS ########################### # define file names here so when dump/exo is called # it can use the defines to get the files correct later define FS_1 exo_1_top_left.faceset define FS_2 exo_2_top_right.faceset # extract the outside surface of the full mesh # to get outside face to element relationships extract/surfmesh/1 0 0/mo_surf/mo/external cmo/addatt/mo_surf/itetclr0/vint/scalar/nelements////0 finish cmo/addatt/mo_surf/itetclr1/vint/scalar/nelements////0 finish cmo/addatt/mo_surf/idnode0/vint/scalar/nnodes////0 finish cmo/addatt/mo_surf/idelem0/vint/scalar/nelements////0 finish cmo/addatt/mo_surf/idelem1/vint/scalar/nelements////0 finish cmo/addatt/mo_surf/facecol/vint/scalar/nelements////0 finish cmo/addatt/mo_surf/idface0/vint/scalar/nelements////0 finish cmo/addatt/mo_surf/idface1/vint/scalar/nelements////0 finish cmo/addatt/mo_surf/top/VINT/scalar/nnodes/linear/permanent/gxaf finish cmo/addatt/mo_surf/bottom/VINT/scalar/nnodes/linear/permanent/gxaf finish cmo/addatt/mo_surf/left_w/VINT/scalar/nnodes/linear/permanent/gxaf finish cmo/addatt/mo_surf/right_e/VINT/scalar/nnodes/linear/permanent/gxaf finish cmo/addatt/mo_surf/back_n/VINT/scalar/nnodes/linear/permanent/gxaf finish cmo/addatt/mo_surf/front_s/VINT/scalar/nnodes/linear/permanent/gxaf finish cmo/addatt/mo_surf/ncon50/int/scalar/scalar/constant/permanent/x/0 finish cmo/addatt/mo_surf/nconbnd/int/scalar/scalar/constant/permanent/x/0 finish cmo/addatt/mo_surf/icontab/vint/scalar/ncon50/constant/permanent/x/0 finish resetpts/itp geniee finish finish cmo/printatt/mo_surf/-all-/minmax ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH -def- 0.000000000E+00 0.000000000E+00 0.000000000E+00 8306 scalar 1 1 0 1 vector 3 3 0 1 nnodes 8306 8306 0 1 nedges 0 0 0 1 nfaces 0 0 0 1 nelements 16608 16608 0 1 mbndry 16000000 16000000 0 1 ndimensions_topo 2 2 0 1 ndimensions_geom 3 3 0 1 nodes_per_element 3 3 0 1 edges_per_element 3 3 0 1 faces_per_element 3 3 0 1 isetwd 0 0 0 8306 ialias 0 0 0 8306 imt1 1 1 0 8306 itp1 0 0 0 8306 icr1 0 0 0 8306 isn1 0 0 0 8306 xic 0.000000000E+00 1.000000000E+03 1.000000000E+03 8306 yic 0.000000000E+00 1.000000000E+03 1.000000000E+03 8306 zic 0.000000000E+00 6.000000000E+02 6.000000000E+02 8306 xtetwd 0 0 0 16608 itetclr 1 1 0 16608 itettyp 3 3 0 16608 itetoff 0 49821 49821 16608 jtetoff 0 49821 49821 16608 itet 1 8306 8305 16608x3 jtet 1 49824 49823 16608x3 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 8306 8306 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 0 0 0 1 number_of_eltsets 0 0 0 1 number_of_fsets 0 0 0 1 itetclr0 0 0 0 16608 itetclr1 1 2 1 16608 idnode0 1 160583 160582 8306 idelem0 0 0 0 16608 idelem1 1 896606 896605 16608 facecol 1 2 1 16608 idface0 0 0 0 16608 idface1 1 4 3 16608 top -5 1 6 8306 bottom -5 2 7 8306 left_w -5 3 8 8306 right_e -5 5 10 8306 back_n -5 6 11 8306 front_s -5 4 9 8306 ncon50 0 0 0 1 nconbnd 0 0 0 1 0 length attribute: icontab # ## Set surface elements to direction they face based on normals # These values will be written to itetclr, copy to id_side attribute # 1 = bottom 2 = top, 4-6 are sides cmo/addatt/mo_surf/id_side/vint/scalar/nelements cmo/select mo_surf settets/normal assign_color_normal epsilon for zero: 1.0000000E-15 12168 out of 16608 element colors changed. cmo/copyatt/mo_surf mo_surf/id_side itetclr 16608 copied from mo_surf itetclr to -> mo_surf id_side cmo/printatt/mo_surf/id_side/minmax ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH id_side 1 6 5 16608 dump gmv outside_surface.gmv mo_surf cmo/modatt/-def-/-def-/ioflag/x finish # make sure to remove all attributes except idelem1 and idface1 cmo/DELATT/mo_surf/itetclr0 cmo/DELATT/mo_surf/idnode0 cmo/DELATT/mo_surf/idelem0 cmo/DELATT/mo_surf/facecol cmo/DELATT/mo_surf/itetclr1 cmo/DELATT/mo_surf/idface0 ############### SET FACES - top left #################### # TOP LEFT = 1 cmo/copy/mo_tmp/mo_surf cmo/select/mo_tmp # changed to lt so face sets do not overlap # pset / psel / attribute / xic / 1,0,0 / le / 700.0 pset/psel/attribute/xic/1,0,0/lt/700.0 THE PSET psel HAS 3384 POINTS eltset esel inclusive pset,get,psel THE ELTSET esel HAS 7056 ELEMENTS eltset/e_top/itetclr eq 2 THE ELTSET e_top HAS 6320 ELEMENTS eltset/e_keep/inter e_top esel THE ELTSET e_keep HAS 2780 ELEMENTS eltset/e_delete/not/e_keep THE ELTSET e_delete HAS 13828 ELEMENTS rmpoint/element/eltset get e_delete Dudded 6848 points plus parents, and 13828 elements removed. rmpoint/compress 6848 points removed and 0 elements removed. The new point count is 1458. The new element count is 2780. cmo/DELATT/mo_tmp/id_side dump/avs2/FS_1/mo_tmp/0 0 0 2 WARNING: dump/avs iopt_points= 0 WARNING: dump/avs iopt_elements= 0 WARNING: dump/avs iopt_values_node= 0 WARNING: dump/avs iopt_values_elem= 2 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/mo_tmp Released Mesh Object: mo_tmp # TOP RIGHT = 2 cmo/copy/mo_tmp/mo_surf cmo/select/mo_tmp pset/psel/attribute/xic/1,0,0/gt/700.0 THE PSET psel HAS 4632 POINTS eltset esel inclusive pset,get,psel THE ELTSET esel HAS 9552 ELEMENTS eltset/e_top/itetclr eq 2 THE ELTSET e_top HAS 6320 ELEMENTS eltset/e_keep/inter e_top esel THE ELTSET e_keep HAS 3540 ELEMENTS eltset/e_delete/not/e_keep THE ELTSET e_delete HAS 13068 ELEMENTS rmpoint/element/eltset get e_delete Dudded 6489 points plus parents, and 13068 elements removed. rmpoint/compress 6489 points removed and 0 elements removed. The new point count is 1817. The new element count is 3540. cmo/DELATT/mo_tmp/id_side dump/avs2/FS_2/mo_tmp/0 0 0 2 WARNING: dump/avs iopt_points= 0 WARNING: dump/avs iopt_elements= 0 WARNING: dump/avs iopt_values_node= 0 WARNING: dump/avs iopt_values_elem= 2 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/mo_tmp Released Mesh Object: mo_tmp ############################################################ # FINAL FILES # write final output mesh and with faceset file # order of files detirmines integer value of the set # This exo file can be read into GMV, side sets are under Display Surfaces dump/exo/tet_exodus.exo/mo///facesets FS_1 FS_2 ExodusII: Start writing to file: tet_exodus.exo using cmo: mo cmo/set_id/mo/element/e_num_temp cmo/addatt/mo e_num_temp/VINT/scalar/nelements/linear/permanent//0 finish finish sort/mo/index/ascending/ikey_utr/itetclr/e_num_temp cmo/addatt/mo/ikey_utr/vint/scalar/nelements///gax/0 finish SORT: order key written to attribute: ikey_utr finish cmo/DELATT/mo/e_num_temp finish cmo/modatt/mo/itetclr/ioflag/lga finish INITIALIZE exodus Title: LAGRIT TO EXODUS number of dimension: 3 number of nodes: 160583 number of elements: 921085 number of edges: 0 number of edge blocks: 0 number of element blocks: 3 number of face blocks: 0 number of node sets: 0 number of edge sets: 3 number of element sets: 0 number of side sets: 2 number of face sets: 0 number of node maps: 0 number of edge maps: 0 number of face maps: 0 number of element maps: 0 cmo/modatt/mo/itetclr/ioflag/l finish ExodusII: Done writing to file: tet_exodus.exo using cmo: mo dump gmv tet_exodus.gmv mo cmo/modatt/-def-/-def-/ioflag/x finish # DONE finish LaGriT successfully completed