* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * Program: LaGriT V3.300 Linux * * * date_compile: 2017/10/25 Release * * * Run Time: 2018/Nov 6 09:00:49 * * * 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 facesets # facesets 3 bottom=1, top=2, sides=3 define INPUT_MESH stack.inp define CMO_IN momesh # 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 # 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 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 # default # Set surface elements to direction they face based on normals 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 # 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/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/ikey_utr DELATT no action, attribute does not exist: ikey_utr ######## DEFINE FACESET FILES #################### # this version has 3 sets 1=bottom, 2=top, 3=sides # use stack attribute layertyp to set top and bottom # 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 # attribute id_side can add flexibility to colors # when done setting id_side, copy to itetclr 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 ######## WRITE FACESET FILES #################### define/FILENAME/fs1_bot.faceset define/SS_ID/1 cmo/copy/mo_tmp/mo_surf cmo/select/mo_tmp eltset/e_keep/id_side/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 cmo/DELATT/mo_surf/layertyp cmo/DELATT/mo_tmp/id_side 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/id_side/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 cmo/DELATT/mo_surf/layertyp DELATT no action, attribute does not exist: layertyp cmo/DELATT/mo_tmp/id_side 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_sides.faceset define/SS_ID/2 cmo/copy/mo_tmp/mo_surf cmo/select/mo_tmp eltset/e_keep/id_side/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 cmo/DELATT/mo_surf/layertyp DELATT no action, attribute does not exist: layertyp cmo/DELATT/mo_tmp/id_side 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 some informational files and output surfmesh # write the extracted surfmesh with fs colors dump gmv surfmesh_fs3.gmv mo_surf cmo/modatt/-def-/-def-/ioflag/x finish dump avs surfmesh_fs3.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 0 2 2 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 2 2 920 itetclr 1 3 2 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 2 2 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 ncon50 0 0 0 1 nconbnd 0 0 0 1 0 length attribute: icontab id_side 1 3 2 920 cmo delete mo_surf Released Mesh Object: mo_surf 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 ############################ # WRITE final mesh with multi material and all facesets define OUT_EXO mesh_fs3.exo dump/exo/OUT_EXO/CMO_IN///facesets fs1_bot.faceset fs2_top.faceset fs3_sides.fac eset got 1 fs1_bot.faceset got 2 fs2_top.faceset got 3 fs3_sides.faceset ExodusII: Start writing to file: mesh_fs3.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: 3 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: 336 ------------------------------------------ Done ExodusII Side Sets Total: 3 ExodusII: Done writing to ExodusII file: mesh_fs3.exo using cmo: momesh # GMV can display facesets as surfaces dump/gmv/mesh_exo.gmv/CMO_IN cmo/modatt/-def-/-def-/ioflag/x finish finish LaGriT successfully completed