Lower Triangle Feb 9 2018 Triangulated 50m Prism Mesh V1, 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_v1 -rw-r--r-- 1 tamiller sft 2.3M Feb 9 10:06 prism_fs3.exo -rw-r--r-- 1 tamiller sft 2.3M Feb 9 10:06 prism_fs6.exo -rw-r--r-- 1 tamiller sft 5.6K Feb 9 10:06 prism_single_fs3.exo -rw-r--r-- 1 tamiller sft 3.1M Feb 9 10:44 sm_50m_v1.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: 34596 number of elements: 55835 Elements per layer: 11167 Nodes per layer: 5766 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: 11167 2 Side Set tag: top 2 Faces: 11167 3 Side Set tag: right 3 Faces: 515 4 Side Set tag: north 4 Faces: 535 5 Side Set tag: left 5 Faces: 755 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: 1815 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 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 ------------------------------------------------------------------