* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * Program: LaGriT V3.101 Linux m64 * * * date_compile: 2013/04/25 RH gf4.5 * * * Run Time: 2013/Dec 4 13:31:15 * * * 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 # subset to smaller bounding box shape area # # FULL MESH # xic 1.060500000E+05 4.510500000E+05 3.450000000E+05 1102042 # yic 3.869700000E+06 4.442200000E+06 5.725000000E+05 1102042 # zic -4.021226013E+03 1.698531097E+03 5.719757111E+03 1102042 # # NEW MESH Dec 2013 # xic 4.201000000E+05 4.299000000E+05 9.800000000E+03 117450 # yic 4.307100000E+06 4.316900000E+06 9.800000000E+03 117450 # zic -2.208199951E+03 1.081999969E+02 2.316399948E+03 117450 # # xval 1.000000000E+00 9.000000000E+00 8.000000000E+00 117450 # nlayer 1 47 46 117450 # # Note this mesh is numbered from top to bottom # # # 50 km X 50 km domain. Four corners for this domain: # # SW x240000 y4043000 mins # SE x290000 y4043000 # NE x290000 y4093000 maxs # NW x240000 y4093000 # define MO mo # Define for current data file define XMIN 420100. define YMIN 4307100. define ZMIN -2210. define ZMIN -5000. define XMAX 430100. define YMAX 4317100. define ZMAX 110. define ZMAX 5000. define Z_BOT -4000. define Z_BOT -2400. define Z_TOP 1700. define Z_TOP 300. # flat and top layers have next numbers after max xval define Z_BOT_ID 10 define Z_BOT_VAL 10. define Z_TOP_ID 11 define Z_TOP_VAL 11. # tag top and bottom flat points define TOP_LAYER 50 define BOT_LAYER 0 read avs points.inp mo cmo/addatt/mo/xval/VDOUBLE/scalar/nnodes/linear/permanent/gxaf/0.0 finish cmo/addatt/mo/nlayer/VINT/scalar/nnodes/linear/permanent/gxaf/0.0 finish cmo/status/brief The current-mesh-object(CMO) is: mo 1 Mesh Object name: mo number of nodes = 117450 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 cmo/printatt/mo/-xyz-/minmax Warning: nelements = 0 ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH xic 4.201000000E+05 4.299000000E+05 9.800000000E+03 117450 yic 4.307100000E+06 4.316900000E+06 9.800000000E+03 117450 zic -2.208199951E+03 1.081999969E+02 2.316399948E+03 117450 cmo addatt mo imat VINT/scalar/nnodes cmo copyatt/mo mo/imat xval Warning: nelements = 0 in cmo mo Warning: nelements = 0 in cmo mo Warning: Attribute type VDOUBLE written to VINT 117450 copied from mo xval to -> mo imat cmo printatt mo -all- minmax Warning: nelements = 0 ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH -def- 0.000000000E+00 0.000000000E+00 0.000000000E+00 117450 scalar 1 1 0 1 vector 3 3 0 1 nnodes 117450 117450 0 1 nedges 0 0 0 1 nfaces 0 0 0 1 nelements 0 0 0 1 mbndry 16000000 16000000 0 1 ndimensions_topo 3 3 0 1 ndimensions_geom 3 3 0 1 nodes_per_element 4 4 0 1 edges_per_element 6 6 0 1 faces_per_element 4 4 0 1 isetwd 0 0 0 117450 ialias 0 0 0 117450 imt1 1 47 46 117450 itp1 0 0 0 117450 icr1 0 0 0 117450 isn1 0 0 0 117450 xic 4.201000000E+05 4.299000000E+05 9.800000000E+03 117450 yic 4.307100000E+06 4.316900000E+06 9.800000000E+03 117450 zic -2.208199951E+03 1.081999969E+02 2.316399948E+03 117450 0 length attribute: xtetwd 0 length attribute: itetclr 0 length attribute: itettyp 0 length attribute: itetoff 0 length attribute: jtetoff 0 length attribute: itet 0 length attribute: jtet epsilon 1.000000004E-15 1.000000004E-15 0.000000000E+00 1 epsilonl 3.120068113E-09 3.120068113E-09 0.000000000E+00 1 epsilona 4.384175439E-05 4.384175439E-05 0.000000000E+00 1 epsilonv 4.939760845E-02 4.939760845E-02 0.000000000E+00 1 ipointi 1 1 0 1 ipointj 117450 117450 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.201000000E+05 4.201000000E+05 0.000000000E+00 1 ymin 4.307100000E+06 4.307100000E+06 0.000000000E+00 1 zmin -2.208199951E+03 -2.208199951E+03 0.000000000E+00 1 xmax 4.299000000E+05 4.299000000E+05 0.000000000E+00 1 ymax 4.316900000E+06 4.316900000E+06 0.000000000E+00 1 zmax 1.081999969E+02 1.081999969E+02 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 xval 1.000000000E+00 9.000000000E+00 8.000000000E+00 117450 nlayer 1 47 46 117450 imat 1 9 8 117450 ######### # SUBSET ######### pset/pbox/geom/xyz/1,0,0/XMIN YMIN ZMIN/XMAX YMAX ZMAX THE PSET pbox HAS 117450 POINTS pset/pdel/not pbox THE PSET pdel HAS 0 POINTS rmpoint/pset,get,pdel Dudded 0 points plus parents, there are no elements in mesh. rmpoint/compress 0 points removed and 0 elements removed. The new point count is 117450. The new element count is 0. cmo/printatt/mo/-xyz-/minmax Warning: nelements = 0 ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH xic 4.201000000E+05 4.299000000E+05 9.800000000E+03 117450 yic 4.307100000E+06 4.316900000E+06 9.800000000E+03 117450 zic -2.208199951E+03 1.081999969E+02 2.316399948E+03 117450 dump gmv subset.gmv mo cmo/modatt/-def-/-def-/ioflag/x finish cmo printatt mo -all- minmax Warning: nelements = 0 ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH -def- 0.000000000E+00 0.000000000E+00 0.000000000E+00 117450 scalar 1 1 0 1 vector 3 3 0 1 nnodes 117450 117450 0 1 nedges 0 0 0 1 nfaces 0 0 0 1 nelements 0 0 0 1 mbndry 16000000 16000000 0 1 ndimensions_topo 3 3 0 1 ndimensions_geom 3 3 0 1 nodes_per_element 4 4 0 1 edges_per_element 6 6 0 1 faces_per_element 4 4 0 1 isetwd -3 -3 0 117450 ialias 0 0 0 117450 imt1 1 47 46 117450 itp1 0 0 0 117450 icr1 0 0 0 117450 isn1 0 0 0 117450 xic 4.201000000E+05 4.299000000E+05 9.800000000E+03 117450 yic 4.307100000E+06 4.316900000E+06 9.800000000E+03 117450 zic -2.208199951E+03 1.081999969E+02 2.316399948E+03 117450 0 length attribute: xtetwd 0 length attribute: itetclr 0 length attribute: itettyp 0 length attribute: itetoff 0 length attribute: jtetoff 0 length attribute: itet 0 length attribute: jtet epsilon 1.000000004E-15 1.000000004E-15 0.000000000E+00 1 epsilonl 3.120068113E-09 3.120068113E-09 0.000000000E+00 1 epsilona 4.384175439E-05 4.384175439E-05 0.000000000E+00 1 epsilonv 4.939760845E-02 4.939760845E-02 0.000000000E+00 1 ipointi 1 1 0 1 ipointj 117450 117450 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.201000000E+05 4.201000000E+05 0.000000000E+00 1 ymin 4.307100000E+06 4.307100000E+06 0.000000000E+00 1 zmin -2.208199951E+03 -2.208199951E+03 0.000000000E+00 1 xmax 4.299000000E+05 4.299000000E+05 0.000000000E+00 1 ymax 4.316900000E+06 4.316900000E+06 0.000000000E+00 1 zmax 1.081999969E+02 1.081999969E+02 0.000000000E+00 1 kdtree_level 0 0 0 1 max_number_sets 64 64 0 1 number_of_psets 2 2 0 1 number_of_eltsets 0 0 0 1 number_of_fsets 0 0 0 1 xval 1.000000000E+00 9.000000000E+00 8.000000000E+00 117450 nlayer 1 47 46 117450 imat 1 9 8 117450 ######### # ADD FLAT BOTTOM and FLAT TOP # this helps connect and makes it easier to # remove "outside" elements ######### cmo copy motmp mo cmo setatt motmp zic Z_BOT Warning: nelements = 0 117450 values reset for attribute zic cmo select motmp filter 1,0,0 FILTER:Use internal epsilonl value = 0.312006811336E-08 FILTER:Dudding 114950 duplicate points. rmpoint compress 114950 points removed and 0 elements removed. The new point count is 2500. The new element count is 0. cmo setatt motmp imat Z_BOT_ID Warning: nelements = 0 2500 values reset for attribute imat cmo setatt motmp xval Z_BOT_VAL Warning: nelements = 0 2500 values reset for attribute xval cmo setatt motmp nlayer BOT_LAYER Warning: nelements = 0 2500 values reset for attribute nlayer cmo copy motmp2 motmp cmo setatt motmp2 zic Z_TOP Warning: nelements = 0 2500 values reset for attribute zic cmo setatt motmp2 imat Z_TOP_ID Warning: nelements = 0 2500 values reset for attribute imat cmo setatt motmp2 xval Z_TOP_VAL Warning: nelements = 0 2500 values reset for attribute xval cmo setatt motmp2 nlayer TOP_LAYER Warning: nelements = 0 2500 values reset for attribute nlayer #### MERGE ####### cmo printatt mo -all- minmax Warning: nelements = 0 ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH -def- 0.000000000E+00 0.000000000E+00 0.000000000E+00 117450 scalar 1 1 0 1 vector 3 3 0 1 nnodes 117450 117450 0 1 nedges 0 0 0 1 nfaces 0 0 0 1 nelements 0 0 0 1 mbndry 16000000 16000000 0 1 ndimensions_topo 3 3 0 1 ndimensions_geom 3 3 0 1 nodes_per_element 4 4 0 1 edges_per_element 6 6 0 1 faces_per_element 4 4 0 1 isetwd -3 -3 0 117450 ialias 0 0 0 117450 imt1 1 47 46 117450 itp1 0 0 0 117450 icr1 0 0 0 117450 isn1 0 0 0 117450 xic 4.201000000E+05 4.299000000E+05 9.800000000E+03 117450 yic 4.307100000E+06 4.316900000E+06 9.800000000E+03 117450 zic -2.208199951E+03 1.081999969E+02 2.316399948E+03 117450 0 length attribute: xtetwd 0 length attribute: itetclr 0 length attribute: itettyp 0 length attribute: itetoff 0 length attribute: jtetoff 0 length attribute: itet 0 length attribute: jtet epsilon 1.000000004E-15 1.000000004E-15 0.000000000E+00 1 epsilonl 3.120068113E-09 3.120068113E-09 0.000000000E+00 1 epsilona 4.384175439E-05 4.384175439E-05 0.000000000E+00 1 epsilonv 4.939760845E-02 4.939760845E-02 0.000000000E+00 1 ipointi 1 1 0 1 ipointj 117450 117450 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.201000000E+05 4.201000000E+05 0.000000000E+00 1 ymin 4.307100000E+06 4.307100000E+06 0.000000000E+00 1 zmin -2.208199951E+03 -2.208199951E+03 0.000000000E+00 1 xmax 4.299000000E+05 4.299000000E+05 0.000000000E+00 1 ymax 4.316900000E+06 4.316900000E+06 0.000000000E+00 1 zmax 1.081999969E+02 1.081999969E+02 0.000000000E+00 1 kdtree_level 0 0 0 1 max_number_sets 64 64 0 1 number_of_psets 2 2 0 1 number_of_eltsets 0 0 0 1 number_of_fsets 0 0 0 1 xval 1.000000000E+00 9.000000000E+00 8.000000000E+00 117450 nlayer 1 47 46 117450 imat 1 9 8 117450 cmo printatt motmp -all- minmax Warning: nelements = 0 ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH -def- 0.000000000E+00 0.000000000E+00 0.000000000E+00 2500 scalar 1 1 0 1 vector 3 3 0 1 nnodes 2500 2500 0 1 nedges 0 0 0 1 nfaces 0 0 0 1 nelements 0 0 0 1 mbndry 16000000 16000000 0 1 ndimensions_topo 3 3 0 1 ndimensions_geom 3 3 0 1 nodes_per_element 4 4 0 1 edges_per_element 6 6 0 1 faces_per_element 4 4 0 1 isetwd -3 -3 0 2500 ialias 0 0 0 2500 imt1 1 1 0 2500 itp1 0 0 0 2500 icr1 0 0 0 2500 isn1 0 0 0 2500 xic 4.201000000E+05 4.299000000E+05 9.800000000E+03 2500 yic 4.307100000E+06 4.316900000E+06 9.800000000E+03 2500 zic -2.400000000E+03 -2.400000000E+03 0.000000000E+00 2500 0 length attribute: xtetwd 0 length attribute: itetclr 0 length attribute: itettyp 0 length attribute: itetoff 0 length attribute: jtetoff 0 length attribute: itet 0 length attribute: jtet epsilon 1.000000004E-15 1.000000004E-15 0.000000000E+00 1 epsilonl 3.120068113E-09 3.120068113E-09 0.000000000E+00 1 epsilona 4.384175439E-05 4.384175439E-05 0.000000000E+00 1 epsilonv 4.939760845E-02 4.939760845E-02 0.000000000E+00 1 ipointi 1 1 0 1 ipointj 2500 2500 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.201000000E+05 4.201000000E+05 0.000000000E+00 1 ymin 4.307100000E+06 4.307100000E+06 0.000000000E+00 1 zmin -2.208199951E+03 -2.208199951E+03 0.000000000E+00 1 xmax 4.299000000E+05 4.299000000E+05 0.000000000E+00 1 ymax 4.316900000E+06 4.316900000E+06 0.000000000E+00 1 zmax 1.081999969E+02 1.081999969E+02 0.000000000E+00 1 kdtree_level 0 0 0 1 max_number_sets 64 64 0 1 number_of_psets 2 2 0 1 number_of_eltsets 0 0 0 1 number_of_fsets 0 0 0 1 xval 1.000000000E+01 1.000000000E+01 0.000000000E+00 2500 nlayer 0 0 0 2500 imat 10 10 0 2500 cmo printatt motmp2 -all- minmax Warning: nelements = 0 ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH -def- 0.000000000E+00 0.000000000E+00 0.000000000E+00 2500 scalar 1 1 0 1 vector 3 3 0 1 nnodes 2500 2500 0 1 nedges 0 0 0 1 nfaces 0 0 0 1 nelements 0 0 0 1 mbndry 16000000 16000000 0 1 ndimensions_topo 3 3 0 1 ndimensions_geom 3 3 0 1 nodes_per_element 4 4 0 1 edges_per_element 6 6 0 1 faces_per_element 4 4 0 1 isetwd -3 -3 0 2500 ialias 0 0 0 2500 imt1 1 1 0 2500 itp1 0 0 0 2500 icr1 0 0 0 2500 isn1 0 0 0 2500 xic 4.201000000E+05 4.299000000E+05 9.800000000E+03 2500 yic 4.307100000E+06 4.316900000E+06 9.800000000E+03 2500 zic 3.000000000E+02 3.000000000E+02 0.000000000E+00 2500 0 length attribute: xtetwd 0 length attribute: itetclr 0 length attribute: itettyp 0 length attribute: itetoff 0 length attribute: jtetoff 0 length attribute: itet 0 length attribute: jtet epsilon 1.000000004E-15 1.000000004E-15 0.000000000E+00 1 epsilonl 3.120068113E-09 3.120068113E-09 0.000000000E+00 1 epsilona 4.384175439E-05 4.384175439E-05 0.000000000E+00 1 epsilonv 4.939760845E-02 4.939760845E-02 0.000000000E+00 1 ipointi 1 1 0 1 ipointj 2500 2500 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.201000000E+05 4.201000000E+05 0.000000000E+00 1 ymin 4.307100000E+06 4.307100000E+06 0.000000000E+00 1 zmin -2.208199951E+03 -2.208199951E+03 0.000000000E+00 1 xmax 4.299000000E+05 4.299000000E+05 0.000000000E+00 1 ymax 4.316900000E+06 4.316900000E+06 0.000000000E+00 1 zmax 1.081999969E+02 1.081999969E+02 0.000000000E+00 1 kdtree_level 0 0 0 1 max_number_sets 64 64 0 1 number_of_psets 2 2 0 1 number_of_eltsets 0 0 0 1 number_of_fsets 0 0 0 1 xval 1.100000000E+01 1.100000000E+01 0.000000000E+00 2500 nlayer 50 50 0 2500 imat 11 11 0 2500 cmo copy motet motmp copypts/motet/mo/0,0/1,0,0 copypts added nodes 2501 to 119950 copypts/motet/motmp2/0,0/1,0,0 copypts added nodes 119951 to 122450 cmo printatt motet -all- minmax Warning: nelements = 0 ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH -def- 0.000000000E+00 0.000000000E+00 0.000000000E+00 122450 scalar 1 1 0 1 vector 3 3 0 1 nnodes 122450 122450 0 1 nedges 0 0 0 1 nfaces 0 0 0 1 nelements 0 0 0 1 mbndry 16000000 16000000 0 1 ndimensions_topo 3 3 0 1 ndimensions_geom 3 3 0 1 nodes_per_element 4 4 0 1 edges_per_element 6 6 0 1 faces_per_element 4 4 0 1 isetwd -3 -3 0 122450 ialias 0 0 0 122450 imt1 1 47 46 122450 itp1 0 0 0 122450 icr1 0 0 0 122450 isn1 0 0 0 122450 xic 4.201000000E+05 4.299000000E+05 9.800000000E+03 122450 yic 4.307100000E+06 4.316900000E+06 9.800000000E+03 122450 zic -2.400000000E+03 3.000000000E+02 2.700000000E+03 122450 0 length attribute: xtetwd 0 length attribute: itetclr 0 length attribute: itettyp 0 length attribute: itetoff 0 length attribute: jtetoff 0 length attribute: itet 0 length attribute: jtet epsilon 1.000000004E-15 1.000000004E-15 0.000000000E+00 1 epsilonl 3.135235225E-09 3.135235225E-09 0.000000000E+00 1 epsilona 4.426903288E-05 4.426903288E-05 0.000000000E+00 1 epsilonv 5.757794241E-02 5.757794241E-02 0.000000000E+00 1 ipointi 119951 119951 0 1 ipointj 122450 122450 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.201000000E+05 4.201000000E+05 0.000000000E+00 1 ymin 4.307100000E+06 4.307100000E+06 0.000000000E+00 1 zmin -2.400000000E+03 -2.400000000E+03 0.000000000E+00 1 xmax 4.299000000E+05 4.299000000E+05 0.000000000E+00 1 ymax 4.316900000E+06 4.316900000E+06 0.000000000E+00 1 zmax 3.000000000E+02 3.000000000E+02 0.000000000E+00 1 kdtree_level 0 0 0 1 max_number_sets 64 64 0 1 number_of_psets 2 2 0 1 number_of_eltsets 0 0 0 1 number_of_fsets 0 0 0 1 xval 1.000000000E+00 1.100000000E+01 1.000000000E+01 122450 nlayer 0 50 50 122450 imat 1 11 10 122450 dump gmv points_to_connect.gmv motet cmo/modatt/-def-/-def-/ioflag/x finish # STOP and Evaluate subset points finish LaGriT successfully completed