* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * Program: LaGriT V3.200 Linux m64 * * * date_compile: 2016/10/25 Ubu gf * * * Run Time: 2016/Nov 16 10:13:30 * * * 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 457 self-degenerate faces ELMTESTD WARNING: mesh has 0 degenerate elements and 508 self-degenerate elements finish cmo/status/brief The current-mesh-object(CMO) is: mo_ice 1 Mesh Object name: mo_ice number of nodes = 1601 number of elements = 688 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 1.023017000E+02 1.275104000E+02 2.520870000E+01 1601 yic 8.497500000E+01 1.096966000E+02 2.472160000E+01 1601 zic 1.724107000E+00 4.792780936E+00 3.068673936E+00 1601 cmo/printatt mo_ice itetclr minmax ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH itetclr 1 1 0 688 cmo/addatt/mo_ice/mat_orig/VINT/scalar/nelements cmo/setatt/mo_ice/itetclr/30 688 values reset for attribute itetclr cmo/setatt/mo_ice/mat_orig/30 688 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 = 1601 number of elements = 688 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 = 413043 number of elements = 798720 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 1.023017000E+02 1.275104000E+02 2.520870000E+01 413043 yic 8.497500000E+01 1.096966000E+02 2.472160000E+01 413043 zic -4.500000000E+01 5.448649372E+00 5.044864937E+01 413043 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 798720 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 798720 interpolate/map/mo node_map/1,0,0/mo_ice itetclr nodes in indexed point set = 413043 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 20653 0 0 5 % 41305 0 0 10 % 61957 0 0 15 % 82609 0 0 20 % 103261 0 0 25 % 123913 0 0 30 % 144565 0 0 35 % 165217 0 0 40 % 185869 0 0 45 % 206521 1230 476 50 % 227173 4353 1977 55 % 247825 8680 4708 60 % 268477 14416 8590 65 % 289129 21340 13220 70 % 309781 29037 18209 75 % 330433 36860 23397 80 % 351085 37990 24099 85 % 371737 37990 24099 90 % 392389 37990 24099 95 % 413041 37990 24099 100 % 413043 37990 24099 Total cmo select mo finish WARNING: Sink points not inside source grid: 388944 Total sink points marked: 388944 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 = 798720 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 15975 0 0 2 % 31949 0 0 4 % 47923 0 0 6 % 63897 0 0 8 % 79871 0 0 10 % 95845 0 0 12 % 111819 0 0 14 % 127793 0 0 16 % 143767 0 0 18 % 159741 0 0 20 % 175715 0 0 22 % 191689 0 0 24 % 207663 0 0 26 % 223637 0 0 28 % 239611 0 0 30 % 255585 0 0 32 % 271559 0 0 34 % 287533 0 0 36 % 303507 0 0 38 % 319481 0 0 40 % 335455 0 0 42 % 351429 0 0 44 % 367403 0 0 46 % 383377 0 0 48 % 399351 136 136 50 % 415325 702 702 52 % 431299 1745 1745 54 % 447273 3074 3074 56 % 463247 4821 4820 58 % 479221 6904 6897 60 % 495195 9379 9328 62 % 511169 12425 12212 64 % 527143 16043 15416 66 % 543117 19866 18722 68 % 559091 23861 22114 70 % 575065 28058 25609 72 % 591039 32494 29271 74 % 607013 37004 32967 76 % 622987 41582 36694 78 % 638961 46327 40526 80 % 654935 48236 42069 82 % 670909 48236 42069 84 % 686883 48236 42069 86 % 702857 48236 42069 88 % 718831 48236 42069 90 % 734805 48236 42069 92 % 750779 48236 42069 94 % 766753 48236 42069 96 % 782727 48236 42069 98 % 798701 48236 42069 100 % 798720 48236 42069 Total cmo select mo finish WARNING: Sink points not inside source grid: 756651 Total sink points marked: 756651 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 413043 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 24099 POINTS eltset/e_ice1/inclusive/pset get p_ice THE ELTSET e_ice1 HAS 53023 ELEMENTS eltset/e_ice2/elem_map/eq 30 THE ELTSET e_ice2 HAS 42069 ELEMENTS eltset/e_ice/union e_ice1 e_ice2 THE ELTSET e_ice HAS 53041 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 798720 ELEMENTS eltset/e_maxl/itetclr/gt/MAX_ICELAYER THE ELTSET e_maxl HAS 645120 ELEMENTS # CHECK bounds of ice over all materials eltset e_all_e30 inter e_ice e_all THE ELTSET e_all_e30 HAS 53041 ELEMENTS # Find elements of ice in permitted layers only eltset/e30/inter/e_ice e_maxl THE ELTSET e30 HAS 51522 ELEMENTS # set itetclr to ice material cmo/setatt/mo/itetclr/eltset,get,e30/30 51522 values reset for attribute itetclr # set node_ice to ice material pset/pts_ice2/eltset/e30 THE PSET pts_ice2 HAS 34577 POINTS cmo setatt mo node_ice 0 413043 values reset for attribute node_ice cmo setatt mo node_ice pset,get,pts_ice2/1 34577 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 7199 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: 413043 number of elements: 798720 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 798720 values reset for attribute itetclr cmo setatt mo_ice itetclr 2 688 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 688 cmo/printatt mo_ice mat_orig minmax ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH mat_orig 30 30 0 688 cmo/printatt mo_ice -xyz- minmax ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH xic 1.023017000E+02 1.275104000E+02 2.520870000E+01 1601 yic 8.497500000E+01 1.096966000E+02 2.472160000E+01 1601 zic 1.724107000E+00 4.792780936E+00 3.068673936E+00 1601 # Full stacked mesh cmo/printatt mo itetclr minmax ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH itetclr 1 30 29 798720 cmo/printatt mo mat_orig minmax ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH mat_orig 1 6 5 798720 cmo/printatt mo -xyz- minmax ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH xic 1.023017000E+02 1.275104000E+02 2.520870000E+01 413043 yic 8.497500000E+01 1.096966000E+02 2.472160000E+01 413043 zic -4.500000000E+01 5.448649372E+00 5.044864937E+01 413043 cmo/printatt mo -all- minmax ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH -def- 0.000000000E+00 0.000000000E+00 0.000000000E+00 413043 scalar 1 1 0 1 vector 3 3 0 1 nnodes 413043 413043 0 1 nedges 0 0 0 1 nfaces 0 0 0 1 nelements 798720 798720 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 413043 ialias 0 0 0 413043 imt1 1 6 5 413043 itp1 0 10 10 413043 icr1 0 0 0 413043 isn1 0 0 0 413043 xic 1.023017000E+02 1.275104000E+02 2.520870000E+01 413043 yic 8.497500000E+01 1.096966000E+02 2.472160000E+01 413043 zic -4.500000000E+01 5.448649372E+00 5.044864937E+01 413043 xtetwd 8 127 119 798720 itetclr 1 30 29 798720 itettyp 7 7 0 798720 itetoff 0 4792314 4792314 798720 jtetoff 0 3993595 3993595 798720 itet 1 413043 413042 798720x6 jtet 2 19993596 19993594 798720x5 epsilon 1.000000004E-15 1.000000004E-15 0.000000000E+00 1 epsilonl 1.367279621E-11 1.367279621E-11 0.000000000E+00 1 epsilona 8.419270360E-10 8.419270360E-10 0.000000000E+00 1 epsilonv 6.980986437E-09 6.980986437E-09 0.000000000E+00 1 ipointi 1 1 0 1 ipointj 413043 413043 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.023017000E+02 1.023017000E+02 0.000000000E+00 1 ymin 8.497500000E+01 8.497500000E+01 0.000000000E+00 1 zmin -4.500000000E+01 -4.500000000E+01 0.000000000E+00 1 xmax 1.275104000E+02 1.275104000E+02 0.000000000E+00 1 ymax 1.096966000E+02 1.096966000E+02 0.000000000E+00 1 zmax 5.448649372E+00 5.448649372E+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 413043 node_ice 0 1 1 413043 mat_orig 1 6 5 798720 ikey_utr 1 798720 798719 798720 finish LaGriT successfully completed