* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * Program: LaGriT V3.101 Linux m64 * * * date_compile: 2013/04/25 RH gf4.5 * * * Run Time: 2013/Nov 6 12:14:48 * * * 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 758 self-degenerate faces ELMTESTD WARNING: mesh has 0 degenerate elements and 876 self-degenerate elements finish Nodes: 2755 Elements: 1122 Done reading GMV binary total bytes: 261544 file: ice_wedge.gmv cmo/printatt mo_ice -xyz- minmax ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH xic 1.030189209E+02 1.445391388E+02 4.152021790E+01 2755 yic 1.234896774E+02 1.485746765E+02 2.508499908E+01 2755 zic 1.081499815E+00 4.203951359E+00 3.122451544E+00 2755 cmo/printatt mo_ice itetclr minmax ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH itetclr 1 1 0 1122 cmo/setatt/mo_ice/itetclr/30 1122 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 1122 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: 698148 Elements: 1358208 Done reading GMV binary total bytes: 86462180 file: full_mesh.gmv cmo/printatt mo -xyz- minmax ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH xic 1.030189209E+02 1.445391388E+02 4.152021790E+01 698148 yic 1.234896774E+02 1.485746765E+02 2.508499908E+01 698148 zic -4.500000000E+01 4.733375072E+00 4.973337507E+01 698148 cmo/addatt/mo/ielemclr/VINT/scalar/nnodes cmo/addatt/mo/mat_orig/VINT/scalar/nelements cmo/copyatt/mo/mo/mat_orig/itetclr 1358208 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 1358208 intrp/map/mo ielemclr/1,0,0/mo_ice itetclr nodes in indexed point set = 698148 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 13964 0 0 2 % 27927 0 0 4 % 41890 0 0 6 % 55853 0 0 8 % 69816 0 0 10 % 83779 0 0 12 % 97742 0 0 14 % 111705 0 0 16 % 125668 0 0 18 % 139631 0 0 20 % 153594 0 0 22 % 167557 0 0 24 % 181520 0 0 26 % 195483 0 0 28 % 209446 0 0 30 % 223409 0 0 32 % 237372 0 0 34 % 251335 0 0 36 % 265298 0 0 38 % 279261 0 0 40 % 293224 0 0 42 % 307187 0 0 44 % 321150 0 0 46 % 335113 740 269 48 % 349076 1870 695 50 % 363039 3175 1219 52 % 377002 5097 2250 54 % 390965 6851 3389 56 % 404928 9315 4981 58 % 418891 11902 6743 60 % 432854 14630 8735 62 % 446817 18114 11143 64 % 460780 21842 13699 66 % 474743 25610 16341 68 % 488706 29685 19049 70 % 502669 33779 21827 72 % 516632 38011 24690 74 % 530595 42628 27686 76 % 544558 47087 30694 78 % 558521 51575 33646 80 % 572484 52147 34021 82 % 586447 52147 34021 84 % 600410 52147 34021 86 % 614373 52147 34021 88 % 628336 52147 34021 90 % 642299 52147 34021 92 % 656262 52147 34021 94 % 670225 52147 34021 96 % 684188 52147 34021 98 % 698148 52147 34021 Total cmo select mo finish WARNING: Sink points not inside source grid: 664127 Total sink points marked: 664127 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 698148 # 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 34021 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 76859 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 1358208 ELEMENTS eltset/e_maxl/itetclr/gt/MAX_ICELAYER THE ELTSET e_maxl HAS 1088640 ELEMENTS # CHECK bounds of ice over all materials eltset e_all_e30 inter e_ice e_all THE ELTSET e_all_e30 HAS 76859 ELEMENTS # Find elements of ice in permitted layers only eltset/e30/inter/e_ice e_maxl THE ELTSET e30 HAS 74575 ELEMENTS # set itetclr to ice material cmo/setatt/mo/itetclr/eltset,get,e30/30 74575 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: 698148 number of elements: 1358208 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 1358208 cmo/printatt mo mat_orig minmax ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH mat_orig 1 6 5 1358208 cmo/printatt mo -xyz- minmax ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH xic 1.030189209E+02 1.445391388E+02 4.152021790E+01 698148 yic 1.234896774E+02 1.485746765E+02 2.508499908E+01 698148 zic -4.500000000E+01 4.733375072E+00 4.973337507E+01 698148 # ice wedge mesh object cmo/printatt mo_ice itetclr minmax ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH itetclr 30 30 0 1122 cmo/printatt mo_ice mat_orig minmax ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH mat_orig 30 30 0 1122 cmo/printatt mo_ice -xyz- minmax ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH xic 1.030189209E+02 1.445391388E+02 4.152021790E+01 2755 yic 1.234896774E+02 1.485746765E+02 2.508499908E+01 2755 zic 1.081499815E+00 4.203951359E+00 3.122451544E+00 2755 # merge for viewing in same space cmo setatt mo itetclr 1 1358208 values reset for attribute itetclr cmo setatt mo_ice itetclr 2 1122 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