kmaterna
commented
9 months ago Hi, thanks for reporting this! I just tested it in a newly-created environment using your config file, and surprisingly I didn't reproduce the bug. Mine worked all the way through outputs. I did reproduce the pgk_resources DeprecationWarning, so I'll work on that in a bit. But I think the IEEE-float/string-to-float error and the DeprecationWarning are separate issues.
The main issue could be deeper. What operating system are you working on, and what did you use to compile visr.f? We might even have different versions of visr.exe or visr.f -- do you think that source code changed recently? If it's an issue with visr fortran code, I don't know how useful I'll be at debugging....
One way to debug: based on the screen output, it looks like you have the main visr output file, Output/visr/strain_output.txt. The rest of my code is just reading/re-packaging the contents of this file. On my system, this file has a ton of strain values, and doesn't have any occurrences of the string '** ' which your system is reading. This suggests to me that visr didn't finish properly on your system, and produced broken output. Could you also upload that output file here?
I wonder also if it might work on a different dataset, just breaking here due to a floating point issue on one specific computer architecture...? I hope that's not the answer.
fengnian-chang
Hi, Great work! But I encountered an error when applying the vsir method to the test dataset, it says:
Hello! We are... Computing strain using : visr Input data from : ../test/testing_data/NorCal_stationvels.txt Calculation range : [-125.0, -120.0, 38.0, 42.0] Putting the outputs : Output/visr/
Reading file ../test/testing_data/NorCal_stationvels.txt 284 stations before applying cleaning. 233 stations after imposing bounding box.
Computing strain via Visr method. VISR executable found at ../contrib/visr/visr.exe Calling visr.exe fortran code to compute strain: ../contrib/visr/visr.exe < visr_strain.drv Note: The following floating-point exceptions are signalling: IEEE_INVALID_FLAG IEEE_DIVIDE_BY_ZERO IEEE_UNDERFLOW_FLAG IEEE_DENORMAL Traceback (most recent call last): File "/home/gmtsar-user/anaconda3/envs/Strain_2D/bin/strain_rate_compute.py", line 4, in
import('pkg_resources').run_script('Strain-2D==1.1.1', 'strain_rate_compute.py')
File "/home/gmtsar-user/anaconda3/envs/Strain_2D/lib/python3.11/site-packages/pkg_resources/init.py", line 722, in run_script
self.require(requires)[0].run_script(script_name, ns)
File "/home/gmtsar-user/anaconda3/envs/Strain_2D/lib/python3.11/site-packages/pkg_resources/init.py", line 1561, in run_script
exec(code, namespace, namespace)
File "/home/gmtsar-user/anaconda3/envs/Strain_2D/lib/python3.11/site-packages/Strain_2D-1.1.1-py3.11.egg/EGG-INFO/scripts/strain_rate_compute.py", line 11, in
internal_coordinator.strain_coordinator(MyParams);
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
File "/home/gmtsar-user/anaconda3/envs/Strain_2D/lib/python3.11/site-packages/Strain_2D-1.1.1-py3.11.egg/strain/internal_coordinator.py", line 26, in strain_coordinator
[Ve, Vn, rot, exx, exy, eyy, vels, resids] = constructed_object.compute(velField); # computing strain
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
File "/home/gmtsar-user/anaconda3/envs/Strain_2D/lib/python3.11/site-packages/Strain_2D-1.1.1-py3.11.egg/strain/models/strain_visr.py", line 27, in compute
[Ve, Vn, rot_grd, exx_grd, exy_grd, eyy_grd] = compute_visr(myVelfield, self._strain_range, self._grid_inc,
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
File "/home/gmtsar-user/anaconda3/envs/Strain_2D/lib/python3.11/site-packages/Strain_2D-1.1.1-py3.11.egg/strain/models/strain_visr.py", line 80, in compute_visr
[Ve, Vn, rot, exx, exy, eyy] = make_output_grids_from_strain_out(strain_output_file, xdata, ydata);
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
File "/home/gmtsar-user/anaconda3/envs/Strain_2D/lib/python3.11/site-packages/Strain_2D-1.1.1-py3.11.egg/strain/models/strain_visr.py", line 178, in make_output_grids_from_strain_out
rotation.append(float(line[53:60]));
^^^^^^^^^^^^^^^^^^
ValueError: could not convert string to float: '** '
The config file was attached. 05_strain_config_visr.txt
BTW, the DeprecationWarning error was reported when I ran 'python setup.py install'.