* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * Program: LaGriT V3.101 Linux m64 * * * date_compile: 2013/04/25 RH gf4.5 * * * Run Time: 2014/Mar 17 13:40:38 * * * 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: 381 Elements: 126 Done reading GMV binary total bytes: 40704 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 = 381 number of elements = 126 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.810E+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.260E+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.110E-11 35 epsilona REAL scalar scalar const perma l 5.551E-10 36 epsilonv REAL scalar scalar const perma l 1.110E-11 37 ipointi INT scalar scalar const perma l 1.000E+00 38 ipointj INT scalar scalar const perma l 3.810E+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 1.000E+00 56 zmax REAL scalar scalar const perma l 4.980E+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 381 scalar 1 1 0 1 vector 3 3 0 1 nnodes 381 381 0 1 nedges 0 0 0 1 nfaces 0 0 0 1 nelements 380 380 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 381 ialias 0 0 0 381 imt1 1 1 0 381 itp1 0 2 2 381 icr1 0 0 0 381 isn1 0 0 0 381 xic 0.000000000E+00 1.000000000E+00 1.000000000E+00 381 yic 0.000000000E+00 1.000000000E+00 1.000000000E+00 381 zic -4.500000000E+01 4.980000019E+00 4.998000002E+01 381 xtetwd 0 0 0 380 itetclr 1 1 0 380 itettyp 3 4 1 380 itetoff 0 1514 1514 380 jtetoff 0 1514 1514 380 itet 0 381 381 380x* jtet 0 3794 3794 380x* 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 381 381 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 380 itetclr1 2 6 4 380 idnode0 1 381 380 381 idelem0 0 0 0 380 idelem1 1 126 125 380 facecol 1 5 4 380 idface0 0 0 0 380 idface1 1 5 4 380 layertyp -2.000000000E+00 2.000000000E+00 4.000000000E+00 381 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 379 out of 380 element colors changed. cmo/copyatt/mo_surf mo_surf/id_side itetclr 380 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 11 10 380 # 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 0 ELEMENTS eltset/e_delete/not/e_keep THE ELTSET e_delete HAS 380 ELEMENTS rmpoint/element/eltset get e_delete Dudded 381 points plus parents, there are no elements in mesh. rmpoint/compress 381 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_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 0 ELEMENTS eltset/e_delete/not/e_keep THE ELTSET e_delete HAS 380 ELEMENTS rmpoint/element/eltset get e_delete Dudded 381 points plus parents, there are no elements in mesh. rmpoint/compress 381 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_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 126 ELEMENTS eltset/e_delete/not/e_keep THE ELTSET e_delete HAS 254 ELEMENTS rmpoint/element/eltset get e_delete Dudded 127 points plus parents, and 254 elements removed. rmpoint/compress 127 points removed and 0 elements removed. The new point count is 254. The new element count is 126. # 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_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 126 ELEMENTS eltset/e_delete/not/e_keep THE ELTSET e_delete HAS 254 ELEMENTS rmpoint/element/eltset get e_delete Dudded 127 points plus parents, and 254 elements removed. rmpoint/compress 127 points removed and 0 elements removed. The new point count is 254. The new element count is 126. # 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 379 ELEMENTS rmpoint/element/eltset get e_delete Dudded 378 points plus parents, and 379 elements removed. rmpoint/compress 378 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 379 ELEMENTS rmpoint/element/eltset get e_delete Dudded 378 points plus parents, and 379 elements removed. rmpoint/compress 378 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 381 scalar 1 1 0 1 vector 3 3 0 1 nnodes 381 381 0 1 nedges 0 0 0 1 nfaces 0 0 0 1 nelements 380 380 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 381 ialias 0 0 0 381 imt1 1 1 0 381 itp1 0 2 2 381 icr1 0 0 0 381 isn1 0 0 0 381 xic 0.000000000E+00 1.000000000E+00 1.000000000E+00 381 yic 0.000000000E+00 1.000000000E+00 1.000000000E+00 381 zic -4.500000000E+01 4.980000019E+00 4.998000002E+01 381 xtetwd 0 0 0 380 itetclr 1 11 10 380 itettyp 3 4 1 380 itetoff 0 1514 1514 380 jtetoff 0 1514 1514 380 itet 0 381 381 380x* jtet 0 3794 3794 380x* 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 381 381 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 380 itetclr1 2 6 4 380 idnode0 1 381 380 381 idelem0 0 0 0 380 idelem1 1 126 125 380 facecol 1 5 4 380 idface0 0 0 0 380 idface1 1 5 4 380 layertyp -2.000000000E+00 2.000000000E+00 4.000000000E+00 381 ncon50 0 0 0 1 nconbnd 0 0 0 1 0 length attribute: icontab id_side 1 11 10 380 cmo printatt mo_pri -all- minmax ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH -def- 0.000000000E+00 0.000000000E+00 0.000000000E+00 381 scalar 1 1 0 1 vector 3 3 0 1 nnodes 381 381 0 1 nedges 0 0 0 1 nfaces 0 0 0 1 nelements 126 126 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 381 ialias 0 0 0 381 imt1 2 6 4 381 itp1 10 10 0 381 icr1 0 0 0 381 isn1 0 0 0 381 xic 0.000000000E+00 1.000000000E+00 1.000000000E+00 381 yic 0.000000000E+00 1.000000000E+00 1.000000000E+00 381 zic -4.500000000E+01 4.980000019E+00 4.998000002E+01 381 xtetwd 0 0 0 126 itetclr 2 6 4 126 itettyp 7 7 0 126 itetoff 0 750 750 126 jtetoff 0 625 625 126 itet 1 381 380 126x6 jtet 2 16000627 16000625 126x5 epsilon 1.000000004E-15 1.000000004E-15 0.000000000E+00 1 epsilonl 1.110223114E-11 1.110223114E-11 0.000000000E+00 1 epsilona 5.551116016E-10 5.551116016E-10 0.000000000E+00 1 epsilonv 1.109778936E-11 1.109778936E-11 0.000000000E+00 1 ipointi 1 1 0 1 ipointj 381 381 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 1.000000000E+00 1.000000000E+00 0.000000000E+00 1 zmax 4.980000019E+00 4.980000019E+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 381 cmo printatt mo_pri -xyz- minmax ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH xic 0.000000000E+00 1.000000000E+00 1.000000000E+00 381 yic 0.000000000E+00 1.000000000E+00 1.000000000E+00 381 zic -4.500000000E+01 4.980000019E+00 4.998000002E+01 381 cmo select mo_pri quality volume itetclr QUALITY TOTALS for 5 MATERIALS: ........................................................... eltset/-etmp-/itetclr/eq/2/ THE ELTSET -etmp- HAS 6 ELEMENTS finish epsilonvol: 1.1097789E-11 --------------------------------------- VOLUMES for ELEMENT COLOR 2 ........................... All elements have volume 9.9999905E-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.1000E-01 and 0.1000E-01: 8 element volumes b/w 0.1000E-01 and 0.1000E-01: 0 element volumes b/w 0.1000E-01 and 0.1000E-01: 0 element volumes b/w 0.1000E-01 and 0.1000E-01: 0 element volumes b/w 0.1000E-01 and 0.1000E-01: 1 min volume = 9.9999905E-03 max volume = 1.0000229E-02 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 9.9999905E-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.1250E-01 and 0.1250E-01: 4 element volumes b/w 0.1250E-01 and 0.1250E-01: 0 element volumes b/w 0.1250E-01 and 0.1250E-01: 0 element volumes b/w 0.1250E-01 and 0.1250E-01: 0 element volumes b/w 0.1250E-01 and 0.1250E-01: 16 min volume = 1.2499809E-02 max volume = 1.2500048E-02 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.5000E-01 and 0.9106E-01: 50 element volumes b/w 0.9106E-01 and 0.1659E+00: 0 element volumes b/w 0.1659E+00 and 0.3021E+00: 10 element volumes b/w 0.3021E+00 and 0.5502E+00: 10 element volumes b/w 0.5502E+00 and 0.1002E+01: 15 min volume = 4.9999952E-02 max volume = 1.0020018E+00 eltset/-etmp-/delete/ ELTSET -etmp- DELETED finish --------------------------------------- SUMMARY VOLUME/AREA of ELEMENT COLORS for mo_pri 1 5 ........................................................... Color Num. Elements Volume Fractional Volume 2 6 0.5999994E-01 0.002400958 3 9 0.9000015E-01 0.003601447 4 6 0.5999994E-01 0.002400958 5 20 0.2500000E+00 0.010004002 6 85 0.2453000E+02 0.981592636 Total elements: 126 Total Volume: 0.2499000E+02 ----------------------------------------------------------- 126 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_5 _lef.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: 381 number of elements: 126 number of edges: 0 number of edge blocks: 0 number of element blocks: 5 number of face blocks: 0 number of node sets: 0 number of edge sets: 3 number of element sets: 0 number of side sets: 4 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 126 values reset for attribute itetclr cmo setatt mo_pri imt 1 381 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_5_lef.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: 381 number of elements: 126 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: 4 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 381 scalar 1 1 0 1 vector 3 3 0 1 nnodes 381 381 0 1 nedges 0 0 0 1 nfaces 0 0 0 1 nelements 126 126 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 381 ialias 0 0 0 381 imt1 1 1 0 381 itp1 10 10 0 381 icr1 0 0 0 381 isn1 0 0 0 381 xic 0.000000000E+00 1.000000000E+00 1.000000000E+00 381 yic 0.000000000E+00 1.000000000E+00 1.000000000E+00 381 zic -4.500000000E+01 4.980000019E+00 4.998000002E+01 381 xtetwd 0 0 0 126 itetclr 1 1 0 126 itettyp 7 7 0 126 itetoff 0 750 750 126 jtetoff 0 625 625 126 itet 1 381 380 126x6 jtet 1 16000000 15999999 126x5 epsilon 1.000000004E-15 1.000000004E-15 0.000000000E+00 1 epsilonl 1.110223114E-11 1.110223114E-11 0.000000000E+00 1 epsilona 5.551116016E-10 5.551116016E-10 0.000000000E+00 1 epsilonv 1.109778936E-11 1.109778936E-11 0.000000000E+00 1 ipointi 1 1 0 1 ipointj 381 381 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 1.000000000E+00 1.000000000E+00 0.000000000E+00 1 zmax 4.980000019E+00 4.980000019E+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 381 ikey_utr 1 126 125 126 finish LaGriT successfully completed