ValueError: cannot convert float NaN to integer #3

New Issue

asia · 2025-05-05T16:33:19+02:00

asia commented

2025-05-05 16:33:19 +02:00

Stack trace:

File "/net/afscra/people/plgtestisl/insilicolab-dev/wpm-1745839867994-m4tFSz0e/app-sandbox/680798cba646a24438a96aec-10/src/shf_wrapper.py", line 57, in
main(sys.argv)
File "/net/afscra/people/plgtestisl/insilicolab-dev/wpm-1745839867994-m4tFSz0e/app-sandbox/680798cba646a24438a96aec-10/src/shf_wrapper.py", line 52, in main
shf(**vars(args))
File "/net/afscra/people/plgtestisl/insilicolab-dev/wpm-1745839867994-m4tFSz0e/app-sandbox/680798cba646a24438a96aec-10/src/seismic_hazard_forecasting.py", line 360, in main
act_rate, bin_counts, bin_edges, out, pprs, rt, idx, u_e = calc_bins(np.array(datenum_data), time_unit,
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
File "/net/pr2/projects/plgrid/plggeposapps/pythonenvs/sha-env/lib/python3.11/site-packages/igfash/rate.py", line 170, in calc_bins
out, ppr = apply_beast(act_rate, **kwargs) # Apply BEAST
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
File "/net/pr2/projects/plgrid/plggeposapps/pythonenvs/sha-env/lib/python3.11/site-packages/igfash/rate.py", line 42, in apply_beast
ncp_median = int(np.max(result.trend.ncp_median)) #Ensure it's integer
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
ValueError: cannot convert float NaN to integer

Script call (with parameters):
python shf_wrapper.py LGCD_Catalog.mat --mc_file=mccalc_output.mat --pdf_file=PDF.mat --m_file=m.mat --m_select=True --mag_label=Mw --mc=10.85 --m_max=0.1 --m_kde_method=adaptiveKDE --xy_select=True --grid_dim=5000 --xy_win_method=False --rate_select=False --time_win_duration=100.0 --forecast_select=True --forecast_len=100 --time_unit=hours --model=Lasocki2013 --products_string=PGA --verbose=True

Stack trace: > File "/net/afscra/people/plgtestisl/insilicolab-dev/wpm-1745839867994-m4tFSz0e/app-sandbox/680798cba646a24438a96aec-10/src/shf_wrapper.py", line 57, in <module> > main(sys.argv) > File "/net/afscra/people/plgtestisl/insilicolab-dev/wpm-1745839867994-m4tFSz0e/app-sandbox/680798cba646a24438a96aec-10/src/shf_wrapper.py", line 52, in main > shf(**vars(args)) > File "/net/afscra/people/plgtestisl/insilicolab-dev/wpm-1745839867994-m4tFSz0e/app-sandbox/680798cba646a24438a96aec-10/src/seismic_hazard_forecasting.py", line 360, in main > act_rate, bin_counts, bin_edges, out, pprs, rt, idx, u_e = calc_bins(np.array(datenum_data), time_unit, > ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^ > File "/net/pr2/projects/plgrid/plggeposapps/pythonenvs/sha-env/lib/python3.11/site-packages/igfash/rate.py", line 170, in calc_bins > out, ppr = apply_beast(act_rate, **kwargs) # Apply BEAST > ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^ > File "/net/pr2/projects/plgrid/plggeposapps/pythonenvs/sha-env/lib/python3.11/site-packages/igfash/rate.py", line 42, in apply_beast > ncp_median = int(np.max(result.trend.ncp_median)) #Ensure it's integer > ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^ > ValueError: cannot convert float NaN to integer Script call (with parameters): `python shf_wrapper.py LGCD_Catalog.mat --mc_file=mccalc_output.mat --pdf_file=PDF.mat --m_file=m.mat --m_select=True --mag_label=Mw --mc=10.85 --m_max=0.1 --m_kde_method=adaptiveKDE --xy_select=True --grid_dim=5000 --xy_win_method=False --rate_select=False --time_win_duration=100.0 --forecast_select=True --forecast_len=100 --time_unit=hours --model=Lasocki2013 --products_string=PGA --verbose=True`

inputs.zip

205 KiB

Login to join this conversation.

No Label

No Milestone

No Assignees

1 Participants

Notifications

Due Date

The due date is invalid or out of range. Please use the format 'yyyy-mm-dd'.

No due date set.

Dependencies

No dependencies set.

Reference: official-apps/SeismicHazardForecasting#3