Write stor file for FEHM mesh files written in 2002 /scratch/er/stauffer/TA53/Grid tet2_ascii.stor (re-connect Delaunay) AMatbld3d_stor: Matrix compress_eps: 0.1000000E-07 AMatbld3d_stor: Local epsilon: 0.1000000E-14 AMatbld3d_stor: *****Negative Coefficients ****** AMatbld3d_stor: Total Number of Negative Coefficients 3 AMatbld3d_stor: Number of Significant Negative Coefs 3 AMatbld3d_stor: Number of 'zero' (< 10e-8 *max) coefs 13 Negative coef at row 1809 Column 2171 with value -2.9478856E-05 Negative coef at row 12580 Column 12921 with value -2.9478856E-05 Negative coef at row 13212 Column 13936 with value -1.4848462E-04 AMatbld3d_stor: npoints = 23025 ncoefs = 323605 AMatbld3d_stor: Number of unique coefs = 80915 AMatbld3d_stor: Maximum num. connections to a node = 27 AMatbld3d_stor: Volume min = 1.7846156E+01 AMatbld3d_stor: Volume max = 6.0627587E+03 AMatbld3d_stor: Total Volume: 2.5931491E+07 AMatbld3d_stor: abs(Aij/xij) min = 2.1319242E-06 AMatbld3d_stor: abs(Aij/xij) max = 7.6769150E+02 AMatbld3d_stor: (Aij/xij) max = 1.4848462E-04 AMatbld3d_stor: (Aij/xij) min = -7.6769150E+02 AMatbld3d_stor Matrix coefficient values stored as scalar area/distance AMatbld3d_stor Matrix compression used for graph and coefficient values ascii STOR file written with name tet2_ascii.stor FEHM files from tet2 reconnected mesh: -rw-r--r-- 1 tamiller sft 9405897 Jun 7 09:50 tet2_ascii.stor -rw-r--r-- 1 tamiller sft 7048515 Jun 7 09:50 tet2.fehmn -rw-r--r-- 1 tamiller sft 253683 Jun 7 09:50 tet2_material.zone -rw-r--r-- 1 tamiller sft 59789 Jun 7 09:50 tet2_outside.zone -rw-r--r-- 1 tamiller sft 344317 Jun 7 09:50 tet2_outside_vor.area *********dump_material_lists******** Minimum material ID value = 1 Maximum material ID value = 12 Total possible materials = 12 Material Puye 1 has 722 nodes. #nodes/nnodes is 0.313572213054E-01 Material Puye 2 has 1265 nodes. #nodes/nnodes is 0.549402832985E-01 Material Puye 3 has 1805 nodes. #nodes/nnodes is 0.783930495381E-01 Material Basalt 4 has 1805 nodes. #nodes/nnodes is 0.783930495381E-01 Material Guaje 5 has 3249 nodes. #nodes/nnodes is 0.141107484698 Material Qbo 6 has 2527 nodes. #nodes/nnodes is 0.109750270844 Material Qct 7 has 1800 nodes. #nodes/nnodes is 0.781758949161E-01 Material Qtsk 8 has 726 nodes. #nodes/nnodes is 0.315309427679E-01 Material Qbt1g 9 has 2476 nodes. #nodes/nnodes is 0.107535287738 Material Qbt1v 10 has 2693 nodes. #nodes/nnodes is 0.116959825158 Material Qbt2 11 has 2155 nodes. #nodes/nnodes is 0.935939177871E-01 Material Air 12 has 1802 nodes. #nodes/nnodes is 0.782627612352E-01 tet_ascii.stor (original tet mesh) -rw-r--r-- 1 tamiller sft 9405646 Jun 7 09:49 tet_ascii.stor AMatbld3d_stor: Matrix compress_eps: 0.1000000E-07 AMatbld3d_stor: Local epsilon: 0.1000000E-14 AMatbld3d_stor: *****Negative Coefficients ****** AMatbld3d_stor: Total Number of Negative Coefficients 29 AMatbld3d_stor: Number of Significant Negative Coefs 29 AMatbld3d_stor: Number of 'zero' (< 10e-8 *max) coefs 12 AMatbld3d_stor: npoints = 23025 ncoefs = 323563 AMatbld3d_stor: Number of unique coefs = 80945 AMatbld3d_stor: Maximum num. connections to a node = 27 AMatbld3d_stor: Volume min = 1.7846156E+01 AMatbld3d_stor: Volume max = 6.0627587E+03 AMatbld3d_stor: Total Volume: 2.5931490E+07 AMatbld3d_stor: abs(Aij/xij) min = 1.6787165E-06 AMatbld3d_stor: abs(Aij/xij) max = 7.6769150E+02 AMatbld3d_stor: (Aij/xij) max = 2.6412472E-03 AMatbld3d_stor: (Aij/xij) min = -7.6769150E+02 AMatbld3d_stor Matrix coefficient values stored as scalar area/distance AMatbld3d_stor Matrix compression used for graph and coefficient values ascii STOR file written with name tet_ascii.stor NAME MIN MAX DIFFERENCE LENGTH xic 4.996510200E+05 4.999253400E+05 2.743200000E+02 23025 yic 5.396407800E+05 5.399151000E+05 2.743200000E+02 23025 zic 1.779858264E+03 2.124456000E+03 3.445977360E+02 23025 -------------------------------------------- elements with aspect ratio < .01: 8369 elements with aspect ratio b/w .01 and .02: 616 elements with aspect ratio b/w .02 and .05: 4895 elements with aspect ratio b/w .05 and .1 : 5311 elements with aspect ratio b/w .1 and .2 : 13332 elements with aspect ratio b/w .2 and .5 : 71993 elements with aspect ratio b/w .5 and 1. : 27525 min aspect ratio = 0.2014E-06 max aspect ratio = 0.9514E+00 epsilonvol: 5.7579460E-06 --------------------------------------- element volumes b/w 0.2865E-02 and 0.3768E-01: 801 element volumes b/w 0.3768E-01 and 0.4956E+00: 3098 element volumes b/w 0.4956E+00 and 0.6519E+01: 4619 element volumes b/w 0.6519E+01 and 0.8574E+02: 29264 element volumes b/w 0.8574E+02 and 0.1128E+04: 94259 min volume = 2.8645104E-03 max volume = 1.1278047E+03 ----------------------------------------------------------- 132041 total elements evaluated. Source 2002 Mesh: 23025 nodes 131476 tets Tet Mesh in FEHM format: -rwxr-xr-x 1 stauffer sft 7025345 Mar 5 2002 lansce_feet.fehmn (coords in feet) -rwxr-xr-x 1 stauffer sft 7025345 Mar 5 2002 lansce_qbt2new.fehmn (coords in meters) Tet Mesh in AVS format (feet): ==> TA53_FEET_geo <== 0000023025 131476 00 00 00 0000000001 0.163925510E+07 0.177045337E+07 0.583935903E+04 ==> TA53_FEET_geo_NOAIR <== 0000021223 0000119728 000000000 000000000 0000000000 0000000001 1.639255103325E+06 1.770453372600E+06 5.839359033105E+03 We could not find the original mesh directory but this directory appears to have the mesh work: /scratch/er/stauffer/TA53/Grid/Sheets * Surfaces and z extents. "bot" and "top" are flat surfaces created by us. * "bot" clips the last layer before the water table so that the bottom surface * is perpendicular to the Z-axis. "top" is a flat surface created above the * topographical surface. This is because Delaunay connect on points of the * topo surface does not recover the surface well and creates -ve coefficients. * So, we make a layer of air on top of the topography and it seems to do ok. Stack Layers and Write FEHM and stor files: TA53_qbt2moved_mkgrid.in:dump/fehm/lansce_qbt2new Missing File: /scratch/ees6/er/stauffer/TA53/Grid/lansce_qbt2new.gmv as used in GMV attr file Grid_try2.attr Geology defined with shape files stacked into mesh (not used?): MDAH_geology_SPfeet.inp is areah_met.gmv converted to feet FEHM dumped stor file: ==> lansce_qbt2new_FEHM.stor <== FEHM V3.3.0linUbuntu 17-02-23 06/05/20 14:38:36 title:initial state for TA53 model w/qbt2_bot at 6850ft 159524 23025 365099 3 0 0 -1 6.415827861760E+01 1.293016628016E+02 1.302695299608E+02 1.312805975112E+02 1.323089487854E+02 12 Materials in lansce_qbt2new_material.zone probably from stacked surfaces in Sheets directory: stack/layers/avs/ & bot.inp 1 1 & w_tab.inp 2 1 & tpf_bot.inp 3 2 & tb4_bot.inp 4 2 & tb4_top.inp 5 2 & qbog_bot.inp 5 2 & qbof_bot.inp 6 4 & qct_bot.inp 7 4 & qbtt_bot.inp 8 4 & qbt1g_bot.inp 9 6 & qbt1v_bot.inp 10 6 & qbt2_bot_new.inp 11 5 & surface.inp 12 5 & top.inp 13 2 & Source surface files: header.readme ncols 19 nrows 19 xllcorner 1639275 yllcorner 1770475 cellsize 50 NODATA_value -9999 byteorder MSBFIRST -rw-r--r-- 1 stauffer sft 1444 Oct 30 2001 qbof_bot -rw-r--r-- 1 stauffer sft 1444 Oct 30 2001 qbog_bot -rw-r--r-- 1 stauffer sft 1444 Oct 30 2001 qbt1g_bot -rw-r--r-- 1 stauffer sft 1444 Oct 30 2001 qbt1v_bot -rw-r--r-- 1 stauffer sft 1444 Oct 30 2001 qbt2_bot -rw-r--r-- 1 stauffer sft 1444 Oct 30 2001 qbt3_bot -rw-r--r-- 1 stauffer sft 1444 Oct 30 2001 qbtt_bot -rw-r--r-- 1 stauffer sft 1444 Oct 30 2001 qct_bot -rw-r--r-- 1 stauffer sft 1444 Oct 30 2001 surface -rw-r--r-- 1 stauffer sft 1444 Oct 30 2001 tb4_bot -rw-r--r-- 1 stauffer sft 1444 Oct 30 2001 tb4_top -rw-r--r-- 1 stauffer sft 1444 Oct 30 2001 tpf_bot -rw-r--r-- 1 stauffer sft 1444 Oct 30 2001 tpt_bot -rw-r--r-- 1 stauffer sft 1444 Oct 30 2001 w_tab ==================================================================================================== Use pylagrit to read fehm file and write avs file ==================================================================================================== es04:Grid:bsh% cat /home/tamiller/scripts/read_fehm.txt susfty-brine:mesh_2002_files:bsh% ipython In [1]: import pylagrit In [2]: l=pylagrit.PyLaGriT() AttributeError: 'PyLaGriT' object has no attribute 'lagrit_exe' In [3]: !which lagrit /n/swdev/mesh_tools/lagrit/install-Ubuntu-14.04-x86_64-gcc4.8.4/lagrit In [4]: l=pylagrit.PyLaGriT(lagrit_exe='/n/swdev/mesh_tools/lagrit/install-Ubuntu-14.04-x86_64-gcc4.8.4/lagrit') * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * Program: LaGriT V3.200 Linux m64 * * * date_compile: 2016/10/25 Ubu gf * * * Run Time: 2017/Jun 6 12:21:36 * * * Manual: http://lagrit.lanl.gov * * * * * * * * * * * * * * * * * * * * * * * * * * * * * -----oOo----- LaGriT V3 LA-CC-15-069 LaGriT Copyright: This program was prepared by Los Alamos National Security, LLC at Los Alamos National Laboratory (LANL) under contract No. DE-AC52-06NA25396 with the U.S. Department of Energy (DOE). All rights in the program are reserved by the DOE and Los Alamos National Security, LLC. Permission is granted to the public to copy and use this software without charge, provided that this Notice and any statement of authorship are reproduced on all copies. Neither the U.S. Government nor LANS makes any warranty, express or implied, or assumes any liability or responsibility for the use of this software. -----oOo----- Output log file: outx3dgen Command log file: logx3dgen In [5]: mo=l.read l.read l.read_fehm l.read_nonblocking l.read_script l.readline l.readlines In [5]: mo=l.read_fehm('lansce_qbt2new.fehmn') read/temp.inp/mo1 read/temp.inp/mo1 read/avs/temp.inp/mo1 0 0 geniee geniee: mesh has 107494 jtet loops, max cycle length= 11 finish cmo/status/brief The current-mesh-object(CMO) is: mo1 1 Mesh Object name: mo1 number of nodes = 23025 number of elements = 131476 dimensions geometry = 3 element type = qua dimensions topology = 2 4 nodes 4 faces 4 edges boundary flag = 16000000 status = active finish finish In [6]: exit -rw-r----- 1 tamiller tamiller 5692271 Jun 6 12:23 temp.inp