Lower Triangle Feb 9 2018 Triangulated 50m Prism Mesh V2, with sublayers Top surface triangulation created in this directory tri_top.inp Edge lengths arranged around 50m length Mesh directory /scratch/sft/tam/IDEAS/East_River/Lower_Triangle/lagrit_tri_50m_v2 -rw-r--r-- 1 tamiller sft 6.6M Feb 9 11:11 prism_fs3.exo -rw-r--r-- 1 tamiller sft 6.6M Feb 9 11:12 prism_fs6.exo -rw-r--r-- 1 tamiller sft 6.9K Feb 9 11:10 prism_single_fs3.exo -rw-r--r-- 1 tamiller sft 7.5M Feb 9 11:42 sm_50m_v2.tar.gz prism_fs3.exo prism_fs6.exo prism_single_fs3.exo README.txt images/s*png stack_tri_exc.out.txt See README.txt Title: LAGRIT TO EXODUSII number of dimension: 3 number of nodes: 98022 number of elements: 178672 Elements per layer: 11167 Nodes per layer: 5766 Total stacked surfaces: 17 Total material layers: 5 Material layers are ordered from top = 1 to bottom = 5. Layer 1: dz = 0.1, sublayers = 0, cells 1 Layer 2: dz = 0.1, sublayers = 2, cells 3 Layer 3: dz = 0.6, sublayers = 0, cells 1 Layer 4: dz = 1.0, sublayers = 7, cells 8 Layer 5: dz = 7.0, sublayers = 2, cells 3 Total depth = 30 meters 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: 11167 2 Side Set tag: top 2 Faces: 11167 3 Side Set tag: right 3 Faces: 1659 4 Side Set tag: north 4 Faces: 1712 5 Side Set tag: left 5 Faces: 2427 6 Side Set tag: outlet extend 6 Faces: 10 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: 11167 2 Side Set tag: top 2 Faces: 11167 3 Side Set tag: all sides 3 Faces: 5808 4) prism_single_fs3.exo Single Column Prism Title: LAGRIT TO EXODUSII number of dimension: 3 number of nodes: 51 number of elements: 16 1 Side Set tag: bottom 1 Faces: 1 2 Side Set tag: top 2 Faces: 1 3 Side Set tag: sides 3 Faces: 48 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 TOP surface interpolated values from DEM: NAME MIN MAX DIFFERENCE xic 3.259128385E+05 3.311760044E+05 5.263165900E+03 yic 4.314867400E+06 4.320370711E+06 5.503310400E+03 zic 2.729530188E+03 3.775072021E+03 1.045541834E+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 ------------------------------------------------------------------