dfnWorks icon indicating copy to clipboard operation
dfnWorks copied to clipboard

FEHM area stor file with nan values

Open millerta opened this issue 7 months ago • 1 comments

I encountered a problem with full_mesh_vol_area.stor for FEHM. It was mostly -nan- values. As a workaround I have been using full_mesh.stor

On servers at /home/tamiller/src/dfnWorks/work/dfn_wells

% tail fehm_output_20m/full_mesh_vol_area.stor -nan -nan -nan -nan -nan -nan -nan -nan -nan -nan -nan -nan -nan -nan -nan

See driver.py for the commands in case this is a user error.

millerta avatar Jun 24 '25 14:06 millerta

In order to better understand where code may be going wrong, I created and ran 3 FEHM examples. All 3 full_mesh_vol_area.stor files have bad values. I am running Program: pydfnWorks v2.9.2

My local directory: /home/tamiller/src/dfnWorks/examples fehm_example = master version from github with log_normal stochastic fractures fehm_3rect_nondarcy = user rectangles on normal vectors fehm_example_user_rects = user rectangles from 4_user_rects

tail fehm*/*/*area.stor

==> fehm_3rect_nondarcy/output/full_mesh_vol_area.stor <==
           -nan            -nan            -nan            -nan            -nan
           -nan            -nan            -nan            -nan            -nan
           -nan            -nan            -nan            -nan            -nan

==> fehm_example/test_output/full_mesh_vol_area.stor <==
-4.583627560665e-314 -2.031039468596e-315 0.000000000000e+00 -4.263999204049e-314 -2.992371760261e-313
-5.091534412096e-314 0.000000000000e+00 -4.557089359571e-314 -3.503138378274e-313 -4.263999204049e-314
0.000000000000e+00 -3.972736173935e-314 -1.419631080854e-313 -1.516866271414e-314 -1.011332731999e-314
-4.557089359571e-314 -3.972736173935e-314 0.000000000000e+00

==> fehm_example_user_rects/output/full_mesh_vol_area.stor <==
0.000000000000e+00 -1.889143344219e+247 -2.789221678699e+247 -7.516370983288e+247 -1.889143344219e+247
0.000000000000e+00 -1.563517904820e+247 -2.638501323365e+247 -6.468498517581e+247 -1.563517904820e+247
0.000000000000e+00 -3.047890740362e+247 -3.217799814529e+246 -3.926367068489e+247 -3.047890740362e+247
0.000000000000e+00

Also there is a Warning written during user defined fractures:

4_user_rects/output.log:[2025-05-15 10:08:09] WARNING: Warning: You have defined stopCondition = 1 (P32 program stopping condition) but have no stochastic shape families defined. Automatically setting stopCondition to 0 for use with user defined polygons and nPoly.

millerta avatar Jul 23 '25 18:07 millerta