On Nov 25, 2013, at 4:23 PM, Zunsheng Jiao wrote:


Hi Terry, Phil,
 
I have a 10 km by 10 km Petrel model (Box) for the CERC project, and tried to use your Lagrit scripts to create a computation mesh for FEHM.
 
First I reformat the Petrel output file to avoid DOS chars.: awk -f reformat.awk Petrel_grid.dat > points.dat;
 
Then use read points.dat into avs (points.inp): lagrit < points.lgi;
 
Then use connect.lgi to collect cells: lagrit < connect.lgi.
At this step, I changed the xmin, xmax, ymin, ymax, zmin, and zmax in the script connect.lgi according the model domain. Also I added  2 zones in the script color_elements.mlgi from 7 zones to 9 zones.    The script is run through, and output file tet_trunc.inp. But the file tet_trunc.inp only show 1 zone.
 
Domain size:
 
Xmin: 420100, Xmax 430100 (meter)
Ymin: 4307100, Ymax 4317100 (meter)
Zmin: -2208, Zmax 109 (elevation meter)
 
Model has 9 zones (formations).
 
Could you please help me to find what I am missing?
 
Petrel grid file petrel_grid.dat is uploaded on your transfer site.
 
Have a great holidays,
 
John
 
 
From: Phil Stauffer [mailto:pheeel.gud@gmail.com] 
Sent: Friday, September 06, 2013 12:25 PM
To: Zunsheng Jiao
Subject: you should have the mesh
 
The one you need to run with is in the SUBSET subdirectory.

The workflow should be easy to follow from the lagrit files.

P

Terry had to fix some issues with material types that were outside of the 1-7 range.
 
Terry A. Miller
LANL EES-16 505-667-8009
MESH Software https://lagrit.lanl.gov/
MESH Projects   https://meshing.lanl.gov

September 5, 2013 1:51:29 PM MDT
/scratch/fwo/tam/wy_petrel_v01/subset_box/

-rw-r--r-- 1 tamiller sftus 5107726 Sep  5 13:11 tet.fehmn
-rw-r--r-- 1 tamiller sftus  189443 Sep  5 13:11 tet_material.zone
-rw-r--r-- 1 tamiller sftus  293106 Sep  5 13:11 tet_multi_mat.zone
-rw-r--r-- 1 tamiller sftus   46897 Sep  5 13:11 tet_outside.zone
-rw-r--r-- 1 tamiller sftus  269893 Sep  5 13:11 tet_outside_vor.area
-rw-r--r-- 1 tamiller sftus 4376108 Sep  5 13:11 tet_fehm.gmv
-rw-r--r-- 1 tamiller sftus 6135860 Sep  5 13:11 tet.stor

/scratch/fwo/tam/wy_petrel_v01/subset_box/images
 

-------- Original Message --------
Subject:	RE: mesh from gslib points
Date:	Wed, 4 Sep 2013 22:39:58 +0000
From:	Zunsheng Jiao <JJiao@uwyo.edu>
To:	Stauffer, Philip <stauffer@lanl.gov>

Phil,
 
We could start with a 50 km X 50 km domain. Four corners for this domain:
 
SW x240000 y4043000
SE  x290000  y4043000
NE x290000 y4093000
NW x240000 y4093000
 
Thanks,
 
John
 
From: Stauffer, Philip [mailto:stauffer@lanl.gov] 
Sent: Wednesday, September 04, 2013 1:56 PM
To: Zunsheng Jiao
Subject: Fwd: mesh from gslib points
 
John,

Terry is working on it. there are some issues. How much of the domain do you need ?  If we trim some of the edges it will be much better computationally.  Not sure how eclipse computes on such a bad mesh.  Not voronoi at all in some places.

P



-------- Original Message --------
Subject:
mesh from gslib points
Date:
Wed, 4 Sep 2013 13:51:30 -0600
From:
Miller, Terry A <tamiller@lanl.gov>
To:
Stauffer, Philip H <stauffer@lanl.gov>
 

Hi Phil,
 
I have created awk and lagrit scripts to correct point values, subset,
then connect into an FEHM mesh.
 
This is a NS slice, 2 elements wide. 30x exaggeration in Z.
I had to cut the north end off where layers go up to the surface and can not connect delaunay.
 
This mesh has a few neg ccoef, but they are on the outside of the mesh so should be ok.
 
I can make a big box shape subset of this in the morning.
 
/scratch/fwo/tam/wy_gslib_v01/subset_NS