work folder /home/lucia/sandbox/barrow_site3/ input files /home/lucia/sandbox/barrow_site3/input/ documents /home/lucia/sandbox/barrow_site3/docs/ facesets /home/lucia/sandbox/barrow_site3/face/
ID LAGRIT MESH (EXO eb_prop1) Material Description 30 material 30 (300000) ice 5 material 5 (50000) moss layer (active layer 4) 4 material 4 (40000) upper organic layer (active layer 3) 3 material 3 (30000) lower organic layer (active layer 2) 2 material 2 (20000) active layer 1 (could be organic or mineral) 1 material 1 (10000) frozen mineral soil ID Name Cell Thick Num Cells Mat Thick 5 moss active layer 4 .2 1 0.20 4 upper organic layer 3 .3 1 0.30 3 lower organic layer 2 .1 1 0.10 2 layer 1 (organic or mineral) .05 1 0.05 1 frozen mineral soil .34545 2 - 11 0.69 - 3.799 30 ice .34545 0 - 9 0.00 - 3.730
* * * * * * * * * * * * * * * * * * * * * * * * *
* * *
* * Program: LaGriT V3.101 Linux m64 *
* * date_compile: 2013/04/25 RH gf4.5 *
* * Run Time: 2013/May 7 08:34:54 *
* * 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/proto_full_poly.inp/mo_pri
read/avs/proto_full_poly.inp/mo_pri
cmo/addatt/mo_pri/layertyp/VDOUBLE/scalar/nnodes/linear/permanent/gxaf/0.0
finish
cmo/addatt/mo_pri/elemclr/VDOUBLE/scalar/nnodes/linear/permanent/gxaf/0.0
finish
cmo/addatt/mo_pri/itetclr_/VDOUBLE/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 = 51344 number of elements = 94080
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 = 51344 number of elements = 94080
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 5.134E+04
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 9.408E+04
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 8.220E-12
35 epsilona REAL scalar scalar const perma l 3.043E-10
36 epsilonv REAL scalar scalar const perma l 7.375E-10
37 ipointi INT scalar scalar const perma l 1.000E+00
38 ipointj INT scalar scalar const perma l 5.134E+04
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 4.635E+01
52 ymin REAL scalar scalar const perma l 1.062E+02
53 zmin REAL scalar scalar const perma l 5.430E-02
54 xmax REAL scalar scalar const perma l 7.184E+01
55 ymax REAL scalar scalar const perma l 1.326E+02
56 zmax REAL scalar scalar const perma l 4.992E+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
67 elemclr VDOU scalar nnodes linea perma gxa 0.000E+00
68 itetclr_ VDOU 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/VDOUBLE/scalar/nnodes/linear/permanent/gxaf
finish
cmo/addatt/mo_surf/elemclr/VDOUBLE/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 8434
scalar 1 1 0 1
vector 3 3 0 1
nnodes 8434 8434 0 1
nedges 0 0 0 1
nfaces 0 0 0 1
nelements 14704 14704 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 8434
ialias 0 0 0 8434
imt1 1 1 0 8434
itp1 0 2 2 8434
icr1 0 0 0 8434
isn1 0 0 0 8434
xic 4.634709930E+01 7.183999634E+01 2.549289703E+01 8434
yic 1.062200012E+02 1.326065216E+02 2.638652039E+01 8434
zic 5.430011824E-02 4.992260933E+00 4.937960815E+00 8434
xtetwd 0 0 0 14704
itetclr 1 1 0 14704
itettyp 3 4 1 14704
itetoff 0 46268 46268 14704
jtetoff 0 46268 46268 14704
itet 0 8434 8434 14704x*
jtet 0 147034 147034 14704x*
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 8434 8434 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 14704
itetclr1 1 30 29 14704
idnode0 1 51344 51343 8434
idelem0 0 0 0 14704
idelem1 1 94080 94079 14704
facecol 1 6 5 14704
idface0 0 0 0 14704
idface1 1 5 4 14704
layertyp -2.000000000E+00 2.000000000E+00 4.000000000E+00 8434
elemclr 3.000000000E+01 3.100000000E+01 1.000000000E+00 8434
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
8432 out of 14704 element colors changed.
cmo/copyatt/mo_surf mo_surf/id_side itetclr
14704 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 14704
# 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 960 ELEMENTS
eltset/e_delete/not/e_keep
THE ELTSET e_delete HAS 13744 ELEMENTS
rmpoint/element/eltset get e_delete
Dudded 7394 points plus parents, and 13744 elements removed.
rmpoint/compress
7394 points removed and 0 elements removed.
The new point count is 1040.
The new element count is 960.
# 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 480 ELEMENTS
eltset/e_delete/not/e_keep
THE ELTSET e_delete HAS 14224 ELEMENTS
rmpoint/element/eltset get e_delete
Dudded 7906 points plus parents, and 14224 elements removed.
rmpoint/compress
7906 points removed and 0 elements removed.
The new point count is 528.
The new element count is 480.
# 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 240 ELEMENTS
eltset/e_delete/not/e_keep
THE ELTSET e_delete HAS 14464 ELEMENTS
rmpoint/element/eltset get e_delete
Dudded 8162 points plus parents, and 14464 elements removed.
rmpoint/compress
8162 points removed and 0 elements removed.
The new point count is 272.
The new element count is 240.
# 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 480 ELEMENTS
eltset/e_delete/not/e_keep
THE ELTSET e_delete HAS 14224 ELEMENTS
rmpoint/element/eltset get e_delete
Dudded 7906 points plus parents, and 14224 elements removed.
rmpoint/compress
7906 points removed and 0 elements removed.
The new point count is 528.
The new element count is 480.
# 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 6272 ELEMENTS
eltset/e_delete/not/e_keep
THE ELTSET e_delete HAS 8432 ELEMENTS
rmpoint/element/eltset get e_delete
Dudded 5225 points plus parents, and 8432 elements removed.
rmpoint/compress
5225 points removed and 0 elements removed.
The new point count is 3209.
The new element count is 6272.
# 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 6272 ELEMENTS
eltset/e_delete/not/e_keep
THE ELTSET e_delete HAS 8432 ELEMENTS
rmpoint/element/eltset get e_delete
Dudded 5225 points plus parents, and 8432 elements removed.
rmpoint/compress
5225 points removed and 0 elements removed.
The new point count is 3209.
The new element count is 6272.
# 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
############################
# 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 8434
scalar 1 1 0 1
vector 3 3 0 1
nnodes 8434 8434 0 1
nedges 0 0 0 1
nfaces 0 0 0 1
nelements 14704 14704 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 8434
ialias 0 0 0 8434
imt1 1 1 0 8434
itp1 0 2 2 8434
icr1 0 0 0 8434
isn1 0 0 0 8434
xic 4.634709930E+01 7.183999634E+01 2.549289703E+01 8434
yic 1.062200012E+02 1.326065216E+02 2.638652039E+01 8434
zic 5.430011824E-02 4.992260933E+00 4.937960815E+00 8434
xtetwd 0 0 0 14704
itetclr 1 6 5 14704
itettyp 3 4 1 14704
itetoff 0 46268 46268 14704
jtetoff 0 46268 46268 14704
itet 0 8434 8434 14704x*
jtet 0 147034 147034 14704x*
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 8434 8434 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 14704
itetclr1 1 30 29 14704
idnode0 1 51344 51343 8434
idelem0 0 0 0 14704
idelem1 1 94080 94079 14704
facecol 1 6 5 14704
idface0 0 0 0 14704
idface1 1 5 4 14704
layertyp -2.000000000E+00 2.000000000E+00 4.000000000E+00 8434
elemclr 3.000000000E+01 3.100000000E+01 1.000000000E+00 8434
ncon50 0 0 0 1
nconbnd 0 0 0 1
0 length attribute: icontab
id_side 1 6 5 14704
cmo printatt mo_pri -all- minmax
ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH
-def- 0.000000000E+00 0.000000000E+00 0.000000000E+00 51344
scalar 1 1 0 1
vector 3 3 0 1
nnodes 51344 51344 0 1
nedges 0 0 0 1
nfaces 0 0 0 1
nelements 94080 94080 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 51344
ialias 0 0 0 51344
imt1 1 5 4 51344
itp1 0 10 10 51344
icr1 0 0 0 51344
isn1 0 0 0 51344
xic 4.634709930E+01 7.183999634E+01 2.549289703E+01 51344
yic 1.062200012E+02 1.326065216E+02 2.638652039E+01 51344
zic 5.430011824E-02 4.992260933E+00 4.937960815E+00 51344
xtetwd 0 0 0 94080
itetclr 1 30 29 94080
itettyp 7 7 0 94080
itetoff 0 564474 564474 94080
jtetoff 0 470395 470395 94080
itet 1 51344 51343 94080x6
jtet 2 16470397 16470395 94080x5
epsilon 1.000000004E-15 1.000000004E-15 0.000000000E+00 1
epsilonl 8.220209176E-12 8.220209176E-12 0.000000000E+00 1
epsilona 3.043165085E-10 3.043165085E-10 0.000000000E+00 1
epsilonv 7.375461151E-10 7.375461151E-10 0.000000000E+00 1
ipointi 1 1 0 1
ipointj 51344 51344 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 4.634709930E+01 4.634709930E+01 0.000000000E+00 1
ymin 1.062200012E+02 1.062200012E+02 0.000000000E+00 1
zmin 5.430011824E-02 5.430011824E-02 0.000000000E+00 1
xmax 7.183999634E+01 7.183999634E+01 0.000000000E+00 1
ymax 1.326065216E+02 1.326065216E+02 0.000000000E+00 1
zmax 4.992260933E+00 4.992260933E+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 51344
elemclr 3.000000000E+01 3.100000000E+01 1.000000000E+00 51344
itetclr_ 1.000000000E+00 5.000000000E+00 4.000000000E+00 94080
ikey_utr 1 94080 94079 94080
cmo printatt mo_pri -xyz- minmax
ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH
xic 4.634709930E+01 7.183999634E+01 2.549289703E+01 51344
yic 1.062200012E+02 1.326065216E+02 2.638652039E+01 51344
zic 5.430011824E-02 4.992260933E+00 4.937960815E+00 51344
cmo select mo_pri
quality volume itetclr
QUALITY TOTALS for 6 MATERIALS:
...........................................................
eltset/-etmp-/itetclr/eq/1/
THE ELTSET -etmp- HAS 60864 ELEMENTS
finish
epsilonvol: 7.3754612E-10
---------------------------------------
VOLUMES for ELEMENT COLOR 1
...........................
element volumes b/w 0.7890E-02 and 0.1085E-01: 271
element volumes b/w 0.1085E-01 and 0.1492E-01: 1395
element volumes b/w 0.1492E-01 and 0.2052E-01: 18833
element volumes b/w 0.2052E-01 and 0.2822E-01: 39444
element volumes b/w 0.2822E-01 and 0.3880E-01: 921
min volume = 7.8900892E-03 max volume = 3.8803718E-02
eltset/-etmp-/delete/
ELTSET -etmp- DELETED
finish
eltset/-etmp-/itetclr/eq/2/
THE ELTSET -etmp- HAS 6272 ELEMENTS
finish
---------------------------------------
VOLUMES for ELEMENT COLOR 2
...........................
element volumes b/w 0.1142E-02 and 0.1570E-02: 25
element volumes b/w 0.1570E-02 and 0.2160E-02: 130
element volumes b/w 0.2160E-02 and 0.2970E-02: 1833
element volumes b/w 0.2970E-02 and 0.4084E-02: 4079
element volumes b/w 0.4084E-02 and 0.5616E-02: 205
min volume = 1.1419920E-03 max volume = 5.6163476E-03
eltset/-etmp-/delete/
ELTSET -etmp- DELETED
finish
eltset/-etmp-/itetclr/eq/3/
THE ELTSET -etmp- HAS 6272 ELEMENTS
finish
---------------------------------------
VOLUMES for ELEMENT COLOR 3
...........................
element volumes b/w 0.2284E-02 and 0.3141E-02: 25
element volumes b/w 0.3141E-02 and 0.4319E-02: 130
element volumes b/w 0.4319E-02 and 0.5940E-02: 1833
element volumes b/w 0.5940E-02 and 0.8168E-02: 4079
element volumes b/w 0.8168E-02 and 0.1123E-01: 205
min volume = 2.2839732E-03 max volume = 1.1232642E-02
eltset/-etmp-/delete/
ELTSET -etmp- DELETED
finish
eltset/-etmp-/itetclr/eq/4/
THE ELTSET -etmp- HAS 6272 ELEMENTS
finish
---------------------------------------
VOLUMES for ELEMENT COLOR 4
...........................
element volumes b/w 0.6852E-02 and 0.9423E-02: 25
element volumes b/w 0.9423E-02 and 0.1296E-01: 130
element volumes b/w 0.1296E-01 and 0.1782E-01: 1833
element volumes b/w 0.1782E-01 and 0.2450E-01: 4079
element volumes b/w 0.2450E-01 and 0.3370E-01: 205
min volume = 6.8519305E-03 max volume = 3.3697978E-02
eltset/-etmp-/delete/
ELTSET -etmp- DELETED
finish
eltset/-etmp-/itetclr/eq/5/
THE ELTSET -etmp- HAS 6272 ELEMENTS
finish
---------------------------------------
VOLUMES for ELEMENT COLOR 5
...........................
element volumes b/w 0.4568E-02 and 0.6282E-02: 25
element volumes b/w 0.6282E-02 and 0.8638E-02: 130
element volumes b/w 0.8638E-02 and 0.1188E-01: 1833
element volumes b/w 0.1188E-01 and 0.1634E-01: 4079
element volumes b/w 0.1634E-01 and 0.2247E-01: 205
min volume = 4.5679464E-03 max volume = 2.2465283E-02
eltset/-etmp-/delete/
ELTSET -etmp- DELETED
finish
eltset/-etmp-/itetclr/eq/30/
THE ELTSET -etmp- HAS 8128 ELEMENTS
finish
---------------------------------------
VOLUMES for ELEMENT COLOR 30
...........................
element volumes b/w 0.9353E-02 and 0.1243E-01: 7
element volumes b/w 0.1243E-01 and 0.1652E-01: 207
element volumes b/w 0.1652E-01 and 0.2196E-01: 1704
element volumes b/w 0.2196E-01 and 0.2919E-01: 5265
element volumes b/w 0.2919E-01 and 0.3880E-01: 945
min volume = 9.3529829E-03 max volume = 3.8803751E-02
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 60864 0.1308482E+04 0.740990154
2 6272 0.1984116E+02 0.011235998
3 6272 0.3968212E+02 0.022471889
4 6272 0.1190466E+03 0.067415773
5 6272 0.7936425E+02 0.044943777
30 8128 0.1994401E+03 0.112942409
Total elements: 94080 Total Volume: 0.1765856E+04
-----------------------------------------------------------
94080 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
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: 51344
number of elements: 94080
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: 6
number of face sets: 0
number of node maps: 0
number of edge maps: 0
number of face maps: 0
number of element maps: 0
cmo/modatt/mo_pri/itetclr/ioflag/l
finish
ExodusII: Done writing to 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
############################
# version with single material and top faceset only
cmo setatt mo_pri itetclr 1
94080 values reset for attribute itetclr
cmo setatt mo_pri imt 1
51344 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_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: 51344
number of elements: 94080
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: 6
number of face sets: 0
number of node maps: 0
number of edge maps: 0
number of face maps: 0
number of element maps: 0
cmo/modatt/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 51344
scalar 1 1 0 1
vector 3 3 0 1
nnodes 51344 51344 0 1
nedges 0 0 0 1
nfaces 0 0 0 1
nelements 94080 94080 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 51344
ialias 0 0 0 51344
imt1 1 1 0 51344
itp1 0 10 10 51344
icr1 0 0 0 51344
isn1 0 0 0 51344
xic 4.634709930E+01 7.183999634E+01 2.549289703E+01 51344
yic 1.062200012E+02 1.326065216E+02 2.638652039E+01 51344
zic 5.430011824E-02 4.992260933E+00 4.937960815E+00 51344
xtetwd 0 0 0 94080
itetclr 1 1 0 94080
itettyp 7 7 0 94080
itetoff 0 564474 564474 94080
jtetoff 0 470395 470395 94080
itet 1 51344 51343 94080x6
jtet 2 16000000 15999998 94080x5
epsilon 1.000000004E-15 1.000000004E-15 0.000000000E+00 1
epsilonl 8.220209176E-12 8.220209176E-12 0.000000000E+00 1
epsilona 3.043165085E-10 3.043165085E-10 0.000000000E+00 1
epsilonv 7.375461151E-10 7.375461151E-10 0.000000000E+00 1
ipointi 1 1 0 1
ipointj 51344 51344 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 4.634709930E+01 4.634709930E+01 0.000000000E+00 1
ymin 1.062200012E+02 1.062200012E+02 0.000000000E+00 1
zmin 5.430011824E-02 5.430011824E-02 0.000000000E+00 1
xmax 7.183999634E+01 7.183999634E+01 0.000000000E+00 1
ymax 1.326065216E+02 1.326065216E+02 0.000000000E+00 1
zmax 4.992260933E+00 4.992260933E+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 51344
elemclr 3.000000000E+01 3.100000000E+01 1.000000000E+00 51344
itetclr_ 1.000000000E+00 5.000000000E+00 4.000000000E+00 94080
ikey_utr 1 94080 94079 94080
finish
LaGriT successfully completed