Dear MRST developers,
I wanted to leave a comment on your website about the 3D EDFM model's workflow, specifically regarding the step where a 3D fracture surface is transformed into a 2D set of points. I noticed that the operation only involves rotating the points about the origin to align with the xy-plane.
Why only the xy-plane?
Additionally, when my fracture surface is perpendicular to the x-axis, the code encounters the following error: Error using horzcat Dimensions of arrays being concatenated are not consistent. Error in fractureplanegeneralgrid (line 67) F.faces.neighbors(i,:) = [cellNo(neigh),0]; Error in Example_1_manualEDFM (line 164) Fgrid(i).grid = fractureplanegeneralgrid(V,C,points,planenormal,aperture,tol);
Could you please advise on how to resolve this issue?
Thank you for your time and assistance!