========================================================== Version 00 mesh Surfaces in box shape using relative elevation (depth from DEM) use ev_surfs_rel/data_rel_mat*zero*inp interpolate each ev_grid data onto new template box to make a new set of elevations grids Surface 08 is the resampled DEM slightly smaller than original to use as top surface. FUll TOP ==> vers00/08_surf_zero_m.inp <== copied to /interp_iso_surfs_full 17871 17600 4 0 0 000001 9.207226614650E+00 2.750035658156E+00 2.315332795154E+03 SUBSET TOP ==> /interp_iso_surfs_box/08_surf_zero_m.inp <== 8645 8448 4 0 0 00001 1.323923184000E+01 1.660855200000E+01 2.315845438369E+03 SUBSET TOP buffer ==> /interp_iso_surfs_box/08_surf_zero_m_p05.inp <== 8645 8448 4 0 0 00001 1.323923184000E+01 1.660855200000E+01 2.316674000000E+03 Full Surfaces at DEM resolution are: # Grid_size: 301 x 141 # Grid_space: 1613800.000000,1614100.000000,1771230.000000,1771370.000000 Full Surfaces at MESH resolution are: translated and resampled near to zero relative to full DEM # Grid_size: 161 x 111 # Grid_space: 1613900.000000,1614060.000000,1771240.000000,1771350.000000 min x, y 9.207226615 57.97522661 max x, y 2.750035658 36.27803566 length x = 48.768 length y = 33.528 DEM Original Coordinates min x, y are 1.613869792564E+06 1.771230977573E+06 (491907.513 539871.202 meters) max x, y are 1.614098792564E+06 1.771372977573E+06 (491977.312 539914.484 meters) length x = 229 feet (69.97992 meters) length y = 142 feet (43.2816 meters) DEM Translated Coordinates have min x,y at 0,0 max x, y are 229.00 142.00 (69.79 43.28 meters) DIRECTORIES: EARTHVISION surfaces interpolated from relative elevations /ev_surfs_rel relative elevations computed for each surface using ev in /ev_surfs_rel Data mat07 is bottom surface, mat01 is top surface below the ground Input Data hole relative elevations # hole data points about 12 feet apart # bounding box 1613980, 1614042 x 1771266, 1771312 # layers are between .32 and 3.96 feet thick ==> data_rel_mat01.pts <== 1614042.2655 1771265.46392 0.664219119849 -rw-r--r-- 1 tam radiusd 532 Sep 12 13:00 data_rel_mat01.pts Earthvision calculated surface at DEM resolution -rw-r--r-- 1 tam radiusd 1607318 Sep 12 13:16 data_rel_mat01_surf.dat -rw-r--r-- 1 tam radiusd 4818043 Sep 18 14:11 data_rel_mat01_surf.inp # full dem size grids # run ev utilities on linux or sun for earthvision license # ev_2grid creates earthvision binary file # Calculated by ev_2grid with normal minimum tension algorithm # convert to ascii with ev_export # use lagrit to write AVS inp files # 301 x 141 = 42441 ev_2grid -o "tmp.2grd" -r 1613800,1614100,1771230,1771370 -s 301,141 -d $f -G $f:r.out ev_export -o $f:r_surf.dat tmp.2grd lagrit < mksurf_ev.lgi # Grid_size: 301 x 141 # Grid_space: 1613800.000000,1614100.000000,1771230.000000,1771370.000000 ==> data_rel_mat01_surf.dat <== -613800 1771230 .4198848605156 1 1 Earthvision calculated surface at MESH resolution Negative depths to subtract corrected to zero -rw-r--r-- 1 tam radiusd 667992 Sep 18 14:12 data_rel_mat01_surf_sub.dat -rw-r--r-- 1 tam radiusd 2024883 Sep 21 11:00 data_rel_mat01_surf_sub.inp # Gaid_size: 161 x 111 # Grid_space: 1613900.000000,1614060.000000,1771240.000000,1771350.000000 ==> data_rel_mat01_surf_sub.dat <== 1613900 1771240 -.288172125816 1 1 ==> head data_rel_mat01_surf_sub.inp <== 17871 17600 0 0 0 000001 1.613900000000E+06 1.771240000000E+06 0.000000000000E+00 Quad template with Interpolated DEM ==> dem_sub_quad.inp <== 17871 17600 000001 1.613900000000E+06 1.771240000000E+06 7.596236204574E+03 Original coordinates in feet Quad elev subtracted from DEM quad ==> data_rel_mat01_surf_sub_elev.inp <== 17871 17600 000001 1.613900000000E+06 1.771240000000E+06 7.596236204574E+03 FULL Surface Stack into Mesh /interp_iso_surfs_full copied from work in /grid/vers00 Use /ev_surfs_rel/data_rel_mat0*_surf_sub_elev.inp Rename so mat01 is top surf 08 translate to 0,0 convert to meters from data_rel_mat01_surf_sub_elev.out cmo printatt cmoiso -xyz- minmax ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH xic 1.613900000E+06 1.614060000E+06 1.600000000E+02 17871 yic 1.771240000E+06 1.771350000E+06 1.100000000E+02 17871 zic 7.589187858E+03 7.601478891E+03 1.229103329E+01 17871 depth 0.000000000E+00 1.301148772E+00 1.301148772E+00 17871 SUBSET Surfaces /interp_iso_surfs_box (with same names as full) Numerical Mesh in /grid/vers00_box FEHM files in fehm_vers00_box Top surface of model in fehm_vers00_top2d ========================================================== Data Full size ev surfaces from EV are in /ev_surfs_rel These are using template larger than DEM and relative elevations Original DEM and scattered hole data are in /data The DEM is used to create surfaces based on relative elevations of hole data The hole data are represented as scattered x,y,z points where z is actual or relative The surfaces use the same template x,y as used in the original DEM data/dem_1191.inp are quads with the elevations from sheetij in feet data/dem_1191_tri.inp are triangles in feet min x, y are 1.613869792564E+06 1.771230977573E+06 (491907.513 539871.202 meters) max x, y are 1.614098792564E+06 1.771372977573E+06 (491977.312 539914.484 meters) length x = 229 feet (69.97992 meters) length y = 142 feet (43.2816 meters) number of nodes = 32890 number of elements = 65036 xic 1.613869793E+06 1.614098793E+06 2.290000000E+02 32890 yic 1.771230978E+06 1.771372978E+06 1.420000000E+02 32890 zic 7.588093000E+03 7.604049000E+03 1.595600000E+01 32890 scattered data elevations (these are depth in feet calculated from thickness in cm): NOTE: mat01 are the materials at highest elevation, mat07 are lowest (the interpolated surfaces are numbered from bottom to top 1 to 8) layer_data_rel* are the point data with relative (depth) from ground layer_data_abs* are the point data with estimated absolute elevations dem_1191* are the dem elevations on 1 foot spacing dem_box* are the dem elevations at subset box around the data The x and y coordinates are translated to zero to avoid round-off errors. The x, y translation is based on DEM coordinates feet 1613869.7925636, 1771230.9775733 meters 491907.513, 539871.202 DEM in data directory header >>> ncols 230 nrows 143 xllcorner 1613869.2925636 yllcorner 1771230.4775733 cellsize 1 NODATA_value -9999 ******************************************** * LAGRIT read DEM * instead of lower corner use cell center define XTRANS 1613869.2925636 + .5 = 1613869.7925636 define YTRANS 1771230.4775733 + .5 = 1771230.9775733 ******************************************** define nx 230 define ny 143 define minx 1613869.2925636 define miny 1771230.4775733 * add .5 to x and y define minx 1613869.7925636 define miny 1771230.9775733 define NL 6 define dx 1.0 define dy 1.0 read/sheetij/ARCFile/nx,ny/minx,miny/dx,dy/skip NL yflip ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH xic 1.613869793E+06 1.614098793E+06 2.290000000E+02 32890 yic 1.771230978E+06 1.771372978E+06 1.420000000E+02 32890 zic 7.588093000E+03 7.604049000E+03 1.595600000E+01 32890 trans 1,0,0 XTRANS YTRANS 0./0. 0. 0. cmo/addatt/cmo/transauto/VDOUBLE/scalar/vector/user/permanent/x/0.0 ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH xic 0.000000000E+00 2.290000000E+02 2.290000000E+02 32890 yic 0.000000000E+00 1.420000000E+02 1.420000000E+02 32890 zic 7.588093000E+03 7.604049000E+03 1.595600000E+01 32890 * convert to meters ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH xic 0.000000000E+00 6.979920000E+01 6.979920000E+01 32890 yic 0.000000000E+00 4.328160000E+01 4.328160000E+01 32890 zic 2.312850746E+03 2.317714135E+03 4.863388800E+00 32890 ========================================================== PROJECT INFO: Date: Thu, 14 Aug 2008 08:50:42 -0600 From: David Moulton Reply-To: David Moulton To: gaz@lanl.gov Hi George, I'm attaching the paper "Runoff from a semiarid ponderosa pine hillslope in New Mexico", B.P. Wilcox et al (1997). Figure 5 shows the layer structure at 11 of the neutron probe locations. I've extracted the thickness of these layers and entered them in a table along with the relative locations in (x,y) of the probes. Hopefully this is enough to give you a starting point to play around with. I'm going to meet Gary Langhorst next Tuesday to get some help with the DEM data, and we are going to check for data that Brent had on the system out there. If you wanted to touch base with Brent to see if he can send us data, or can tell us where to look that would be cool. In the next e-mail I'll send some other papers that are part of this study and important for the lateral subsurface modeling. Data from these may be useful too. Cheers, David Tube X (in) Y (in) X (m) Y(m) A A+E A+Bw Bw Bt CB+Btb CB 1601 393.71 1027.22 10 26.09 20.25 0.00 0.00 0.00 22.70 0.00 50.31 1602 -52.86 1181.34 -1.34 30.01 20.25 0.00 0.00 0.00 40.49 0.00 60.12 1603 -207.07 1214.83 -5.26 30.86 25.31 0.00 0.00 0.00 24.69 0.00 59.88 1604 -232.56 1620.91 -5.91 41.17 0.00 0.00 30.06 0.00 30.06 0.00 60.12 1605 62.47 1519.06 1.59 38.58 9.76 0.00 0.00 20.12 29.27 0.00 50.00 1606 350.8 1421.23 8.91 36.1 9.76 0.00 0.00 19.51 28.66 0.00 51.22 1607 224.74 1466.79 5.71 37.26 0.00 0.00 40.24 0.00 34.76 14.63 0.00 1608 -78.34 1574.01 -1.99 39.98 0.00 0.00 39.88 0.00 40.49 0.00 9.82 1609 -362.64 1280.51 -9.21 32.52 0.00 0.00 30.06 0.00 44.79 0.00 19.63 1610 96 1122.37 2.44 28.51 0.00 30.67 0.00 0.00 39.26 0.00 49.69 1611 223.4 1084.84 5.67 27.56 0.00 0.00 30.49 0.00 29.88 0.00 0.00 ========================================================== MESH INFO: FEHM 3D mesh /scratch/nts/tam/ponderosa/fehm_vers00_box Mesh Object name: cmotet number of nodes = 143124 number of elements = 824160 dimensions geometry = 3 element type = tet dimensions topology = 3 4 nodes 4 faces 6 edges -rw-r--r-- 1 tam pf 43952418 Sep 24 06:54 tet_asciic.fehmn -rw-r--r-- 1 tam pf 57677083 Sep 24 06:57 tet_asciic.stor -rw-r--r-- 1 tam pf 1574554 Sep 24 06:54 tet_asciic_material.zone -rw-r--r-- 1 tam pf 8793754 Sep 24 06:57 tet_asciic_multi_mat.zone -rw-r--r-- 1 tam pf 250663 Sep 24 06:57 tet_asciic_outside.zone -rw-r--r-- 1 tam pf 1446196 Sep 24 06:57 tet_asciic_outside.area -rw-r--r-- 1 tam pf 250663 Sep 24 06:57 tet_outside.zone -rw-r--r-- 1 tam pf 1446196 Sep 24 06:57 tet_outside.area *********dump_material_lists******** Minimum material ID value = 1 Maximum material ID value = 9 Total possible materials = 9 Material R 1 has 17290 nodes. #nodes/nnodes is 0.120804337498 Material CB 2 has 18235 nodes. #nodes/nnodes is 0.127407003717 Material CB+Btb 3 has 8243 nodes. #nodes/nnodes is 0.575934154998E-01 Material Bt 4 has 25388 nodes. #nodes/nnodes is 0.177384645482 Material BW 5 has 15354 nodes. #nodes/nnodes is 0.107277605433 Material A+BW 6 has 17856 nodes. #nodes/nnodes is 0.124758950281 Material A+E 7 has 3582 nodes. #nodes/nnodes is 0.250272490987E-01 Material A 8 has 11242 nodes. #nodes/nnodes is 0.785472736927E-01 Material top 9 has 25934 nodes. #nodes/nnodes is 0.181199519298 Face Area Area_x Area_y Area_z Face top 1 has 8645 nodes. 0.8590038E+03 0.4381603E+02 0.6617173E+02 0.7848446E+03 Face bottom 2 has 8645 nodes. 0.9109310E+03 0.7540813E+02 0.1019969E+03 0.7848446E+03 Face left_w 3 has 985 nodes. 0.9586601E+02 0.8929169E+02 0.2743154E+02 0.5945772E+01 Face right_e 5 has 815 nodes. 0.9427740E+02 0.8929242E+02 0.2743183E+02 0.5945758E+01 Face back_n 6 has 1719 nodes. 0.1942835E+03 0.5553294E+02 0.1841655E+03 0.1226311E+02 Face front_s 4 has 1963 nodes. 0.1940093E+03 0.5553418E+02 0.1841655E+03 0.1226317E+02 NAME MIN MAX DIFFERENCE LENGTH xic 1.323923184E+01 5.735634960E+01 4.411711776E+01 143124 yic 5.181600000E+00 3.531229920E+01 3.013069920E+01 143124 zic 2.311900000E+03 2.316674000E+03 4.774000000E+00 143124 epsilonl, epsilonaspect: 1.1909894E-11 1.6893659E-33 -------------------------------------------- elements with aspect ratio < .01: 23038 elements with aspect ratio b/w .01 and .02: 4583 elements with aspect ratio b/w .02 and .05: 43196 elements with aspect ratio b/w .05 and .1 : 469143 elements with aspect ratio b/w .1 and .2 : 42344 elements with aspect ratio b/w .2 and .5 : 96688 elements with aspect ratio b/w .5 and 1. : 145168 min aspect ratio = 0.2408E-05 max aspect ratio = 0.9974E+00 epsilonvol: 1.4090908E-09 --------------------------------------- element volumes b/w 0.2715E-07 and 0.5130E-06: 48 element volumes b/w 0.5130E-06 and 0.9694E-05: 7940 element volumes b/w 0.9694E-05 and 0.1832E-03: 521936 element volumes b/w 0.1832E-03 and 0.3462E-02: 103818 element volumes b/w 0.3462E-02 and 0.6542E-01: 190418 element volumes b/w 0.3462E-02 and 0.6542E-01: 190418 min volume = 2.7147114E-08 max volume = 6.5420364E-02 ----------------------------------------------------------- 824160 total elements evaluated. ***Construct and Compress Sparse Matrix:3D*** Matbld3d_stor: *****Negative Coefficients ****** Matbld3d_stor: Total Number of Negative Coefficients 436 Matbld3d_stor: Number of Significant Negative Coefs 436 Matbld3d_stor: Number of 'zero' (< 10e-8*max) coefs 0 Negative coef at row 932 Column 9577 with value -5.9319523E-05 Negative coef at row 1065 Column 9577 with value -4.9899054E-04 Negative coef at row 1065 Column 9710 with value -1.4448471E-04 Negative coef at row 124484 Column 132995 with value -1.7818098E-05 Negative coef at row 125306 Column 125439 with value -1.1733053E-05 Negative coef at row 132221 Column 140733 with value -4.4194789E-06 Matbld3d_stor: npoints = 143124 n connections = 2079350 Matbld3d_stor: Number of written coefs = 395418 Matbld3d_stor: Maximum num. connections to a node = 27 Matbld3d_stor: Volume min = 2.3214395E-04 Matbld3d_stor: Volume max = 2.0745765E-01 Matbld3d_stor: Total Volume: 3.7468481E+03 Matbld3d_stor: abs(Aij/xij) min = 0.0000000E+00 Matbld3d_stor: abs(Aij/xij) max = 9.5458649E+00 Matbld3d_stor: (Aij/xij) max = 4.7193920E-03 Matbld3d_stor: (Aij/xij) min = -9.5458649E+00 "mesh_info.txt" 100 lines, 5085 characters voronoi% cat mesh_info.txt /scratch/nts/tam/ponderosa/fehm_vers00_box -rw-r--r-- 1 tam pf 43952418 Sep 24 06:54 tet_asciic.fehmn -rw-r--r-- 1 tam pf 57677083 Sep 24 06:57 tet_asciic.stor -rw-r--r-- 1 tam pf 1574554 Sep 24 06:54 tet_asciic_material.zone -rw-r--r-- 1 tam pf 8793754 Sep 24 06:57 tet_asciic_multi_mat.zone -rw-r--r-- 1 tam pf 250663 Sep 24 06:57 tet_asciic_outside.zone -rw-r--r-- 1 tam pf 1446196 Sep 24 06:57 tet_asciic_outside.area -rw-r--r-- 1 tam pf 250663 Sep 24 06:57 tet_outside.zone -rw-r--r-- 1 tam pf 1446196 Sep 24 06:57 tet_outside.area Sep 23 17:59 tet_color.gmv -> /ponderosa/grid/vers00_box/tet_color.gmv Sep 24 06:51 tet_color.inp -> /ponderosa/grid/vers00_box/tet_color.inp Mesh Object name: cmotet number of nodes = 143124 number of elements = 824160 dimensions geometry = 3 element type = tet dimensions topology = 3 4 nodes 4 faces 6 edges *********dump_material_lists******** Minimum material ID value = 1 Maximum material ID value = 9 Total possible materials = 9 Material R 1 has 17290 nodes. #nodes/nnodes is 0.120804337498 Material CB 2 has 18235 nodes. #nodes/nnodes is 0.127407003717 Material CB+Btb 3 has 8243 nodes. #nodes/nnodes is 0.575934154998E-01 Material Bt 4 has 25388 nodes. #nodes/nnodes is 0.177384645482 Material BW 5 has 15354 nodes. #nodes/nnodes is 0.107277605433 Material A+BW 6 has 17856 nodes. #nodes/nnodes is 0.124758950281 Material A+E 7 has 3582 nodes. #nodes/nnodes is 0.250272490987E-01 Material A 8 has 11242 nodes. #nodes/nnodes is 0.785472736927E-01 Material top 9 has 25934 nodes. #nodes/nnodes is 0.181199519298 Face Area Area_x Area_y Area_z Face top 1 has 8645 nodes. 0.8590038E+03 0.4381603E+02 0.6617173E+02 0.7848446E+03 Face bottom 2 has 8645 nodes. 0.9109310E+03 0.7540813E+02 0.1019969E+03 0.7848446E+03 Face left_w 3 has 985 nodes. 0.9586601E+02 0.8929169E+02 0.2743154E+02 0.5945772E+01 Face right_e 5 has 815 nodes. 0.9427740E+02 0.8929242E+02 0.2743183E+02 0.5945758E+01 Face back_n 6 has 1719 nodes. 0.1942835E+03 0.5553294E+02 0.1841655E+03 0.1226311E+02 Face front_s 4 has 1963 nodes. 0.1940093E+03 0.5553418E+02 0.1841655E+03 0.1226317E+02 NAME MIN MAX DIFFERENCE LENGTH xic 1.323923184E+01 5.735634960E+01 4.411711776E+01 143124 yic 5.181600000E+00 3.531229920E+01 3.013069920E+01 143124 zic 2.311900000E+03 2.316674000E+03 4.774000000E+00 143124 epsilonl, epsilonaspect: 1.1909894E-11 1.6893659E-33 -------------------------------------------- elements with aspect ratio < .01: 23038 elements with aspect ratio b/w .01 and .02: 4583 elements with aspect ratio b/w .02 and .05: 43196 elements with aspect ratio b/w .05 and .1 : 469143 elements with aspect ratio b/w .1 and .2 : 42344 elements with aspect ratio b/w .2 and .5 : 96688 elements with aspect ratio b/w .5 and 1. : 145168 min aspect ratio = 0.2408E-05 max aspect ratio = 0.9974E+00 epsilonvol: 1.4090908E-09 --------------------------------------- element volumes b/w 0.2715E-07 and 0.5130E-06: 48 element volumes b/w 0.5130E-06 and 0.9694E-05: 7940 element volumes b/w 0.9694E-05 and 0.1832E-03: 521936 element volumes b/w 0.1832E-03 and 0.3462E-02: 103818 element volumes b/w 0.3462E-02 and 0.6542E-01: 190418 element volumes b/w 0.3462E-02 and 0.6542E-01: 190418 min volume = 2.7147114E-08 max volume = 6.5420364E-02 ----------------------------------------------------------- 824160 total elements evaluated. ***Construct and Compress Sparse Matrix:3D*** Matbld3d_stor: *****Negative Coefficients ****** Matbld3d_stor: Total Number of Negative Coefficients 436 Matbld3d_stor: Number of Significant Negative Coefs 436 Matbld3d_stor: Number of 'zero' (< 10e-8*max) coefs 0 Negative coef at row 932 Column 9577 with value -5.9319523E-05 Negative coef at row 1065 Column 9577 with value -4.9899054E-04 Negative coef at row 1065 Column 9710 with value -1.4448471E-04 Negative coef at row 124484 Column 132995 with value -1.7818098E-05 Negative coef at row 125306 Column 125439 with value -1.1733053E-05 Negative coef at row 132221 Column 140733 with value -4.4194789E-06 Matbld3d_stor: npoints = 143124 n connections = 2079350 Matbld3d_stor: Number of written coefs = 395418 Matbld3d_stor: Maximum num. connections to a node = 27 Matbld3d_stor: Volume min = 2.3214395E-04 Matbld3d_stor: Volume max = 2.0745765E-01 Matbld3d_stor: Total Volume: 3.7468481E+03 Matbld3d_stor: abs(Aij/xij) min = 0.0000000E+00 Matbld3d_stor: abs(Aij/xij) max = 9.5458649E+00 Matbld3d_stor: (Aij/xij) max = 4.7193920E-03 Matbld3d_stor: (Aij/xij) min = -9.5458649E+00 FEHM 2D top /scratch/nts/tam/ponderosa/fehm_vers00_top2d -rw-r--r-- 1 tam pf 4533 Sep 29 11:55 mesh_info.txt -rw-r--r-- 1 tam pf 2864203 Sep 29 11:55 tri_asciic.stor -rw-r--r-- 1 tam pf 1171417 Sep 29 11:55 tri.fehmn -rw-r--r-- 1 tam pf 2955 Sep 29 11:55 tri_interface.zone -rw-r--r-- 1 tam pf 95182 Sep 29 11:55 tri_material.zone -rw-r--r-- 1 tam pf 27132 Sep 29 11:55 tri_multi_mat.zone -rw-r--r-- 1 tam pf 25264 Sep 29 11:55 tri_outside.area -rw-r--r-- 1 tam pf 4472 Sep 29 11:55 tri_outside.zone In mesh_info.txt >> The current-mesh-object(CMO) is: cmotet number of nodes = 8645 number of elements = 16896 dimensions geometry = 3 element type = tri dimensions topology = 2 3 nodes 3 faces 3 edges *********dump_material_lists******** Minimum material ID value = 5 Maximum material ID value = 8 Total possible materials = 4 Material Bw 5 (orange) has 1235 nodes. #nodes/nnodes is 0.142857142857 Material A+Bw 6 (brown) has 2622 nodes. #nodes/nnodes is 0.303296703297 Material A+E 7 (pink) has 734 nodes. #nodes/nnodes is 0.849045691151E-01 Material A 8 (tan) has 4054 nodes. #nodes/nnodes is 0.468941584731 outside face num_nodes left_w Area Area_x Area_y Area_z Face left_w 3 has 65 nodes. 0.1967340E+02 0.1872383E+02 0.5702240E+01 0.1078411E+01 Face right_e 5 has 65 nodes. 0.1964057E+02 0.1870137E+02 0.5731645E+01 0.6671539E+00 Face back_n 6 has 133 nodes. 0.4045200E+02 0.1163733E+02 0.3865222E+02 0.9159316E+00 Face front_s 4 has 133 nodes. 0.4042500E+02 0.1171434E+02 0.3858012E+02 0.1583097E+01 epsilonl, epsilonaspect: 1.1879474E-11 1.6764541E-33 -------------------------------------------- elements with aspect ratio b/w .5 and 1. : 16896 min aspect ratio = 0.8234E+00 max aspect ratio = 0.8368E+00 epsilonvol: 8.4072781E-10 --------------------------------------- element volumes b/w 0.4642E-01 and 0.4656E-01: 9435 element volumes b/w 0.4656E-01 and 0.4670E-01: 6478 element volumes b/w 0.4670E-01 and 0.4684E-01: 749 element volumes b/w 0.4684E-01 and 0.4698E-01: 193 element volumes b/w 0.4698E-01 and 0.4713E-01: 41 element volumes b/w 0.4698E-01 and 0.4713E-01: 41 min volume = 4.6419610E-02 max volume = 4.7126727E-02 ----------------------------------------------------------- 16896 total elements evaluated. xic 1.323923184E+01 5.735634960E+01 4.411711776E+01 8645 yic 5.181600000E+00 3.531229920E+01 3.013069920E+01 8645 zic 2.313765510E+03 2.316613896E+03 2.848386022E+00 8645 Matbldtri: total neg 1816 mincoef= -8.1402466E-03 maxcoef= 1.0189416E+00 Negative coeff2d: sort edge vor_coeff row_node col_node Negative coeff2d: 1 51806 -8.1402466E-03 7471 7605 Negative coeff2d: 2 52734 -8.1402466E-03 7605 7471 Negative coeff2d: 3 52719 -7.4859908E-03 7602 7736 Negative coeff2d: 28 53639 -5.6396757E-03 7734 7868 Negative coeff2d: 29 59137 -5.4503064E-03 8528 8394 Negative coeff2d: 30 58239 -5.4503064E-03 8394 8528 ------------ Matbldtri Number of Nodes = 8645 Max. Connection to a Node = 7 Number of Area Coefficients = 59725 Minimum Voronoi area = 2.3234120E-02 Maximum Voronoi area = 9.4166063E-02 Minimum Voronoi edge length = -8.1402466E-03 Maximum Voronoi edge length = 1.0189416E+00 Total Voronoi area = 7.8683064E+02 Total Triangle area = 7.8683064E+02 I was able to find mesh nodes associated with your 2d tri mesh. but still need to do same for your 3d mesh. Here is the directory and files. /scratch/nts/tam/ponderosa/fehm_vers00_top2d layer_data_xy.inp ===> 11 wells x, y coordinates in feet wellpts_xy.inp ===> 11 wells x, y translated to mesh, meters tri_well_pts.dat ===> list of mesh points nearest to 11 wells tri_well_material.zone ===> zone format for well id 1 thru 11 NOTES ===================================================== For boundary work: /scratch/er/gable/grid_gen/hillslope_runoff/data/well_logs_02 /home/cwg/src/lagrit/work/boundary_components/lagrit_g_linux Create an attribute:id_numb < dotask_command = 'cmo/addatt/' // < > cmo(1:icharlnf(cmo)) // < > '/' // < > 'id_numb' // < > '/vint/scalar/nnodes/linear/permanent/afgx/0/' // Note: this is only implemented for identification of node based boundary component searching. It will not work with the 'element' argument. Hi Carl and George, In my home directory on the EES unix systems you will see a subdirectory ponderosa: /home/moulton/ponderosa In this subdirectory there is dem data in ascii (Arc info format: header with information + a long list of Z values). dem_1991.asc This is a rectangle that surrounds the ponderosa site. The data is in State Plane coordinates with 1ft grid spacing. The layer data from the 11 neutron tubes that I extracted from the paper has been processed into comma delimated columns of the form x,y,z,index Files are layer_data_abs.dat - absolute z in ft (surface was read from probe GPS data) layer_data_rel.dat - relative z in ft (surface is treated as 0) Not sure which is more useful. In locations where layers pinch out, the depth is simply the same as it was in the layer above. This may or may not be the best way to handle this, but if you are fitting surfaces it may make the most sense. Let me know if you have additional questions, and/or would like some further changes in the way the data is tabulated. Cheers, David PS. Once we are rolling with this I will generate the data from the additional probes that we found data on.