* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * Program: LaGriT V3.101 Linux m64 * * * date_compile: 2013/04/25 RH gf4.5 * * * Run Time: 2013/Jun 19 15:45:44 * * * 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 = 634148 number of elements = 3686223 dimensions geometry = 3 element type = tet dimensions topology = 3 4 nodes 4 faces 6 edges boundary flag = 22000000 status = active finish cmo printatt mo -all- minmax ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH -def- 0.000000000E+00 0.000000000E+00 0.000000000E+00 634148 scalar 1 1 0 1 vector 3 3 0 1 nnodes 634148 634148 0 1 nedges 0 0 0 1 nfaces 0 0 0 1 nelements 3686223 3686223 0 1 mbndry 22000000 22000000 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 634148 ialias 0 0 0 634148 imt1 1 5 4 634148 itp1 0 12 12 634148 icr1 0 0 0 634148 isn1 0 0 0 634148 xic 0.000000000E+00 1.000000000E+03 1.000000000E+03 634148 yic 0.000000000E+00 1.000000000E+03 1.000000000E+03 634148 zic 0.000000000E+00 6.000000000E+02 6.000000000E+02 634148 xtetwd 0 0 0 3686223 itetclr 1 5 4 3686223 itettyp 5 5 0 3686223 itetoff 0 14744888 14744888 3686223 jtetoff 0 14744888 14744888 3686223 itet 1 634148 634147 3686223x4 jtet 2 36744749 36744747 3686223x4 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 634148 634148 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 634148 bottom -5 2 7 634148 left_w -5 3 8 634148 right_e -5 5 10 634148 back_n -5 6 11 634148 front_s -5 4 9 634148 ### 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 3686223 values reset for attribute itetclr cmo printatt mo imt minmax ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH imt1 1 5 4 634148 # 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 29575 POINTS eltset/e3/exclusive pset,get,p3 THE ELTSET e3 HAS 146894 ELEMENTS cmo setatt mo itetclr eltset,get,e3 3 146894 values reset for attribute itetclr # fault find imt 4, set to 2 # watch that it does not get too wide pset/p4/attribute imt/1,0,0/eq 4 THE PSET p4 HAS 478130 POINTS eltset/e4/exclusive pset,get,p4 THE ELTSET e4 HAS 2458156 ELEMENTS cmo setatt mo itetclr eltset,get,e4 2 2458156 values reset for attribute itetclr # write element sets for each material eltset/e_mat3/itetclr/eq/3 THE ELTSET e_mat3 HAS 146894 ELEMENTS eltset/e_mat2/itetclr/eq/2 THE ELTSET e_mat2 HAS 2458156 ELEMENTS eltset/e_mat1/itetclr/eq/1 THE ELTSET e_mat1 HAS 1081173 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 18141 scalar 1 1 0 1 vector 3 3 0 1 nnodes 18141 18141 0 1 nedges 0 0 0 1 nfaces 0 0 0 1 nelements 36278 36278 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 18141 ialias 0 0 0 18141 imt1 1 1 0 18141 itp1 0 0 0 18141 icr1 0 0 0 18141 isn1 0 0 0 18141 xic 0.000000000E+00 1.000000000E+03 1.000000000E+03 18141 yic 0.000000000E+00 1.000000000E+03 1.000000000E+03 18141 zic 0.000000000E+00 6.000000000E+02 6.000000000E+02 18141 xtetwd 0 0 0 36278 itetclr 1 1 0 36278 itettyp 3 3 0 36278 itetoff 0 108831 108831 36278 jtetoff 0 108831 108831 36278 itet 1 18141 18140 36278x3 jtet 1 108834 108833 36278x3 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 18141 18141 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 36278 itetclr1 1 2 1 36278 idnode0 1 634148 634147 18141 idelem0 0 0 0 36278 idelem1 1 3539329 3539328 36278 facecol 1 2 1 36278 idface0 0 0 0 36278 idface1 1 4 3 36278 top -5 1 6 18141 bottom -5 2 7 18141 left_w -5 3 8 18141 right_e -5 5 10 18141 back_n -5 6 11 18141 front_s -5 4 9 18141 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 26758 out of 36278 element colors changed. cmo/copyatt/mo_surf mo_surf/id_side itetclr 36278 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 36278 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 7033 POINTS eltset esel inclusive pset,get,psel THE ELTSET esel HAS 14548 ELEMENTS eltset/e_top/itetclr eq 2 THE ELTSET e_top HAS 11600 ELEMENTS eltset/e_keep/inter e_top esel THE ELTSET e_keep HAS 3360 ELEMENTS eltset/e_delete/not/e_keep THE ELTSET e_delete HAS 32918 ELEMENTS rmpoint/element/eltset get e_delete Dudded 16378 points plus parents, and 32918 elements removed. rmpoint/compress 16378 points removed and 0 elements removed. The new point count is 1763. The new element count is 3360. 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 10624 POINTS eltset esel inclusive pset,get,psel THE ELTSET esel HAS 21730 ELEMENTS eltset/e_top/itetclr eq 2 THE ELTSET e_top HAS 11600 ELEMENTS eltset/e_keep/inter e_top esel THE ELTSET e_keep HAS 8240 ELEMENTS eltset/e_delete/not/e_keep THE ELTSET e_delete HAS 28038 ELEMENTS rmpoint/element/eltset get e_delete Dudded 13956 points plus parents, and 28038 elements removed. rmpoint/compress 13956 points removed and 0 elements removed. The new point count is 4185. The new element count is 8240. 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: 634148 number of elements: 3686223 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