Skip to content

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

@kralovcova

Description

@kralovcova

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

Metadata

Metadata

Assignees

No one assigned

    Labels

    No labels
    No labels

    Type

    No type
    No fields configured for issues without a type.

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions