* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * Program: LaGriT V3.101 Linux m64 * * * date_compile: 2013/04/25 RH gf4.5 * * * Run Time: 2013/Nov 6 12:13:36 * * * 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 # Creates a 3D 3D 3D ice wedge # read the outer polygon points (no elements) # connect points into flat shape of triangles cmo/create/mo///triplane read/avs/poly_sim1.inp/mo cmo/status/brief The current-mesh-object(CMO) is: mo 1 Mesh Object name: mo number of nodes = 18 number of elements = 0 dimensions geometry = 2 element type = tri dimensions topology = 2 3 nodes 3 faces 3 edges boundary flag = 16000000 status = active finish triangulate/counterclockwise filter/1,0,0 FILTER:Use internal epsilonl value = 0.107713032103E-10 FILTER:Dudding 2 duplicate points. finish rmpoint/compress 2 points removed and 0 elements removed. The new point count is 16. The new element count is 16. finish resetpts/itp geniee finish finish cmo/setatt/mo/imt/1,0,0/1 16 values reset for attribute imt1 cmo/setatt/mo/itetclr/1,0,0/1 16 values reset for attribute itetclr resetpts/itp geniee finish dump/tmp_poly_tri1.gmv/mo dump/gmv/tmp_poly_tri1.gmv/mo cmo/modatt/-def-/-def-/ioflag/x finish finish # check triangulation quality epsilonl, epsilonaspect: 1.0771303E-11 1.2496971E-33 -------------------------------------------- elements with aspect ratio < .01: 3 elements with aspect ratio b/w .01 and .02: 1 elements with aspect ratio b/w .02 and .05: 0 elements with aspect ratio b/w .05 and .1 : 0 elements with aspect ratio b/w .1 and .2 : 0 elements with aspect ratio b/w .2 and .5 : 4 elements with aspect ratio b/w .5 and 1. : 8 min aspect ratio = 0.1531E-04 max aspect ratio = 0.2610E+00 epsilonvol: 5.2251201E-10 --------------------------------------- element volumes b/w 0.2534E-01 and 0.8997E-01: 2 element volumes b/w 0.8997E-01 and 0.3195E+00: 0 element volumes b/w 0.3195E+00 and 0.1134E+01: 0 element volumes b/w 0.1134E+01 and 0.4029E+01: 0 element volumes b/w 0.4029E+01 and 0.1431E+02: 14 min volume = 2.5336035E-02 max volume = 1.4305248E+01 ----------------------------------------------------------- 16 total elements evaluated. # read the inner polygon points (no elements) # connect points into flat shape of triangles cmo/create/mo_inner///triplane read/avs/poly_sim2.inp/mo_inner cmo/status/brief The current-mesh-object(CMO) is: mo_inner 1 Mesh Object name: mo number of nodes = 16 number of elements = 16 dimensions geometry = 2 element type = tri dimensions topology = 2 3 nodes 3 faces 3 edges boundary flag = 16000000 status = inactive 2 Mesh Object name: mo_inner number of nodes = 9 number of elements = 0 dimensions geometry = 2 element type = tri dimensions topology = 2 3 nodes 3 faces 3 edges boundary flag = 16000000 status = active finish triangulate/counterclockwise filter/1,0,0 FILTER:Use internal epsilonl value = 0.978291118521E-11 FILTER:Dudding 0 duplicate points. finish rmpoint/compress 0 points removed and 0 elements removed. The new point count is 9. The new element count is 7. finish resetpts/itp geniee finish finish cmo/setatt/mo_inner/imt/1,0,0/2 9 values reset for attribute imt1 cmo/setatt/mo_inner/itetclr/1,0,0/2 7 values reset for attribute itetclr resetpts/itp geniee finish dump/gmv/tmp_poly_tri2.gmv/mo_inner cmo/modatt/-def-/-def-/ioflag/x finish ############### CHECK tri surfs ############# # the number of triangles is related to number of poly points cmo select mo quality epsilonl, epsilonaspect: 1.0771303E-11 1.2496971E-33 -------------------------------------------- elements with aspect ratio < .01: 3 elements with aspect ratio b/w .01 and .02: 1 elements with aspect ratio b/w .02 and .05: 0 elements with aspect ratio b/w .05 and .1 : 0 elements with aspect ratio b/w .1 and .2 : 0 elements with aspect ratio b/w .2 and .5 : 4 elements with aspect ratio b/w .5 and 1. : 8 min aspect ratio = 0.1531E-04 max aspect ratio = 0.2610E+00 epsilonvol: 5.2251201E-10 --------------------------------------- element volumes b/w 0.2534E-01 and 0.8997E-01: 2 element volumes b/w 0.8997E-01 and 0.3195E+00: 0 element volumes b/w 0.3195E+00 and 0.1134E+01: 0 element volumes b/w 0.1134E+01 and 0.4029E+01: 0 element volumes b/w 0.4029E+01 and 0.1431E+02: 14 min volume = 2.5336035E-02 max volume = 1.4305248E+01 ----------------------------------------------------------- 16 total elements evaluated. cmo select mo_inner quality epsilonl, epsilonaspect: 9.7829112E-12 9.3627695E-34 -------------------------------------------- elements with aspect ratio < .01: 2 elements with aspect ratio b/w .01 and .02: 0 elements with aspect ratio b/w .02 and .05: 0 elements with aspect ratio b/w .05 and .1 : 0 elements with aspect ratio b/w .1 and .2 : 0 elements with aspect ratio b/w .2 and .5 : 0 elements with aspect ratio b/w .5 and 1. : 5 min aspect ratio = 0.1076E-02 max aspect ratio = 0.9446E+00 epsilonvol: 4.3101858E-10 --------------------------------------- element volumes b/w 0.2540E+00 and 0.8900E+00: 2 element volumes b/w 0.8900E+00 and 0.3118E+01: 0 element volumes b/w 0.3118E+01 and 0.1093E+02: 0 element volumes b/w 0.1093E+02 and 0.3829E+02: 0 element volumes b/w 0.3829E+02 and 0.1342E+03: 5 min volume = 2.5399350E-01 max volume = 1.3415328E+02 ----------------------------------------------------------- 7 total elements evaluated. # Both shapes should be within same area, 2nd poly inside cmo printatt mo -xyz- minmax ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH xic 1.030189210E+02 1.445391350E+02 4.152021400E+01 16 yic 1.234896800E+02 1.485746700E+02 2.508499000E+01 16 zic 0.000000000E+00 0.000000000E+00 0.000000000E+00 16 cmo printatt mo_inner -xyz- minmax ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH xic 1.048481294E+02 1.431229863E+02 3.827485692E+01 9 yic 1.249985574E+02 1.468198879E+02 2.182133046E+01 9 zic 0.000000000E+00 0.000000000E+00 0.000000000E+00 9 # Write merged poly shapes addmesh/merge/mo_merge/mo/mo_inner dump/gmv/tmp_poly.gmv/mo_merge cmo/modatt/-def-/-def-/ioflag/x finish # There should be 2 materials cmo/printatt/mo_merge/imt/minmax ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH imt1 1 2 1 25 cmo/printatt/mo_merge/itetclr/minmax ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH itetclr 1 2 1 23 ############################################# # Break big triangles into smaller for better distribution cmo select mo_merge massage/5/1.e-5/1.e-5/strictmergelength nodes in point set = 25 recon/0/0.10000000E-04/ RECON2D reconnection loop2d: 1 Recon2d: iteration number= 1 number of flips= 6 Recon2d: iteration number= 2 number of flips= 2 Recon2d: iteration number= 3 number of flips= 2 RECON: done. finish nodes in point set = 25 SGD: 0 total relaxations on point set of 25 indept. nodes. nodes in point set = 25 Face iteration, number of faces: 1 6 Face iteration, number of faces: 1 10 Face iteration, number of faces: 1 19 Face iteration, number of faces: 1 15 Face iteration, number of faces: 1 18 Face iteration, number of faces: 1 17 Face iteration, number of faces: 1 15 Face iteration, number of faces: 1 12 Face iteration, number of faces: 1 12 Face iteration, number of faces: 1 4 Face iteration, number of faces: 1 1 recon/0/0.10000000E-04/ RECON2D reconnection loop2d: 1 Recon2d: iteration number= 1 number of flips= 33 Recon2d: iteration number= 2 number of flips= 14 RECON: done. finish smooth/position/esug/1,0,0 SMOOTH: position option nodes in point set = 154 SMOOTH: Smooth 2D triangular mesh using ESUG Min. area BEFORE smoothing= 0.8873791E-04 Max. area BEFORE smoothing= 0.9470015E+01 Iteration= 1; root mean square error= 0.3807578E+00 Iteration= 2; root mean square error= 0.2705322E+00 Iteration= 3; root mean square error= 0.1458594E+00 Iteration= 4; root mean square error= 0.7747996E-01 Iteration= 5; root mean square error= 0.4680546E-01 Iteration= 6; root mean square error= 0.3083506E-01 Iteration= 7; root mean square error= 0.2381777E-01 Iteration= 8; root mean square error= 0.2072364E-01 Iteration= 9; root mean square error= 0.1963048E-01 Iteration= 10; root mean square error= 0.1727875E-01 Iteration= 11; root mean square error= 0.1355319E-01 Iteration= 12; root mean square error= 0.1130126E-01 Iteration= 13; root mean square error= 0.8328037E-02 Iteration= 14; root mean square error= 0.2963465E-02 Min. area AFTER smoothing= 0.4480048E-01 Max. area AFTER smoothing= 0.8392478E+01 recon 0 RECON2D reconnection loop2d: 1 Recon2d: iteration number= 1 number of flips= 18 Recon2d: iteration number= 2 number of flips= 6 Recon2d: iteration number= 3 number of flips= 2 RECON: done. smooth/position/esug/1,0,0 SMOOTH: position option nodes in point set = 154 SMOOTH: Smooth 2D triangular mesh using ESUG Min. area BEFORE smoothing= 0.4100831E-01 Max. area BEFORE smoothing= 0.8392478E+01 Iteration= 1; root mean square error= 0.2604119E+00 Iteration= 2; root mean square error= 0.1395853E+00 Iteration= 3; root mean square error= 0.5818236E-01 Iteration= 4; root mean square error= 0.3025865E-01 Iteration= 5; root mean square error= 0.2298531E-01 Iteration= 6; root mean square error= 0.1887356E-01 Iteration= 7; root mean square error= 0.1581417E-01 Iteration= 8; root mean square error= 0.1336645E-01 Iteration= 9; root mean square error= 0.1134880E-01 Iteration= 10; root mean square error= 0.9660507E-02 Iteration= 11; root mean square error= 0.8235859E-02 Iteration= 12; root mean square error= 0.7027705E-02 Iteration= 13; root mean square error= 0.6000089E-02 Iteration= 14; root mean square error= 0.5124448E-02 Iteration= 15; root mean square error= 0.4377475E-02 Iteration= 16; root mean square error= 0.3739825E-02 Iteration= 17; root mean square error= 0.3195265E-02 Iteration= 18; root mean square error= 0.2730085E-02 Iteration= 19; root mean square error= 0.2332652E-02 Min. area AFTER smoothing= 0.4480048E-01 Max. area AFTER smoothing= 0.7818767E+01 recon 0 RECON2D reconnection loop2d: 1 Recon2d: iteration number= 1 number of flips= 8 RECON: done. smooth/position/esug/1,0,0 SMOOTH: position option nodes in point set = 154 SMOOTH: Smooth 2D triangular mesh using ESUG Min. area BEFORE smoothing= 0.4480048E-01 Max. area BEFORE smoothing= 0.7818767E+01 Iteration= 1; root mean square error= 0.1769684E+00 Iteration= 2; root mean square error= 0.7603170E-01 Iteration= 3; root mean square error= 0.3553703E-01 Iteration= 4; root mean square error= 0.2238394E-01 Iteration= 5; root mean square error= 0.1617508E-01 Iteration= 6; root mean square error= 0.1236252E-01 Iteration= 7; root mean square error= 0.9726600E-02 Iteration= 8; root mean square error= 0.7809633E-02 Iteration= 9; root mean square error= 0.6373162E-02 Iteration= 10; root mean square error= 0.5270093E-02 Iteration= 11; root mean square error= 0.4403881E-02 Iteration= 12; root mean square error= 0.3709841E-02 Iteration= 13; root mean square error= 0.3144101E-02 Iteration= 14; root mean square error= 0.2676445E-02 Iteration= 15; root mean square error= 0.2285630E-02 Iteration= 16; root mean square error= 0.1956331E-02 Iteration= 17; root mean square error= 0.1677178E-02 Min. area AFTER smoothing= 0.4480048E-01 Max. area AFTER smoothing= 0.7570071E+01 recon 0 RECON2D reconnection loop2d: 1 Recon2d: iteration number= 1 number of flips= 1 RECON: done. smooth/position/esug/1,0,0 SMOOTH: position option nodes in point set = 154 SMOOTH: Smooth 2D triangular mesh using ESUG Min. area BEFORE smoothing= 0.4480048E-01 Max. area BEFORE smoothing= 0.7570071E+01 Iteration= 1; root mean square error= 0.6610979E-01 Iteration= 2; root mean square error= 0.1807925E-01 Iteration= 3; root mean square error= 0.9076955E-02 Iteration= 4; root mean square error= 0.5608533E-02 Iteration= 5; root mean square error= 0.3920864E-02 Iteration= 6; root mean square error= 0.2872525E-02 Iteration= 7; root mean square error= 0.2157107E-02 Iteration= 8; root mean square error= 0.1648499E-02 Iteration= 9; root mean square error= 0.1277148E-02 Iteration= 10; root mean square error= 0.1000810E-02 Iteration= 11; root mean square error= 0.7922049E-03 Iteration= 12; root mean square error= 0.6328947E-03 Min. area AFTER smoothing= 0.4480048E-01 Max. area AFTER smoothing= 0.7546913E+01 recon 0 RECON2D reconnection loop2d: 1 Recon2d: iteration number= 1 number of flips= 1 RECON: done. smooth/position/esug/1,0,0 SMOOTH: position option nodes in point set = 154 SMOOTH: Smooth 2D triangular mesh using ESUG Min. area BEFORE smoothing= 0.4480048E-01 Max. area BEFORE smoothing= 0.7546913E+01 Iteration= 1; root mean square error= 0.8796955E-02 Iteration= 2; root mean square error= 0.4045150E-02 Iteration= 3; root mean square error= 0.2128849E-02 Iteration= 4; root mean square error= 0.1395891E-02 Iteration= 5; root mean square error= 0.9747832E-03 Iteration= 6; root mean square error= 0.7089844E-03 Iteration= 7; root mean square error= 0.5306092E-03 Iteration= 8; root mean square error= 0.4056227E-03 Iteration= 9; root mean square error= 0.3157035E-03 Iteration= 10; root mean square error= 0.2497183E-03 Iteration= 11; root mean square error= 0.2004419E-03 Iteration= 12; root mean square error= 0.1630202E-03 Iteration= 13; root mean square error= 0.1341271E-03 Iteration= 14; root mean square error= 0.1114541E-03 Iteration= 15; root mean square error= 0.9338473E-04 Iteration= 16; root mean square error= 0.7877731E-04 Min. area AFTER smoothing= 0.4480048E-01 Max. area AFTER smoothing= 0.7539480E+01 recon 0 RECON2D reconnection loop2d: 1 RECON: done. smooth/position/esug/1,0,0 SMOOTH: position option nodes in point set = 154 SMOOTH: Smooth 2D triangular mesh using ESUG Min. area BEFORE smoothing= 0.4480048E-01 Max. area BEFORE smoothing= 0.7539480E+01 Iteration= 1; root mean square error= 0.6681758E-04 Iteration= 2; root mean square error= 0.5691796E-04 Iteration= 3; root mean square error= 0.4864832E-04 Iteration= 4; root mean square error= 0.4168850E-04 Iteration= 5; root mean square error= 0.3579591E-04 Iteration= 6; root mean square error= 0.3078332E-04 Iteration= 7; root mean square error= 0.2650357E-04 Iteration= 8; root mean square error= 0.2283912E-04 Iteration= 9; root mean square error= 0.1969462E-04 Iteration= 10; root mean square error= 0.1699176E-04 Iteration= 11; root mean square error= 0.1466553E-04 Iteration= 12; root mean square error= 0.1266150E-04 Iteration= 13; root mean square error= 0.1093376E-04 Iteration= 14; root mean square error= 0.9443362E-05 Iteration= 15; root mean square error= 0.8157160E-05 Iteration= 16; root mean square error= 0.7046811E-05 Iteration= 17; root mean square error= 0.6088036E-05 Iteration= 18; root mean square error= 0.5259988E-05 Iteration= 19; root mean square error= 0.4544741E-05 Iteration= 20; root mean square error= 0.3926863E-05 Iteration= 21; root mean square error= 0.3393057E-05 Iteration= 22; root mean square error= 0.2931856E-05 Iteration= 23; root mean square error= 0.2533367E-05 Iteration= 24; root mean square error= 0.2189053E-05 Iteration= 25; root mean square error= 0.1891542E-05 Min. area AFTER smoothing= 0.4480048E-01 Max. area AFTER smoothing= 0.7538605E+01 recon 0 RECON2D reconnection loop2d: 1 RECON: done. massage/1/1.e-5/1.e-5/strictmergelength nodes in point set = 154 recon/0/0.10000000E-04/ RECON2D reconnection loop2d: 1 RECON: done. finish nodes in point set = 154 SGD: 9 total relaxations on point set of 154 indept. nodes. rmpoint/compress 0 points removed and 0 elements removed. The new point count is 154. The new element count is 228. finish nodes in point set = 154 recon/0/0.10000000E-04/ RECON2D reconnection loop2d: 1 Recon2d: iteration number= 1 number of flips= 6 RECON: done. finish nodes in point set = 154 SGD: 5 total relaxations on point set of 154 indept. nodes. rmpoint/compress 0 points removed and 0 elements removed. The new point count is 154. The new element count is 228. finish nodes in point set = 154 recon/0/0.10000000E-04/ RECON2D reconnection loop2d: 1 RECON: done. finish nodes in point set = 154 SGD: 0 total relaxations on point set of 154 indept. nodes. nodes in point set = 154 Face iteration, number of faces: 1 74 Face iteration, number of faces: 1 72 Face iteration, number of faces: 1 111 Face iteration, number of faces: 1 145 Face iteration, number of faces: 1 169 Face iteration, number of faces: 1 215 Face iteration, number of faces: 1 225 Face iteration, number of faces: 1 237 Face iteration, number of faces: 1 237 Face iteration, number of faces: 1 243 Face iteration, number of faces: 1 225 Face iteration, number of faces: 1 156 Face iteration, number of faces: 1 63 Face iteration, number of faces: 1 11 recon/0/0.10000000E-04/ RECON2D reconnection loop2d: 1 Recon2d: iteration number= 1 number of flips= 351 Recon2d: iteration number= 2 number of flips= 26 Recon2d: iteration number= 3 number of flips= 4 Recon2d: iteration number= 4 number of flips= 1 RECON: done. finish smooth/position/esug/1,0,0 SMOOTH: position option nodes in point set = 2337 SMOOTH: Smooth 2D triangular mesh using ESUG Min. area BEFORE smoothing= 0.3178610E-02 Max. area BEFORE smoothing= 0.3221297E+00 Iteration= 1; root mean square error= 0.7216973E-01 Iteration= 2; root mean square error= 0.2875830E-01 Iteration= 3; root mean square error= 0.1469273E-01 Iteration= 4; root mean square error= 0.1015737E-01 Iteration= 5; root mean square error= 0.8017529E-02 Iteration= 6; root mean square error= 0.6733688E-02 Iteration= 7; root mean square error= 0.5852517E-02 Iteration= 8; root mean square error= 0.5200298E-02 Iteration= 9; root mean square error= 0.4693340E-02 Iteration= 10; root mean square error= 0.4284973E-02 Iteration= 11; root mean square error= 0.3946819E-02 Iteration= 12; root mean square error= 0.3660587E-02 Iteration= 13; root mean square error= 0.3413970E-02 Iteration= 14; root mean square error= 0.3198388E-02 Iteration= 15; root mean square error= 0.3007680E-02 Iteration= 16; root mean square error= 0.2837303E-02 Iteration= 17; root mean square error= 0.2683827E-02 Iteration= 18; root mean square error= 0.2544611E-02 Iteration= 19; root mean square error= 0.2417576E-02 Iteration= 20; root mean square error= 0.2301066E-02 Iteration= 21; root mean square error= 0.2193734E-02 Iteration= 22; root mean square error= 0.2094475E-02 Iteration= 23; root mean square error= 0.2002368E-02 Iteration= 24; root mean square error= 0.1916637E-02 Iteration= 25; root mean square error= 0.1836625E-02 Min. area AFTER smoothing= 0.2202948E-01 Max. area AFTER smoothing= 0.3506410E+00 recon 0 RECON2D reconnection loop2d: 1 Recon2d: iteration number= 1 number of flips= 141 Recon2d: iteration number= 2 number of flips= 2 Recon2d: iteration number= 3 number of flips= 2 RECON: done. smooth/position/esug/1,0,0 SMOOTH: position option nodes in point set = 2337 SMOOTH: Smooth 2D triangular mesh using ESUG Min. area BEFORE smoothing= 0.2202948E-01 Max. area BEFORE smoothing= 0.3506410E+00 Iteration= 1; root mean square error= 0.4212969E-01 Iteration= 2; root mean square error= 0.1602277E-01 Iteration= 3; root mean square error= 0.6465802E-02 Iteration= 4; root mean square error= 0.4018310E-02 Iteration= 5; root mean square error= 0.3019440E-02 Iteration= 6; root mean square error= 0.2496622E-02 Iteration= 7; root mean square error= 0.2178013E-02 Iteration= 8; root mean square error= 0.1962512E-02 Iteration= 9; root mean square error= 0.1804676E-02 Iteration= 10; root mean square error= 0.1681802E-02 Iteration= 11; root mean square error= 0.1581621E-02 Iteration= 12; root mean square error= 0.1497062E-02 Iteration= 13; root mean square error= 0.1423811E-02 Iteration= 14; root mean square error= 0.1359097E-02 Iteration= 15; root mean square error= 0.1301057E-02 Iteration= 16; root mean square error= 0.1248388E-02 Iteration= 17; root mean square error= 0.1200145E-02 Iteration= 18; root mean square error= 0.1155620E-02 Iteration= 19; root mean square error= 0.1114273E-02 Iteration= 20; root mean square error= 0.1075681E-02 Iteration= 21; root mean square error= 0.1039503E-02 Iteration= 22; root mean square error= 0.1005466E-02 Iteration= 23; root mean square error= 0.9733438E-03 Iteration= 24; root mean square error= 0.9429461E-03 Iteration= 25; root mean square error= 0.9141137E-03 Min. area AFTER smoothing= 0.2135722E-01 Max. area AFTER smoothing= 0.3370724E+00 recon 0 RECON2D reconnection loop2d: 1 Recon2d: iteration number= 1 number of flips= 47 RECON: done. smooth/position/esug/1,0,0 SMOOTH: position option nodes in point set = 2337 SMOOTH: Smooth 2D triangular mesh using ESUG Min. area BEFORE smoothing= 0.2135722E-01 Max. area BEFORE smoothing= 0.3370724E+00 Iteration= 1; root mean square error= 0.2483761E-01 Iteration= 2; root mean square error= 0.1013069E-01 Iteration= 3; root mean square error= 0.4471117E-02 Iteration= 4; root mean square error= 0.2982078E-02 Iteration= 5; root mean square error= 0.2300648E-02 Iteration= 6; root mean square error= 0.1908735E-02 Iteration= 7; root mean square error= 0.1658022E-02 Iteration= 8; root mean square error= 0.1485672E-02 Iteration= 9; root mean square error= 0.1360137E-02 Iteration= 10; root mean square error= 0.1264145E-02 Iteration= 11; root mean square error= 0.1187670E-02 Iteration= 12; root mean square error= 0.1124628E-02 Iteration= 13; root mean square error= 0.1071180E-02 Iteration= 14; root mean square error= 0.1024815E-02 Iteration= 15; root mean square error= 0.9838388E-03 Iteration= 16; root mean square error= 0.9470735E-03 Iteration= 17; root mean square error= 0.9136776E-03 Iteration= 18; root mean square error= 0.8830357E-03 Iteration= 19; root mean square error= 0.8546880E-03 Iteration= 20; root mean square error= 0.8282839E-03 Iteration= 21; root mean square error= 0.8035512E-03 Iteration= 22; root mean square error= 0.7802751E-03 Iteration= 23; root mean square error= 0.7582834E-03 Iteration= 24; root mean square error= 0.7374362E-03 Iteration= 25; root mean square error= 0.7176180E-03 Min. area AFTER smoothing= 0.2121985E-01 Max. area AFTER smoothing= 0.3285175E+00 recon 0 RECON2D reconnection loop2d: 1 Recon2d: iteration number= 1 number of flips= 19 RECON: done. smooth/position/esug/1,0,0 SMOOTH: position option nodes in point set = 2337 SMOOTH: Smooth 2D triangular mesh using ESUG Min. area BEFORE smoothing= 0.2121985E-01 Max. area BEFORE smoothing= 0.3285175E+00 Iteration= 1; root mean square error= 0.1544333E-01 Iteration= 2; root mean square error= 0.6060531E-02 Iteration= 3; root mean square error= 0.2828812E-02 Iteration= 4; root mean square error= 0.1914662E-02 Iteration= 5; root mean square error= 0.1500378E-02 Iteration= 6; root mean square error= 0.1264076E-02 Iteration= 7; root mean square error= 0.1113287E-02 Iteration= 8; root mean square error= 0.1008603E-02 Iteration= 9; root mean square error= 0.9310210E-03 Iteration= 10; root mean square error= 0.8705363E-03 Iteration= 11; root mean square error= 0.8215017E-03 Iteration= 12; root mean square error= 0.7805313E-03 Iteration= 13; root mean square error= 0.7454817E-03 Iteration= 14; root mean square error= 0.7149297E-03 Iteration= 15; root mean square error= 0.6878901E-03 Iteration= 16; root mean square error= 0.6636571E-03 Iteration= 17; root mean square error= 0.6417106E-03 Iteration= 18; root mean square error= 0.6216574E-03 Iteration= 19; root mean square error= 0.6031948E-03 Iteration= 20; root mean square error= 0.5860854E-03 Iteration= 21; root mean square error= 0.5701403E-03 Iteration= 22; root mean square error= 0.5552073E-03 Iteration= 23; root mean square error= 0.5411622E-03 Iteration= 24; root mean square error= 0.5279026E-03 Iteration= 25; root mean square error= 0.5153431E-03 Min. area AFTER smoothing= 0.2123605E-01 Max. area AFTER smoothing= 0.3303549E+00 recon 0 RECON2D reconnection loop2d: 1 Recon2d: iteration number= 1 number of flips= 2 RECON: done. smooth/position/esug/1,0,0 SMOOTH: position option nodes in point set = 2337 SMOOTH: Smooth 2D triangular mesh using ESUG Min. area BEFORE smoothing= 0.2123605E-01 Max. area BEFORE smoothing= 0.3303549E+00 Iteration= 1; root mean square error= 0.5217928E-02 Iteration= 2; root mean square error= 0.2099332E-02 Iteration= 3; root mean square error= 0.1037197E-02 Iteration= 4; root mean square error= 0.7936706E-03 Iteration= 5; root mean square error= 0.6802550E-03 Iteration= 6; root mean square error= 0.6130688E-03 Iteration= 7; root mean square error= 0.5678881E-03 Iteration= 8; root mean square error= 0.5346321E-03 Iteration= 9; root mean square error= 0.5085416E-03 Iteration= 10; root mean square error= 0.4871028E-03 Iteration= 11; root mean square error= 0.4688760E-03 Iteration= 12; root mean square error= 0.4529826E-03 Iteration= 13; root mean square error= 0.4388570E-03 Iteration= 14; root mean square error= 0.4261180E-03 Iteration= 15; root mean square error= 0.4144981E-03 Iteration= 16; root mean square error= 0.4038022E-03 Iteration= 17; root mean square error= 0.3938835E-03 Iteration= 18; root mean square error= 0.3846285E-03 Iteration= 19; root mean square error= 0.3759475E-03 Iteration= 20; root mean square error= 0.3677676E-03 Iteration= 21; root mean square error= 0.3600295E-03 Iteration= 22; root mean square error= 0.3526835E-03 Iteration= 23; root mean square error= 0.3456879E-03 Iteration= 24; root mean square error= 0.3390072E-03 Iteration= 25; root mean square error= 0.3326112E-03 Min. area AFTER smoothing= 0.2124897E-01 Max. area AFTER smoothing= 0.3306224E+00 recon 0 RECON2D reconnection loop2d: 1 Recon2d: iteration number= 1 number of flips= 2 RECON: done. smooth/position/esug/1,0,0 SMOOTH: position option nodes in point set = 2337 SMOOTH: Smooth 2D triangular mesh using ESUG Min. area BEFORE smoothing= 0.2124897E-01 Max. area BEFORE smoothing= 0.3306224E+00 Iteration= 1; root mean square error= 0.5892357E-02 Iteration= 2; root mean square error= 0.2918728E-02 Iteration= 3; root mean square error= 0.1398954E-02 Iteration= 4; root mean square error= 0.9323396E-03 Iteration= 5; root mean square error= 0.7206417E-03 Iteration= 6; root mean square error= 0.5975444E-03 Iteration= 7; root mean square error= 0.5181494E-03 Iteration= 8; root mean square error= 0.4638257E-03 Iteration= 9; root mean square error= 0.4249875E-03 Iteration= 10; root mean square error= 0.3961114E-03 Iteration= 11; root mean square error= 0.3738362E-03 Iteration= 12; root mean square error= 0.3560494E-03 Iteration= 13; root mean square error= 0.3413910E-03 Iteration= 14; root mean square error= 0.3289686E-03 Iteration= 15; root mean square error= 0.3181852E-03 Iteration= 16; root mean square error= 0.3086338E-03 Iteration= 17; root mean square error= 0.3000321E-03 Iteration= 18; root mean square error= 0.2921801E-03 Iteration= 19; root mean square error= 0.2849337E-03 Iteration= 20; root mean square error= 0.2781868E-03 Iteration= 21; root mean square error= 0.2718601E-03 Iteration= 22; root mean square error= 0.2658929E-03 Iteration= 23; root mean square error= 0.2602381E-03 Iteration= 24; root mean square error= 0.2548582E-03 Iteration= 25; root mean square error= 0.2497234E-03 Min. area AFTER smoothing= 0.2125940E-01 Max. area AFTER smoothing= 0.3305441E+00 recon 0 RECON2D reconnection loop2d: 1 Recon2d: iteration number= 1 number of flips= 1 RECON: done. filter/1,0,0 FILTER:Use internal epsilonl value = 0.107713032103E-10 FILTER:Dudding 154 duplicate points. rmpoint/compress 154 points removed and 0 elements removed. The new point count is 2183. The new element count is 4188. # Flatten input elevations and interpolate to poly shape read/gmv/lidar_site3.gmv/mo_lidar Reading GMV binary file gmvinputieee no geometry for mesh object geniee finish Nodes: 1000000 Elements: 1996002 Done reading GMV binary total bytes: 542181544 file: lidar_site3.gmv #read / gmv / poly40_topo.gmv / mo_lidar cmo/copyatt/mo_lidar/mo_lidar/z_val/zic CMO_COPYATT WARNING: attribute does not exist: attribute now created: z_val cmo/addatt/mo_lidar/z_val/VDOUBLE/scalar/nnodes/linear/permanent/l finish 1000000 copied from mo_lidar zic to -> mo_lidar z_val cmo/setatt/mo_lidar/zic/1,0,0/0. 1000000 values reset for attribute zic cmo/addatt/mo_merge/z_val/vdouble/scalar/nnodes intrp/continuous/mo_merge z_val/1,0,0/mo_lidar z_val/linear nodes in indexed point set = 2183 INTRP METHOD: continuous FUNCTION: linear TIEBREAKER: maxtie FLAG: plus1 0.610480022430E+01 cmo select mo_lidar finish cmo kdtree build cmo/addatt//v2/INT/scalar/scalar/constant/temporary//2.0 finish cmo/addatt//linkt/VINT/v2/nelements//temporary/x/0.0 finish cmo/addatt//v12/INT/scalar/scalar/constant/temporary//12.0 finish cmo/addatt//sbox/VDOUBLE/v12/nelements/linear/temporary/x/0.0 finish finish Sink point Elems Searched Elements Found Percent Done 45 45 45 2 % 89 89 89 4 % 133 133 133 6 % 177 177 177 8 % 221 221 221 10 % 265 265 265 12 % 309 309 309 14 % 353 353 353 16 % 397 397 397 18 % 441 441 441 20 % 485 485 485 22 % 529 529 529 24 % 573 573 573 26 % 617 617 617 28 % 661 661 661 30 % 705 705 705 32 % 749 749 749 34 % 793 793 793 36 % 837 837 837 38 % 881 881 881 40 % 925 925 925 42 % 969 969 969 44 % 1013 1013 1013 46 % 1057 1057 1057 48 % 1101 1101 1101 50 % 1145 1145 1145 52 % 1189 1189 1189 54 % 1233 1233 1233 56 % 1277 1277 1277 58 % 1321 1321 1321 60 % 1365 1365 1365 62 % 1409 1409 1409 64 % 1453 1453 1453 66 % 1497 1497 1497 68 % 1541 1541 1541 70 % 1585 1585 1585 72 % 1629 1629 1629 74 % 1673 1673 1673 76 % 1717 1717 1717 78 % 1761 1761 1761 80 % 1805 1805 1805 82 % 1849 1849 1849 84 % 1893 1893 1893 86 % 1937 1937 1937 88 % 1981 1981 1981 90 % 2025 2025 2025 92 % 2069 2069 2069 94 % 2113 2113 2113 96 % 2157 2157 2157 98 % 2183 2183 2183 Total cmo select mo_merge finish intrp/continuous done. cmo kdtree release cmo/DELATT/mo_lidar/sbox finish cmo/DELATT/mo_lidar/v12 finish cmo/DELATT/mo_lidar/linkt finish cmo/DELATT/mo_lidar/v2 finish finish cmo/copyatt/mo_merge/mo_merge/zic/z_val 2183 copied from mo_merge z_val to -> mo_merge zic dump/gmv/tmp_poly_massage.gmv/mo_merge cmo/modatt/-def-/-def-/ioflag/x finish cmo/printatt/mo_merge/zic/minmax ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH zic 4.381499767E+00 4.734310412E+00 3.528106447E-01 2183 # REFINE - Depends on Parameter for horizontal spacing (edge lengths) cmo/copy/mo_poly_tri/mo_merge cmo/select/mo_poly_tri cmo/setatt/mo_poly_tri/imt/1 2183 values reset for attribute imt1 cmo/setatt/mo_poly_tri/itetclr/1 4188 values reset for attribute itetclr cmo/setatt/mo_poly_tri/itp/0 2183 values reset for attribute itp1 resetpts/itp geniee finish # each call to refine2d will split edges into smaller triangles refine2d Edge iteration refine2db: 0 0 Edge-refined tris: before= 4188 after= 16752 refine2d Edge iteration refine2db: 0 0 Edge-refined tris: before= 16752 after= 67008 #refine2d dump/tmp_poly_refine.inp/mo_poly_tri dump/avs/tmp_poly_refine.inp/mo_poly_tri cmo/modatt/-def-/-def-/ioflag/x finish finish resetpts itp geniee finish cmo delete mo_poly_tri Released Mesh Object: mo_poly_tri # Extract the edge lines around the polygon read/avs/tmp_poly_refine.inp/mo_poly_tri cmo/addatt/mo_poly_tri/z_val/VDOUBLE/scalar/nnodes/linear/permanent/gxaf/0.0 finish geniee finish cmo/status/brief The current-mesh-object(CMO) is: mo_poly_tri 1 Mesh Object name: mo number of nodes = 16 number of elements = 16 dimensions geometry = 2 element type = tri dimensions topology = 2 3 nodes 3 faces 3 edges boundary flag = 16000000 status = inactive 2 Mesh Object name: mo_inner number of nodes = 9 number of elements = 7 dimensions geometry = 2 element type = tri dimensions topology = 2 3 nodes 3 faces 3 edges boundary flag = 16000000 status = inactive 3 Mesh Object name: mo_merge number of nodes = 2183 number of elements = 4188 dimensions geometry = 2 element type = tri dimensions topology = 2 3 nodes 3 faces 3 edges boundary flag = 16000000 status = inactive 4 Mesh Object name: mo_lidar number of nodes = 1000000 number of elements = 1996002 dimensions geometry = 3 element type = tri dimensions topology = 2 3 nodes 3 faces 3 edges boundary flag = 36000000 status = inactive 5 Mesh Object name: mo_poly_tri number of nodes = 33857 number of elements = 67008 dimensions geometry = 3 element type = tri dimensions topology = 2 3 nodes 3 faces 3 edges boundary flag = 16000000 status = active finish extract/surfmesh/1,0,0/mo_line_work/mo_poly_tri cmo/addatt/mo_line_work/itetclr0/vint/scalar/nelements////0 finish cmo/addatt/mo_line_work/itetclr1/vint/scalar/nelements////0 finish cmo/addatt/mo_line_work/idnode0/vint/scalar/nnodes////0 finish cmo/addatt/mo_line_work/idelem0/vint/scalar/nelements////0 finish cmo/addatt/mo_line_work/idelem1/vint/scalar/nelements////0 finish cmo/addatt/mo_line_work/facecol/vint/scalar/nelements////0 finish cmo/addatt/mo_line_work/idface0/vint/scalar/nelements////0 finish cmo/addatt/mo_line_work/idface1/vint/scalar/nelements////0 finish cmo/addatt/mo_line_work/z_val/VDOUBLE/scalar/nnodes/linear/permanent/gxaf finish cmo/addatt/mo_line_work/ncon50/int/scalar/scalar/constant/permanent/x/0 finish cmo/addatt/mo_line_work/nconbnd/int/scalar/scalar/constant/permanent/x/0 finish cmo/addatt/mo_line_work/icontab/vint/scalar/ncon50/constant/permanent/x/0 finish resetpts/itp geniee finish finish cmo/status/brief The current-mesh-object(CMO) is: mo_line_work 1 Mesh Object name: mo number of nodes = 16 number of elements = 16 dimensions geometry = 2 element type = tri dimensions topology = 2 3 nodes 3 faces 3 edges boundary flag = 16000000 status = inactive 2 Mesh Object name: mo_inner number of nodes = 9 number of elements = 7 dimensions geometry = 2 element type = tri dimensions topology = 2 3 nodes 3 faces 3 edges boundary flag = 16000000 status = inactive 3 Mesh Object name: mo_merge number of nodes = 2183 number of elements = 4188 dimensions geometry = 2 element type = tri dimensions topology = 2 3 nodes 3 faces 3 edges boundary flag = 16000000 status = inactive 4 Mesh Object name: mo_lidar number of nodes = 1000000 number of elements = 1996002 dimensions geometry = 3 element type = tri dimensions topology = 2 3 nodes 3 faces 3 edges boundary flag = 36000000 status = inactive 5 Mesh Object name: mo_poly_tri number of nodes = 33857 number of elements = 67008 dimensions geometry = 3 element type = tri dimensions topology = 2 3 nodes 3 faces 3 edges boundary flag = 16000000 status = inactive 6 Mesh Object name: mo_line_work number of nodes = 704 number of elements = 704 dimensions geometry = 3 element type = lin dimensions topology = 1 2 nodes 2 faces 1 edges boundary flag = 16000000 status = active dump/gmv/tmp_line.gmv/mo_line_work cmo/modatt/-def-/-def-/ioflag/x finish # Create distance attibute from outer polygon edge compute/distance_field/mo_merge/mo_line_work/dfield cmo/addatt/mo_merge/temp_tmp_var/VDOUBLE/scalar/nnodes/linear/temporary//-1.0 finish cmo/addatt/mo_merge/xobj_tmp_var/VDOUBLE/scalar/nnodes/linear/temporary//-1.0 finish cmo/addatt/mo_merge/yobj_tmp_var/VDOUBLE/scalar/nnodes/linear/temporary//-1.0 finish cmo/addatt/mo_merge/zobj_tmp_var/VDOUBLE/scalar/nnodes/linear/temporary//-1.0 finish cmo/addatt/mo_merge/dfield/VDOUBLE/scalar/nnodes/linear/permanent//-1.0 finish interpolate/voronoi/mo_merge/xobj_tmp_var/1 0 0/mo_line_work/xic/keepatt nodes in indexed point set = 2183 cmo/addatt/mo_merge/pt_gtg/VINT/scalar/nnodes//permanent/agfx/0 finish INTRP METHOD: voronoi FUNCTION: linear TIEBREAKER: maxtie FLAG: plus1 0.145539135000E+03 pt_gtg attribute will be added and kept. Sink point Points searched Points Found Percent Done 547 547 547 25 % 1093 1093 1093 50 % 1639 1639 1639 75 % 2183 2183 2183 Total cmo select mo_merge finish intrp/voronoi done. cmo kdtree release finish finish interpolate/voronoi/mo_merge/yobj_tmp_var/1 0 0/mo_line_work/yic/keepatt nodes in indexed point set = 2183 INTRP METHOD: voronoi FUNCTION: linear TIEBREAKER: maxtie FLAG: plus1 0.149574670000E+03 SKIPPING POINT SEARCH... using lookup attribute pt_gtg Sink point Points searched Points Found Percent Done 547 547 547 25 % 1093 1093 1093 50 % 1639 1639 1639 75 % 2183 2183 2183 Total cmo select mo_merge finish intrp/voronoi done. cmo kdtree release finish finish interpolate/voronoi/mo_merge/zobj_tmp_var/1 0 0/mo_line_work/zic nodes in indexed point set = 2183 INTRP METHOD: voronoi FUNCTION: linear TIEBREAKER: maxtie FLAG: plus1 0.565469289831E+01 SKIPPING POINT SEARCH... using lookup attribute pt_gtg Sink point Points searched Points Found Percent Done 547 547 547 25 % 1093 1093 1093 50 % 1639 1639 1639 75 % 2183 2183 2183 Total cmo DELATT mo_merge/pt_gtg finish cmo select mo_merge finish intrp/voronoi done. cmo kdtree release finish finish math/sub/mo_merge/temp_tmp_var/1,0,0/mo_merge/xic/mo_merge/xobj_tmp_var finish math/power/mo_merge/dfield/1,0,0/mo_merge/temp_tmp_var/2.0 finish math/sub/mo_merge/temp_tmp_var/1,0,0/mo_merge/yic/mo_merge/yobj_tmp_var finish math/power/mo_merge/temp_tmp_var/1,0,0/mo_merge/temp_tmp_var/2.0 finish math/add/mo_merge/dfield/1 0 0/mo_merge/dfield/mo_merge/temp_tmp_var finish math/sub/mo_merge/temp_tmp_var/1,0,0/mo_merge/zic/mo_merge/zobj_tmp_var finish math/power/mo_merge/temp_tmp_var/1,0,0/mo_merge/temp_tmp_var/2.0 finish math/add/mo_merge/dfield/1 0 0/mo_merge/dfield/mo_merge/temp_tmp_var finish math/power/mo_merge/dfield/1,0,0/mo_merge/dfield/0.5 finish cmo/delatt/mo_merge/temp_tmp_var finish cmo/delatt/mo_merge/xobj_tmp_var finish cmo/delatt/mo_merge/yobj_tmp_var finish cmo/delatt/mo_merge/zobj_tmp_var finish COMPUTE distance_field: Done. cmo printatt mo_merge dfield minmax ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH dfield 1.445849715E-14 9.911063668E+00 9.911063668E+00 2183 dump/tmp_poly_refine.gmv/mo_merge dump/gmv/tmp_poly_refine.gmv/mo_merge cmo/modatt/-def-/-def-/ioflag/x finish finish # Sizes for ice wedge used in math routines to define shape define/HALF_WIDTH/1.4 define/DEPTH/2.8 define/SLOPE/2.0 define/L1/0.02 define/L2/0.24 define/L3/0.1 define/L4/0.05 # Create top and bottom layer of ice wedge mesh object # top of ice wedge is translated to requested layer # to put top under a specified layer # These distances are same as defined for full mesh define/L1/0.02 # this is all material 2 define/L2/0.24 # this is all material 3 define/L3/0.24 # this is all material 4 define/L4/0.5 # this is all material 5 define/L5/5 # subtract each previous layer until ice top is positioned correctly # CHECK Starting elevation is same as top elevation cmo printatt mo_merge zic minmax ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH zic 4.381499767E+00 4.734310412E+00 3.528106447E-01 2183 math/sub/mo_merge/zic/1,0,0/mo_merge/zic/L1 #dump / avs / tmp_poly40_1.inp / mo_merge math/sub/mo_merge/zic/1,0,0/mo_merge/zic/L2 #dump / avs / tmp_poly40_2.inp / mo_merge math/sub/mo_merge/zic/1,0,0/mo_merge/zic/L3 #dump / avs / tmp_poly40_3.inp / mo_merge #math / sub / mo_merge / zic / 1,0,0 / mo_merge / zic / L4 #math / sub / mo_merge / zic / 1,0,0 / mo_merge / zic / L5 #dump / avs / tmp_poly40_4.inp / mo_merge # top file for stack routine dump/avs/tmp_top.inp/mo_merge cmo/modatt/-def-/-def-/ioflag/x finish # CHECK top elevation should be located at intended layer cmo printatt mo_merge zic minmax ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH zic 3.881499767E+00 4.234310412E+00 3.528106447E-01 2183 # math routines to pinchout ice wedge at the bottom cmo/addatt/mo_merge/z_ele/vdouble/scalar/nnodes cmo/copyatt/mo_merge/mo_merge/z_save/zic CMO_COPYATT WARNING: attribute does not exist: attribute now created: z_save cmo/addatt/mo_merge/z_save/VDOUBLE/scalar/nnodes/linear/permanent/l finish 2183 copied from mo_merge zic to -> mo_merge z_save cmo/printatt/mo_merge/zic/minmax ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH zic 3.881499767E+00 4.234310412E+00 3.528106447E-01 2183 cmo/addatt/mo_merge/depth/vdouble/scalar/nnodes cmo/setatt/mo_merge/depth/1,0,0/DEPTH 2183 values reset for attribute depth cmo/printatt/mo_merge/imt/minmax ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH imt1 1 2 1 2183 #eltset / e2 / itetclr / eq / 2 #pset / p2 / eltset / e2 #cmo / setatt / mo_merge / imt / pset get p2 / 2 pset/pdown/attribute/dfield/1,0,0/lt/HALF_WIDTH THE PSET pdown HAS 572 POINTS #pset / pdown / attribute / imt / 1,0,0 / eq / 1 #math / divide / mo_merge / slope / pset get pdown / mo_merge / depth / mo_merge math/multiply/mo_merge/z_ele/pset get pdown/mo_merge/dfield/SLOPE math/sub/mo_merge/z_ele/pset get pdown/mo_merge/z_ele/DEPTH math/add/mo_merge/z_ele/1,0,0/mo_merge/z_ele/mo_merge/z_save #math / ceiling / mo_merge / z_save / 1,0,0 / mo_merge / z_save / 0. cmo/printatt/mo_merge/z_ele/minmax ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH z_ele 1.081499767E+00 4.234310412E+00 3.152810645E+00 2183 cmo/copyatt/mo_merge/mo_merge/zic/z_ele 2183 copied from mo_merge z_ele to -> mo_merge zic # bottom layer for stack routine dump/avs/tmp_bot.inp/mo_merge cmo/modatt/-def-/-def-/ioflag/x finish # CHECK bottom elevation should be DEPTH below top layer cmo printatt mo_merge zic minmax ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH zic 1.081499767E+00 4.234310412E+00 3.152810645E+00 2183 # stack top and bottom into a prism mesh ice wedge object cmo/create/mo_stack stack/layers/avs/tmp_bot.inp 1 tmp_top.inp 1/flip Layers to create: 2 Max material number: 1 Reading 2 surface files... cmo/addatt/mo_stack/layertyp/VINT/scalar/nnodes//permanent/agfx finish Read first surface: tmp_bot.inp read/avs/tmp_bot.inp/def1 cmo/addatt/def1/z_val/VDOUBLE/scalar/nnodes/linear/permanent/gxaf/0.0 finish cmo/addatt/def1/dfield/VDOUBLE/scalar/nnodes/linear/permanent/gxaf/0.0 finish cmo/addatt/def1/z_ele/VDOUBLE/scalar/nnodes/linear/permanent/gxaf/0.0 finish cmo/addatt/def1/depth/VDOUBLE/scalar/nnodes/linear/permanent/gxaf/0.0 finish geniee finish cmo/status/brief The current-mesh-object(CMO) is: def1 1 Mesh Object name: mo number of nodes = 16 number of elements = 16 dimensions geometry = 2 element type = tri dimensions topology = 2 3 nodes 3 faces 3 edges boundary flag = 16000000 status = inactive 2 Mesh Object name: mo_inner number of nodes = 9 number of elements = 7 dimensions geometry = 2 element type = tri dimensions topology = 2 3 nodes 3 faces 3 edges boundary flag = 16000000 status = inactive 3 Mesh Object name: mo_merge number of nodes = 2183 number of elements = 4188 dimensions geometry = 2 element type = tri dimensions topology = 2 3 nodes 3 faces 3 edges boundary flag = 16000000 status = inactive 4 Mesh Object name: mo_lidar number of nodes = 1000000 number of elements = 1996002 dimensions geometry = 3 element type = tri dimensions topology = 2 3 nodes 3 faces 3 edges boundary flag = 36000000 status = inactive 5 Mesh Object name: mo_poly_tri number of nodes = 33857 number of elements = 67008 dimensions geometry = 3 element type = tri dimensions topology = 2 3 nodes 3 faces 3 edges boundary flag = 16000000 status = inactive 6 Mesh Object name: mo_line_work number of nodes = 704 number of elements = 704 dimensions geometry = 3 element type = lin dimensions topology = 1 2 nodes 2 faces 1 edges boundary flag = 16000000 status = inactive 7 Mesh Object name: mo_stack number of nodes = 0 number of elements = 0 dimensions geometry = 3 element type = tet dimensions topology = 3 4 nodes 4 faces 6 edges boundary flag = 16000000 status = inactive 8 Mesh Object name: def1 number of nodes = 2183 number of elements = 4188 dimensions geometry = 3 element type = tri dimensions topology = 2 3 nodes 3 faces 3 edges boundary flag = 16000000 status = active finish finish 4188 of 4188 downward elements flipped up Read surface: tmp_top.inp read/avs/tmp_top.inp/def2 cmo/addatt/def2/z_val/VDOUBLE/scalar/nnodes/linear/permanent/gxaf/0.0 finish cmo/addatt/def2/dfield/VDOUBLE/scalar/nnodes/linear/permanent/gxaf/0.0 finish geniee finish cmo/status/brief The current-mesh-object(CMO) is: def2 1 Mesh Object name: mo number of nodes = 16 number of elements = 16 dimensions geometry = 2 element type = tri dimensions topology = 2 3 nodes 3 faces 3 edges boundary flag = 16000000 status = inactive 2 Mesh Object name: mo_inner number of nodes = 9 number of elements = 7 dimensions geometry = 2 element type = tri dimensions topology = 2 3 nodes 3 faces 3 edges boundary flag = 16000000 status = inactive 3 Mesh Object name: mo_merge number of nodes = 2183 number of elements = 4188 dimensions geometry = 2 element type = tri dimensions topology = 2 3 nodes 3 faces 3 edges boundary flag = 16000000 status = inactive 4 Mesh Object name: mo_lidar number of nodes = 1000000 number of elements = 1996002 dimensions geometry = 3 element type = tri dimensions topology = 2 3 nodes 3 faces 3 edges boundary flag = 36000000 status = inactive 5 Mesh Object name: mo_poly_tri number of nodes = 33857 number of elements = 67008 dimensions geometry = 3 element type = tri dimensions topology = 2 3 nodes 3 faces 3 edges boundary flag = 16000000 status = inactive 6 Mesh Object name: mo_line_work number of nodes = 704 number of elements = 704 dimensions geometry = 3 element type = lin dimensions topology = 1 2 nodes 2 faces 1 edges boundary flag = 16000000 status = inactive 7 Mesh Object name: mo_stack number of nodes = 0 number of elements = 0 dimensions geometry = 3 element type = tri dimensions topology = 2 3 nodes 3 faces 3 edges boundary flag = 16000000 status = inactive 8 Mesh Object name: def1 number of nodes = 2183 number of elements = 4188 dimensions geometry = 3 element type = tri dimensions topology = 2 3 nodes 3 faces 3 edges boundary flag = 16000000 status = inactive 9 Mesh Object name: cmoprev number of nodes = 2183 number of elements = 4188 dimensions geometry = 3 element type = tri dimensions topology = 2 3 nodes 3 faces 3 edges boundary flag = 16000000 status = inactive 10 Mesh Object name: def2 number of nodes = 2183 number of elements = 4188 dimensions geometry = 3 element type = tri dimensions topology = 2 3 nodes 3 faces 3 edges boundary flag = 16000000 status = active finish finish cmo/delete/def2/ Released Mesh Object: def2 finish 4188 of 4188 downward elements flipped up Total equal elevations: 1611 of 2183 Layers compared: tmp_bot.inp and tmp_top.inp cmo/delete/cmoprev/ Released Mesh Object: cmoprev finish cmo/copy/cmoprev/cmonxt/ finish addmesh/merge/mo_stack/def1/cmonxt/ searching for mo_stack attributes in cmonxt attribute not found: 66 layertyp ADDMESH: 1 attributes not copied from cmonxt to mo_stack finish warning: Column has 0 height. column, node: 8 2191 Warning: Column has 0 height: 8 warning: Column has 0 height. column, node: 9 2192 Warning: Column has 0 height: 9 warning: Column has 0 height. column, node: 10 2193 Warning: Column has 0 height: 10 warning: Column has 0 height. column, node: 11 2194 Warning: Column has 0 height: 11 warning: Column has 0 height. column, node: 12 2195 Warning: Column has 0 height: 12 warning: Column has 0 height. column, node: 13 2196 Warning: Column has 0 height: 13 warning: Column has 0 height. column, node: 14 2197 Warning: Column has 0 height: 14 warning: Column has 0 height. column, node: 15 2198 Warning: Column has 0 height: 15 warning: Column has 0 height. column, node: 16 2199 Warning: Column has 0 height: 16 warning: Column has 0 height. column, node: 17 2200 Total warnings for flat columns: 3222 ................................................................ surface name layer color type llcorner zic tmp_bot.inp 1 1 -1 1) 1.114700E+00 tmp_top.inp 2 1 -2 2184) 3.914700E+00 Elements per layer: 4188 stacked total: 8376 Nodes per layer: 2183 stacked total: 4366 No Truncating layer specified. files read: 2 from total: 2 ................................................................ cmo/addatt/mo_stack/nlayers/INT/scalar/scalar/constant// finish cmo/addatt/mo_stack/nnperlayer/INT/scalar/scalar/constant// finish cmo/addatt/mo_stack/neperlayer/INT/scalar/scalar/constant// finish cmo/setatt/mo_stack/nlayers 2 1 values reset for attribute nlayers finish cmo/setatt/mo_stack/nnperlayer 2183 1 values reset for attribute nnperlayer finish cmo/setatt/mo_stack/neperlayer 4188 1 values reset for attribute neperlayer finish Released Mesh Object: def1 Released Mesh Object: cmonxt Released Mesh Object: cmoprev stack done. stack/fill/mo_prism/mo_stack cmo/copy/mo_prism/mo_stack/ finish geniee finish filter/1,0,0 FILTER:Use internal epsilonl value = 0.999999993923E-08 FILTER:Dudding 1611 duplicate points. rmpoint/compress 1611 points removed and 0 elements removed. The new point count is 2755. The new element count is 4188. rmpoint/zero_volume/1.e-12 Dudded 0 points plus parents, plus 3066 elements removed. # Final output file to be used later dump/gmv/tmp_ice_wedge.gmv/mo_prism cmo/modatt/-def-/-def-/ioflag/x finish ############# CHECK FINAL WORK ##################### # write surface triangulation sizes cmo select mo_poly_tri quality epsilonl, epsilonaspect: 1.0771588E-11 1.2497962E-33 -------------------------------------------- elements with aspect ratio < .01: 0 elements with aspect ratio b/w .01 and .02: 0 elements with aspect ratio b/w .02 and .05: 0 elements with aspect ratio b/w .05 and .1 : 0 elements with aspect ratio b/w .1 and .2 : 32 elements with aspect ratio b/w .2 and .5 : 64 elements with aspect ratio b/w .5 and 1. : 66912 min aspect ratio = 0.1428E+00 max aspect ratio = 0.1000E+01 epsilonvol: 8.1593473E-11 --------------------------------------- element volumes b/w 0.1330E-02 and 0.2303E-02: 64 element volumes b/w 0.2303E-02 and 0.3988E-02: 32 element volumes b/w 0.3988E-02 and 0.6907E-02: 8080 element volumes b/w 0.6907E-02 and 0.1196E-01: 48768 element volumes b/w 0.1196E-01 and 0.2072E-01: 10064 min volume = 1.3295367E-03 max volume = 2.0718972E-02 ----------------------------------------------------------- 67008 total elements evaluated. quality/edge_min epsilonl, epsilonaspect: 1.0771588E-11 1.2497962E-33 -------------------------------------------- element norm min edge length < .01: 0 element norm edge length b/w .01 and .02: 0 element norm edge length b/w .02 and .05: 0 element norm edge length b/w .05 and .1 : 32 element norm edge length b/w .1 and .2 : 32 element norm edge length b/w .2 and .5 : 6592 element norm edge length b/w .5 and 1. : 60352 min min edge length = 0.1600E-01 max min edge length = 0.1995E+00 ----------------------------------------------------------- 67008 total elements evaluated. quality/edge_max epsilonl, epsilonaspect: 1.0771588E-11 1.2497962E-33 -------------------------------------------- element norm max edge length < .01: 0 element norm edge length b/w .01 and .02: 0 element norm edge length b/w .02 and .05: 0 element norm edge length b/w .05 and .1 : 0 element norm edge length b/w .1 and .2 : 0 element norm edge length b/w .2 and .5 : 528 element norm edge length b/w .5 and 1. : 66480 min max edge length = 0.9991E-01 max max edge length = 0.2502E+00 ----------------------------------------------------------- 67008 total elements evaluated. # coordinates for polygon surface (top elevations) cmo printatt mo_poly_tri -xyz- minmax ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH xic 1.030189210E+02 1.445391350E+02 4.152021400E+01 33857 yic 1.234896800E+02 1.485746700E+02 2.508499000E+01 33857 zic 4.381499767E+00 4.734310412E+00 3.528106447E-01 33857 # coordinates for stacked ice mesh object cmo printatt mo_prism -xyz- minmax ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH xic 1.030189210E+02 1.445391350E+02 4.152021400E+01 2755 yic 1.234896800E+02 1.485746700E+02 2.508499000E+01 2755 zic 1.081499767E+00 4.203951535E+00 3.122451767E+00 2755 # there should be one material for ice mesh object cmo printatt mo_prism imt minmax ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH imt1 1 1 0 2755 cmo printatt mo_prism itetclr minmax ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH itetclr 1 1 0 1122 ################################################## # DONE # file tmp_poly40_refine.inp surface for full mesh layers # file tmp_ice_wedge.gmv 3D ice wedge mesh for interpolation finish LaGriT successfully completed