* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * Program: LaGriT V3.101 Linux m64 * * * date_compile: 2013/04/25 RH gf4.5 * * * Run Time: 2013/Jul 11 15:12:57 * * * 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 # Interpolate ice wedge mesh on to full stacked mesh # This will overwrite material layers with ice 30 # Input the ice wedge mesh object read/gmv/ice_wedge.gmv/mo_ice Reading GMV binary file gmvinputieee cmo/addatt/mo_ice/layertyp/VDOUBLE/scalar/nnodes/linear/permanent/gxa/0.0 finish no geometry for mesh object geniee geniee: mesh has 8770 self-degenerate faces ELMTESTD WARNING: mesh has 0 degenerate elements and 9836 self-degenerate elements finish Nodes: 28388 Elements: 12714 Done reading GMV binary total bytes: 2855620 file: ice_wedge.gmv cmo/printatt mo_ice -xyz- minmax ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH xic 7.933999634E+01 1.909821930E+02 1.116421967E+02 28388 yic 4.647000122E+01 1.617087402E+02 1.152387390E+02 28388 zic 9.624001384E-01 4.356940746E+00 3.394540608E+00 28388 cmo/printatt mo_ice itetclr minmax ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH itetclr 1 17 16 12714 cmo/setatt/mo_ice/itetclr/30 12714 values reset for attribute itetclr cmo/copyatt/mo_ice/mo_ice/mat_orig/itetclr CMO_COPYATT WARNING: attribute does not exist: attribute now created: mat_orig cmo/addatt/mo_ice/mat_orig/VINT/scalar/nelements/user/permanent/l finish 12714 copied from mo_ice itetclr to -> mo_ice mat_orig # Input the full stacked mesh read/gmv/full_mesh.gmv/mo Reading GMV binary file gmvinputieee cmo/addatt/mo/layertyp/VDOUBLE/scalar/nnodes/linear/permanent/gxa/0.0 finish no geometry for mesh object geniee finish Nodes: 5242753 Elements: 10285056 Done reading GMV binary total bytes: 661870456 file: full_mesh.gmv cmo/printatt mo -xyz- minmax ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH xic 7.933999634E+01 1.909821930E+02 1.116421967E+02 5242753 yic 4.647000122E+01 1.617087402E+02 1.152387390E+02 5242753 zic -1.000000000E+00 4.908469200E+00 5.908469200E+00 5242753 cmo/addatt/mo/ielemclr/VINT/scalar/nnodes cmo/addatt/mo/mat_orig/VINT/scalar/nelements cmo/copyatt/mo/mo/mat_orig/itetclr 10285056 copied from mo itetclr to -> mo mat_orig # Interpolate ice wedge mesh on to full mesh node attribute ielemclr # Save itetclr to mat_orig before overwriting cmo/printatt mo itetclr minmax ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH itetclr 1 6 5 10285056 intrp/map/mo ielemclr/1,0,0/mo_ice itetclr nodes in indexed point set = 5242753 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 104856 0 0 2 % 209711 0 0 4 % 314566 0 0 6 % 419421 0 0 8 % 524276 0 0 10 % 629131 0 0 12 % 733986 0 0 14 % 838841 0 0 16 % 943696 0 0 18 % 1048551 0 0 20 % 1153406 0 0 22 % 1258261 0 0 24 % 1363116 2927 834 26 % 1467971 12450 4031 28 % 1572826 22783 7528 30 % 1677681 32711 11294 32 % 1782536 44936 17612 34 % 1887391 59927 26264 36 % 1992246 76559 36143 38 % 2097101 94996 47756 40 % 2201956 115653 61548 42 % 2306811 137974 76906 44 % 2411666 161919 93749 46 % 2516521 187622 111988 48 % 2621376 215548 131680 50 % 2726231 246591 153145 52 % 2831086 279633 175621 54 % 2935941 313789 198809 56 % 3040796 348601 222228 58 % 3145651 383297 246073 60 % 3250506 420126 270597 62 % 3355361 457529 295368 64 % 3460216 494978 320648 66 % 3565071 532935 346204 68 % 3669926 572428 372327 70 % 3774781 612982 398990 72 % 3879636 633492 412122 74 % 3984491 633520 412144 76 % 4089346 633520 412144 78 % 4194201 633520 412144 80 % 4299056 633520 412144 82 % 4403911 633520 412144 84 % 4508766 633520 412144 86 % 4613621 633520 412144 88 % 4718476 633520 412144 90 % 4823331 633520 412144 92 % 4928186 633520 412144 94 % 5033041 633520 412144 96 % 5137896 633520 412144 98 % 5242751 633520 412144 100 % 5242753 633520 412144 Total cmo select mo finish WARNING: Sink points not inside source grid: 4830609 Total sink points marked: 4830609 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/ielemclr/minmax ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH ielemclr 30 31 1 5242753 # Find nodes set with ice value 30 cmo/select/mo pset/p_ice/attribute/ielemclr/1,0,0/eq/30 THE PSET p_ice HAS 412144 POINTS # Find elements with nodes set to ice # note inclusion or exclusion options select more or less eltset/e_ice/inclusive/pset get p_ice THE ELTSET e_ice HAS 982057 ELEMENTS #eltset / e_ice / exclusive / pset get p_ice # 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 10285056 ELEMENTS eltset/e_maxl/itetclr/gt/MAX_ICELAYER THE ELTSET e_maxl HAS 7499520 ELEMENTS # CHECK bounds of ice over all materials eltset e_all_e30 inter e_ice e_all THE ELTSET e_all_e30 HAS 982057 ELEMENTS # Find elements of ice in permitted layers only eltset/e30/inter/e_ice e_maxl THE ELTSET e30 HAS 954300 ELEMENTS # set itetclr to ice material cmo/setatt/mo/itetclr/eltset,get,e30/30 954300 values reset for attribute itetclr 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: 5242753 number of elements: 10285056 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 file: proto_full_poly.exo using cmo: mo finish dump/proto_full_poly.gmv/mo dump/gmv/proto_full_poly.gmv/mo cmo/modatt/-def-/-def-/ioflag/x finish 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 ####################### # Full stacked mesh cmo/printatt mo itetclr minmax ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH itetclr 1 30 29 10285056 cmo/printatt mo mat_orig minmax ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH mat_orig 1 6 5 10285056 cmo/printatt mo -xyz- minmax ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH xic 7.933999634E+01 1.909821930E+02 1.116421967E+02 5242753 yic 4.647000122E+01 1.617087402E+02 1.152387390E+02 5242753 zic -1.000000000E+00 4.908469200E+00 5.908469200E+00 5242753 # ice wedge mesh object cmo/printatt mo_ice itetclr minmax ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH itetclr 30 30 0 12714 cmo/printatt mo_ice mat_orig minmax ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH mat_orig 30 30 0 12714 cmo/printatt mo_ice -xyz- minmax ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH xic 7.933999634E+01 1.909821930E+02 1.116421967E+02 28388 yic 4.647000122E+01 1.617087402E+02 1.152387390E+02 28388 zic 9.624001384E-01 4.356940746E+00 3.394540608E+00 28388 # merge for viewing in same space cmo setatt mo itetclr 1 10285056 values reset for attribute itetclr cmo setatt mo_ice itetclr 2 12714 values reset for attribute itetclr addmesh/merge/mo_merge/mo/mo_ice searching for mo_merge attributes in mo_ice attribute not found: 67 ielemclr attribute not found: 69 ikey_utr ADDMESH: 2 attributes not copied from mo_ice to mo_merge dump gmv ck_merge_full_ice.gmv mo_merge cmo/modatt/-def-/-def-/ioflag/x finish # DONE # Final Files: # proto_full_poly.exo # proto_full_poly.gmv # proto_full_poly.inp finish LaGriT successfully completed