* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * Program: LaGriT V3.300 Linux * * * date_compile: 2017/10/25 Release * * * Run Time: 2018/Nov 6 10:04:57 * * * Manual: http://lagrit.lanl.gov * * * * * * * * * * * * * * * * * * * * * * * * * * * * * -----oOo----- LaGriT V3 LA-CC-15-069 https://github.com/lanl/LaGriT Copyright 2015. Los Alamos National Security, LLC. This software was produced under U.S. Government contract DE-AC52-06NA25396 for Los Alamos National Laboratory (LANL), which is operated by Los Alamos National Security, LLC for the U.S. Department of Energy. The U.S. Government has rights to use, reproduce, and distribute this software. Neither the government nor Los Alamos National Security, LLC makes any warranty, express or implied, or assumes any liability for the use of this software. If software is modified to produce derivative works, such modified software should be clearly marked, so as not to confuse it with the version available from LANL. -----oOo----- Output log file: outx3dgen Command log file: logx3dgen # write exodus file with side facesets # interpolated from boundary attributes itetclr and lay1_out # overwrite side values with the top layer outlet value # # There are 2 easy ways to find the top layer # 1. make attribute and assign 1 - nlayers (in the stack mesh) # 2. use layertyp = -2 to pset top and inclusive connected sides # then make element set of faces on sides but not top # # This example includes both methods to find and set elements by layer id # But uses method 2 which is more general # # See surfmesh_fs_all.inp and .gmv to check for correct fs assignments # define INPUT_MESH stack.inp define INPUT_BNDRY bndry_line_colors.inp define CMO_IN momesh define CMO_BNDRY moline # READ mesh read/INPUT_MESH/CMO_IN read/avs/stack.inp/momesh cmo/addatt/momesh/layertyp/VINT/scalar/nnodes/linear/permanent/gxaf/0.0 finish geniee finish cmo/status/brief The current-mesh-object(CMO) is: momesh 1 Mesh Object name: momesh number of nodes = 1000 number of elements = 1344 dimensions geometry = 3 element type = pri dimensions topology = 3 6 nodes 5 faces 9 edges boundary flag = 16000000 status = active finish finish resetpts/itp geniee finish ################### color mesh by nlayers ################# # METHOD 1. FIND AND SET ilayer 1-N starting at top # THIS IS NOT GENERAL # this might take awhile but needed if doing layers in addition to top layer # this assumes each layer is neperlayer=336 with nelements=1344 # this assumes element order from top to bottom cmo/addatt/CMO_IN/ilayer/VINT/scalar/nelements/linear/permanent/gxaf/0.0 cmo/setatt/CMO_IN ilayer 1 1344 values reset for attribute ilayer cmo/setatt/CMO_IN ilayer/673,1008 2 336 values reset for attribute ilayer cmo/setatt/CMO_IN ilayer/337,672 3 336 values reset for attribute ilayer cmo/setatt/CMO_IN ilayer/1,336,1 4 336 values reset for attribute ilayer # check ilayer colors cmo printatt CMO_IN ilayer minmax ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH ilayer 1 4 3 1344 dump gmv tmp_ilayers.gmv CMO_IN cmo/modatt/-def-/-def-/ioflag/x finish ############################################################ # ASSUME INPUT BNDRY LINE ATTRIBUTES itetclr lay1_out read/INPUT_BNDRY/CMO_BNDRY read/avs/bndry_line_colors.inp/moline cmo/addatt/moline/idnode0/VINT/scalar/nnodes/linear/permanent/gxaf/0.0 finish cmo/addatt/moline/elev/VDOUBLE/scalar/nnodes/linear/permanent/gxaf/0.0 finish cmo/addatt/moline/var1/VDOUBLE/scalar/nnodes/linear/permanent/gxaf/0.0 finish cmo/addatt/moline/var2/VDOUBLE/scalar/nnodes/linear/permanent/gxaf/0.0 finish cmo/addatt/moline/nkey/VINT/scalar/nnodes/linear/permanent/gxaf/0.0 finish cmo/addatt/moline/itetclr0/VINT/scalar/nelements/linear/permanent/gxaf/0.0 finish cmo/addatt/moline/itetclr1/VINT/scalar/nelements/linear/permanent/gxaf/0.0 finish cmo/addatt/moline/idelem0/VINT/scalar/nelements/linear/permanent/gxaf/0.0 finish cmo/addatt/moline/idelem1/VINT/scalar/nelements/linear/permanent/gxaf/0.0 finish cmo/addatt/moline/facecol/VINT/scalar/nelements/linear/permanent/gxaf/0.0 finish cmo/addatt/moline/idface0/VINT/scalar/nelements/linear/permanent/gxaf/0.0 finish cmo/addatt/moline/idface1/VINT/scalar/nelements/linear/permanent/gxaf/0.0 finish cmo/addatt/moline/ikey/VINT/scalar/nelements/linear/permanent/gxaf/0.0 finish cmo/addatt/moline/cid/VINT/scalar/nelements/linear/permanent/gxaf/0.0 finish cmo/addatt/moline/ctype/VINT/scalar/nelements/linear/permanent/gxaf/0.0 finish cmo/addatt/moline/loopid/VINT/scalar/nelements/linear/permanent/gxaf/0.0 finish cmo/addatt/moline/lay1_out/VINT/scalar/nelements/linear/permanent/gxaf/0.0 finish geniee finish cmo/status/brief The current-mesh-object(CMO) is: moline 1 Mesh Object name: momesh number of nodes = 1000 number of elements = 1344 dimensions geometry = 3 element type = pri dimensions topology = 3 6 nodes 5 faces 9 edges boundary flag = 16000000 status = inactive 2 Mesh Object name: moline number of nodes = 62 number of elements = 62 dimensions geometry = 3 element type = lin dimensions topology = 1 2 nodes 2 faces 1 edges boundary flag = 16000000 status = active finish finish resetpts/itp geniee finish # CHECK overlap of mesh and bndry cmo/printatt/CMO_IN/-xyz- minmax ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH xic 0.000000000E+00 1.200000000E+01 1.200000000E+01 1000 yic 0.000000000E+00 3.500000000E+00 3.500000000E+00 1000 zic 0.000000000E+00 7.398997995E+00 7.398997995E+00 1000 cmo/printatt/CMO_BNDRY/-xyz- minmax ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH xic 0.000000000E+00 1.200000000E+01 1.200000000E+01 62 yic 0.000000000E+00 3.500000000E+00 3.500000000E+00 62 zic 0.000000000E+00 0.000000000E+00 0.000000000E+00 62 # CHECK ATTRIBUTES # these will be used to set facesets on sides and top layer cmo/printatt/CMO_BNDRY/itetclr minmax ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH itetclr 1 5 4 62 cmo/printatt/CMO_BNDRY/lay1_out minmax ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH lay1_out 0 1 1 62 # PREPARE stack mesh for Exodus # sort based on element itetclr number and median location # Note Exodus expects materials to appear in mesh in order 1 to N # this is also done so numbering for Exodus matches other formats cmo/select/CMO_IN 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/CMO_IN/index/ascending/ikey/itetclr xmed ymed zmed cmo/addatt/momesh/ikey/vint/scalar/nelements///gax/0 finish SORT: order key written to attribute: ikey reorder/CMO_IN/ikey geniee finish cmo/DELATT/CMO_IN/xmed cmo/DELATT/CMO_IN/ymed cmo/DELATT/CMO_IN/zmed cmo/DELATT/CMO_IN/ikey # EXTRACT surface with cell and face attributes # to get the outside face to element relationships # this will be used to write faceset files for dump/exo extract/surfmesh/1 0 0/mo_surf/CMO_IN/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/layertyp/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 586 scalar 1 1 0 1 vector 3 3 0 1 nnodes 586 586 0 1 nedges 0 0 0 1 nfaces 0 0 0 1 nelements 920 920 0 1 mbndry 16000000 16000000 0 1 ndimensions_topo 2 2 0 1 ndimensions_geom 3 3 0 1 nodes_per_element 10 10 0 1 edges_per_element 12 12 0 1 faces_per_element 10 10 0 1 isetwd 0 0 0 586 ialias 0 0 0 586 imt1 1 1 0 586 itp1 0 2 2 586 icr1 0 0 0 586 isn1 0 0 0 586 xic 0.000000000E+00 1.200000000E+01 1.200000000E+01 586 yic 0.000000000E+00 3.500000000E+00 3.500000000E+00 586 zic 0.000000000E+00 7.398997995E+00 7.398997995E+00 586 xtetwd 0 0 0 920 itetclr 1 1 0 920 itettyp 3 4 1 920 itetoff 0 3004 3004 920 jtetoff 0 3004 3004 920 itet 0 586 586 920x* jtet 0 9194 9194 920x* 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 586 586 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 920 itetclr1 1 3 2 920 idnode0 1 1000 999 586 idelem0 0 0 0 920 idelem1 1 1344 1343 920 facecol 1 3 2 920 idface0 0 0 0 920 idface1 1 5 4 920 layertyp -2 2 4 586 ncon50 0 0 0 1 nconbnd 0 0 0 1 0 length attribute: icontab # This is the default for finding mesh facesets 1-6 # Set surface elements to direction they face based on normals # This can be skipped or just overwrite cmo/addatt/mo_surf/id_side/vint/scalar/nelements cmo/select mo_surf settets/normal assign_color_normal epsilon for zero: 1.0000000E-15 584 out of 920 element colors changed. cmo/copyatt/mo_surf mo_surf/id_side itetclr 920 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 920 dump gmv tmp_fs_normals.gmv mo_surf cmo/modatt/-def-/-def-/ioflag/x finish # This selection is more accurate for top and bottom # use stack attribute to avoid surface with multiple normal directions # set all sides to default 3 all cmo/select mo_surf cmo setatt mo_surf/id_side/3 920 values reset for attribute id_side pset/ptop/attribute layertyp/1,0,0/eq -2 THE PSET ptop HAS 200 POINTS pset/pbot/attribute layertyp/1,0,0/eq -1 THE PSET pbot HAS 200 POINTS eltset/etop/exclusive pset,get,ptop THE ELTSET etop HAS 336 ELEMENTS eltset/ebot/exclusive pset,get,pbot THE ELTSET ebot HAS 336 ELEMENTS cmo setatt mo_surf/id_side/eltset,get,etop 2 336 values reset for attribute id_side cmo setatt mo_surf/id_side/eltset,get,ebot 1 336 values reset for attribute id_side cmo printatt mo_surf/id_side/minmax ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH id_side 1 3 2 920 cmo/copyatt/mo_surf mo_surf/itetclr id_side 920 copied from mo_surf id_side to -> mo_surf itetclr dump gmv tmp_surf_sides.gmv mo_surf cmo/modatt/-def-/-def-/ioflag/x finish ############## NEW BNDRY SET SIDES ############################ # based on boundary lines, set all sides same as boundary lines # ELEMENT SET of all sides cmo select mo_surf eltset/esides/id_side eq 3 THE ELTSET esides HAS 248 ELEMENTS pset/psides/eltset/esides THE PSET psides HAS 310 POINTS # POINT SET of sides, top, and bottom cmo/setatt/mo_surf/imt pset,get,ptop 2 200 values reset for attribute imt1 cmo/setatt/mo_surf/imt pset,get,pbot 1 200 values reset for attribute imt1 cmo/setatt/mo_surf/imt pset,get,psides 3 310 values reset for attribute imt1 dump gmv tmp_surf_sides.gmv mo_surf cmo/modatt/-def-/-def-/ioflag/x finish # PREPARE FOR INTERPOLATE # Make 3D mesh flat on same plane as boundary lines cmo select mo_surf cmo/addatt/mo_surf/zsave/vdouble/scalar/nnodes cmo copyatt mo_surf mo_surf zsave zic 586 copied from mo_surf zic to -> mo_surf zsave cmo setatt mo_surf zic 0. 586 values reset for attribute zic cmo setatt CMO_BNDRY zic 0. 62 values reset for attribute zic # INTERPOLATE boundary faces to side faces # FIX numbering so mesh itetclr 1 and 2 remain top and bottom # ADD 2 to the boundary line materials assuming they start with 1 math/add/CMO_BNDRY/itetclr/1,0,0/CMO_BNDRY itetclr/2 cmo printatt CMO_BNDRY itetclr minmax ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH itetclr 3 7 4 62 # SET ATTRIBUTE id_side and lay1_out from boundary lines cmo/addatt/mo_surf/lay1_out/VINT/scalar/nelements/linear/permanent/gxaf/0.0 cmo/addatt/mo_surf/ilayer/VINT/scalar/nelements/linear/permanent/gxaf/0.0 interpolate/map/mo_surf id_side/eltset,get,esides/CMO_BNDRY itetclr/keepatt elements in indexed set = 248 cmo/addatt/mo_surf/el_gtg/VINT/scalar/nelements//permanent/agfx/0 finish INTRP METHOD: map FUNCTION: user TIEBREAKER: maxtie FLAG: plus1 0.800000000000E+01 el_gtg attribute will be added and kept. cmo select moline finish cmo kdtree build cmo/addatt//v2/INT/scalar/scalar/constant/temporary//2.0 finish cmo/addatt//linkt/VINT/v2/nelements//temporary/x/0.0 finish cmo/addatt//v12/INT/scalar/scalar/constant/temporary//12.0 finish cmo/addatt//sbox/VDOUBLE/v12/nelements/linear/temporary/x/0.0 finish finish Sink point Elems Searched Elements Found Percent Done 63 63 63 25 % 125 125 125 50 % 187 187 187 75 % 248 248 248 Total cmo select mo_surf finish Indexed sink points: 248 Total unchanged: 672 intrp/map done. cmo kdtree release cmo/DELATT/moline/sbox finish cmo/DELATT/moline/v12 finish cmo/DELATT/moline/linkt finish cmo/DELATT/moline/v2 finish finish interpolate/map/mo_surf lay1_out/eltset,get,esides/CMO_BNDRY lay1_out elements in indexed set = 248 INTRP METHOD: map FUNCTION: linear TIEBREAKER: maxtie FLAG: plus1 0.200000000000E+01 SKIPPING ELEMENT SEARCH... using lookup attribute el_gtg Sink point Elems Searched Elements Found Percent Done 63 63 63 25 % 125 125 125 50 % 187 187 187 75 % 248 248 248 Total cmo DELATT mo_surf/el_gtg finish cmo select mo_surf finish Indexed sink points: 248 Total unchanged: 672 intrp/map done. cmo kdtree release finish cmo printatt CMO_BNDRY id_sides minmax cannot attribute id_sides in moline cmo printatt mo_surf id_sides minmax cannot attribute id_sides in mo_surf cmo printatt mo_surf lay1_out minmax ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH lay1_out 0 1 1 920 cmo copyatt mo_surf mo_surf zic zsave 586 copied from mo_surf zsave to -> mo_surf zic cmo DELATT mo_surf zsave dump gmv tmp_interpolate_sides.gmv mo_surf cmo/modatt/-def-/-def-/ioflag/x finish # CHECK facesets have 2 plus bndry num side values # total = bndry + top and bottom = 5 + 2 =7 cmo printatt mo_surf itetclr minmax ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH itetclr 1 3 2 920 # SET FINAL FACE VALUES for SIDES (outlet will overwrite these) # TOP and Bottom were defined earlier, do the sides eltset/e3/id_side/eq 3 THE ELTSET e3 HAS 72 ELEMENTS eltset/e4/id_side/eq 4 THE ELTSET e4 HAS 24 ELEMENTS eltset/e5/id_side/eq 5 THE ELTSET e5 HAS 96 ELEMENTS eltset/e6/id_side/eq 6 THE ELTSET e6 HAS 28 ELEMENTS eltset/e7/id_side/eq 7 THE ELTSET e7 HAS 28 ELEMENTS cmo setatt mo_surf/id_side/eltset,get,etop 2 336 values reset for attribute id_side cmo setatt mo_surf/id_side/eltset,get,ebot 1 336 values reset for attribute id_side cmo printatt mo_surf/id_side/minmax ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH id_side 1 7 6 920 # check material numbers, must be greater than 0 # copy id_side to itetclr which is done except for outlets cmo/copyatt/mo_surf mo_surf/itetclr id_side 920 copied from mo_surf id_side to -> mo_surf itetclr cmo printatt CMO_BNDRY imt minmax ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH imt1 1 5 4 62 cmo printatt mo_surf imt minmax ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH imt1 1 3 2 586 cmo printatt mo_surf itetclr minmax ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH itetclr 1 7 6 920 dump gmv surfmesh_fs_sides.gmv mo_surf cmo/modatt/-def-/-def-/ioflag/x finish dump avs surfmesh_fs_sides.inp mo_surf cmo/modatt/-def-/-def-/ioflag/x finish ############## NEW SET OUTLET LAYER 1 ########################### # based on boundary lines, set all sides same as boundary lines # There are 2 easy ways to find the top layer # - make attribute and assign 1 - nlayers (in the stack mesh) # - use layertyp = -2 to pset top and inclusive connected sides # then make element set of faces on sides but not top # # FOR THIS WORKFLOW USE top surface points to find layer 1 top layer cmo select mo_surf cmo setatt mo_surf ilayer 0 920 values reset for attribute ilayer # METHOD 2. FIND AND SET top ilayer 1 # SET ATTRIBUTE ilayer pset/ptop/attribute layertyp/1,0,0/eq -2 THE PSET ptop HAS 200 POINTS eltset/elay_inc/inclusive pset,get,ptop THE ELTSET elay_inc HAS 398 ELEMENTS cmo setatt mo_surf ilayer eltset,get,elay_inc 1 398 values reset for attribute ilayer cmo setatt mo_surf ilayer eltset,get,etop 0 336 values reset for attribute ilayer # INTERSECT ilayer with lay1_out and overwrite itetclr # This faceset will be value of itetclr + 1 = 8 eltset/e1/ilayer eq 1 THE ELTSET e1 HAS 62 ELEMENTS eltset/e2/lay1_out eq 1 THE ELTSET e2 HAS 24 ELEMENTS eltset/e_out1/inter e1 e2 THE ELTSET e_out1 HAS 6 ELEMENTS cmo/setatt/mo_surf/itetclr/eltset,get,e_out1 8 6 values reset for attribute itetclr dump gmv surfmesh_fs_all.gmv mo_surf cmo/modatt/-def-/-def-/ioflag/x finish dump avs surfmesh_fs_all.inp mo_surf cmo/modatt/-def-/-def-/ioflag/x finish cmo/printatt/mo_surf/-all- minmax ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH -def- 0.000000000E+00 0.000000000E+00 0.000000000E+00 586 scalar 1 1 0 1 vector 3 3 0 1 nnodes 586 586 0 1 nedges 0 0 0 1 nfaces 0 0 0 1 nelements 920 920 0 1 mbndry 16000000 16000000 0 1 ndimensions_topo 2 2 0 1 ndimensions_geom 3 3 0 1 nodes_per_element 10 10 0 1 edges_per_element 12 12 0 1 faces_per_element 10 10 0 1 isetwd 1 6 5 586 ialias 0 0 0 586 imt1 1 3 2 586 itp1 0 2 2 586 icr1 0 0 0 586 isn1 0 0 0 586 xic 0.000000000E+00 1.200000000E+01 1.200000000E+01 586 yic 0.000000000E+00 3.500000000E+00 3.500000000E+00 586 zic 0.000000000E+00 7.398997995E+00 7.398997995E+00 586 xtetwd 2 3860 3858 920 itetclr 1 8 7 920 itettyp 3 4 1 920 itetoff 0 3004 3004 920 jtetoff 0 3004 3004 920 itet 0 586 586 920x* jtet 0 9194 9194 920x* epsilon 1.000000004E-15 1.000000004E-15 0.000000000E+00 1 epsilonl 2.775557562E-12 2.775557562E-12 0.000000000E+00 1 epsilona 3.469446952E-11 3.469446952E-11 0.000000000E+00 1 epsilonv 3.469446952E-11 3.469446952E-11 0.000000000E+00 1 ipointi 1 1 0 1 ipointj 586 586 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.200000000E+01 1.200000000E+01 0.000000000E+00 1 ymax 3.500000000E+00 3.500000000E+00 0.000000000E+00 1 zmax 0.000000000E+00 0.000000000E+00 0.000000000E+00 1 kdtree_level 0 0 0 1 max_number_sets 64 64 0 1 number_of_psets 3 3 0 1 number_of_eltsets 0 0 0 1 number_of_fsets 0 0 0 1 itetclr0 0 0 0 920 itetclr1 1 3 2 920 idnode0 1 1000 999 586 idelem0 0 0 0 920 idelem1 1 1344 1343 920 facecol 1 3 2 920 idface0 0 0 0 920 idface1 1 5 4 920 layertyp -2 2 4 586 ncon50 0 0 0 1 nconbnd 0 0 0 1 0 length attribute: icontab id_side 1 7 6 920 lay1_out 0 1 1 920 ilayer 0 1 1 920 # 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 cmo/DELATT/mo_surf/ikey_utr DELATT no action, attribute does not exist: ikey_utr cmo/DELATT/mo_surf/nlayers DELATT no action, attribute does not exist: nlayers cmo/DELATT/mo_surf/nnperlayer DELATT no action, attribute does not exist: nnperlayer cmo/DELATT/mo_surf/neperlayer DELATT no action, attribute does not exist: neperlayer cmo/DELATT/mo_surf/layertyp cmo/DELATT/mo_surf/ncon50 cmo/DELATT/mo_surf/nconbnd cmo/DELATT/mo_surf/icontab DELATT no action, attribute does not exist: icontab cmo/DELATT/mo_surf/id_side cmo/DELATT/mo_surf/ilayer cmo/DELATT/mo_surf/lay1_out # CHECK only user attributes should be idelem1 and idface1 cmo/printatt/mo_surf/-all- minmax ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH -def- 0.000000000E+00 0.000000000E+00 0.000000000E+00 586 scalar 1 1 0 1 vector 3 3 0 1 nnodes 586 586 0 1 nedges 0 0 0 1 nfaces 0 0 0 1 nelements 920 920 0 1 mbndry 16000000 16000000 0 1 ndimensions_topo 2 2 0 1 ndimensions_geom 3 3 0 1 nodes_per_element 10 10 0 1 edges_per_element 12 12 0 1 faces_per_element 10 10 0 1 isetwd 1 6 5 586 ialias 0 0 0 586 imt1 1 3 2 586 itp1 0 2 2 586 icr1 0 0 0 586 isn1 0 0 0 586 xic 0.000000000E+00 1.200000000E+01 1.200000000E+01 586 yic 0.000000000E+00 3.500000000E+00 3.500000000E+00 586 zic 0.000000000E+00 7.398997995E+00 7.398997995E+00 586 xtetwd 2 3860 3858 920 itetclr 1 8 7 920 itettyp 3 4 1 920 itetoff 0 3004 3004 920 jtetoff 0 3004 3004 920 itet 0 586 586 920x* jtet 0 9194 9194 920x* epsilon 1.000000004E-15 1.000000004E-15 0.000000000E+00 1 epsilonl 2.775557562E-12 2.775557562E-12 0.000000000E+00 1 epsilona 3.469446952E-11 3.469446952E-11 0.000000000E+00 1 epsilonv 3.469446952E-11 3.469446952E-11 0.000000000E+00 1 ipointi 1 1 0 1 ipointj 586 586 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.200000000E+01 1.200000000E+01 0.000000000E+00 1 ymax 3.500000000E+00 3.500000000E+00 0.000000000E+00 1 zmax 0.000000000E+00 0.000000000E+00 0.000000000E+00 1 kdtree_level 0 0 0 1 max_number_sets 64 64 0 1 number_of_psets 3 3 0 1 number_of_eltsets 0 0 0 1 number_of_fsets 0 0 0 1 idelem1 1 1344 1343 920 idface1 1 5 4 920 ######## WRITE FACESET FILES #################### # itetclr now has all faces colored from 1 to N=8 # LOOP through 1 to N and write fsID.faceset files # we assume 1 is bottom and 2 is top define/FILENAME/fs1_bot.faceset define/SS_ID/1 cmo/copy/mo_tmp/mo_surf cmo/select/mo_tmp eltset/e_keep/itetclr/eq/SS_ID THE ELTSET e_keep HAS 336 ELEMENTS eltset/e_delete/not/e_keep THE ELTSET e_delete HAS 584 ELEMENTS rmpoint/element/eltset get e_delete Dudded 386 points plus parents, and 584 elements removed. rmpoint/compress 386 points removed and 0 elements removed. RMPOINT: new point count is 200 RMPOINT: new element count is 336 dump/avs2/FILENAME/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 define/FILENAME/fs2_top.faceset define/SS_ID/2 cmo/copy/mo_tmp/mo_surf cmo/select/mo_tmp eltset/e_keep/itetclr/eq/SS_ID THE ELTSET e_keep HAS 336 ELEMENTS eltset/e_delete/not/e_keep THE ELTSET e_delete HAS 584 ELEMENTS rmpoint/element/eltset get e_delete Dudded 386 points plus parents, and 584 elements removed. rmpoint/compress 386 points removed and 0 elements removed. RMPOINT: new point count is 200 RMPOINT: new element count is 336 dump/avs2/FILENAME/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 ## CHECK that top number of faces = bottom = layer define/FILENAME/fs3.faceset define/SS_ID/3 cmo/copy/mo_tmp/mo_surf cmo/select/mo_tmp eltset/e_keep/itetclr/eq/SS_ID THE ELTSET e_keep HAS 72 ELEMENTS eltset/e_delete/not/e_keep THE ELTSET e_delete HAS 848 ELEMENTS rmpoint/element/eltset get e_delete Dudded 486 points plus parents, and 848 elements removed. rmpoint/compress 486 points removed and 0 elements removed. RMPOINT: new point count is 100 RMPOINT: new element count is 72 dump/avs2/FILENAME/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 define/FILENAME/fs3.faceset define/SS_ID/3 cmo/copy/mo_tmp/mo_surf cmo/select/mo_tmp eltset/e_keep/itetclr/eq/SS_ID THE ELTSET e_keep HAS 72 ELEMENTS eltset/e_delete/not/e_keep THE ELTSET e_delete HAS 848 ELEMENTS rmpoint/element/eltset get e_delete Dudded 486 points plus parents, and 848 elements removed. rmpoint/compress 486 points removed and 0 elements removed. RMPOINT: new point count is 100 RMPOINT: new element count is 72 dump/avs2/FILENAME/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 define/FILENAME/fs4.faceset define/SS_ID/4 cmo/copy/mo_tmp/mo_surf cmo/select/mo_tmp eltset/e_keep/itetclr/eq/SS_ID THE ELTSET e_keep HAS 18 ELEMENTS eltset/e_delete/not/e_keep THE ELTSET e_delete HAS 902 ELEMENTS rmpoint/element/eltset get e_delete Dudded 558 points plus parents, and 902 elements removed. rmpoint/compress 558 points removed and 0 elements removed. RMPOINT: new point count is 28 RMPOINT: new element count is 18 dump/avs2/FILENAME/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 define/FILENAME/fs5.faceset define/SS_ID/5 cmo/copy/mo_tmp/mo_surf cmo/select/mo_tmp eltset/e_keep/itetclr/eq/SS_ID THE ELTSET e_keep HAS 96 ELEMENTS eltset/e_delete/not/e_keep THE ELTSET e_delete HAS 824 ELEMENTS rmpoint/element/eltset get e_delete Dudded 461 points plus parents, and 824 elements removed. rmpoint/compress 461 points removed and 0 elements removed. RMPOINT: new point count is 125 RMPOINT: new element count is 96 dump/avs2/FILENAME/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 define/FILENAME/fs6.faceset define/SS_ID/6 cmo/copy/mo_tmp/mo_surf cmo/select/mo_tmp eltset/e_keep/itetclr/eq/SS_ID THE ELTSET e_keep HAS 28 ELEMENTS eltset/e_delete/not/e_keep THE ELTSET e_delete HAS 892 ELEMENTS rmpoint/element/eltset get e_delete Dudded 546 points plus parents, and 892 elements removed. rmpoint/compress 546 points removed and 0 elements removed. RMPOINT: new point count is 40 RMPOINT: new element count is 28 dump/avs2/FILENAME/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 define/FILENAME/fs7.faceset define/SS_ID/7 cmo/copy/mo_tmp/mo_surf cmo/select/mo_tmp eltset/e_keep/itetclr/eq/SS_ID THE ELTSET e_keep HAS 28 ELEMENTS eltset/e_delete/not/e_keep THE ELTSET e_delete HAS 892 ELEMENTS rmpoint/element/eltset get e_delete Dudded 546 points plus parents, and 892 elements removed. rmpoint/compress 546 points removed and 0 elements removed. RMPOINT: new point count is 40 RMPOINT: new element count is 28 dump/avs2/FILENAME/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 define/FILENAME/fs8.faceset define/SS_ID/8 cmo/copy/mo_tmp/mo_surf cmo/select/mo_tmp eltset/e_keep/itetclr/eq/SS_ID THE ELTSET e_keep HAS 6 ELEMENTS eltset/e_delete/not/e_keep THE ELTSET e_delete HAS 914 ELEMENTS rmpoint/element/eltset get e_delete Dudded 572 points plus parents, and 914 elements removed. rmpoint/compress 572 points removed and 0 elements removed. RMPOINT: new point count is 14 RMPOINT: new element count is 6 dump/avs2/FILENAME/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 ############################ # WRITE final mesh with multi material and all facesets define OUT_EXO mesh_fs8.exo dump/exo/OUT_EXO/CMO_IN///facesets fs1_bot.faceset fs2_top.faceset fs3.faceset f s4.faceset fs5.faceset fs6.faceset fs7.faceset fs8.faceset got 1 fs1_bot.faceset got 2 fs2_top.faceset got 3 fs3.faceset got 4 fs4.faceset got 5 fs5.faceset got 6 fs6.faceset got 7 fs7.faceset got 8 fs8.faceset ExodusII: Start writing to file: mesh_fs8.exo using cmo: momesh cmo/set_id/momesh/element/e_num_temp cmo/addatt/momesh e_num_temp/VINT/scalar/nelements/linear/permanent//0 finish finish sort/momesh/index/ascending/ikey_utr/itetclr/e_num_temp cmo/addatt/momesh/ikey_utr/vint/scalar/nelements///gax/0 finish SORT: order key written to attribute: ikey_utr finish cmo/DELATT/momesh/e_num_temp finish Title: LAGRIT TO EXODUSII number of dimension: 3 number of nodes: 1000 number of elements: 1344 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: 0 number of element sets: 0 number of side sets: 8 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 ------------------------------------------ EXPSS loop: 1 Side Set tag: 1 Faces: 336 2 Side Set tag: 2 Faces: 336 3 Side Set tag: 3 Faces: 72 4 Side Set tag: 4 Faces: 18 5 Side Set tag: 5 Faces: 96 6 Side Set tag: 6 Faces: 28 7 Side Set tag: 7 Faces: 28 8 Side Set tag: 8 Faces: 6 ------------------------------------------ Done ExodusII Side Sets Total: 8 ExodusII: Done writing to ExodusII file: mesh_fs8.exo using cmo: momesh # GMV can display facesets as surfaces dump/gmv/mesh_exo.gmv/CMO_IN cmo/modatt/-def-/-def-/ioflag/x finish # REPORT ATTRIBUTE VALUES of MESH cmo/printatt/CMO_IN/-all- minmax ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH -def- 0.000000000E+00 0.000000000E+00 0.000000000E+00 1000 scalar 1 1 0 1 vector 3 3 0 1 nnodes 1000 1000 0 1 nedges 0 0 0 1 nfaces 0 0 0 1 nelements 1344 1344 0 1 mbndry 16000000 16000000 0 1 ndimensions_topo 3 3 0 1 ndimensions_geom 3 3 0 1 nodes_per_element 6 6 0 1 edges_per_element 9 9 0 1 faces_per_element 5 5 0 1 isetwd 0 0 0 1000 ialias 0 0 0 1000 imt1 1 3 2 1000 itp1 0 12 12 1000 icr1 0 0 0 1000 isn1 0 0 0 1000 xic 0.000000000E+00 1.200000000E+01 1.200000000E+01 1000 yic 0.000000000E+00 3.500000000E+00 3.500000000E+00 1000 zic 0.000000000E+00 7.398997995E+00 7.398997995E+00 1000 xtetwd 0 0 0 1344 itetclr 1 3 2 1344 itettyp 7 7 0 1344 itetoff 0 8058 8058 1344 jtetoff 0 6715 6715 1344 itet 1 1000 999 1344x6 jtet 4 16006711 16006707 1344x5 epsilon 1.000000004E-15 1.000000004E-15 0.000000000E+00 1 epsilonl 3.225347286E-12 3.225347286E-12 0.000000000E+00 1 epsilona 4.685033946E-11 4.685033946E-11 0.000000000E+00 1 epsilonv 6.900211864E-11 6.900211864E-11 0.000000000E+00 1 ipointi 1 1 0 1 ipointj 1000 1000 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.200000000E+01 1.200000000E+01 0.000000000E+00 1 ymax 3.500000000E+00 3.500000000E+00 0.000000000E+00 1 zmax 7.398997995E+00 7.398997995E+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 layertyp -2 2 4 1000 ilayer 1 4 3 1344 ikey_utr 1 1344 1343 1344 finish LaGriT successfully completed