* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * Program: LaGriT V3.101 Linux m64 * * * date_compile: 2013/04/25 RH gf4.5 * * * Run Time: 2013/Nov 14 13:17:27 * * * 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 the file with the 3D polygon mesh # this is a triangle mesh stacked into prism elements # such that tri faces show top and bottom, quads on sides # itetclr materials have been defined, facesets have not read/full_mesh.gmv/mo_pri read/gmv/full_mesh.gmv/mo_pri Reading GMV binary file gmvinputieee cmo/addatt/mo_pri/layertyp/VDOUBLE/scalar/nnodes/linear/permanent/gxa/0.0 finish no geometry for mesh object geniee finish Nodes: 396 Elements: 131 Done reading GMV binary total bytes: 42372 file: full_mesh.gmv finish cmo status mo_pri The current-mesh-object(CMO) is: mo_pri 1 Mesh Object name: mo_pri number of nodes = 396 number of elements = 131 dimensions geometry = 3 element type = pri dimensions topology = 3 6 nodes 5 faces 9 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 3.960E+02 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.310E+02 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 6.000E+00 12 edges_per_eleme INT scalar scalar const perma l 9.000E+00 13 faces_per_eleme INT scalar scalar const perma l 5.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 1.111E-11 35 epsilona REAL scalar scalar const perma l 5.554E-10 36 epsilonv REAL scalar scalar const perma l 7.772E-12 37 ipointi INT scalar scalar const perma l 1.000E+00 38 ipointj INT scalar scalar const perma l 3.960E+02 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 -4.500E+01 54 xmax REAL scalar scalar const perma l 1.000E+00 55 ymax REAL scalar scalar const perma l 7.000E-01 56 zmax REAL scalar scalar const perma l 5.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 66 layertyp VDOU scalar nnodes linea perma gxa 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/mo_pri/index/ascending/ikey/itetclr xmed ymed zmed cmo/addatt/mo_pri/ikey/vint/scalar/nelements///gax/0 finish SORT: order key written to attribute: ikey reorder/mo_pri/ikey geniee finish dump/gmv/out_tmp_sort1.gmv/mo_pri cmo/modatt/-def-/-def-/ioflag/x finish cmo/DELATT/mo_pri/ikey # sort nodes based on mesh convention z, y, then x sort/mo_pri/index/ascending/ikey/zic yic xic cmo/addatt/mo_pri/ikey/vint/scalar/nnodes///gax/0 finish SORT: order key written to attribute: ikey reorder/mo_pri/ikey dump gmv tmp_sort.gmv mo_pri cmo/modatt/-def-/-def-/ioflag/x finish cmo DELATT mo_pri xmed cmo DELATT mo_pri ymed cmo DELATT mo_pri zmed cmo DELATT mo_pri ikey #---------------------------------------------------------------- extract/surfmesh/1 0 0/mo_surf/mo_pri/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/VDOUBLE/scalar/nnodes/linear/permanent/gxa 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 396 scalar 1 1 0 1 vector 3 3 0 1 nnodes 396 396 0 1 nedges 0 0 0 1 nfaces 0 0 0 1 nelements 395 395 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 396 ialias 0 0 0 396 imt1 1 1 0 396 itp1 0 2 2 396 icr1 0 0 0 396 isn1 0 0 0 396 xic 0.000000000E+00 1.000000000E+00 1.000000000E+00 396 yic 0.000000000E+00 6.999999881E-01 6.999999881E-01 396 zic -4.500000000E+01 5.000000000E+00 5.000000000E+01 396 xtetwd 0 0 0 395 itetclr 1 1 0 395 itettyp 3 4 1 395 itetoff 0 1574 1574 395 jtetoff 0 1574 1574 395 itet 0 396 396 395x* jtet 0 3944 3944 395x* 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 396 396 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 395 itetclr1 1 6 5 395 idnode0 1 396 395 396 idelem0 0 0 0 395 idelem1 1 131 130 395 facecol 1 6 5 395 idface0 0 0 0 395 idface1 1 5 4 395 layertyp -2.000000000E+00 2.000000000E+00 4.000000000E+00 396 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 # 3 = east right # 4 = north back # 5 = west left # 6 = south front cmo/addatt/mo_surf/id_side/vint/scalar/nelements cmo/select mo_surf settets/normal assign_color_normal epsilon for zero: 1.0000000E-15 394 out of 395 element colors changed. cmo/copyatt/mo_surf mo_surf/id_side itetclr 395 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 395 # this was previous convention, use settets convention instead. # so that top is now 2 and 1 is bottom, etc. # the settets convention is more widely used # cmo / setatt / mo_surf / id_side / eltset get etop / 1 # cmo / setatt / mo_surf / id_side / eltset get ebot / 2 # cmo / setatt / mo_surf / id_side / eltset get efro / 3 # cmo / setatt / mo_surf / id_side / eltset get ebac / 4 # cmo / setatt / mo_surf / id_side / eltset get elef / 5 # cmo / setatt / mo_surf / id_side / eltset get erig / 6 define/FILENAME/fs_3_rig.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 131 ELEMENTS eltset/e_delete/not/e_keep THE ELTSET e_delete HAS 264 ELEMENTS rmpoint/element/eltset get e_delete Dudded 132 points plus parents, and 264 elements removed. rmpoint/compress 132 points removed and 0 elements removed. The new point count is 264. The new element count is 131. # 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 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_bac.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 131 ELEMENTS eltset/e_delete/not/e_keep THE ELTSET e_delete HAS 264 ELEMENTS rmpoint/element/eltset get e_delete Dudded 132 points plus parents, and 264 elements removed. rmpoint/compress 132 points removed and 0 elements removed. The new point count is 264. The new element count is 131. # 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 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_lef.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 0 ELEMENTS eltset/e_delete/not/e_keep THE ELTSET e_delete HAS 395 ELEMENTS rmpoint/element/eltset get e_delete Dudded 396 points plus parents, there are no elements in mesh. rmpoint/compress 396 points removed and 0 elements removed. The new point count is 0. The new element count is 0. # 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 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. WARNING: dumpavs WARNING: nnodes=0 nelements = 0 WARNING: No output dump/gmv/tmp_faceset.gmv mo_tmp DUMP/GMV early exit: no nodes in cmo: mo_tmp DUMP/GMV ERROR: could not open file for: mo_tmp cmo/delete/mo_tmp Released Mesh Object: mo_tmp finish define/FILENAME/fs_6_fro.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 131 ELEMENTS eltset/e_delete/not/e_keep THE ELTSET e_delete HAS 264 ELEMENTS rmpoint/element/eltset get e_delete Dudded 132 points plus parents, and 264 elements removed. rmpoint/compress 132 points removed and 0 elements removed. The new point count is 264. The new element count is 131. # 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 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_bot.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 1 ELEMENTS eltset/e_delete/not/e_keep THE ELTSET e_delete HAS 394 ELEMENTS rmpoint/element/eltset get e_delete Dudded 393 points plus parents, and 394 elements removed. rmpoint/compress 393 points removed and 0 elements removed. The new point count is 3. The new element count is 1. # 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 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_top.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 1 ELEMENTS eltset/e_delete/not/e_keep THE ELTSET e_delete HAS 394 ELEMENTS rmpoint/element/eltset get e_delete Dudded 393 points plus parents, and 394 elements removed. rmpoint/compress 393 points removed and 0 elements removed. The new point count is 3. The new element count is 1. # 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 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 # this version layer has number of elements = 6272 #---------------------------------------------------------------- ######################## Add Top Corners Section ################ # Find facesets at top corners for Boundary Conditions # Using x,y points of the input polygons (before triangulation) # Points are found on outside surface of full mesh, no internal nodes # Write exodus file to check corner point facesets # fs_7_Csides.faceset and fs_8_Ctop.faceset # (exodus assigns side set id based on file count) #define FPOLY_PTS boundary_points.inp #infile bc_corner_faceset.mlgi #dump / exo / mesh_bc.exo / mo_pri / / / facesets & # fs_7_Csides.faceset fs_8_Ctop.faceset ############################ # 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 396 scalar 1 1 0 1 vector 3 3 0 1 nnodes 396 396 0 1 nedges 0 0 0 1 nfaces 0 0 0 1 nelements 395 395 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 396 ialias 0 0 0 396 imt1 1 1 0 396 itp1 0 2 2 396 icr1 0 0 0 396 isn1 0 0 0 396 xic 0.000000000E+00 1.000000000E+00 1.000000000E+00 396 yic 0.000000000E+00 6.999999881E-01 6.999999881E-01 396 zic -4.500000000E+01 5.000000000E+00 5.000000000E+01 396 xtetwd 0 0 0 395 itetclr 1 6 5 395 itettyp 3 4 1 395 itetoff 0 1574 1574 395 jtetoff 0 1574 1574 395 itet 0 396 396 395x* jtet 0 3944 3944 395x* 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 396 396 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 395 itetclr1 1 6 5 395 idnode0 1 396 395 396 idelem0 0 0 0 395 idelem1 1 131 130 395 facecol 1 6 5 395 idface0 0 0 0 395 idface1 1 5 4 395 layertyp -2.000000000E+00 2.000000000E+00 4.000000000E+00 396 ncon50 0 0 0 1 nconbnd 0 0 0 1 0 length attribute: icontab id_side 1 6 5 395 cmo printatt mo_pri -all- minmax ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH -def- 0.000000000E+00 0.000000000E+00 0.000000000E+00 396 scalar 1 1 0 1 vector 3 3 0 1 nnodes 396 396 0 1 nedges 0 0 0 1 nfaces 0 0 0 1 nelements 131 131 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 396 ialias 0 0 0 396 imt1 1 6 5 396 itp1 10 10 0 396 icr1 0 0 0 396 isn1 0 0 0 396 xic 0.000000000E+00 1.000000000E+00 1.000000000E+00 396 yic 0.000000000E+00 6.999999881E-01 6.999999881E-01 396 zic -4.500000000E+01 5.000000000E+00 5.000000000E+01 396 xtetwd 0 0 0 131 itetclr 1 6 5 131 itettyp 7 7 0 131 itetoff 0 780 780 131 jtetoff 0 650 650 131 itet 1 396 395 131x6 jtet 2 16000652 16000650 131x5 epsilon 1.000000004E-15 1.000000004E-15 0.000000000E+00 1 epsilonl 1.110553822E-11 1.110553822E-11 0.000000000E+00 1 epsilona 5.554423588E-10 5.554423588E-10 0.000000000E+00 1 epsilonv 7.771561040E-12 7.771561040E-12 0.000000000E+00 1 ipointi 1 1 0 1 ipointj 396 396 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 -4.500000000E+01 -4.500000000E+01 0.000000000E+00 1 xmax 1.000000000E+00 1.000000000E+00 0.000000000E+00 1 ymax 6.999999881E-01 6.999999881E-01 0.000000000E+00 1 zmax 5.000000000E+00 5.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 layertyp -2.000000000E+00 2.000000000E+00 4.000000000E+00 396 cmo printatt mo_pri -xyz- minmax ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH xic 0.000000000E+00 1.000000000E+00 1.000000000E+00 396 yic 0.000000000E+00 6.999999881E-01 6.999999881E-01 396 zic -4.500000000E+01 5.000000000E+00 5.000000000E+01 396 cmo select mo_pri quality volume itetclr QUALITY TOTALS for 6 MATERIALS: ........................................................... eltset/-etmp-/itetclr/eq/1/ THE ELTSET -etmp- HAS 2 ELEMENTS finish epsilonvol: 7.7715610E-12 --------------------------------------- VOLUMES for ELEMENT COLOR 1 ........................... element volumes b/w 0.3500E-02 and 0.3500E-02: 1 element volumes b/w 0.3500E-02 and 0.3500E-02: 0 element volumes b/w 0.3500E-02 and 0.3500E-02: 0 element volumes b/w 0.3500E-02 and 0.3500E-02: 0 element volumes b/w 0.3500E-02 and 0.3500E-02: 1 min volume = 3.4999132E-03 max volume = 3.5000800E-03 eltset/-etmp-/delete/ ELTSET -etmp- DELETED finish eltset/-etmp-/itetclr/eq/2/ THE ELTSET -etmp- HAS 9 ELEMENTS finish --------------------------------------- VOLUMES for ELEMENT COLOR 2 ........................... All elements have volume 3.4999132E-03 eltset/-etmp-/delete/ ELTSET -etmp- DELETED finish eltset/-etmp-/itetclr/eq/3/ THE ELTSET -etmp- HAS 9 ELEMENTS finish --------------------------------------- VOLUMES for ELEMENT COLOR 3 ........................... element volumes b/w 0.7000E-02 and 0.7000E-02: 8 element volumes b/w 0.7000E-02 and 0.7000E-02: 0 element volumes b/w 0.7000E-02 and 0.7000E-02: 0 element volumes b/w 0.7000E-02 and 0.7000E-02: 0 element volumes b/w 0.7000E-02 and 0.7000E-02: 1 min volume = 6.9999932E-03 max volume = 7.0001601E-03 eltset/-etmp-/delete/ ELTSET -etmp- DELETED finish eltset/-etmp-/itetclr/eq/4/ THE ELTSET -etmp- HAS 6 ELEMENTS finish --------------------------------------- VOLUMES for ELEMENT COLOR 4 ........................... All elements have volume 3.4999132E-03 eltset/-etmp-/delete/ ELTSET -etmp- DELETED finish eltset/-etmp-/itetclr/eq/5/ THE ELTSET -etmp- HAS 20 ELEMENTS finish --------------------------------------- VOLUMES for ELEMENT COLOR 5 ........................... element volumes b/w 0.8750E-02 and 0.8750E-02: 4 element volumes b/w 0.8750E-02 and 0.8750E-02: 0 element volumes b/w 0.8750E-02 and 0.8750E-02: 0 element volumes b/w 0.8750E-02 and 0.8750E-02: 0 element volumes b/w 0.8750E-02 and 0.8750E-02: 16 min volume = 8.7498663E-03 max volume = 8.7500332E-03 eltset/-etmp-/delete/ ELTSET -etmp- DELETED finish eltset/-etmp-/itetclr/eq/6/ THE ELTSET -etmp- HAS 85 ELEMENTS finish --------------------------------------- VOLUMES for ELEMENT COLOR 6 ........................... element volumes b/w 0.3500E-01 and 0.6372E-01: 50 element volumes b/w 0.6372E-01 and 0.1160E+00: 0 element volumes b/w 0.1160E+00 and 0.2112E+00: 10 element volumes b/w 0.2112E+00 and 0.3845E+00: 10 element volumes b/w 0.3845E+00 and 0.7000E+00: 15 min volume = 3.4999966E-02 max volume = 7.0000199E-01 eltset/-etmp-/delete/ ELTSET -etmp- DELETED finish --------------------------------------- SUMMARY VOLUME/AREA of ELEMENT COLORS for mo_pri 1 6 ........................................................... Color Num. Elements Volume Fractional Volume 1 2 0.6999993E-02 0.000400000 2 9 0.6299994E-01 0.003599997 3 9 0.6300011E-01 0.003600006 4 6 0.4199996E-01 0.002399998 5 20 0.1750000E+00 0.010000000 6 85 0.1715000E+02 0.980000000 Total elements: 131 Total Volume: 0.1750000E+02 ----------------------------------------------------------- 131 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/mo_pri///facesets fs_1_bot.faceset fs_2_top.faceset fs_3 _rig.faceset fs_4_bac.faceset fs_6_fro.faceset ExodusII: Start writing to file: mesh_fsets.exo using cmo: mo_pri cmo/set_id/mo_pri/element/e_num_temp cmo/addatt/mo_pri e_num_temp/VINT/scalar/nelements/linear/permanent//0 finish finish sort/mo_pri/index/ascending/ikey_utr/itetclr/e_num_temp cmo/addatt/mo_pri/ikey_utr/vint/scalar/nelements///gax/0 finish SORT: order key written to attribute: ikey_utr finish cmo/DELATT/mo_pri/e_num_temp finish cmo/modatt/mo_pri/itetclr/ioflag/lga finish INITIALIZE exodus Title: LAGRIT TO EXODUS number of dimension: 3 number of nodes: 396 number of elements: 131 number of edges: 0 number of edge blocks: 0 number of element blocks: 6 number of face blocks: 0 number of node sets: 0 number of edge sets: 3 number of element sets: 0 number of side sets: 5 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_pri/itetclr/ioflag/l finish ExodusII: Done writing to file: mesh_fsets.exo using cmo: mo_pri # fs_7_Csides.faceset fs_8_Ctop.faceset dump/avs2/mesh_fsets.inp/mo_pri cmo/modatt/-def-/-def-/ioflag/x finish dump/gmv/mesh_fsets.gmv/mo_pri cmo/modatt/-def-/-def-/ioflag/x finish ############################ # version with single material and top faceset only cmo setatt mo_pri itetclr 1 131 values reset for attribute itetclr cmo setatt mo_pri imt 1 396 values reset for attribute imt1 resetpts itp geniee finish dump/exo/mesh_fsets_1mat.exo/mo_pri///facesets fs_1_bot.faceset fs_2_top.faceset fs_3_rig.faceset fs_4_bac.faceset fs_6_fro.faceset ExodusII: Start writing to file: mesh_fsets_1mat.exo using cmo: mo_pri cmo/set_id/mo_pri/element/e_num_temp cmo/addatt/mo_pri e_num_temp/VINT/scalar/nelements/linear/permanent//0 finish finish sort/mo_pri/index/ascending/ikey_utr/itetclr/e_num_temp cmo/DELATT/mo_pri/ikey_utr finish cmo/addatt/mo_pri/ikey_utr/vint/scalar/nelements///gax/0 finish SORT: order key written to attribute: ikey_utr finish cmo/DELATT/mo_pri/e_num_temp finish cmo/modatt/mo_pri/itetclr/ioflag/lga finish INITIALIZE exodus Title: LAGRIT TO EXODUS number of dimension: 3 number of nodes: 396 number of elements: 131 number of edges: 0 number of edge blocks: 0 number of element blocks: 1 number of face blocks: 0 number of node sets: 0 number of edge sets: 3 number of element sets: 0 number of side sets: 5 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_pri/itetclr/ioflag/l finish ExodusII: Done writing to file: mesh_fsets_1mat.exo using cmo: mo_pri dump/avs2/mesh_fsets_1mat.inp/mo_pri cmo/modatt/-def-/-def-/ioflag/x finish dump/gmv/mesh_fsets_1mat.gmv/mo_pri cmo/modatt/-def-/-def-/ioflag/x finish cmo printatt mo_pri -all- minmax ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH -def- 0.000000000E+00 0.000000000E+00 0.000000000E+00 396 scalar 1 1 0 1 vector 3 3 0 1 nnodes 396 396 0 1 nedges 0 0 0 1 nfaces 0 0 0 1 nelements 131 131 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 396 ialias 0 0 0 396 imt1 1 1 0 396 itp1 10 10 0 396 icr1 0 0 0 396 isn1 0 0 0 396 xic 0.000000000E+00 1.000000000E+00 1.000000000E+00 396 yic 0.000000000E+00 6.999999881E-01 6.999999881E-01 396 zic -4.500000000E+01 5.000000000E+00 5.000000000E+01 396 xtetwd 0 0 0 131 itetclr 1 1 0 131 itettyp 7 7 0 131 itetoff 0 780 780 131 jtetoff 0 650 650 131 itet 1 396 395 131x6 jtet 1 16000000 15999999 131x5 epsilon 1.000000004E-15 1.000000004E-15 0.000000000E+00 1 epsilonl 1.110553822E-11 1.110553822E-11 0.000000000E+00 1 epsilona 5.554423588E-10 5.554423588E-10 0.000000000E+00 1 epsilonv 7.771561040E-12 7.771561040E-12 0.000000000E+00 1 ipointi 1 1 0 1 ipointj 396 396 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 -4.500000000E+01 -4.500000000E+01 0.000000000E+00 1 xmax 1.000000000E+00 1.000000000E+00 0.000000000E+00 1 ymax 6.999999881E-01 6.999999881E-01 0.000000000E+00 1 zmax 5.000000000E+00 5.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 layertyp -2.000000000E+00 2.000000000E+00 4.000000000E+00 396 ikey_utr 1 131 130 131 finish LaGriT successfully completed