* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * Program: LaGriT V3.202 Linux m64 * * * date_compile: 2017/04/20 mtools * * * Run Time: 2017/Apr 30 13:32:51 * * * 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 # hex mesh # with itetclr and imt set to material colors # write facesets for box shape define INPUT_MESH hex.gmv read/INPUT_MESH/moin read/gmv/hex.gmv/moin Reading GMV binary file gmvinputieee no geometry for mesh object geniee finish Nodes: 2704 Elements: 1875 Done reading GMV binary total bytes: 282424 file: hex.gmv finish cmo printatt moin -all- minmax ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH -def- 0.000000000E+00 0.000000000E+00 0.000000000E+00 2704 scalar 1 1 0 1 vector 3 3 0 1 nnodes 2704 2704 0 1 nedges 0 0 0 1 nfaces 0 0 0 1 nelements 1875 1875 0 1 mbndry 16000000 16000000 0 1 ndimensions_topo 3 3 0 1 ndimensions_geom 3 3 0 1 nodes_per_element 8 8 0 1 edges_per_element 12 12 0 1 faces_per_element 6 6 0 1 isetwd 0 0 0 2704 ialias 0 0 0 2704 imt1 1 3 2 2704 itp1 0 12 12 2704 icr1 0 1 1 2704 isn1 0 0 0 2704 xic 0.000000000E+00 1.000000000E+02 1.000000000E+02 2704 yic 0.000000000E+00 1.000000000E+02 1.000000000E+02 2704 zic 0.000000000E+00 3.000000000E+00 3.000000000E+00 2704 xtetwd 0 0 0 1875 itetclr 1 3 2 1875 itettyp 8 8 0 1875 itetoff 0 14992 14992 1875 jtetoff 0 11244 11244 1875 itet 1 2704 2703 1875x8 jtet 2 16011226 16011224 1875x6 epsilon 1.000000004E-15 1.000000004E-15 0.000000000E+00 1 epsilonl 3.140891380E-11 3.140891380E-11 0.000000000E+00 1 epsilona 4.442890500E-09 4.442890500E-09 0.000000000E+00 1 epsilonv 6.661338148E-09 6.661338148E-09 0.000000000E+00 1 ipointi 1 1 0 1 ipointj 2704 2704 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+02 1.000000000E+02 0.000000000E+00 1 ymax 1.000000000E+02 1.000000000E+02 0.000000000E+00 1 zmax 3.000000000E+00 3.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 cmo setatt moin imt 1 2704 values reset for attribute imt1 resetpts itp geniee finish cmo status moin The current-mesh-object(CMO) is: moin 1 Mesh Object name: moin number of nodes = 2704 number of elements = 1875 dimensions geometry = 3 element type = hex dimensions topology = 3 8 nodes 6 faces 12 edges boundary flag = 16000000 status = active NAME TYPE RANK LENGTH INTER PERSI IO VALUE 1 -def- VDOU scalar nnodes linea tempo agl 0.000E+00 2 scalar INT scalar scalar const perma l 1.000E+00 3 vector INT scalar scalar const perma l 3.000E+00 4 nnodes INT scalar scalar const perma l 2.704E+03 5 nedges INT scalar scalar const perma l 0.000E+00 6 nfaces INT scalar scalar const perma l 0.000E+00 7 nelements INT scalar scalar const perma l 1.875E+03 8 mbndry INT scalar scalar const perma l 1.600E+07 9 ndimensions_top INT scalar scalar const perma l 3.000E+00 10 ndimensions_geo INT scalar scalar const perma l 3.000E+00 11 nodes_per_eleme INT scalar scalar const perma l 8.000E+00 12 edges_per_eleme INT scalar scalar const perma l 1.200E+01 13 faces_per_eleme INT scalar scalar const perma l 6.000E+00 14 isetwd VDOU scalar nnodes or perma l 0.000E+00 15 ialias VINT scalar nnodes seque perma l 0.000E+00 16 imt1 VINT scalar nnodes max perma gal 0.000E+00 17 itp1 VINT scalar nnodes min perma gal 0.000E+00 18 icr1 VINT scalar nnodes min perma gal 0.000E+00 19 isn1 VINT scalar nnodes user perma gal 0.000E+00 20 xic VDOU scalar nnodes linea perma l 0.000E+00 21 yic VDOU scalar nnodes linea perma l 0.000E+00 22 zic VDOU scalar nnodes linea perma l 0.000E+00 23 xtetwd VDOU scalar nelements or perma l 0.000E+00 24 itetclr VINT scalar nelements user perma l 0.000E+00 25 itettyp VINT scalar nelements user perma l 0.000E+00 26 itetoff VINT scalar nelements user perma l 0.000E+00 27 jtetoff VINT scalar nelements user perma l 0.000E+00 28 itet VINT nodes_per_ nelements user perma l 0.000E+00 29 jtet VINT faces_per_ nelements user perma l 0.000E+00 30 ipolydat CHAR scalar scalar const perma l yes 31 vor2d CHAR scalar scalar const perma l yes 32 vor3d CHAR scalar scalar const perma l no 33 epsilon REAL scalar scalar const perma l 1.000E-15 34 epsilonl REAL scalar scalar const perma l 3.141E-11 35 epsilona REAL scalar scalar const perma l 4.443E-09 36 epsilonv REAL scalar scalar const perma l 6.661E-09 37 ipointi INT scalar scalar const perma l 1.000E+00 38 ipointj INT scalar scalar const perma l 2.704E+03 39 idebug INT scalar scalar const perma l 0.000E+00 40 itypconv_sm INT scalar scalar const perma l 1.000E+00 41 maxiter_sm INT scalar scalar const perma l 2.500E+01 42 tolconv_sm REAL scalar scalar const perma l 1.000E+00 43 nnfreq INT scalar scalar const perma l 1.000E+00 44 ivoronoi INT scalar scalar const perma l 1.000E+00 45 iopt2to2 INT scalar scalar const perma l 2.000E+00 46 dumptype CHAR scalar scalar const perma l binary 47 velname CHAR scalar scalar const perma l vels 48 densname CHAR scalar scalar const perma l ric 49 presname CHAR scalar scalar const perma l pic 50 enername CHAR scalar scalar const perma l eic 51 xmin REAL scalar scalar const perma l 0.000E+00 52 ymin REAL scalar scalar const perma l 0.000E+00 53 zmin REAL scalar scalar const perma l 0.000E+00 54 xmax REAL scalar scalar const perma l 1.000E+02 55 ymax REAL scalar scalar const perma l 1.000E+02 56 zmax REAL scalar scalar const perma l 3.000E+00 57 kdtree_level INT scalar scalar const perma l 0.000E+00 58 max_number_sets INT scalar scalar const perma l 6.400E+01 59 number_of_psets INT scalar scalar const perma l 0.000E+00 60 number_of_eltse INT scalar scalar const perma l 0.000E+00 61 psetnames VCHA scalar max_number_sets const perma l 62 eltsetnames VCHA scalar max_number_sets const perma l 63 geom_name CHAR scalar max_number_sets const perma l -defaultge 64 fsetnames VCHA scalar max_number_sets const perma l 65 number_of_fsets INT scalar scalar const perma l 0.000E+00 ## Extract the outside surface of the polygon mesh # The extracted surface will have the element and element face # that were extracted from the 3D mesh elements # attributes, idelem1 and idface1 are written to faceset files # with each file grouped into seperate faceset files #### IMPORTANT ################### # The cells MUST be pre-sorted by itetclr material values # Exodus will reorder elements internally, we do not want the # cells to be re-ordered. If they are, faceset correlation to # the cell numbers will no longer be correct. # Check that facesets are correct by reading the exo file into GMV # and selecting Surfaces under Display. #--------- BEFORE EXTRACT -------------------------------------- # sort based on itetclr values and cell location # secondary sort conventions are up to user # xmed, ymed, zmed will arrange into columns (after itetclr sort) 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/moin/index/ascending/ikey/itetclr xmed ymed zmed cmo/addatt/moin/ikey/vint/scalar/nelements///gax/0 finish SORT: order key written to attribute: ikey reorder/moin/ikey geniee finish dump/gmv/out_tmp_sort1.gmv/moin cmo/modatt/-def-/-def-/ioflag/x finish cmo/DELATT/moin/ikey # sort nodes based on mesh convention z, y, then x sort/moin/index/ascending/ikey/zic yic xic cmo/addatt/moin/ikey/vint/scalar/nnodes///gax/0 finish SORT: order key written to attribute: ikey reorder/moin/ikey dump gmv tmp_sort.gmv moin cmo/modatt/-def-/-def-/ioflag/x finish cmo DELATT moin xmed cmo DELATT moin ymed cmo DELATT moin zmed cmo DELATT moin ikey cmo printatt moin itetclr minmax ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH itetclr 1 3 2 1875 cmo printatt moin -xyz- minmax ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH xic 0.000000000E+00 1.000000000E+02 1.000000000E+02 2704 yic 0.000000000E+00 1.000000000E+02 1.000000000E+02 2704 zic 0.000000000E+00 3.000000000E+00 3.000000000E+00 2704 #---------------------------------------------------------------- extract/surfmesh/1 0 0/mo_surf/moin/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/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 1552 scalar 1 1 0 1 vector 3 3 0 1 nnodes 1552 1552 0 1 nedges 0 0 0 1 nfaces 0 0 0 1 nelements 1550 1550 0 1 mbndry 16000000 16000000 0 1 ndimensions_topo 2 2 0 1 ndimensions_geom 3 3 0 1 nodes_per_element 4 4 0 1 edges_per_element 4 4 0 1 faces_per_element 4 4 0 1 isetwd 0 0 0 1552 ialias 0 0 0 1552 imt1 1 1 0 1552 itp1 0 0 0 1552 icr1 0 1 1 1552 isn1 0 0 0 1552 xic 0.000000000E+00 1.000000000E+02 1.000000000E+02 1552 yic 0.000000000E+00 1.000000000E+02 1.000000000E+02 1552 zic 0.000000000E+00 3.000000000E+00 3.000000000E+00 1552 xtetwd 0 0 0 1550 itetclr 1 1 0 1550 itettyp 4 4 0 1550 itetoff 0 6196 6196 1550 jtetoff 0 6196 6196 1550 itet 1 1552 1551 1550x4 jtet 1 6200 6199 1550x4 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 1552 1552 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 1550 itetclr1 1 3 2 1550 idnode0 1 2704 2703 1552 idelem0 0 0 0 1550 idelem1 1 1875 1874 1550 facecol 1 3 2 1550 idface0 0 0 0 1550 idface1 1 6 5 1550 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 # Note settets are 1=bot,2=top,3=right,4=back,5=left,6=front # Top faces will face top,right,and left # Set to FEHM convention: # 1 = top # 2 = bottom # 3 = west left # 4 = south front # 5 = east right # 6 = north back cmo/addatt/mo_surf/id_side/vint/scalar/nelements cmo/select mo_surf settets/normal assign_color_normal epsilon for zero: 1.0000000E-15 925 out of 1550 element colors changed. cmo/copyatt/mo_surf mo_surf/id_side itetclr 1550 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 1550 dump gmv tmp_settets.gmv mo_surf cmo/modatt/-def-/-def-/ioflag/x finish # copy final colors into attribute for facesets cmo/copyatt/mo_surf/mo_surf/id_side/itetclr 1550 copied from mo_surf itetclr to -> mo_surf id_side define/FILENAME/fs_3.faceset define/SS_ID/3 infile output_faceset.mlgi # assume attribute id_side has faceset id numbers # SS_ID is defined in driver command file cmo/copy/mo_tmp/mo_surf cmo/select/mo_tmp eltset/e_keep/id_side/eq/SS_ID THE ELTSET e_keep HAS 75 ELEMENTS eltset/e_delete/not/e_keep THE ELTSET e_delete HAS 1475 ELEMENTS rmpoint/element/eltset get e_delete Dudded 1448 points plus parents, and 1475 elements removed. rmpoint/compress 1448 points removed and 0 elements removed. The new point count is 104. The new element count is 75. # make sure to remove all attributes except idelem1 and idface1 cmo/DELATT/mo_tmp/itetclr0 cmo/DELATT/mo_tmp/idnode0 cmo/DELATT/mo_tmp/idelem0 cmo/DELATT/mo_tmp/facecol cmo/DELATT/mo_tmp/itetclr1 cmo/DELATT/mo_tmp/idface0 cmo/DELATT/mo_tmp/id_side cmo/DELATT/mo_tmp/layertyp DELATT no action, attribute does not exist: layertyp 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 dump/gmv/tmp_faceset.gmv mo_tmp cmo/modatt/-def-/-def-/ioflag/x finish cmo/delete/mo_tmp Released Mesh Object: mo_tmp finish define/FILENAME/fs_4.faceset define/SS_ID/4 infile output_faceset.mlgi # assume attribute id_side has faceset id numbers # SS_ID is defined in driver command file cmo/copy/mo_tmp/mo_surf cmo/select/mo_tmp eltset/e_keep/id_side/eq/SS_ID THE ELTSET e_keep HAS 75 ELEMENTS eltset/e_delete/not/e_keep THE ELTSET e_delete HAS 1475 ELEMENTS rmpoint/element/eltset get e_delete Dudded 1448 points plus parents, and 1475 elements removed. rmpoint/compress 1448 points removed and 0 elements removed. The new point count is 104. The new element count is 75. # make sure to remove all attributes except idelem1 and idface1 cmo/DELATT/mo_tmp/itetclr0 cmo/DELATT/mo_tmp/idnode0 cmo/DELATT/mo_tmp/idelem0 cmo/DELATT/mo_tmp/facecol cmo/DELATT/mo_tmp/itetclr1 cmo/DELATT/mo_tmp/idface0 cmo/DELATT/mo_tmp/id_side cmo/DELATT/mo_tmp/layertyp DELATT no action, attribute does not exist: layertyp 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 dump/gmv/tmp_faceset.gmv mo_tmp cmo/modatt/-def-/-def-/ioflag/x finish cmo/delete/mo_tmp Released Mesh Object: mo_tmp finish define/FILENAME/fs_5.faceset define/SS_ID/5 infile output_faceset.mlgi # assume attribute id_side has faceset id numbers # SS_ID is defined in driver command file cmo/copy/mo_tmp/mo_surf cmo/select/mo_tmp eltset/e_keep/id_side/eq/SS_ID THE ELTSET e_keep HAS 75 ELEMENTS eltset/e_delete/not/e_keep THE ELTSET e_delete HAS 1475 ELEMENTS rmpoint/element/eltset get e_delete Dudded 1448 points plus parents, and 1475 elements removed. rmpoint/compress 1448 points removed and 0 elements removed. The new point count is 104. The new element count is 75. # make sure to remove all attributes except idelem1 and idface1 cmo/DELATT/mo_tmp/itetclr0 cmo/DELATT/mo_tmp/idnode0 cmo/DELATT/mo_tmp/idelem0 cmo/DELATT/mo_tmp/facecol cmo/DELATT/mo_tmp/itetclr1 cmo/DELATT/mo_tmp/idface0 cmo/DELATT/mo_tmp/id_side cmo/DELATT/mo_tmp/layertyp DELATT no action, attribute does not exist: layertyp 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 dump/gmv/tmp_faceset.gmv mo_tmp cmo/modatt/-def-/-def-/ioflag/x finish cmo/delete/mo_tmp Released Mesh Object: mo_tmp finish define/FILENAME/fs_6.faceset define/SS_ID/6 infile output_faceset.mlgi # assume attribute id_side has faceset id numbers # SS_ID is defined in driver command file cmo/copy/mo_tmp/mo_surf cmo/select/mo_tmp eltset/e_keep/id_side/eq/SS_ID THE ELTSET e_keep HAS 75 ELEMENTS eltset/e_delete/not/e_keep THE ELTSET e_delete HAS 1475 ELEMENTS rmpoint/element/eltset get e_delete Dudded 1448 points plus parents, and 1475 elements removed. rmpoint/compress 1448 points removed and 0 elements removed. The new point count is 104. The new element count is 75. # make sure to remove all attributes except idelem1 and idface1 cmo/DELATT/mo_tmp/itetclr0 cmo/DELATT/mo_tmp/idnode0 cmo/DELATT/mo_tmp/idelem0 cmo/DELATT/mo_tmp/facecol cmo/DELATT/mo_tmp/itetclr1 cmo/DELATT/mo_tmp/idface0 cmo/DELATT/mo_tmp/id_side cmo/DELATT/mo_tmp/layertyp DELATT no action, attribute does not exist: layertyp 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 dump/gmv/tmp_faceset.gmv mo_tmp cmo/modatt/-def-/-def-/ioflag/x finish cmo/delete/mo_tmp Released Mesh Object: mo_tmp finish define/FILENAME/fs_1.faceset define/SS_ID/1 infile output_faceset.mlgi # assume attribute id_side has faceset id numbers # SS_ID is defined in driver command file cmo/copy/mo_tmp/mo_surf cmo/select/mo_tmp eltset/e_keep/id_side/eq/SS_ID THE ELTSET e_keep HAS 625 ELEMENTS eltset/e_delete/not/e_keep THE ELTSET e_delete HAS 925 ELEMENTS rmpoint/element/eltset get e_delete Dudded 876 points plus parents, and 925 elements removed. rmpoint/compress 876 points removed and 0 elements removed. The new point count is 676. The new element count is 625. # make sure to remove all attributes except idelem1 and idface1 cmo/DELATT/mo_tmp/itetclr0 cmo/DELATT/mo_tmp/idnode0 cmo/DELATT/mo_tmp/idelem0 cmo/DELATT/mo_tmp/facecol cmo/DELATT/mo_tmp/itetclr1 cmo/DELATT/mo_tmp/idface0 cmo/DELATT/mo_tmp/id_side cmo/DELATT/mo_tmp/layertyp DELATT no action, attribute does not exist: layertyp 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 dump/gmv/tmp_faceset.gmv mo_tmp cmo/modatt/-def-/-def-/ioflag/x finish cmo/delete/mo_tmp Released Mesh Object: mo_tmp finish define/FILENAME/fs_2.faceset define/SS_ID/2 infile output_faceset.mlgi # assume attribute id_side has faceset id numbers # SS_ID is defined in driver command file cmo/copy/mo_tmp/mo_surf cmo/select/mo_tmp eltset/e_keep/id_side/eq/SS_ID THE ELTSET e_keep HAS 625 ELEMENTS eltset/e_delete/not/e_keep THE ELTSET e_delete HAS 925 ELEMENTS rmpoint/element/eltset get e_delete Dudded 876 points plus parents, and 925 elements removed. rmpoint/compress 876 points removed and 0 elements removed. The new point count is 676. The new element count is 625. # make sure to remove all attributes except idelem1 and idface1 cmo/DELATT/mo_tmp/itetclr0 cmo/DELATT/mo_tmp/idnode0 cmo/DELATT/mo_tmp/idelem0 cmo/DELATT/mo_tmp/facecol cmo/DELATT/mo_tmp/itetclr1 cmo/DELATT/mo_tmp/idface0 cmo/DELATT/mo_tmp/id_side cmo/DELATT/mo_tmp/layertyp DELATT no action, attribute does not exist: layertyp 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 dump/gmv/tmp_faceset.gmv mo_tmp cmo/modatt/-def-/-def-/ioflag/x finish cmo/delete/mo_tmp Released Mesh Object: mo_tmp finish #---------------------------------------------------------------- ## CHECK that top number of faces = bottom = layer #---------------------------------------------------------------- ############################ # write some informational files and output dump/gmv/surf_sides.gmv/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 1552 scalar 1 1 0 1 vector 3 3 0 1 nnodes 1552 1552 0 1 nedges 0 0 0 1 nfaces 0 0 0 1 nelements 1550 1550 0 1 mbndry 16000000 16000000 0 1 ndimensions_topo 2 2 0 1 ndimensions_geom 3 3 0 1 nodes_per_element 4 4 0 1 edges_per_element 4 4 0 1 faces_per_element 4 4 0 1 isetwd 0 0 0 1552 ialias 0 0 0 1552 imt1 1 1 0 1552 itp1 0 0 0 1552 icr1 0 1 1 1552 isn1 0 0 0 1552 xic 0.000000000E+00 1.000000000E+02 1.000000000E+02 1552 yic 0.000000000E+00 1.000000000E+02 1.000000000E+02 1552 zic 0.000000000E+00 3.000000000E+00 3.000000000E+00 1552 xtetwd 0 0 0 1550 itetclr 1 6 5 1550 itettyp 4 4 0 1550 itetoff 0 6196 6196 1550 jtetoff 0 6196 6196 1550 itet 1 1552 1551 1550x4 jtet 1 6200 6199 1550x4 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 1552 1552 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 1550 itetclr1 1 3 2 1550 idnode0 1 2704 2703 1552 idelem0 0 0 0 1550 idelem1 1 1875 1874 1550 facecol 1 3 2 1550 idface0 0 0 0 1550 idface1 1 6 5 1550 ncon50 0 0 0 1 nconbnd 0 0 0 1 0 length attribute: icontab id_side 1 6 5 1550 cmo printatt moin -all- minmax ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH -def- 0.000000000E+00 0.000000000E+00 0.000000000E+00 2704 scalar 1 1 0 1 vector 3 3 0 1 nnodes 2704 2704 0 1 nedges 0 0 0 1 nfaces 0 0 0 1 nelements 1875 1875 0 1 mbndry 16000000 16000000 0 1 ndimensions_topo 3 3 0 1 ndimensions_geom 3 3 0 1 nodes_per_element 8 8 0 1 edges_per_element 12 12 0 1 faces_per_element 6 6 0 1 isetwd 0 0 0 2704 ialias 0 0 0 2704 imt1 1 1 0 2704 itp1 0 12 12 2704 icr1 0 1 1 2704 isn1 0 0 0 2704 xic 0.000000000E+00 1.000000000E+02 1.000000000E+02 2704 yic 0.000000000E+00 1.000000000E+02 1.000000000E+02 2704 zic 0.000000000E+00 3.000000000E+00 3.000000000E+00 2704 xtetwd 0 0 0 1875 itetclr 1 3 2 1875 itettyp 8 8 0 1875 itetoff 0 14992 14992 1875 jtetoff 0 11244 11244 1875 itet 1 2704 2703 1875x8 jtet 2 16011248 16011246 1875x6 epsilon 1.000000004E-15 1.000000004E-15 0.000000000E+00 1 epsilonl 3.140891380E-11 3.140891380E-11 0.000000000E+00 1 epsilona 4.442890500E-09 4.442890500E-09 0.000000000E+00 1 epsilonv 6.661338148E-09 6.661338148E-09 0.000000000E+00 1 ipointi 1 1 0 1 ipointj 2704 2704 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+02 1.000000000E+02 0.000000000E+00 1 ymax 1.000000000E+02 1.000000000E+02 0.000000000E+00 1 zmax 3.000000000E+00 3.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 cmo printatt moin -xyz- minmax ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH xic 0.000000000E+00 1.000000000E+02 1.000000000E+02 2704 yic 0.000000000E+00 1.000000000E+02 1.000000000E+02 2704 zic 0.000000000E+00 3.000000000E+00 3.000000000E+00 2704 cmo select moin quality volume itetclr QUALITY TOTALS for 3 MATERIALS: ........................................................... eltset/-etmp-/itetclr/eq/1/ THE ELTSET -etmp- HAS 1650 ELEMENTS finish epsilonvol: 6.6613381E-09 --------------------------------------- VOLUMES for ELEMENT COLOR 1 ........................... All elements have volume 1.6000000E+01 eltset/-etmp-/delete/ ELTSET -etmp- DELETED finish eltset/-etmp-/itetclr/eq/2/ THE ELTSET -etmp- HAS 135 ELEMENTS finish --------------------------------------- VOLUMES for ELEMENT COLOR 2 ........................... All elements have volume 1.6000000E+01 eltset/-etmp-/delete/ ELTSET -etmp- DELETED finish eltset/-etmp-/itetclr/eq/3/ THE ELTSET -etmp- HAS 90 ELEMENTS finish --------------------------------------- VOLUMES for ELEMENT COLOR 3 ........................... All elements have volume 1.6000000E+01 eltset/-etmp-/delete/ ELTSET -etmp- DELETED finish --------------------------------------- SUMMARY VOLUME/AREA of ELEMENT COLORS for moin 1 3 ........................................................... Color Num. Elements Volume Fractional Volume 1 1650 0.2640000E+05 0.880000000 2 135 0.2160000E+04 0.072000000 3 90 0.1440000E+04 0.048000000 Total elements: 1875 Total Volume: 0.3000000E+05 ----------------------------------------------------------- 1875 total elements evaluated. ############################ # write final mesh with multi material and all facesets # node imt does not seem to be included in exo files dump/exo/mesh_fsets.exo/moin///facesets fs_1.faceset fs_2.faceset fs_3.faceset f s_4.faceset fs_5.faceset fs_6.faceset ExodusII: Start writing to file: mesh_fsets.exo using cmo: moin cmo/set_id/moin/element/e_num_temp cmo/addatt/moin e_num_temp/VINT/scalar/nelements/linear/permanent//0 finish finish sort/moin/index/ascending/ikey_utr/itetclr/e_num_temp cmo/addatt/moin/ikey_utr/vint/scalar/nelements///gax/0 finish SORT: order key written to attribute: ikey_utr finish cmo/DELATT/moin/e_num_temp finish cmo/modatt/moin/itetclr/ioflag/lga finish INITIALIZE exodus Title: LAGRIT TO EXODUS number of dimension: 3 number of nodes: 2704 number of elements: 1875 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: 6 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/moin/itetclr/ioflag/l finish ExodusII: Done writing to ExodusII file: mesh_fsets.exo using cmo: moin dump/avs2/mesh_fsets.inp/moin cmo/modatt/-def-/-def-/ioflag/x finish dump/gmv/mesh_fsets.gmv/moin cmo/modatt/-def-/-def-/ioflag/x finish finish LaGriT successfully completed