July 28 2011 Version 01 Unstructured Mesh with basalt interface Using triangulated template and stacked layers this is the version with the well screens added Note: all wells look visually connected vertically except perhaps R-53 It seems to be added 7 meters from the original mesh column Note: This mesh has a large number of negative ccoefs mostly along the top along NW slope. But also along some of the pinch out areas. number of nodes = 452618 number of elements = 2582936 Material 1 has 435574 nodes. #nodes/nnodes is 0.962343510000 Material 2 has 17044 nodes. #nodes/nnodes is 0.376564800000E-01 FILES: -rw-r--r-- 1 tam pf 138036308 Jul 29 13:53 tet.fehmn -rw-r--r-- 1 tam pf 135675729 Jul 27 15:34 tet_zero.stor -rw-r--r-- 1 tam pf 4978851 Jul 27 15:32 tet_material.zone -rw-r--r-- 1 tam pf 587362 Jul 27 15:32 tet_multi_mat.zone -rw-r--r-- 1 tam pf 191916 Jul 28 08:32 outside_top.zone_material.zone -rw-r--r-- 1 tam pf 191916 Jul 28 08:32 outside_bottom.zone_material.zone -rw-r--r-- 1 tam pf 12175 Jul 28 08:32 outside_west.zone_material.zone -rw-r--r-- 1 tam pf 23769 Jul 28 08:32 outside_south.zone_material.zone -rw-r--r-- 1 tam pf 37728 Jul 28 08:32 outside_east.zone_material.zone -rw-r--r-- 1 tam pf 30391 Jul 28 08:32 outside_north.zone_material.zone well screen mesh values and all attribute fields this is a list of each top and bottom screen node and associated attributes # lineno imt1 vorvol ccoef dcenter top_zone top_screen dwells id_node z x y where imt1 is material 1 or 2 (basalt) vorvol is voronoi volume for that node ccoef as computed in lagrit for fehm dcenter tags nodes with 1 if they are exact screen center of well zone and screen numbers are defined based on mysql dwells is the distance the well column is from exact screen x,y id_node is the mesh node number as used by fehm -rw-r--r-- 1 tam pf 22315 Jul 28 14:17 well_screens_top.dat -rw-r--r-- 1 tam pf 27040 Jul 28 14:17 well_screens_bot.dat There are scripts and lagrit input files in directory well_work These are used to query data and find zones for all the wells For instance, these are awked from the above full field files -rw-r--r-- 1 tam pf 4568 Jul 29 13:29 tops_mat_elev_node.dat -rw-r--r-- 1 tam pf 5540 Jul 29 13:29 bots_mat_elev_node.dat -rw-r--r-- 1 tam pf 10910 Jul 29 13:39 tops_vorvol_coords.dat -rw-r--r-- 1 tam pf 13232 Jul 29 13:41 bots_vorvol_coords.dat These files simply give zone, node, tag and elev and are used to list the screen tops and bottoms as individual sets of top and bot nodes. A few screens have nodes very close (in *_dbl.dat) ==> bots_node.dat <== bot_zone id_node dcenter elev 2010001 450754 1.00 1194.5100 2010004 450771 1.00 1228.6500 ==> bots_node_dbl.dat <== bot_zone id_node dcenter elev 2010001 450754 1.00 1194.5100 2010001 50521 0.00 1194.0000 ==> tops_node.dat <== top_zone id_node dcenter elev 2010001 448784 1.00 1639.5200 2010004 448801 1.00 1680.0600 ==> tops_node_dbl.dat <== top_zone id_node dcenter elev 2010001 448784 1.00 1639.5200 2010004 448801 1.00 1680.0600 Visualization files -rw-r--r-- 1 tam pf 422588431 Jul 28 08:34 tet_attrib.inp -rw-r--r-- 1 tam pf 134626976 Jul 28 08:34 tet_attrib.gmv -rw-r--r-- 1 tam pf 31819 Jul 28 14:17 well_screens_top.inp -rw-r--r-- 1 tam pf 9384 Jul 28 14:17 well_screens_top.gmv -rw-r--r-- 1 tam pf 38569 Jul 28 14:17 well_screens_bot.inp -rw-r--r-- 1 tam pf 11328 Jul 28 14:17 well_screens_bot.gmv -rw-r--r-- 1 tam pf 26420875 Jul 28 08:48 surfmesh_external.inp -rw-r--r-- 1 tam pf 34586196 Jul 28 08:59 surfmesh_external.gmv -rw-r--r-- 1 tam pf 148340156 Jul 27 15:32 tet_fehm.gmv -rw-r--r-- 1 tam pf 467850327 Jul 27 15:33 tet_fehm.inp -rw-r--r-- 1 tam pf 18946792 Jul 27 15:42 subset_outside_tets.gmv -rw-r--r-- 1 tam pf 155821676 Jul 28 08:32 outside_by_psets.gmv -rw-r--r-- 1 tam pf 318854609 Jul 28 08:33 outside_by_psets.inp -rw-r--r-- 1 tam pf 4079376 Jul 28 08:34 subset_outside_pts.gmv Mesh Info files -rw-r--r-- 1 tam pf 340404 Jul 27 15:37 write_fehm.out -rw-r--r-- 1 tam pf 52627 Jul 28 08:39 outside_by_psets.out -rw-r--r-- 1 tam pf 25482 Jul 28 10:36 surfmesh.out screen zone names described # 7 digit zone code constructed as folows: # ignore letters tagged on to screen numbers # ignore PZ or PZ-E or Piezometer screens # ignore wells labeled with CH (also noted as PZ) # # [1] First digit indicates group: # 2 - water supply wells O,PM,G,LA (all single screen) # 3 - CdV # 4 - R # 5 - TW or Test Well # 6 - DT or Test Well DT # # [2:3] Next 2 digits are from well name except for (O=01,PM=02,G=03,LA=04) # # [4:5] Followed by 2 digit flag value found appended to well name # first digit is optional second number in well name # 01 = a or A # 02 = b or B # 03 = i or I (when not CdV) # 04 = r or R # # [6:7] Ending with 2 digit screen number usually 01-09 # For water supply wells this is extracted from well name # # ie O-4 Single Screen = 2-01-00-04 # ie G-1A Single Screen = 2-03-01-01 # ie R-25 Screen 9 = 4-25-00-09 # ie TW-2 Screen 1a = 5-02-00-01 # ie CdV-R-15-3 Screen 6 = 3-15-30-06 # ie CdV-16-2(i) Screen 1 = 3-16-20-01 # ie CdV-16-2(i)r Screen 1 = 3-16-24-01 # ln -s /scratch/er/tam/er_monitoring_2011/grid_2mat_v01/tet_addwells/tet_clr_wells.inp tet_wells.inp ln -s /scratch/er/tam/er_monitoring_2011/grid_2mat_v01/tet_addwells/tet_clr_wells.gmv tet_wells.gmv *** Compress Area Coefficient Values *** AMatbld3d_stor: *****Negative Coefficients ****** AMatbld3d_stor: Total Number of Negative Coefficients 4151 AMatbld3d_stor: Number of Significant Negative Coefs 3879 AMatbld3d_stor: Number of 'zero' (< 10e-8*max) coefs 272 AMatbld3d_stor: npoints = 452618 n connections = 5026988 AMatbld3d_stor: Number of written coefs = 501583 AMatbld3d_stor: Maximum num. connections to a node = 29 AMatbld3d_stor: Volume min = -1.0457233E+07 AMatbld3d_stor: Volume max = 7.2075968E+06 AMatbld3d_stor: Total Volume: 1.5283947E+11 AMatbld3d_stor: abs(Aij/xij) min = 1.1780068E-04 AMatbld3d_stor: abs(Aij/xij) max = 1.3476710E+04 AMatbld3d_stor: (Aij/xij) max = 2.7407515E+03 AMatbld3d_stor: (Aij/xij) min = -1.3476710E+04 Compression used for graph and coefficient values The area coefficient values were written as scalar area/distance ascii STOR file written with name tet_zero.stor *** SPARSE COEFFICIENT MATRIX _astor SUCCESSFUL ***