Lower Triangle Aug 13 2018 Triangulated 25m Prism Mesh V1, with no sublayers Top surface triangulation created in this directory tri_top.inp Edge lengths arranged around 25m length (DEM resolution is 10m) Mesh directory /scratch/sft/tam/IDEAS/East_River/Lower_Triangle/lagrit_tri_25m_v1 -rw-r--r-- 1 tamiller sft 19663008 Aug 9 17:21 prism_fs3.exo -rw-r--r-- 1 tamiller sft 19664124 Aug 9 17:21 prism_fs6.exo -rw-r--r-- 1 tamiller sft 5728 Aug 9 17:21 prism_single_fs3.exo Title: LAGRIT TO EXODUSII number of dimension: 3 number of nodes: 284874 number of elements: 469935 Elements per layer: 93987 Nodes per layer: 47479 Total stacked surfaces: 6 Total material layers: 5 Each layer is 1 cell thickness (thickness = dz). The second layer should have a thickness of 0.3m (so it reaches a depth of 0.4m) Material layers are ordered from top = 1 to bottom = 5. Layer 1: dz = 0.1, sublayers = 0 Layer 2: dz = 0.3, sublayers = 0 Layer 3: dz = 0.6, sublayers = 0 Layer 4: dz = 8.0, sublayers = 0 Layer 5: dz = 21.0, sublayers = 0 Exodus Mesh Files --------------------------------------------------------- 1) prism_fs6.exo Prism with 6 facesets, exclusive (each set is defined once and is unique) This mesh needs user input for left-right, north, and river outlet. layer1 outlet faces selected at elevation less than 2771 meters. EXPSS loop: 1 Side Set tag: bottom 1 Faces: 93987 2 Side Set tag: top 2 Faces: 93987 3 Side Set tag: right 3 Faces: 1408 4 Side Set tag: north 4 Faces: 1390 5 Side Set tag: left 5 Faces: 2022 6 Side Set tag: outlet extend 6 Faces: 25 3) prism_fs3.exo Prism with 3 facesets, exclusive (each set is defined once and is unique) This is the default for mesh facesets and does not need user input. EXPSS loop: 1 Side Set tag: bottom 1 Faces: 93987 2 Side Set tag: top 2 Faces: 93987 3 Side Set tag: all sides 3 Faces: 4845 4) prism_single_fs3.exo Single Column Prism Title: LAGRIT TO EXODUSII number of dimension: 3 number of nodes: 18 number of elements: 5 1 Side Set tag: bottom 1 Faces: 1 2 Side Set tag: top 2 Faces: 1 3 Side Set tag: sides 3 Faces: 15 Source data ---------------------------------------------------- Source Data From LBNL Zexuan Xu zexuanxu@lbl.gov USGS 10 meter DEM product, not the Lidar DEM which has small differences /scratch/sft/tam/IDEAS/East_River/Lower_Triangle/data USGS_DEM_LT10_XYZ.xyz 296173 points 10m DEM connected into triangle surface cropped at boundary. dem_tri_crop.inp 147920 nodes 294150 triangles NAME MIN MAX DIFFERENCE xic 3.259128385E+05 3.311760044E+05 5.263165929E+03 yic 4.314867400E+06 4.320370711E+06 5.503310380E+03 zic 2.759333252E+03 3.775072021E+03 1.015738770E+03 Input river_line.inp and north_line.inp created by hand. Left is west side of river, right is east side of river, back is the north end. METHOD NOTES ---------------------------------------------------- In General, the facesets can be defined for bottom, top, and sides (no direction). All other facesets need user input. Within pylagrit or other script: Directions for faces can be detirmined using mesh center point. Then use 3 points to form planes across center and into 4 quadrants. User inputs for left-right-front-back directions: User inputs for River elements and outlets: - Edge faces selected at elevation at the river outlet - Lines input from file or click - Coordinate min max for below or above - Geometry Box shapes ------------------------------------------------------------------