************************************************** Fine deep mesh for Area B polygon #1 Version 4.2.0.B.1 Jan 11 2016 replaces Nov 16 2016 data not at cell centers Terry Miller tamiller@lanl.gov Arctic Projects listed at: https://meshing.lanl.gov *************************************************** Based on V3.2.C.1 Area C fine deep mesh in /scratch/sft2/lucia/grid_gen/arctic/arctic_fine_deep_poly06/ Modified so moss has 1 cell layer instead of 2 cells. The moss layer thickness remains .02 meters. Updated scripts in /scratch/sft/tam/grid_gen/arctic/areaa_fine_deep Layer 1 (top) is moss and should be 2 cm thick with 1 cell. The organic layers combined are 48 cm thick with 24 cells (2 cm each). Top organic and mineral layer are the same. Thickness is 18 cm with 9 cells. Lowest organic layer is about 3/4 the thickness of upper 2 layers. Set thickness to 12 cm with 6 cells. The bottom Materials 5 and 6 are both mineral soils with higher resolution at the top of the mineral soil. The bottom material 6 has variable thickness with higher resolution layers at the top. The domain bottom should be flat at Z = -45. The top of the ice wedge should be at 50 cm depth (corresponds to the top of the inorganic layer). Top of wedge starts at top of mineral layer (replaces cells only in materials 5 and 6) Top wedge width (distance from polygon) is aprox 1.2 Depth from top of ice to tip is aprox 2.8 *************************************************** Input Data *************************************************** Sent by Lauren Charsley-Groffman lcg@lanl.gov Earth & Environmental Sciences Division, EES-14 Los Alamos National Laboratory Replaces AreaB_QGIS_clip.tif (Nov 3 2016) 2013DEM_AreaB.png image of data 2013DEM_AreaB_about.txt data descriptin with extents 2013DEM_AreaB_xyz.txt ARCtools raster-> x,y,z points 2013DEM_AreaB.txt reformat point data, precision fix DEM Raster: 869 * 874 = 759506 points at .25 spacing Left/Bot (llcorner): 585725.75 7910149.25 Right/Top (urcorner): 585943.00 7910367.75 DEM ELEVATION POINTS: min XY: 585725.875 7910149.375 (DEM extents +.125) max XY: 585942.875 7910367.625 (DEM extents -.125) Polygon_HCP_areaB.zip (Oct 27 2016) 1 585828.3017 7910249.9348 5.0791 2 585840.6016 7910258.9466 5.0466 3 585853.5104 7910247.6210 5.0969 4 585847.1778 7910235.4429 5.1341 5 585833.7819 7910234.2250 5.1177 Translate polygon to match DEM data at 0,0 define XMIN 585725.75 define YMIN 7910149.25 trans/1,0,0/XMIN YMIN 0. / 0. 0. 0. *************************************************** Mesh Description *************************************************** Mesh file: mesh_fsets.exo Mesh Element type: prism Software: LaGriT, ExodusII, offset_polygon Total nodes: 413043 Total elements: 798720 Total Volume: 19973.59 Element type: prism MATERIALS: ID Name Mat Thickness Num Cells 1 moss .02 1 2 upper organic .18 9 3 upper mineral .18 9 4 lower organic .12 6 5 lower mineral .50 20 6 deep mineral ~49 85 7 ice ~2.8 -- from file write_exo.out.txt (length=nodes) NODES MIN MAX DIFFERENCE LENGTH xic 1.025517000E+02 1.277604000E+02 2.520870000E+01 413043 yic 8.497500000E+01 1.096966000E+02 2.472160000E+01 413043 zic -4.500000000E+01 5.441364405E+00 5.044136441E+01 413043 from file stack_full_poly.out.txt awk -f ../lagrit_input/report_tops.awk stack_full_poly.out.txt LAYER TOP Z MIN Z MAX DEPTH 1 5.0200 5.4414 0.02 2 5.0000 5.4214 0.18 3 4.8200 5.2414 0.18 4 4.6400 5.0614 0.12 5 4.5200 4.9414 0.50 6 4.0200 4.4414 49.44 from file write_exo.out.txt Color Num. Elements Volume Fractional Volume 1 6144 0.7944030E+01 0.000397727 2 55296 0.7149627E+02 0.003579540 3 55296 0.7149627E+02 0.003579540 4 36864 0.4766418E+02 0.002386360 5 91186 0.1465585E+03 0.007337616 6 502355 0.1944353E+05 0.973462043 30 51579 0.1848990E+03 0.009257175 Total elements: 798720 Total Volume: 0.1997359E+05 from file write_exo.out.txt grep SS_ write_exo.out.txt | grep ELEM EXODUSII side sets: 8 SET 1 SS_bottom HAS 6144 Faces SET 2 SS_top HAS 6144 Faces SET 3 SS_east_right HAS 4160 Faces SET 4 SS_north_back HAS 8320 Faces SET 5 SS_west_left HAS 4160 Faces SET 6 SS_south_front HAS 4160 Faces SET 7 SS_corner_side HAS 10 Faces SET 8 SS_corner_top HAS 10 Faces *************************************************** Archive Location Work Directory: /scratch/sft/tam/grid_gen/arctic/areab_fine_deep/ Archive Directory: areab_fine_deep_poly01 ***************************************************