lanl / dfnWorks

dfnWorks is a parallelized computational suite to generate three-dimensional discrete fracture networks (DFN) and simulate flow and transport. If you download the software please fill out our interest form to stay up to date on releases https://goo.gl/forms/VE39oKsyp4LVC6Gj2 and join our google group https://groups.google.com/d/forum/dfnworks-users . Precompiled Docker Container https://hub.docker.com/r/ees16/dfnworks
https://dfnworks.lanl.gov/
Other
73 stars 42 forks source link

Error using well-package during "Gathering points on fractures" #82

Open kralovcova opened 3 months ago

kralovcova commented 3 months ago

Hello, using well package, sometimes an error occurs during "Gathering points on fractures".

I run the "well example" model. The original model generates 400 fractures. I also tried to increase the fracture density . The problem arises at a higher fracture density. If I increase nPoly to 1000, then processing the well data ends with an error: "IndexError: index 4 is out of bounds for axis 0 with size 4":

--> Working on well inject
--> Interpolating well coordinates into a polyline
well_inject.dat
--> Interpolating well coordinates into a polyline: Complete
--> Writing polyline into avs file: well_inject_line.inp
--> Writing polyline into avs file : well_inject_line.inp : Complete
--> Running: /dfnWorks//bin/lagrit < find_well_inject_segment.lgi -log find_well_inject_segment.log -out find_well_inject_segment.out > find_well_inject_segment.dump
--> LaGriT script find_well_inject_segment.lgi was run successfully
--> Finding well points on DFN for injection
--> Gathering points on fractures
Traceback (most recent call last):
   File "/dfnWorks/work/02_well_example/[driver.py](https://urldefense.com/v3/__http:/driver.py__;!!Bt8fGhp8LhKGRg!CpprKkKV6y8RqxFKcXZT5GypourQmJ2GPGDKfUnoKE2FMiXiQ_YPnnuJ_FMpX-7tdzYr7cxwnCEIT1UQBQnYbdgBIpg$)", line 78, in <module>
     DFN.find_well_intersection_points(wells)
   File "/usr/local/lib/python3.10/dist-packages/pydfnworks/dfnGen/well_package/[wells.py](https://urldefense.com/v3/__http:/wells.py__;!!Bt8fGhp8LhKGRg!CpprKkKV6y8RqxFKcXZT5GypourQmJ2GPGDKfUnoKE2FMiXiQ_YPnnuJ_FMpX-7tdzYr7cxwnCEIT1UQBQnY7mzUkxc$)", line 591, in find_well_intersection_points
     self.run_find_well_intersection_points(well)
   File "/usr/local/lib/python3.10/dist-packages/pydfnworks/dfnGen/well_package/[wells.py](https://urldefense.com/v3/__http:/wells.py__;!!Bt8fGhp8LhKGRg!CpprKkKV6y8RqxFKcXZT5GypourQmJ2GPGDKfUnoKE2FMiXiQ_YPnnuJ_FMpX-7tdzYr7cxwnCEIT1UQBQnY7mzUkxc$)", line 640, in run_find_well_intersection_points
     self.well_point_of_intersection(well)
   File "/usr/local/lib/python3.10/dist-packages/pydfnworks/dfnGen/well_package/[wells.py](https://urldefense.com/v3/__http:/wells.py__;!!Bt8fGhp8LhKGRg!CpprKkKV6y8RqxFKcXZT5GypourQmJ2GPGDKfUnoKE2FMiXiQ_YPnnuJ_FMpX-7tdzYr7cxwnCEIT1UQBQnY7mzUkxc$)", line 850, in well_point_of_intersection
     points = self.gather_points()
   File "/usr/local/lib/python3.10/dist-packages/pydfnworks/dfnGen/meshing/udfm/[map2continuum_helper.py](https://urldefense.com/v3/__http:/map2continuum_helper.py__;!!Bt8fGhp8LhKGRg!CpprKkKV6y8RqxFKcXZT5GypourQmJ2GPGDKfUnoKE2FMiXiQ_YPnnuJ_FMpX-7tdzYr7cxwnCEIT1UQBQnYP7w_mOo$)", line 72, in gather_points
     point = self.polygons[name][j]
IndexError: index 4 is out of bounds for axis 0 with size 4