* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * Program: LaGriT V3.200 Linux m64 * * * date_compile: 2016/10/25 Ubu gf * * * Run Time: 2016/Nov 30 11:21:40 * * * Manual: http://lagrit.lanl.gov * * * * * * * * * * * * * * * * * * * * * * * * * * * * * -----oOo----- LaGriT V3 LA-CC-15-069 LaGriT Copyright: This program was prepared by Los Alamos National Security, LLC at Los Alamos National Laboratory (LANL) under contract No. DE-AC52-06NA25396 with the U.S. Department of Energy (DOE). All rights in the program are reserved by the DOE and Los Alamos National Security, LLC. Permission is granted to the public to copy and use this software without charge, provided that this Notice and any statement of authorship are reproduced on all copies. Neither the U.S. Government nor LANS makes any warranty, express or implied, or assumes any liability or responsibility for the use of this software. -----oOo----- Output log file: outx3dgen Command log file: logx3dgen # Interpolate ice wedge mesh on to full stacked mesh # This will overwrite material layers with ice 30 # Input the ice wedge mesh object read/avs/ice_wedge.inp/mo_ice cmo/addatt/mo_ice/layertyp/VINT/scalar/nnodes/linear/permanent/gxaf/0.0 finish geniee geniee: mesh has 469 self-degenerate faces ELMTESTD WARNING: mesh has 0 degenerate elements and 540 self-degenerate elements finish cmo/status/brief The current-mesh-object(CMO) is: mo_ice 1 Mesh Object name: mo_ice number of nodes = 1421 number of elements = 636 dimensions geometry = 3 element type = pri dimensions topology = 3 6 nodes 5 faces 9 edges boundary flag = 16000000 status = active finish cmo/printatt mo_ice -xyz- minmax ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH xic 9.219250000E+01 1.163908000E+02 2.419830000E+01 1421 yic 2.731400000E+00 2.140300000E+01 1.867160000E+01 1421 zic 1.835156625E+00 4.938156068E+00 3.102999443E+00 1421 cmo/printatt mo_ice itetclr minmax ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH itetclr 1 1 0 636 cmo/addatt/mo_ice/mat_orig/VINT/scalar/nelements cmo/setatt/mo_ice/itetclr/30 636 values reset for attribute itetclr cmo/setatt/mo_ice/mat_orig/30 636 values reset for attribute mat_orig # Input the full stacked mesh read/avs/full_mesh.inp/mo cmo/addatt/mo/layertyp/VINT/scalar/nnodes/linear/permanent/gxaf/0.0 finish geniee finish cmo/status/brief The current-mesh-object(CMO) is: mo 1 Mesh Object name: mo_ice number of nodes = 1421 number of elements = 636 dimensions geometry = 3 element type = pri dimensions topology = 3 6 nodes 5 faces 9 edges boundary flag = 16000000 status = inactive 2 Mesh Object name: mo number of nodes = 381603 number of elements = 732160 dimensions geometry = 3 element type = pri dimensions topology = 3 6 nodes 5 faces 9 edges boundary flag = 16000000 status = active finish cmo/printatt mo -xyz- minmax ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH xic 9.219250000E+01 1.163908000E+02 2.419830000E+01 381603 yic 2.731400000E+00 2.140300000E+01 1.867160000E+01 381603 zic -4.500000000E+01 5.504421902E+00 5.050442190E+01 381603 cmo/addatt/mo/node_map/VINT/scalar/nnodes cmo/addatt/mo/node_ice/VINT/scalar/nnodes cmo/addatt/mo/elem_map/VINT/scalar/nelements cmo/addatt/mo/mat_orig/VINT/scalar/nelements cmo/copyatt/mo/mo/mat_orig/itetclr 732160 copied from mo itetclr to -> mo mat_orig # Interpolate ice wedge mesh on to full mesh node attribute node_map # Save itetclr to mat_orig before overwriting cmo/printatt mo itetclr minmax ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH itetclr 1 6 5 732160 interpolate/map/mo node_map/1,0,0/mo_ice itetclr nodes in indexed point set = 381603 INTRP METHOD: map FUNCTION: user TIEBREAKER: maxtie FLAG: plus1 0.310000000000E+02 cmo select mo_ice 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 19081 0 0 5 % 38161 0 0 10 % 57241 0 0 15 % 76321 0 0 20 % 95401 0 0 25 % 114481 0 0 30 % 133561 0 0 35 % 152641 0 0 40 % 171721 0 0 45 % 190801 1317 525 50 % 209881 4315 2101 55 % 228961 8647 4949 60 % 248041 14370 8922 65 % 267121 21601 13746 70 % 286201 29468 18831 75 % 305281 37709 24168 80 % 324361 38746 24842 85 % 343441 38746 24842 90 % 362521 38746 24842 95 % 381601 38746 24842 100 % 381603 38746 24842 Total cmo select mo finish WARNING: Sink points not inside source grid: 356761 Total sink points marked: 356761 with flag: 31.00 intrp/map done. cmo kdtree release cmo/DELATT/mo_ice/sbox finish cmo/DELATT/mo_ice/v12 finish cmo/DELATT/mo_ice/linkt finish cmo/DELATT/mo_ice/v2 finish finish interpolate/map/mo elem_map/1,0,0/mo_ice itetclr elements in indexed set = 732160 INTRP METHOD: map FUNCTION: user TIEBREAKER: maxtie FLAG: plus1 0.310000000000E+02 cmo select mo_ice 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 14644 0 0 2 % 29287 0 0 4 % 43930 0 0 6 % 58573 0 0 8 % 73216 0 0 10 % 87859 0 0 12 % 102502 0 0 14 % 117145 0 0 16 % 131788 0 0 18 % 146431 0 0 20 % 161074 0 0 22 % 175717 0 0 24 % 190360 0 0 26 % 205003 0 0 28 % 219646 0 0 30 % 234289 0 0 32 % 248932 0 0 34 % 263575 0 0 36 % 278218 0 0 38 % 292861 0 0 40 % 307504 0 0 42 % 322147 0 0 44 % 336790 0 0 46 % 351433 0 0 48 % 366076 173 173 50 % 380719 834 834 52 % 395362 1925 1925 54 % 410005 3350 3350 56 % 424648 5221 5219 58 % 439291 7425 7407 60 % 453934 10062 9974 62 % 468577 13274 12989 64 % 483220 17124 16330 66 % 497863 21218 19771 68 % 512506 25503 23308 70 % 527149 30001 26934 72 % 541792 34697 30696 74 % 556435 39639 34575 76 % 571078 44610 38447 78 % 585721 49863 42474 80 % 600364 51918 44054 82 % 615007 51918 44054 84 % 629650 51918 44054 86 % 644293 51918 44054 88 % 658936 51918 44054 90 % 673579 51918 44054 92 % 688222 51918 44054 94 % 702865 51918 44054 96 % 717508 51918 44054 98 % 732151 51918 44054 100 % 732160 51918 44054 Total cmo select mo finish WARNING: Sink points not inside source grid: 688106 Total sink points marked: 688106 with flag: 31.00 intrp/map done. cmo kdtree release cmo/DELATT/mo_ice/sbox finish cmo/DELATT/mo_ice/v12 finish cmo/DELATT/mo_ice/linkt finish cmo/DELATT/mo_ice/v2 finish finish # CHECK nodes should be values 30 and 31 (outside ice) cmo/printatt/mo/node_map/minmax ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH node_map 30 31 1 381603 cmo/select/mo # Find nodes and elems set with ice value 30 # note inclusion or exclusion options select more or less # eltset / e_ice / exclusive / pset get p_ice pset/p_ice/attribute/node_map/1,0,0/eq/30 THE PSET p_ice HAS 24842 POINTS eltset/e_ice1/inclusive/pset get p_ice THE ELTSET e_ice1 HAS 54436 ELEMENTS eltset/e_ice2/elem_map/eq 30 THE ELTSET e_ice2 HAS 44054 ELEMENTS eltset/e_ice/union e_ice1 e_ice2 THE ELTSET e_ice HAS 54487 ELEMENTS # Find elements in layers where ice exists # Note Material order 1 at top to bottom n define MAX_ICELAYER 4 eltset e_all itetclr ge 1 THE ELTSET e_all HAS 732160 ELEMENTS eltset/e_maxl/itetclr/gt/MAX_ICELAYER THE ELTSET e_maxl HAS 591360 ELEMENTS # CHECK bounds of ice over all materials eltset e_all_e30 inter e_ice e_all THE ELTSET e_all_e30 HAS 54487 ELEMENTS # Find elements of ice in permitted layers only eltset/e30/inter/e_ice e_maxl THE ELTSET e30 HAS 52987 ELEMENTS # set itetclr to ice material cmo/setatt/mo/itetclr/eltset,get,e30/30 52987 values reset for attribute itetclr # set node_ice to ice material pset/pts_ice2/eltset/e30 THE PSET pts_ice2 HAS 35890 POINTS cmo setatt mo node_ice 0 381603 values reset for attribute node_ice cmo setatt mo node_ice pset,get,pts_ice2/1 35890 values reset for attribute node_ice # make sure all nodes on outside faces are captured pset/pout/attribute itp/pset,get,pts_ice2/ge 10 THE PSET pout HAS 8639 POINTS pset/pfix/attribute node_map/pset,get,pout/lt 10 THE PSET pfix HAS 0 POINTS cmo setatt mo node_map pset,get,pfix 3 No values acquired for indexed set of node_map dump/proto_full_poly.gmv/mo dump/gmv/proto_full_poly.gmv/mo cmo/modatt/-def-/-def-/ioflag/x finish finish # remove extra attributes cmo/DELATT/mo/ikey_utr DELATT no action, attribute does not exist: ikey_utr cmo/DELATT/mo/node_map cmo/DELATT/mo/elem_map # Final Mesh with ice: dump/proto_full_poly.exo/mo dump/exo/proto_full_poly.exo/mo ExodusII: Start writing to file: proto_full_poly.exo using cmo: mo cmo/set_id/mo/element/e_num_temp cmo/addatt/mo e_num_temp/VINT/scalar/nelements/linear/permanent//0 finish finish sort/mo/index/ascending/ikey_utr/itetclr/e_num_temp cmo/addatt/mo/ikey_utr/vint/scalar/nelements///gax/0 finish SORT: order key written to attribute: ikey_utr finish cmo/DELATT/mo/e_num_temp finish cmo/modatt/mo/itetclr/ioflag/lga finish INITIALIZE exodus Title: LAGRIT TO EXODUS number of dimension: 3 number of nodes: 381603 number of elements: 732160 number of edges: 0 number of edge blocks: 0 number of element blocks: 7 number of face blocks: 0 number of node sets: 0 number of edge sets: 3 number of element sets: 0 number of side sets: 0 number of face sets: 0 number of node maps: 0 number of edge maps: 0 number of face maps: 0 number of element maps: 0 cmo/modatt/mo/itetclr/ioflag/l finish ExodusII: Done writing to ExodusII file: proto_full_poly.exo using cmo: mo finish dump/proto_full_poly.inp/mo dump/avs/proto_full_poly.inp/mo cmo/modatt/-def-/-def-/ioflag/x finish finish ############# CHECK FINAL MESH ####################### # merge for viewing in same space cmo copy mo_merge mo cmo setatt mo_merge itetclr 1 732160 values reset for attribute itetclr cmo setatt mo_ice itetclr 2 636 values reset for attribute itetclr addmesh/merge/mo_merge/mo_merge/mo_ice cmo/move/-tmp_source_internal-/mo_merge finish searching for mo_merge attributes in mo_ice attribute not found: 67 node_ice attribute not found: 69 ikey_utr ADDMESH: 2 attributes not copied from mo_ice to mo_merge cmo/delete/-tmp_source_internal-/ Released Mesh Object: -tmp_source_internal- finish dump gmv ck_merge_full_ice.gmv mo_merge cmo/modatt/-def-/-def-/ioflag/x finish cmo delete mo_merge Released Mesh Object: mo_merge # ice wedge mesh object cmo/printatt mo_ice itetclr minmax ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH itetclr 2 2 0 636 cmo/printatt mo_ice mat_orig minmax ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH mat_orig 30 30 0 636 cmo/printatt mo_ice -xyz- minmax ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH xic 9.219250000E+01 1.163908000E+02 2.419830000E+01 1421 yic 2.731400000E+00 2.140300000E+01 1.867160000E+01 1421 zic 1.835156625E+00 4.938156068E+00 3.102999443E+00 1421 # Full stacked mesh cmo/printatt mo itetclr minmax ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH itetclr 1 30 29 732160 cmo/printatt mo mat_orig minmax ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH mat_orig 1 6 5 732160 cmo/printatt mo -xyz- minmax ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH xic 9.219250000E+01 1.163908000E+02 2.419830000E+01 381603 yic 2.731400000E+00 2.140300000E+01 1.867160000E+01 381603 zic -4.500000000E+01 5.504421902E+00 5.050442190E+01 381603 cmo/printatt mo -all- minmax ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH -def- 0.000000000E+00 0.000000000E+00 0.000000000E+00 381603 scalar 1 1 0 1 vector 3 3 0 1 nnodes 381603 381603 0 1 nedges 0 0 0 1 nfaces 0 0 0 1 nelements 732160 732160 0 1 mbndry 16000000 16000000 0 1 ndimensions_topo 3 3 0 1 ndimensions_geom 3 3 0 1 nodes_per_element 6 6 0 1 edges_per_element 9 9 0 1 faces_per_element 5 5 0 1 isetwd 0 7 7 381603 ialias 0 0 0 381603 imt1 1 6 5 381603 itp1 0 10 10 381603 icr1 0 0 0 381603 isn1 0 0 0 381603 xic 9.219250000E+01 1.163908000E+02 2.419830000E+01 381603 yic 2.731400000E+00 2.140300000E+01 1.867160000E+01 381603 zic -4.500000000E+01 5.504421902E+00 5.050442190E+01 381603 xtetwd 8 127 119 732160 itetclr 1 30 29 732160 itettyp 7 7 0 732160 itetoff 0 4392954 4392954 732160 jtetoff 0 3660795 3660795 732160 itet 1 381603 381602 732160x6 jtet 2 19660796 19660794 732160x5 epsilon 1.000000004E-15 1.000000004E-15 0.000000000E+00 1 epsilonl 1.310793643E-11 1.310793643E-11 0.000000000E+00 1 epsilona 7.737994692E-10 7.737994692E-10 0.000000000E+00 1 epsilonv 5.066826361E-09 5.066826361E-09 0.000000000E+00 1 ipointi 1 1 0 1 ipointj 381603 381603 0 1 idebug 0 0 0 1 itypconv_sm 1 1 0 1 maxiter_sm 25 25 0 1 tolconv_sm 1.000000000E+00 1.000000000E+00 0.000000000E+00 1 nnfreq 1 1 0 1 ivoronoi 1 1 0 1 iopt2to2 2 2 0 1 xmin 9.219250000E+01 9.219250000E+01 0.000000000E+00 1 ymin 2.731400000E+00 2.731400000E+00 0.000000000E+00 1 zmin -4.500000000E+01 -4.500000000E+01 0.000000000E+00 1 xmax 1.163908000E+02 1.163908000E+02 0.000000000E+00 1 ymax 2.140300000E+01 2.140300000E+01 0.000000000E+00 1 zmax 5.504421902E+00 5.504421902E+00 0.000000000E+00 1 kdtree_level 0 0 0 1 max_number_sets 64 64 0 1 number_of_psets 4 4 0 1 number_of_eltsets 0 0 0 1 number_of_fsets 0 0 0 1 layertyp -2 2 4 381603 node_ice 0 1 1 381603 mat_orig 1 6 5 732160 ikey_utr 1 732160 732159 732160 finish LaGriT successfully completed