* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * Program: LaGriT V3.200 Linux m64 * * * date_compile: 2016/10/25 Ubu gf * * * Run Time: 2016/Nov 30 11:22:28 * * * Manual: http://lagrit.lanl.gov * * * * * * * * * * * * * * * * * * * * * * * * * * * * * -----oOo----- LaGriT V3 LA-CC-15-069 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/proto_full_poly.inp/mo_pri read/avs/proto_full_poly.inp/mo_pri cmo/addatt/mo_pri/layertyp/VINT/scalar/nnodes/linear/permanent/gxaf/0.0 finish cmo/addatt/mo_pri/node_ice/VINT/scalar/nnodes/linear/permanent/gxaf/0.0 finish cmo/addatt/mo_pri/mat_orig/VINT/scalar/nelements/linear/permanent/gxaf/0.0 finish cmo/addatt/mo_pri/ikey_utr/VINT/scalar/nelements/linear/permanent/gxaf/0.0 finish geniee finish cmo/status/brief The current-mesh-object(CMO) is: mo_pri 1 Mesh Object name: mo_pri number of nodes = 381603 number of elements = 732160 dimensions geometry = 3 element type = pri dimensions topology = 3 6 nodes 5 faces 9 edges boundary flag = 16000000 status = active finish finish cmo status mo_pri The current-mesh-object(CMO) is: mo_pri 1 Mesh Object name: mo_pri number of nodes = 381603 number of elements = 732160 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.816E+05 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 7.322E+05 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.311E-11 35 epsilona REAL scalar scalar const perma l 7.738E-10 36 epsilonv REAL scalar scalar const perma l 5.067E-09 37 ipointi INT scalar scalar const perma l 1.000E+00 38 ipointj INT scalar scalar const perma l 3.816E+05 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 9.219E+01 52 ymin REAL scalar scalar const perma l 2.731E+00 53 zmin REAL scalar scalar const perma l -4.500E+01 54 xmax REAL scalar scalar const perma l 1.164E+02 55 ymax REAL scalar scalar const perma l 2.140E+01 56 zmax REAL scalar scalar const perma l 5.504E+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 VINT scalar nnodes linea perma gxa 0.000E+00 67 node_ice VINT scalar nnodes linea perma gxa 0.000E+00 68 mat_orig VINT scalar nelements linea perma gxa 0.000E+00 69 ikey_utr VINT scalar nelements 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/VINT/scalar/nnodes/linear/permanent/gxaf finish cmo/addatt/mo_surf/node_ice/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 30594 scalar 1 1 0 1 vector 3 3 0 1 nnodes 30594 30594 0 1 nedges 0 0 0 1 nfaces 0 0 0 1 nelements 36224 36224 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 30594 ialias 0 0 0 30594 imt1 1 1 0 30594 itp1 0 2 2 30594 icr1 0 0 0 30594 isn1 0 0 0 30594 xic 9.219250000E+01 1.163908000E+02 2.419830000E+01 30594 yic 2.731400000E+00 2.140300000E+01 1.867160000E+01 30594 zic -4.500000000E+01 5.504421902E+00 5.050442190E+01 30594 xtetwd 0 0 0 36224 itetclr 1 1 0 36224 itettyp 3 4 1 36224 itetoff 0 133628 133628 36224 jtetoff 0 133628 133628 36224 itet 0 30594 30594 36224x* jtet 0 362234 362234 36224x* 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 30594 30594 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 36224 itetclr1 1 30 29 36224 idnode0 1 381603 381602 30594 idelem0 0 0 0 36224 idelem1 1 732160 732159 36224 facecol 1 6 5 36224 idface0 0 0 0 36224 idface1 1 5 4 36224 layertyp -2 2 4 30594 node_ice 0 1 1 30594 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 30592 out of 36224 element colors changed. cmo/copyatt/mo_surf mo_surf/id_side itetclr 36224 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 36224 # 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 define/SNAME/SS_east_right 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 # report sideset number of faces eltset/SNAME/id_side/eq/SS_ID THE ELTSET SS_east_right HAS 4160 ELEMENTS eltset/e_keep/id_side/eq/SS_ID THE ELTSET e_keep HAS 4160 ELEMENTS eltset/e_delete/not/e_keep THE ELTSET e_delete HAS 32064 ELEMENTS rmpoint/element/eltset get e_delete Dudded 26271 points plus parents, and 32064 elements removed. rmpoint/compress 26271 points removed and 0 elements removed. The new point count is 4323. The new element count is 4160. # 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 define/SNAME/SS_north_back 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 # report sideset number of faces eltset/SNAME/id_side/eq/SS_ID THE ELTSET SS_north_back HAS 8320 ELEMENTS eltset/e_keep/id_side/eq/SS_ID THE ELTSET e_keep HAS 8320 ELEMENTS eltset/e_delete/not/e_keep THE ELTSET e_delete HAS 27904 ELEMENTS rmpoint/element/eltset get e_delete Dudded 22079 points plus parents, and 27904 elements removed. rmpoint/compress 22079 points removed and 0 elements removed. The new point count is 8515. The new element count is 8320. # 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 define/SNAME/SS_west_left 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 # report sideset number of faces eltset/SNAME/id_side/eq/SS_ID THE ELTSET SS_west_left HAS 4160 ELEMENTS eltset/e_keep/id_side/eq/SS_ID THE ELTSET e_keep HAS 4160 ELEMENTS eltset/e_delete/not/e_keep THE ELTSET e_delete HAS 32064 ELEMENTS rmpoint/element/eltset get e_delete Dudded 26271 points plus parents, and 32064 elements removed. rmpoint/compress 26271 points removed and 0 elements removed. The new point count is 4323. The new element count is 4160. # 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 define/SNAME/SS_south_front 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 # report sideset number of faces eltset/SNAME/id_side/eq/SS_ID THE ELTSET SS_south_front HAS 8320 ELEMENTS eltset/e_keep/id_side/eq/SS_ID THE ELTSET e_keep HAS 8320 ELEMENTS eltset/e_delete/not/e_keep THE ELTSET e_delete HAS 27904 ELEMENTS rmpoint/element/eltset get e_delete Dudded 22079 points plus parents, and 27904 elements removed. rmpoint/compress 22079 points removed and 0 elements removed. The new point count is 8515. The new element count is 8320. # 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 define/SNAME/SS_bottom 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 # report sideset number of faces eltset/SNAME/id_side/eq/SS_ID THE ELTSET SS_bottom HAS 5632 ELEMENTS eltset/e_keep/id_side/eq/SS_ID THE ELTSET e_keep HAS 5632 ELEMENTS eltset/e_delete/not/e_keep THE ELTSET e_delete HAS 30592 ELEMENTS rmpoint/element/eltset get e_delete Dudded 27681 points plus parents, and 30592 elements removed. rmpoint/compress 27681 points removed and 0 elements removed. The new point count is 2913. The new element count is 5632. # 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 define/SNAME/SS_top 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 # report sideset number of faces eltset/SNAME/id_side/eq/SS_ID THE ELTSET SS_top HAS 5632 ELEMENTS eltset/e_keep/id_side/eq/SS_ID THE ELTSET e_keep HAS 5632 ELEMENTS eltset/e_delete/not/e_keep THE ELTSET e_delete HAS 30592 ELEMENTS rmpoint/element/eltset get e_delete Dudded 27681 points plus parents, and 30592 elements removed. rmpoint/compress 27681 points removed and 0 elements removed. The new point count is 2913. The new element count is 5632. # 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 # MUST DEFINE INPUT FILE FPOLY_PTS # output files and SS ID NUMBERS define SS_C_SIDE 7 define SS_NAME_7 SS_corner_side define FILE_C_SIDE fs_7_Csides.faceset define SS_C_TOP 8 define SS_NAME_8 SS_corner_top define FILE_C_TOP fs_8_Ctop.faceset # find and set top corner pts and facesets # use attribute layertyp = -2 from stack command # set top corner point faces (side) = SS_C_SIDE # set top corner point faces (top) = SS_C_TOP read avs FPOLY_PTS mo_cpts cmo/status/brief The current-mesh-object(CMO) is: mo_cpts 1 Mesh Object name: mo_pri number of nodes = 381603 number of elements = 732160 dimensions geometry = 3 element type = pri dimensions topology = 3 6 nodes 5 faces 9 edges boundary flag = 16000000 status = inactive 2 Mesh Object name: mo_surf number of nodes = 30594 number of elements = 36224 dimensions geometry = 3 element type = hyb dimensions topology = 2 10 nodes 10 faces 12 edges boundary flag = 16000000 status = inactive 3 Mesh Object name: mo_cpts number of nodes = 5 number of elements = 0 dimensions geometry = 3 element type = tet dimensions topology = 3 4 nodes 4 faces 6 edges boundary flag = 16000000 status = active finish # CHECK input poly points # zic should be 0, or set to 0 here cmo/printatt/mo_cpts/-xyz- minmax Warning: nelements = 0 ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH xic 9.219250000E+01 1.163908000E+02 2.419830000E+01 5 yic 2.731400000E+00 2.140300000E+01 1.867160000E+01 5 zic 0.000000000E+00 0.000000000E+00 0.000000000E+00 5 cmo/copy/mo_tmp/mo_surf cmo/select/mo_tmp cmo/addatt/mo_tmp/bc_faces/vint/scalar/nelements ########################## # Find the top points 0 distance from poly corner points cmo addatt mo_tmp zicsave VDOUBLE/scalar/nnodes cmo copyatt mo_tmp mo_tmp zicsave zic 30594 copied from mo_tmp zic to -> mo_tmp zicsave cmo setatt mo_tmp zic 0. 30594 values reset for attribute zic compute/distance_field/mo_tmp/mo_cpts/xdist cmo/addatt/mo_tmp/temp_tmp_var/VDOUBLE/scalar/nnodes/linear/temporary//-1.0 finish cmo/addatt/mo_tmp/xobj_tmp_var/VDOUBLE/scalar/nnodes/linear/temporary//-1.0 finish cmo/addatt/mo_tmp/yobj_tmp_var/VDOUBLE/scalar/nnodes/linear/temporary//-1.0 finish cmo/addatt/mo_tmp/zobj_tmp_var/VDOUBLE/scalar/nnodes/linear/temporary//-1.0 finish cmo/addatt/mo_tmp/xdist/VDOUBLE/scalar/nnodes/linear/permanent//-1.0 finish interpolate/voronoi/mo_tmp/xobj_tmp_var/1 0 0/mo_cpts/xic/keepatt nodes in indexed point set = 30594 cmo/addatt/mo_tmp/pt_gtg/VINT/scalar/nnodes//permanent/agfx/0 finish INTRP METHOD: voronoi FUNCTION: linear TIEBREAKER: maxtie FLAG: plus1 0.117390800000E+03 pt_gtg attribute will be added and kept. Sink point Points searched Points Found Percent Done 1531 1535 1531 5 % 3061 3071 3061 10 % 4591 4637 4591 15 % 6121 6212 6121 20 % 7651 7781 7651 25 % 9181 9373 9181 30 % 10711 10903 10711 35 % 12241 12433 12241 40 % 13771 14025 13771 45 % 15301 15617 15301 50 % 16831 17147 16831 55 % 18361 18739 18361 60 % 19891 20269 19891 65 % 21421 21861 21421 70 % 22951 23391 22951 75 % 24481 24964 24481 80 % 26011 26537 26011 85 % 27541 28110 27541 90 % 29071 29683 29071 95 % 30594 31249 30594 Total cmo select mo_tmp finish intrp/voronoi done. cmo kdtree release finish finish interpolate/voronoi/mo_tmp/yobj_tmp_var/1 0 0/mo_cpts/yic/keepatt nodes in indexed point set = 30594 INTRP METHOD: voronoi FUNCTION: linear TIEBREAKER: maxtie FLAG: plus1 0.224029999999E+02 SKIPPING POINT SEARCH... using lookup attribute pt_gtg Sink point Points searched Points Found Percent Done 1531 1531 1531 5 % 3061 3061 3061 10 % 4591 4591 4591 15 % 6121 6121 6121 20 % 7651 7651 7651 25 % 9181 9181 9181 30 % 10711 10711 10711 35 % 12241 12241 12241 40 % 13771 13771 13771 45 % 15301 15301 15301 50 % 16831 16831 16831 55 % 18361 18361 18361 60 % 19891 19891 19891 65 % 21421 21421 21421 70 % 22951 22951 22951 75 % 24481 24481 24481 80 % 26011 26011 26011 85 % 27541 27541 27541 90 % 29071 29071 29071 95 % 30594 30594 30594 Total cmo select mo_tmp finish intrp/voronoi done. cmo kdtree release finish finish interpolate/voronoi/mo_tmp/zobj_tmp_var/1 0 0/mo_cpts/zic nodes in indexed point set = 30594 INTRP METHOD: voronoi FUNCTION: linear TIEBREAKER: maxtie FLAG: plus1 0.100000000000E+01 SKIPPING POINT SEARCH... using lookup attribute pt_gtg Sink point Points searched Points Found Percent Done 1531 1531 1531 5 % 3061 3061 3061 10 % 4591 4591 4591 15 % 6121 6121 6121 20 % 7651 7651 7651 25 % 9181 9181 9181 30 % 10711 10711 10711 35 % 12241 12241 12241 40 % 13771 13771 13771 45 % 15301 15301 15301 50 % 16831 16831 16831 55 % 18361 18361 18361 60 % 19891 19891 19891 65 % 21421 21421 21421 70 % 22951 22951 22951 75 % 24481 24481 24481 80 % 26011 26011 26011 85 % 27541 27541 27541 90 % 29071 29071 29071 95 % 30594 30594 30594 Total cmo DELATT mo_tmp/pt_gtg finish cmo select mo_tmp finish intrp/voronoi done. cmo kdtree release finish finish math/sub/mo_tmp/temp_tmp_var/1,0,0/mo_tmp/xic/mo_tmp/xobj_tmp_var finish math/power/mo_tmp/xdist/1,0,0/mo_tmp/temp_tmp_var/2.0 finish math/sub/mo_tmp/temp_tmp_var/1,0,0/mo_tmp/yic/mo_tmp/yobj_tmp_var finish math/power/mo_tmp/temp_tmp_var/1,0,0/mo_tmp/temp_tmp_var/2.0 finish math/add/mo_tmp/xdist/1 0 0/mo_tmp/xdist/mo_tmp/temp_tmp_var finish math/sub/mo_tmp/temp_tmp_var/1,0,0/mo_tmp/zic/mo_tmp/zobj_tmp_var finish math/power/mo_tmp/temp_tmp_var/1,0,0/mo_tmp/temp_tmp_var/2.0 finish math/add/mo_tmp/xdist/1 0 0/mo_tmp/xdist/mo_tmp/temp_tmp_var finish math/power/mo_tmp/xdist/1,0,0/mo_tmp/xdist/0.5 finish cmo/delatt/mo_tmp/temp_tmp_var finish cmo/delatt/mo_tmp/xobj_tmp_var finish cmo/delatt/mo_tmp/yobj_tmp_var finish cmo/delatt/mo_tmp/zobj_tmp_var finish COMPUTE distance_field: Done. cmo copyatt/mo_tmp mo_tmp/zic zicsave 30594 copied from mo_tmp zicsave to -> mo_tmp zic dump gmv tmp.gmv mo_tmp cmo/modatt/-def-/-def-/ioflag/x finish pset/ptop attribute layertyp/1,0,0/eq -2 THE PSET ptop HAS 2913 POINTS pset p10 attribute xdist/pset,get,ptop lt .001 THE PSET p10 HAS 5 POINTS cmo setatt mo_tmp layertyp pset,get,p10 SS_C_SIDE 5 values reset for attribute layertyp eltset/e_tc/inclusive pset,get,p10 THE ELTSET e_tc HAS 20 ELEMENTS cmo setatt mo_tmp bc_faces eltset,get,e_tc SS_C_SIDE 20 values reset for attribute bc_faces eltset/e_top/id_side/eq 2 THE ELTSET e_top HAS 5632 ELEMENTS eltset e_tc_top inter e_top e_tc THE ELTSET e_tc_top HAS 10 ELEMENTS cmo setatt mo_tmp bc_faces eltset,get,e_tc_top SS_C_TOP 10 values reset for attribute bc_faces dump gmv tmp_bc_sides.gmv mo_tmp cmo/modatt/-def-/-def-/ioflag/x finish ########################## # Write the faceset faces and elements for id SS_C_SIDE cmo copy mo_tmp2 mo_tmp cmo select mo_tmp2 # report element count eltset/SS_NAME_7/bc_faces/eq/SS_C_SIDE THE ELTSET SS_corner_side HAS 10 ELEMENTS eltset/e_keep/bc_faces/eq/SS_C_SIDE THE ELTSET e_keep HAS 10 ELEMENTS eltset/e_delete/not/e_keep THE ELTSET e_delete HAS 36214 ELEMENTS rmpoint/element/eltset get e_delete Dudded 30564 points plus parents, and 36214 elements removed. rmpoint/compress 30564 points removed and 0 elements removed. The new point count is 30. The new element count is 10. # make sure to remove all attributes except idelem1 and idface1 cmo/DELATT/mo_tmp2/itetclr0 cmo/DELATT/mo_tmp2/idnode0 cmo/DELATT/mo_tmp2/idelem0 cmo/DELATT/mo_tmp2/facecol cmo/DELATT/mo_tmp2/itetclr1 cmo/DELATT/mo_tmp2/idface0 cmo/DELATT/mo_tmp2/id_side cmo/DELATT/mo_tmp2/layertyp cmo/DELATT/mo_tmp2/bc_faces dump/avs2/FILE_C_SIDE/mo_tmp2/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 ########################## # Write the faceset faces and elements for id SS_C_TOP cmo select mo_tmp # report element count eltset/SS_NAME_8/bc_faces/eq/SS_C_TOP THE ELTSET SS_corner_top HAS 10 ELEMENTS eltset/e_keep/bc_faces/eq/SS_C_TOP THE ELTSET e_keep HAS 10 ELEMENTS eltset/e_delete/not/e_keep THE ELTSET e_delete HAS 36214 ELEMENTS rmpoint/element/eltset get e_delete Dudded 30574 points plus parents, and 36214 elements removed. rmpoint/compress 30574 points removed and 0 elements removed. The new point count is 20. The new element count is 10. # 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 cmo/DELATT/mo_tmp/bc_faces dump/avs2/FILE_C_TOP/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 # clean up cmo/delete/mo_tmp Released Mesh Object: mo_tmp cmo/delete/mo_tmp2 Released Mesh Object: mo_tmp2 finish dump/exo/mesh_bc.exo/mo_pri///facesets fs_7_Csides.faceset fs_8_Ctop.faceset ExodusII: Start writing to file: mesh_bc.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: 381603 number of elements: 732160 number of edges: 0 number of edge blocks: 0 number of element blocks: 7 number of face blocks: 0 number of node sets: 0 number of edge sets: 3 number of element sets: 0 number of side sets: 2 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 ExodusII file: mesh_bc.exo using cmo: mo_pri ############################ # 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 30594 scalar 1 1 0 1 vector 3 3 0 1 nnodes 30594 30594 0 1 nedges 0 0 0 1 nfaces 0 0 0 1 nelements 36224 36224 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 30594 ialias 0 0 0 30594 imt1 1 1 0 30594 itp1 0 2 2 30594 icr1 0 0 0 30594 isn1 0 0 0 30594 xic 9.219250000E+01 1.163908000E+02 2.419830000E+01 30594 yic 2.731400000E+00 2.140300000E+01 1.867160000E+01 30594 zic -4.500000000E+01 5.504421902E+00 5.050442190E+01 30594 xtetwd 0 0 0 36224 itetclr 1 6 5 36224 itettyp 3 4 1 36224 itetoff 0 133628 133628 36224 jtetoff 0 133628 133628 36224 itet 0 30594 30594 36224x* jtet 0 362234 362234 36224x* 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 30594 30594 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 36224 itetclr1 1 30 29 36224 idnode0 1 381603 381602 30594 idelem0 0 0 0 36224 idelem1 1 732160 732159 36224 facecol 1 6 5 36224 idface0 0 0 0 36224 idface1 1 5 4 36224 layertyp -2 2 4 30594 node_ice 0 1 1 30594 ncon50 0 0 0 1 nconbnd 0 0 0 1 0 length attribute: icontab id_side 1 6 5 36224 cmo printatt mo_pri -all- minmax ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH -def- 0.000000000E+00 0.000000000E+00 0.000000000E+00 381603 scalar 1 1 0 1 vector 3 3 0 1 nnodes 381603 381603 0 1 nedges 0 0 0 1 nfaces 0 0 0 1 nelements 732160 732160 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 381603 ialias 0 0 0 381603 imt1 1 6 5 381603 itp1 0 10 10 381603 icr1 0 0 0 381603 isn1 0 0 0 381603 xic 9.219250000E+01 1.163908000E+02 2.419830000E+01 381603 yic 2.731400000E+00 2.140300000E+01 1.867160000E+01 381603 zic -4.500000000E+01 5.504421902E+00 5.050442190E+01 381603 xtetwd 0 0 0 732160 itetclr 1 30 29 732160 itettyp 7 7 0 732160 itetoff 0 4392954 4392954 732160 jtetoff 0 3660795 3660795 732160 itet 1 381603 381602 732160x6 jtet 3 19660797 19660794 732160x5 epsilon 1.000000004E-15 1.000000004E-15 0.000000000E+00 1 epsilonl 1.310793643E-11 1.310793643E-11 0.000000000E+00 1 epsilona 7.737994692E-10 7.737994692E-10 0.000000000E+00 1 epsilonv 5.066826361E-09 5.066826361E-09 0.000000000E+00 1 ipointi 1 1 0 1 ipointj 381603 381603 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 9.219250000E+01 9.219250000E+01 0.000000000E+00 1 ymin 2.731400000E+00 2.731400000E+00 0.000000000E+00 1 zmin -4.500000000E+01 -4.500000000E+01 0.000000000E+00 1 xmax 1.163908000E+02 1.163908000E+02 0.000000000E+00 1 ymax 2.140300000E+01 2.140300000E+01 0.000000000E+00 1 zmax 5.504421902E+00 5.504421902E+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 381603 node_ice 0 1 1 381603 mat_orig 1 6 5 732160 ikey_utr 1 732160 732159 732160 cmo printatt mo_pri -xyz- minmax ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH xic 9.219250000E+01 1.163908000E+02 2.419830000E+01 381603 yic 2.731400000E+00 2.140300000E+01 1.867160000E+01 381603 zic -4.500000000E+01 5.504421902E+00 5.050442190E+01 381603 cmo select mo_pri quality volume itetclr QUALITY TOTALS for 7 MATERIALS: ........................................................... eltset/-etmp-/itetclr/eq/1/ THE ELTSET -etmp- HAS 5632 ELEMENTS finish epsilonvol: 5.0668264E-09 --------------------------------------- VOLUMES for ELEMENT COLOR 1 ........................... element volumes b/w 0.6213E-03 and 0.7898E-03: 17 element volumes b/w 0.7898E-03 and 0.1004E-02: 409 element volumes b/w 0.1004E-02 and 0.1276E-02: 3384 element volumes b/w 0.1276E-02 and 0.1622E-02: 1774 element volumes b/w 0.1622E-02 and 0.2062E-02: 48 min volume = 6.2127027E-04 max volume = 2.0622457E-03 eltset/-etmp-/delete/ ELTSET -etmp- DELETED finish eltset/-etmp-/itetclr/eq/2/ THE ELTSET -etmp- HAS 50688 ELEMENTS finish --------------------------------------- VOLUMES for ELEMENT COLOR 2 ........................... element volumes b/w 0.6213E-03 and 0.7898E-03: 153 element volumes b/w 0.7898E-03 and 0.1004E-02: 3681 element volumes b/w 0.1004E-02 and 0.1276E-02: 30456 element volumes b/w 0.1276E-02 and 0.1622E-02: 15966 element volumes b/w 0.1622E-02 and 0.2062E-02: 432 min volume = 6.2127021E-04 max volume = 2.0622473E-03 eltset/-etmp-/delete/ ELTSET -etmp- DELETED finish eltset/-etmp-/itetclr/eq/3/ THE ELTSET -etmp- HAS 50688 ELEMENTS finish --------------------------------------- VOLUMES for ELEMENT COLOR 3 ........................... element volumes b/w 0.6213E-03 and 0.7898E-03: 153 element volumes b/w 0.7898E-03 and 0.1004E-02: 3681 element volumes b/w 0.1004E-02 and 0.1276E-02: 30456 element volumes b/w 0.1276E-02 and 0.1622E-02: 15966 element volumes b/w 0.1622E-02 and 0.2062E-02: 432 min volume = 6.2127021E-04 max volume = 2.0622473E-03 eltset/-etmp-/delete/ ELTSET -etmp- DELETED finish eltset/-etmp-/itetclr/eq/4/ THE ELTSET -etmp- HAS 33792 ELEMENTS finish --------------------------------------- VOLUMES for ELEMENT COLOR 4 ........................... element volumes b/w 0.6213E-03 and 0.7898E-03: 102 element volumes b/w 0.7898E-03 and 0.1004E-02: 2454 element volumes b/w 0.1004E-02 and 0.1276E-02: 20304 element volumes b/w 0.1276E-02 and 0.1622E-02: 10644 element volumes b/w 0.1622E-02 and 0.2062E-02: 288 min volume = 6.2127003E-04 max volume = 2.0622461E-03 eltset/-etmp-/delete/ ELTSET -etmp- DELETED finish eltset/-etmp-/itetclr/eq/5/ THE ELTSET -etmp- HAS 80218 ELEMENTS finish --------------------------------------- VOLUMES for ELEMENT COLOR 5 ........................... element volumes b/w 0.7766E-03 and 0.9872E-03: 340 element volumes b/w 0.9872E-03 and 0.1255E-02: 4669 element volumes b/w 0.1255E-02 and 0.1595E-02: 48225 element volumes b/w 0.1595E-02 and 0.2028E-02: 26484 element volumes b/w 0.2028E-02 and 0.2578E-02: 500 min volume = 7.7658784E-04 max volume = 2.5778071E-03 eltset/-etmp-/delete/ ELTSET -etmp- DELETED finish eltset/-etmp-/itetclr/eq/6/ THE ELTSET -etmp- HAS 458155 ELEMENTS finish --------------------------------------- VOLUMES for ELEMENT COLOR 6 ........................... element volumes b/w 0.3106E-02 and 0.7211E-02: 239305 element volumes b/w 0.7211E-02 and 0.1674E-01: 22040 element volumes b/w 0.1674E-01 and 0.3885E-01: 56130 element volumes b/w 0.3885E-01 and 0.9018E-01: 57290 element volumes b/w 0.9018E-01 and 0.2093E+00: 83390 min volume = 3.1063498E-03 max volume = 2.0933958E-01 eltset/-etmp-/delete/ ELTSET -etmp- DELETED finish eltset/-etmp-/itetclr/eq/30/ THE ELTSET -etmp- HAS 52987 ELEMENTS finish --------------------------------------- VOLUMES for ELEMENT COLOR 30 ........................... element volumes b/w 0.1014E-02 and 0.1583E-02: 22553 element volumes b/w 0.1583E-02 and 0.2472E-02: 9869 element volumes b/w 0.2472E-02 and 0.3860E-02: 0 element volumes b/w 0.3860E-02 and 0.6026E-02: 11945 element volumes b/w 0.6026E-02 and 0.9409E-02: 8620 min volume = 1.0141499E-03 max volume = 9.4088156E-03 eltset/-etmp-/delete/ ELTSET -etmp- DELETED finish --------------------------------------- SUMMARY VOLUME/AREA of ELEMENT COLORS for mo_pri 1 7 ........................................................... Color Num. Elements Volume Fractional Volume 1 5632 0.6848947E+01 0.000397497 2 50688 0.6164053E+02 0.003577470 3 50688 0.6164053E+02 0.003577470 4 33792 0.4109368E+02 0.002384980 5 80218 0.1223984E+03 0.007103710 6 458155 0.1676359E+05 0.972918904 30 52987 0.1729907E+03 0.010039970 Total elements: 732160 Total Volume: 0.1723020E+05 ----------------------------------------------------------- 732160 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_5_lef.faceset fs_6_fro.faceset fs_7_Csides.face set fs_8_Ctop.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/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: 381603 number of elements: 732160 number of edges: 0 number of edge blocks: 0 number of element blocks: 7 number of face blocks: 0 number of node sets: 0 number of edge sets: 3 number of element sets: 0 number of side sets: 8 number of face sets: 0 number of node maps: 0 number of edge maps: 0 number of face maps: 0 number of element maps: 0 cmo/modatt/mo_pri/itetclr/ioflag/l finish ExodusII: Done writing to ExodusII file: mesh_fsets.exo using cmo: mo_pri 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 # DONE with multi-material with ice finish LaGriT successfully completed