* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * Program: LaGriT V3.101 Linux m64 * * * date_compile: 2013/04/25 RH gf4.5 * * * Run Time: 2013/Jul 2 16:56:13 * * * 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 30530 self-degenerate faces ELMTESTD WARNING: mesh has 0 degenerate elements and 34552 self-degenerate elements finish Nodes: 95964 Elements: 45315 Done reading GMV binary total bytes: 9993884 file: ice_wedge.gmv cmo/printatt mo_ice -xyz- minmax ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH xic 4.459000015E+01 2.180899963E+02 1.734999962E+02 95964 yic 4.240127182E+01 2.020277863E+02 1.596265144E+02 95964 zic 7.209998965E-01 4.428281307E+00 3.707281411E+00 95964 cmo/printatt mo_ice itetclr minmax ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH itetclr 1 72 71 45315 cmo/setatt/mo_ice/itetclr/30 45315 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 45315 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: 16970053 Elements: 33441792 Done reading GMV binary total bytes: 2140824328 file: full_mesh.gmv cmo/printatt mo -xyz- minmax ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH xic 4.459000015E+01 2.180899963E+02 1.734999962E+02 16970053 yic 4.240127182E+01 2.005397644E+02 1.581384926E+02 16970053 zic -1.000000000E+00 4.997776985E+00 5.997776985E+00 16970053 cmo/addatt/mo/ielemclr/VINT/scalar/nnodes cmo/addatt/mo/mat_orig/VINT/scalar/nelements cmo/copyatt/mo/mo/mat_orig/itetclr 33441792 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 33441792 intrp/map/mo ielemclr/1,0,0/mo_ice itetclr nodes in indexed point set = 16970053 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 339402 0 0 2 % 678803 0 0 4 % 1018204 0 0 6 % 1357605 0 0 8 % 1697006 0 0 10 % 2036407 0 0 12 % 2375808 0 0 14 % 2715209 0 0 16 % 3054610 0 0 18 % 3394011 0 0 20 % 3733412 0 0 22 % 4072813 905 252 24 % 4412214 15458 4270 26 % 4751615 49962 15164 28 % 5091016 87810 27310 30 % 5430417 124956 41128 32 % 5769818 169787 63961 34 % 6109219 224683 94766 36 % 6448620 284990 129846 38 % 6788021 352198 171397 40 % 7127422 426790 220181 42 % 7466823 507469 274713 44 % 7806224 594075 334558 46 % 8145625 686759 399272 48 % 8485026 787481 469313 50 % 8824427 898526 545057 52 % 9163828 1016658 624471 54 % 9503229 1138968 706647 56 % 9842630 1263009 789332 58 % 10182031 1386999 873930 60 % 10521432 1518212 960684 62 % 10860833 1651862 1048332 64 % 11200234 1785831 1138028 66 % 11539635 1921568 1228807 68 % 11879036 2062863 1321393 70 % 12218437 2207532 1415765 72 % 12557838 2279580 1461104 74 % 12897239 2279685 1461192 76 % 13236640 2279685 1461192 78 % 13576041 2279685 1461192 80 % 13915442 2279685 1461192 82 % 14254843 2279685 1461192 84 % 14594244 2279685 1461192 86 % 14933645 2279685 1461192 88 % 15273046 2279685 1461192 90 % 15612447 2279685 1461192 92 % 15951848 2279685 1461192 94 % 16291249 2279685 1461192 96 % 16630650 2279685 1461192 98 % 16970051 2279685 1461192 100 % 16970053 2279685 1461192 Total cmo select mo finish WARNING: Sink points not inside source grid: 15508861 Total sink points marked: 15508861 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 16970053 # 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 1461192 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 3551749 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 33441792 ELEMENTS eltset/e_maxl/itetclr/gt/MAX_ICELAYER THE ELTSET e_maxl HAS 24384640 ELEMENTS # CHECK bounds of ice over all materials eltset e_all_e30 inter e_ice e_all THE ELTSET e_all_e30 HAS 3551749 ELEMENTS # Find elements of ice in permitted layers only eltset/e30/inter/e_ice e_maxl THE ELTSET e30 HAS 3453253 ELEMENTS # set itetclr to ice material cmo/setatt/mo/itetclr/eltset,get,e30/30 3453253 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: 16970053 number of elements: 33441792 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 33441792 cmo/printatt mo mat_orig minmax ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH mat_orig 1 6 5 33441792 cmo/printatt mo -xyz- minmax ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH xic 4.459000015E+01 2.180899963E+02 1.734999962E+02 16970053 yic 4.240127182E+01 2.005397644E+02 1.581384926E+02 16970053 zic -1.000000000E+00 4.997776985E+00 5.997776985E+00 16970053 # ice wedge mesh object cmo/printatt mo_ice itetclr minmax ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH itetclr 30 30 0 45315 cmo/printatt mo_ice mat_orig minmax ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH mat_orig 30 30 0 45315 cmo/printatt mo_ice -xyz- minmax ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH xic 4.459000015E+01 2.180899963E+02 1.734999962E+02 95964 yic 4.240127182E+01 2.020277863E+02 1.596265144E+02 95964 zic 7.209998965E-01 4.428281307E+00 3.707281411E+00 95964 # merge for viewing in same space cmo setatt mo itetclr 1 33441792 values reset for attribute itetclr cmo setatt mo_ice itetclr 2 45315 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