* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * Program: LaGriT V3.101 Linux m64 * * * date_compile: 2013/04/25 RH gf4.5 * * * Run Time: 2013/Oct 15 12:42:32 * * * 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 2D cylinder point set and create Delaunay 2D triangulation define WELL_ID 1 define MAT_ID 2 read avs cylinder_rim.inp mo_rim cmo/addatt/mo_rim/idnode0/VINT/scalar/nnodes/linear/permanent/gxaf/0.0 finish cmo/status/brief The current-mesh-object(CMO) is: mo_rim 1 Mesh Object name: mo_rim number of nodes = 92 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 read avs cylinder_well.inp mo_well cmo/status/brief The current-mesh-object(CMO) is: mo_well 1 Mesh Object name: mo_rim number of nodes = 92 number of elements = 0 dimensions geometry = 3 element type = tet dimensions topology = 3 4 nodes 4 faces 6 edges boundary flag = 16000000 status = inactive 2 Mesh Object name: mo_well number of nodes = 15 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 setatt mo_rim imt MAT_ID Warning: nelements = 0 92 values reset for attribute imt1 cmo setatt mo_well imt WELL_ID Warning: nelements = 0 15 values reset for attribute imt1 # merge the points addmesh/merge/mo/mo_rim/mo_well searching for mo attributes in mo_well attribute not found: 66 idnode0 ADDMESH: 1 attributes not copied from mo_well to mo dump gmv tmp_merge_points.gmv mo cmo/modatt/-def-/-def-/ioflag/x finish # add points to 2D mesh object cmo create cmotri///triplane copypts/cmotri/mo cmo/addatt/cmotri/idnode0/VINT/scalar/nnodes/linear/permanent/gxaf/0.0 finish copypts added nodes 1 to 107 cmo setatt cmotri imt MAT_ID Warning: nelements = 0 107 values reset for attribute imt1 #################################### # CONNECT into triangulation cmo select cmotri cmo/setatt/cmotri/imt 1 Warning: nelements = 0 107 values reset for attribute imt1 connect Coordinates of enclosing triangle are: -0.67500D+02-0.63000D+02-0.00000D+00 0.67500D+02-0.67500D+02-0.00000D+00 0.00000D+00 0.13500D+03 0.00000D+00 Successfully eliminated all multimaterial connections. The mesh is now complete! LaGriT FINISH: connect cmo/setatt/cmotri/itetclr 1 120 values reset for attribute itetclr rmpoint compress 0 points removed and 0 elements removed. The new point count is 107. The new element count is 120. resetpts itp geniee finish # set distance from center for well point selections compute/distance_field/cmotri mo_well/wdist cmo/addatt/cmotri/temp_tmp_var/VDOUBLE/scalar/nnodes/linear/temporary//-1.0 finish cmo/addatt/cmotri/xobj_tmp_var/VDOUBLE/scalar/nnodes/linear/temporary//-1.0 finish cmo/addatt/cmotri/yobj_tmp_var/VDOUBLE/scalar/nnodes/linear/temporary//-1.0 finish cmo/addatt/cmotri/zobj_tmp_var/VDOUBLE/scalar/nnodes/linear/temporary//-1.0 finish cmo/addatt/cmotri/wdist/VDOUBLE/scalar/nnodes/linear/permanent//-1.0 finish interpolate/voronoi/cmotri/xobj_tmp_var/1 0 0/mo_well/xic/keepatt nodes in indexed point set = 107 cmo/addatt/cmotri/pt_gtg/VINT/scalar/nnodes//permanent/agfx/0 finish INTRP METHOD: voronoi FUNCTION: linear TIEBREAKER: maxtie FLAG: plus1 0.300000000000E+01 pt_gtg attribute will be added and kept. Sink point Points searched Points Found Percent Done 28 29 28 25 % 55 56 55 50 % 82 83 82 75 % 107 108 107 Total cmo select cmotri finish intrp/voronoi done. cmo kdtree release finish finish interpolate/voronoi/cmotri/yobj_tmp_var/1 0 0/mo_well/yic/keepatt nodes in indexed point set = 107 INTRP METHOD: voronoi FUNCTION: linear TIEBREAKER: maxtie FLAG: plus1 0.294985582436E+01 SKIPPING POINT SEARCH... using lookup attribute pt_gtg Sink point Points searched Points Found Percent Done 28 28 28 25 % 55 55 55 50 % 82 82 82 75 % 107 107 107 Total cmo select cmotri finish intrp/voronoi done. cmo kdtree release finish finish interpolate/voronoi/cmotri/zobj_tmp_var/1 0 0/mo_well/zic nodes in indexed point set = 107 INTRP METHOD: voronoi FUNCTION: linear TIEBREAKER: maxtie FLAG: plus1 0.100000000000E+01 SKIPPING POINT SEARCH... using lookup attribute pt_gtg Sink point Points searched Points Found Percent Done 28 28 28 25 % 55 55 55 50 % 82 82 82 75 % 107 107 107 Total cmo DELATT cmotri/pt_gtg finish cmo select cmotri finish intrp/voronoi done. cmo kdtree release finish finish math/sub/cmotri/temp_tmp_var/1,0,0/cmotri/xic/cmotri/xobj_tmp_var finish math/power/cmotri/wdist/1,0,0/cmotri/temp_tmp_var/2.0 finish math/sub/cmotri/temp_tmp_var/1,0,0/cmotri/yic/cmotri/yobj_tmp_var finish math/power/cmotri/temp_tmp_var/1,0,0/cmotri/temp_tmp_var/2.0 finish math/add/cmotri/wdist/1 0 0/cmotri/wdist/cmotri/temp_tmp_var finish math/sub/cmotri/temp_tmp_var/1,0,0/cmotri/zic/cmotri/zobj_tmp_var finish math/power/cmotri/temp_tmp_var/1,0,0/cmotri/temp_tmp_var/2.0 finish math/add/cmotri/wdist/1 0 0/cmotri/wdist/cmotri/temp_tmp_var finish math/power/cmotri/wdist/1,0,0/cmotri/wdist/0.5 finish cmo/delatt/cmotri/temp_tmp_var finish cmo/delatt/cmotri/xobj_tmp_var finish cmo/delatt/cmotri/yobj_tmp_var finish cmo/delatt/cmotri/zobj_tmp_var finish COMPUTE distance_field: Done. dump/gmv/tmp_tri1.gmv/cmotri cmo/modatt/-def-/-def-/ioflag/x finish dump/avs/tmp_tri1.inp/cmotri cmo/modatt/-def-/-def-/ioflag/x finish #################################### # Optimize with refine and smooth # Refine by halves, max edge starts near 10. define REFINE1 5.0 refine/rivara///edge/1 0 0/REFINE1///inclusive Face iteration, number of faces: 1 7 Face iteration, number of faces: 1 14 Face iteration, number of faces: 1 21 Face iteration, number of faces: 1 35 Face iteration, number of faces: 1 35 Face iteration, number of faces: 1 46 Face iteration, number of faces: 1 50 Face iteration, number of faces: 1 54 Face iteration, number of faces: 1 44 Face iteration, number of faces: 1 38 Face iteration, number of faces: 1 24 Face iteration, number of faces: 1 21 Face iteration, number of faces: 1 7 recon/0/0.94205548E-11/ RECON2D reconnection loop2d: 1 RECON: done. finish dump/gmv/tmp_tri2.gmv/cmotri cmo/modatt/-def-/-def-/ioflag/x finish dump/avs/tmp_tri2.inp/cmotri cmo/modatt/-def-/-def-/ioflag/x finish define REFINE1 2.0 refine/rivara///edge/1 0 0/REFINE1///inclusive Face iteration, number of faces: 1 170 Face iteration, number of faces: 1 204 Face iteration, number of faces: 1 242 Face iteration, number of faces: 1 108 Face iteration, number of faces: 1 62 Face iteration, number of faces: 1 37 Face iteration, number of faces: 1 28 Face iteration, number of faces: 1 22 Face iteration, number of faces: 1 16 Face iteration, number of faces: 1 9 recon/0/0.94205548E-11/ RECON2D reconnection loop2d: 1 Recon2d: iteration number= 1 number of flips= 14 Recon2d: iteration number= 2 number of flips= 7 Recon2d: iteration number= 3 number of flips= 7 Recon2d: iteration number= 4 number of flips= 7 Recon2d: iteration number= 5 number of flips= 7 Recon2d: iteration number= 6 number of flips= 7 RECON: done. finish dump/gmv/tmp_tri3.gmv/cmotri cmo/modatt/-def-/-def-/ioflag/x finish dump/avs/tmp_tri3.inp/cmotri cmo/modatt/-def-/-def-/ioflag/x finish # Smooth into nicer point distribution # the more times this is repeated, the nicer it is # Select points OUTSIDE the well for movement pset/psmooth/attribute wdist/1,0,0/gt .01 THE PSET psmooth HAS 1386 POINTS smooth/position/esug/pset get psmooth SMOOTH: position option nodes in point set = 1386 SMOOTH: Smooth 2D triangular mesh using ESUG Min. area BEFORE smoothing= 0.5384029E-01 Max. area BEFORE smoothing= 0.8947538E+00 Iteration= 1; root mean square error= 0.1693688E+00 Iteration= 2; root mean square error= 0.1023177E+00 Iteration= 3; root mean square error= 0.8810857E-01 Iteration= 4; root mean square error= 0.8270512E-01 Iteration= 5; root mean square error= 0.7759215E-01 Iteration= 6; root mean square error= 0.7310375E-01 Iteration= 7; root mean square error= 0.6831243E-01 Iteration= 8; root mean square error= 0.6386181E-01 Iteration= 9; root mean square error= 0.5924348E-01 Iteration= 10; root mean square error= 0.5475126E-01 Iteration= 11; root mean square error= 0.5053952E-01 Iteration= 12; root mean square error= 0.4658796E-01 Iteration= 13; root mean square error= 0.4290441E-01 Iteration= 14; root mean square error= 0.3947992E-01 Iteration= 15; root mean square error= 0.3630443E-01 Iteration= 16; root mean square error= 0.3336642E-01 Iteration= 17; root mean square error= 0.3065337E-01 Iteration= 18; root mean square error= 0.2815206E-01 Iteration= 19; root mean square error= 0.2584897E-01 Iteration= 20; root mean square error= 0.2373059E-01 Iteration= 21; root mean square error= 0.2178370E-01 Iteration= 22; root mean square error= 0.1999555E-01 Iteration= 23; root mean square error= 0.1835400E-01 Iteration= 24; root mean square error= 0.1684756E-01 Iteration= 25; root mean square error= 0.1546552E-01 Min. area AFTER smoothing= 0.1023479E+00 Max. area AFTER smoothing= 0.1562668E+01 recon/0 RECON2D reconnection loop2d: 1 Recon2d: iteration number= 1 number of flips= 697 Recon2d: iteration number= 2 number of flips= 522 Recon2d: iteration number= 3 number of flips= 237 Recon2d: iteration number= 4 number of flips= 101 Recon2d: iteration number= 5 number of flips= 33 Recon2d: iteration number= 6 number of flips= 22 Recon2d: iteration number= 7 number of flips= 4 Recon2d: iteration number= 8 number of flips= 2 Recon2d: iteration number= 9 number of flips= 2 RECON: done. rmpoint compress 0 points removed and 0 elements removed. The new point count is 1401. The new element count is 2708. resetpts itp geniee finish dump gmv tmp_smooth1.gmv cmotri cmo/modatt/-def-/-def-/ioflag/x finish pset psmooth delete PSET psmooth DELETED pset/psmooth/attribute wdist/1,0,0/gt .01 THE PSET psmooth HAS 1386 POINTS smooth/position/esug/pset get psmooth SMOOTH: position option nodes in point set = 1386 SMOOTH: Smooth 2D triangular mesh using ESUG Min. area BEFORE smoothing= 0.9857601E-02 Max. area BEFORE smoothing= 0.4095568E+01 Iteration= 1; root mean square error= 0.2252847E+00 Iteration= 2; root mean square error= 0.1294319E+00 Iteration= 3; root mean square error= 0.7591714E-01 Iteration= 4; root mean square error= 0.5312084E-01 Iteration= 5; root mean square error= 0.3789589E-01 Iteration= 6; root mean square error= 0.2875664E-01 Iteration= 7; root mean square error= 0.2307399E-01 Iteration= 8; root mean square error= 0.1912089E-01 Iteration= 9; root mean square error= 0.1614360E-01 Iteration= 10; root mean square error= 0.1379940E-01 Iteration= 11; root mean square error= 0.1190542E-01 Iteration= 12; root mean square error= 0.1035033E-01 Iteration= 13; root mean square error= 0.9058881E-02 Iteration= 14; root mean square error= 0.7977110E-02 Iteration= 15; root mean square error= 0.7064631E-02 Iteration= 16; root mean square error= 0.6290363E-02 Iteration= 17; root mean square error= 0.5629869E-02 Iteration= 18; root mean square error= 0.5063631E-02 Iteration= 19; root mean square error= 0.4575882E-02 Iteration= 20; root mean square error= 0.4153775E-02 Iteration= 21; root mean square error= 0.3786787E-02 Iteration= 22; root mean square error= 0.3466253E-02 Iteration= 23; root mean square error= 0.3185019E-02 Iteration= 24; root mean square error= 0.2937158E-02 Iteration= 25; root mean square error= 0.2717747E-02 Min. area AFTER smoothing= 0.1283839E+00 Max. area AFTER smoothing= 0.7818315E+00 recon/0 RECON2D reconnection loop2d: 1 Recon2d: iteration number= 1 number of flips= 436 Recon2d: iteration number= 2 number of flips= 49 Recon2d: iteration number= 3 number of flips= 10 RECON: done. rmpoint compress 0 points removed and 0 elements removed. The new point count is 1401. The new element count is 2708. resetpts itp geniee finish dump gmv tmp_smooth2.gmv cmotri cmo/modatt/-def-/-def-/ioflag/x finish * smooth entire mesh * do not smooth too close to wells pset psmooth delete PSET psmooth DELETED pset/psmooth/attribute wdist/1,0,0/gt .01 THE PSET psmooth HAS 1386 POINTS smooth/position/esug/pset,get,psmooth SMOOTH: position option nodes in point set = 1386 SMOOTH: Smooth 2D triangular mesh using ESUG Min. area BEFORE smoothing= 0.4847478E-01 Max. area BEFORE smoothing= 0.7710235E+00 Iteration= 1; root mean square error= 0.1089281E+00 Iteration= 2; root mean square error= 0.4605156E-01 Iteration= 3; root mean square error= 0.2838809E-01 Iteration= 4; root mean square error= 0.2252253E-01 Iteration= 5; root mean square error= 0.1948290E-01 Iteration= 6; root mean square error= 0.1728739E-01 Iteration= 7; root mean square error= 0.1557155E-01 Iteration= 8; root mean square error= 0.1416374E-01 Iteration= 9; root mean square error= 0.1298329E-01 Iteration= 10; root mean square error= 0.1197812E-01 Iteration= 11; root mean square error= 0.1111325E-01 Iteration= 12; root mean square error= 0.1036267E-01 Iteration= 13; root mean square error= 0.9706251E-02 Iteration= 14; root mean square error= 0.9127906E-02 Iteration= 15; root mean square error= 0.8614570E-02 Iteration= 16; root mean square error= 0.8155558E-02 Iteration= 17; root mean square error= 0.7742109E-02 Iteration= 18; root mean square error= 0.7367031E-02 Iteration= 19; root mean square error= 0.7024424E-02 Iteration= 20; root mean square error= 0.6709454E-02 Iteration= 21; root mean square error= 0.6418155E-02 Iteration= 22; root mean square error= 0.6147277E-02 Iteration= 23; root mean square error= 0.5894153E-02 Iteration= 24; root mean square error= 0.5656587E-02 Iteration= 25; root mean square error= 0.5432768E-02 Min. area AFTER smoothing= 0.1334647E+00 Max. area AFTER smoothing= 0.6151384E+00 recon/0 RECON2D reconnection loop2d: 1 Recon2d: iteration number= 1 number of flips= 66 RECON: done. rmpoint compress 0 points removed and 0 elements removed. The new point count is 1401. The new element count is 2708. resetpts itp geniee finish dump gmv tmp_smooth3.gmv cmotri cmo/modatt/-def-/-def-/ioflag/x finish # merge thin triangles massage 5. .8/0.1/1,0,0 nodes in point set = 1401 rmpoint/element Dudded 0 points plus parents, and 614 elements removed. finish geniee finish Merged 307 nodes. rmpoint/element Dudded 0 points plus parents, and 468 elements removed. finish geniee finish Merged 234 nodes. rmpoint/element Dudded 0 points plus parents, and 340 elements removed. finish geniee finish Merged 170 nodes. rmpoint/element Dudded 0 points plus parents, and 212 elements removed. finish geniee finish Merged 106 nodes. rmpoint/element Dudded 0 points plus parents, and 120 elements removed. finish geniee finish Merged 60 nodes. rmpoint/element Dudded 0 points plus parents, and 42 elements removed. finish geniee finish Merged 21 nodes. rmpoint/element Dudded 0 points plus parents, and 4 elements removed. finish geniee finish Merged 2 nodes. recon/0/0.10000000E+00/ RECON2D reconnection loop2d: 1 Recon2d: iteration number= 1 number of flips= 163 Recon2d: iteration number= 2 number of flips= 39 Recon2d: iteration number= 3 number of flips= 11 Recon2d: iteration number= 4 number of flips= 2 RECON: done. finish nodes in point set = 501 SGD: 182 total relaxations on point set of 501 indept. nodes. rmpoint/compress 900 points removed and 0 elements removed. The new point count is 501. The new element count is 908. finish nodes in point set = 501 rmpoint/element Dudded 0 points plus parents, and 10 elements removed. finish geniee finish Merged 5 nodes. recon/0/0.10000000E+00/ RECON2D reconnection loop2d: 1 Recon2d: iteration number= 1 number of flips= 63 Recon2d: iteration number= 2 number of flips= 1 RECON: done. finish nodes in point set = 496 SGD: 53 total relaxations on point set of 496 indept. nodes. rmpoint/compress 5 points removed and 0 elements removed. The new point count is 496. The new element count is 898. finish nodes in point set = 496 rmpoint/element Dudded 0 points plus parents, and 4 elements removed. finish geniee finish Merged 2 nodes. recon/0/0.10000000E+00/ RECON2D reconnection loop2d: 1 Recon2d: iteration number= 1 number of flips= 16 Recon2d: iteration number= 2 number of flips= 1 RECON: done. finish nodes in point set = 494 SGD: 13 total relaxations on point set of 494 indept. nodes. rmpoint/compress 2 points removed and 0 elements removed. The new point count is 494. The new element count is 894. finish nodes in point set = 494 recon/0/0.10000000E+00/ RECON2D reconnection loop2d: 1 Recon2d: iteration number= 1 number of flips= 1 RECON: done. finish nodes in point set = 494 SGD: 0 total relaxations on point set of 494 indept. nodes. rmpoint/compress 0 points removed and 0 elements removed. The new point count is 494. The new element count is 894. finish MASSAGE: hit max. iteration limit. nodes in point set = 494 rmpoint compress 0 points removed and 0 elements removed. The new point count is 494. The new element count is 894. resetpts itp geniee finish dump gmv tmp_mass4.gmv cmotri cmo/modatt/-def-/-def-/ioflag/x finish * do not smooth too close to wells pset psmooth delete PSET psmooth DELETED pset/psmooth/attribute wdist/1,0,0/gt .01 THE PSET psmooth HAS 480 POINTS smooth/position/esug/pset,get,psmooth SMOOTH: position option nodes in point set = 480 SMOOTH: Smooth 2D triangular mesh using ESUG Min. area BEFORE smoothing= 0.3516137E+00 Max. area BEFORE smoothing= 0.2171248E+01 Iteration= 1; root mean square error= 0.1792865E+00 Iteration= 2; root mean square error= 0.7124724E-01 Iteration= 3; root mean square error= 0.3737935E-01 Iteration= 4; root mean square error= 0.2641037E-01 Iteration= 5; root mean square error= 0.2077431E-01 Iteration= 6; root mean square error= 0.1725833E-01 Iteration= 7; root mean square error= 0.1477519E-01 Iteration= 8; root mean square error= 0.1288345E-01 Iteration= 9; root mean square error= 0.1137307E-01 Iteration= 10; root mean square error= 0.1012925E-01 Iteration= 11; root mean square error= 0.9082403E-02 Iteration= 12; root mean square error= 0.8187047E-02 Iteration= 13; root mean square error= 0.7411814E-02 Iteration= 14; root mean square error= 0.6734135E-02 Iteration= 15; root mean square error= 0.6137205E-02 Iteration= 16; root mean square error= 0.5608142E-02 Iteration= 17; root mean square error= 0.5136831E-02 Iteration= 18; root mean square error= 0.4715163E-02 Iteration= 19; root mean square error= 0.4336526E-02 Iteration= 20; root mean square error= 0.3995447E-02 Iteration= 21; root mean square error= 0.3687340E-02 Iteration= 22; root mean square error= 0.3408323E-02 Iteration= 23; root mean square error= 0.3155079E-02 Iteration= 24; root mean square error= 0.2924753E-02 Iteration= 25; root mean square error= 0.2714872E-02 Min. area AFTER smoothing= 0.3576634E+00 Max. area AFTER smoothing= 0.1323968E+01 recon/0 RECON2D reconnection loop2d: 1 Recon2d: iteration number= 1 number of flips= 16 RECON: done. rmpoint compress 0 points removed and 0 elements removed. The new point count is 494. The new element count is 894. resetpts itp geniee finish dump gmv tmp_smooth4.gmv cmotri cmo/modatt/-def-/-def-/ioflag/x finish #################### # Final cell and node colors # Select center cylinder set of nodes and cells for well cmo setatt cmotri imt MAT_ID 494 values reset for attribute imt1 cmo setatt cmotri itetclr MAT_ID 894 values reset for attribute itetclr # pset select the center node pset/px/attribute xic/1,0,0/eq 0. THE PSET px HAS 1 POINTS pset/pcenter/attribute yic/pset,get,px/eq 0. THE PSET pcenter HAS 1 POINTS # eltset select cells connected to center node eltset/ewell/inclusive pset,get,pcenter THE ELTSET ewell HAS 7 ELEMENTS # pset select all node vertices of well cells pset/pwell/eltset/ewell THE PSET pwell HAS 8 POINTS # overwrite material MAT_ID with WELL_ID for well cmo setatt cmotri imt pset,get,pwell WELL_ID 8 values reset for attribute imt1 cmo setatt cmotri itetclr eltset,get,ewell WELL_ID 7 values reset for attribute itetclr # reset some lagrit book-keeping arrays resetpts/parent rmpoint/compress 0 points removed and 0 elements removed. The new point count is 494. The new element count is 894. resetpts/itp geniee finish dump/gmv/tmp_smooth_final.gmv/cmotri cmo/modatt/-def-/-def-/ioflag/x finish dump/avs/tmp_smooth_final.inp/cmotri cmo/modatt/-def-/-def-/ioflag/x finish # re-connect with final recon call # to ensure the triangulation is Delaunay recon/1 RECON2D reconnection loop2d: 1 RECON: done. # Write final triangulation dump/gmv/cyl_2D_delaunay.gmv/cmotri cmo/modatt/-def-/-def-/ioflag/x finish dump/avs/cyl_2D_delaunay.inp/cmotri cmo/modatt/-def-/-def-/ioflag/x finish # Write some mesh information cmo/status The current-mesh-object(CMO) is: cmotri 1 Mesh Object name: mo_rim number of nodes = 92 number of elements = 0 dimensions geometry = 3 element type = tet dimensions topology = 3 4 nodes 4 faces 6 edges boundary flag = 16000000 status = inactive 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 9.200E+01 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 0.000E+00 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 4.000E+00 12 edges_per_eleme INT scalar scalar const perma l 6.000E+00 13 faces_per_eleme INT scalar scalar const perma l 4.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 9.421E-12 35 epsilona REAL scalar scalar const perma l 3.997E-10 36 epsilonv REAL scalar scalar const perma l 3.997E-10 37 ipointi INT scalar scalar const perma l 1.000E+00 38 ipointj INT scalar scalar const perma l 9.200E+01 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 -1.500E+01 52 ymin REAL scalar scalar const perma l -1.500E+01 53 zmin REAL scalar scalar const perma l 0.000E+00 54 xmax REAL scalar scalar const perma l 1.500E+01 55 ymax REAL scalar scalar const perma l 1.500E+01 56 zmax REAL scalar scalar const perma l 0.000E+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 idnode0 VINT scalar nnodes linea perma gxa 0.000E+00 2 Mesh Object name: mo_well number of nodes = 15 number of elements = 0 dimensions geometry = 3 element type = tet dimensions topology = 3 4 nodes 4 faces 6 edges boundary flag = 16000000 status = inactive 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 1.500E+01 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 0.000E+00 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 4.000E+00 12 edges_per_eleme INT scalar scalar const perma l 6.000E+00 13 faces_per_eleme INT scalar scalar const perma l 4.000E+00 14 isetwd VDOU scalar nnodes or perma l 0.000E+00 15 ialias VINT scalar nnodes seque perma l 0.000E+00 16 imt1 VINT scalar nnodes max perma gal 0.000E+00 17 itp1 VINT scalar nnodes min perma gal 0.000E+00 18 icr1 VINT scalar nnodes min perma gal 0.000E+00 19 isn1 VINT scalar nnodes user perma gal 0.000E+00 20 xic VDOU scalar nnodes linea perma l 0.000E+00 21 yic VDOU scalar nnodes linea perma l 0.000E+00 22 zic VDOU scalar nnodes linea perma l 0.000E+00 23 xtetwd VDOU scalar nelements or perma l 0.000E+00 24 itetclr VINT scalar nelements user perma l 0.000E+00 25 itettyp VINT scalar nelements user perma l 0.000E+00 26 itetoff VINT scalar nelements user perma l 0.000E+00 27 jtetoff VINT scalar nelements user perma l 0.000E+00 28 itet VINT nodes_per_ nelements user perma l 0.000E+00 29 jtet VINT faces_per_ nelements user perma l 0.000E+00 30 ipolydat CHAR scalar scalar const perma l yes 31 vor2d CHAR scalar scalar const perma l yes 32 vor3d CHAR scalar scalar const perma l no 33 epsilon REAL scalar scalar const perma l 1.000E-15 34 epsilonl REAL scalar scalar const perma l 1.209E-12 35 epsilona REAL scalar scalar const perma l 6.586E-12 36 epsilonv REAL scalar scalar const perma l 6.586E-12 37 ipointi INT scalar scalar const perma l 1.000E+00 38 ipointj INT scalar scalar const perma l 1.500E+01 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 -1.802E+00 52 ymin REAL scalar scalar const perma l -1.950E+00 53 zmin REAL scalar scalar const perma l 0.000E+00 54 xmax REAL scalar scalar const perma l 2.000E+00 55 ymax REAL scalar scalar const perma l 1.950E+00 56 zmax REAL scalar scalar const perma l 0.000E+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 3 Mesh Object name: mo number of nodes = 107 number of elements = 0 dimensions geometry = 3 element type = tet dimensions topology = 3 4 nodes 4 faces 6 edges boundary flag = 16000000 status = inactive NAME TYPE RANK LENGTH INTER PERSI IO VALUE 1 -def- VDOU scalar nnodes linea tempo x 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 1.070E+02 5 nedges INT scalar scalar const perma l 0.000E+00 6 nfaces INT scalar scalar const perma l 0.000E+00 7 nelements INT scalar scalar const perma l 0.000E+00 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 4.000E+00 12 edges_per_eleme INT scalar scalar const perma l 6.000E+00 13 faces_per_eleme INT scalar scalar const perma l 4.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 9.421E-12 35 epsilona REAL scalar scalar const perma l 3.997E-10 36 epsilonv REAL scalar scalar const perma l 3.997E-10 37 ipointi INT scalar scalar const perma l 9.300E+01 38 ipointj INT scalar scalar const perma l 1.070E+02 39 idebug INT scalar scalar const perma l 0.000E+00 40 itypconv_sm INT scalar scalar const perma l 1.000E+00 41 maxiter_sm INT scalar scalar const perma l 2.500E+01 42 tolconv_sm REAL scalar scalar const perma l 1.000E+00 43 nnfreq INT scalar scalar const perma l 1.000E+00 44 ivoronoi INT scalar scalar const perma l 1.000E+00 45 iopt2to2 INT scalar scalar const perma l 2.000E+00 46 dumptype CHAR scalar scalar const perma l binary 47 velname CHAR scalar scalar const perma l vels 48 densname CHAR scalar scalar const perma l ric 49 presname CHAR scalar scalar const perma l pic 50 enername CHAR scalar scalar const perma l eic 51 xmin REAL scalar scalar const perma l -1.500E+01 52 ymin REAL scalar scalar const perma l -1.500E+01 53 zmin REAL scalar scalar const perma l 0.000E+00 54 xmax REAL scalar scalar const perma l 1.500E+01 55 ymax REAL scalar scalar const perma l 1.500E+01 56 zmax REAL scalar scalar const perma l 0.000E+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 idnode0 VINT scalar nnodes linea perma gxa 0.000E+00 4 Mesh Object name: cmotri number of nodes = 494 number of elements = 894 dimensions geometry = 2 element type = tri dimensions topology = 2 3 nodes 3 faces 3 edges boundary flag = 16000000 status = active NAME TYPE RANK LENGTH INTER PERSI IO VALUE 1 -def- VDOU scalar nnodes linea tempo x 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 4.940E+02 5 nedges INT scalar scalar const perma l 0.000E+00 6 nfaces INT scalar scalar const perma l 0.000E+00 7 nelements INT scalar scalar const perma l 8.940E+02 8 mbndry INT scalar scalar const perma l 1.600E+07 9 ndimensions_top INT scalar scalar const perma l 2.000E+00 10 ndimensions_geo INT scalar scalar const perma l 2.000E+00 11 nodes_per_eleme INT scalar scalar const perma l 3.000E+00 12 edges_per_eleme INT scalar scalar const perma l 3.000E+00 13 faces_per_eleme INT scalar scalar const perma l 3.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 9.421E-12 35 epsilona REAL scalar scalar const perma l 3.997E-10 36 epsilonv REAL scalar scalar const perma l 3.997E-10 37 ipointi INT scalar scalar const perma l 1.000E+00 38 ipointj INT scalar scalar const perma l 4.940E+02 39 idebug INT scalar scalar const perma l 0.000E+00 40 itypconv_sm INT scalar scalar const perma l 1.000E+00 41 maxiter_sm INT scalar scalar const perma l 2.500E+01 42 tolconv_sm REAL scalar scalar const perma l 1.000E+00 43 nnfreq INT scalar scalar const perma l 1.000E+00 44 ivoronoi INT scalar scalar const perma l 1.000E+00 45 iopt2to2 INT scalar scalar const perma l 2.000E+00 46 dumptype CHAR scalar scalar const perma l binary 47 velname CHAR scalar scalar const perma l vels 48 densname CHAR scalar scalar const perma l ric 49 presname CHAR scalar scalar const perma l pic 50 enername CHAR scalar scalar const perma l eic 51 xmin REAL scalar scalar const perma l -1.500E+01 52 ymin REAL scalar scalar const perma l -1.500E+01 53 zmin REAL scalar scalar const perma l 0.000E+00 54 xmax REAL scalar scalar const perma l 1.500E+01 55 ymax REAL scalar scalar const perma l 1.500E+01 56 zmax REAL scalar scalar const perma l 0.000E+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 4.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 idnode0 VINT scalar nnodes linea perma gxa 0.000E+00 67 wdist VDOU scalar nnodes linea perma agl -1.000E+00 cmo/printatt//-all-/minmax ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH -def- 0.000000000E+00 0.000000000E+00 0.000000000E+00 494 scalar 1 1 0 1 vector 3 3 0 1 nnodes 494 494 0 1 nedges 0 0 0 1 nfaces 0 0 0 1 nelements 894 894 0 1 mbndry 16000000 16000000 0 1 ndimensions_topo 2 2 0 1 ndimensions_geom 2 2 0 1 nodes_per_element 3 3 0 1 edges_per_element 3 3 0 1 faces_per_element 3 3 0 1 isetwd 0 14 14 494 ialias 0 1401 1401 494 imt1 1 2 1 494 itp1 0 10 10 494 icr1 0 0 0 494 isn1 0 0 0 494 xic -1.500000000E+01 1.500000000E+01 3.000000000E+01 494 yic -1.500000000E+01 1.500000000E+01 3.000000000E+01 494 zic 0.000000000E+00 0.000000000E+00 0.000000000E+00 494 xtetwd 0 1 1 894 itetclr 1 2 1 894 itettyp 3 3 0 894 itetoff 0 2679 2679 894 jtetoff 0 2679 2679 894 itet 1 494 493 894x3 jtet 1 16000152 16000151 894x3 epsilon 1.000000004E-15 1.000000004E-15 0.000000000E+00 1 epsilonl 9.420554752E-12 9.420554752E-12 0.000000000E+00 1 epsilona 3.996802889E-10 3.996802889E-10 0.000000000E+00 1 epsilonv 3.996802889E-10 3.996802889E-10 0.000000000E+00 1 ipointi 1 1 0 1 ipointj 494 494 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.500000000E+01 -1.500000000E+01 0.000000000E+00 1 ymin -1.500000000E+01 -1.500000000E+01 0.000000000E+00 1 zmin 0.000000000E+00 0.000000000E+00 0.000000000E+00 1 xmax 1.500000000E+01 1.500000000E+01 0.000000000E+00 1 ymax 1.500000000E+01 1.500000000E+01 0.000000000E+00 1 zmax 0.000000000E+00 0.000000000E+00 0.000000000E+00 1 kdtree_level 0 0 0 1 max_number_sets 64 64 0 1 number_of_psets 4 4 0 1 number_of_eltsets 0 0 0 1 number_of_fsets 0 0 0 1 idnode0 0 455 455 494 wdist 0.000000000E+00 1.322655919E+01 1.322655919E+01 494 quality epsilonl, epsilonaspect: 9.4205548E-12 8.3604458E-34 -------------------------------------------- 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 : 0 elements with aspect ratio b/w .2 and .5 : 0 elements with aspect ratio b/w .5 and 1. : 894 min aspect ratio = 0.5144E+00 max aspect ratio = 0.9999E+00 epsilonvol: 3.9968029E-10 --------------------------------------- element volumes b/w 0.3577E+00 and 0.4712E+00: 17 element volumes b/w 0.4712E+00 and 0.6209E+00: 103 element volumes b/w 0.6209E+00 and 0.8181E+00: 390 element volumes b/w 0.8181E+00 and 0.1078E+01: 357 element volumes b/w 0.1078E+01 and 0.1420E+01: 27 min volume = 3.5766336E-01 max volume = 1.4201220E+00 ----------------------------------------------------------- 894 total elements evaluated. # end finish LaGriT successfully completed