Lower Triangle Aug 13 2018 Triangulated 10m Prism Mesh V1, with no sublayers Top surface triangulation created in this directory tri_top.inp Edge lengths arranged around 10m length (DEM resolution is 10m) Mesh directory /scratch/sft/tam/IDEAS/East_River/Lower_Triangle/lagrit_tri_10m_v1 -rw-r--r-- 1 tamiller sft 90587392 Aug 13 13:22 prism_fs3.exo -rw-r--r-- 1 tamiller sft 90588508 Aug 13 13:22 prism_fs6.exo -rw-r--r-- 1 tamiller sft 5728 Aug 13 13:22 prism_single_fs3.exo Title: LAGRIT TO EXODUSII number of dimension: 3 number of nodes: 1309554 number of elements: 2171480 Elements per layer: 434296 Nodes per layer: 218259 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: 434296 2 Side Set tag: top 2 Faces: 434296 3 Side Set tag: right 3 Faces: 3147 4 Side Set tag: north 4 Faces: 3175 5 Side Set tag: left 5 Faces: 4725 6 Side Set tag: outlet extend 6 Faces: 53 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: 434296 2 Side Set tag: top 2 Faces: 434296 3 Side Set tag: all sides 3 Faces: 11100 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 ------------------------------------------------------------------