Probabilistic risk issue - vulnerability file

[Lisa Jusufi]

Dear OpenQuake team,

I am having trouble running probability loss assessment. The error I get is regarding the mapping of the exposed assets with the respectable vulnerability curves. The vulnerability functions are the ones from GEM.

I have used the same vulnerability files for scenario losses, and I didn't have a problem running it. 

I would really appreciate if somebody can look at this issue and help me.


PS. I have a problem attaching the hazard file here. I am using the one from ESHM20.

Kindest regards,

Lisa

[Anirudh Rao]

Hi Lisa, 

Could you please report the specific error message you are getting?

Thanks,

Anirudh

[Lisa Jusufi]

2025-10-06T09:19:51.64,INFO,SpawnProcess-1/6340,lisaj@Lisa running C:\Users\lisaj\AppData\Local\Temp\lisaj\calc_464\1 - Probabilistic_risk_annual\job_risk.ini [--hc=427]

2025-10-06T09:19:51.64,INFO,SpawnProcess-1/6340,Using engine version 3.23.2

2025-10-06T09:19:51.91,WARNING,SpawnProcess-1/6340,Using 6 processpool workers

2025-10-06T09:19:52.10,INFO,SpawnProcess-1/6340,Checksum of the inputs: 1940304701 (total size 10.36 MB)

2025-10-06T09:20:03.24,INFO,SpawnProcess-1/6340,Reading C:\Users\lisaj\AppData\Local\Temp\lisaj\calc_464\1 - Probabilistic_risk_annual\Exp_Schools_RNM.xml

2025-10-06T09:20:03.35,INFO,SpawnProcess-1/6340,Read 1_129 assets in 0.11s from C:\Users\lisaj\AppData\Local\Temp\lisaj\calc_464\1 - Probabilistic_risk_annual\Exp_Schools_RNM.csv

2025-10-06T09:20:03.39,INFO,SpawnProcess-1/6340,Inferred exposure mesh in 0.02 seconds

2025-10-06T09:20:03.41,INFO,SpawnProcess-1/6340,Associated 1_129 assets (of 1_129) to 917 sites (of 8_217)

2025-10-06T09:20:03.41,INFO,SpawnProcess-1/6340,Found 43 taxonomies with ~26.3 assets each

2025-10-06T09:20:03.42,INFO,SpawnProcess-1/6340,The most common taxonomy is MUR/LWAL+DNO/H1-R with 233 assets

2025-10-06T09:20:03.46,INFO,SpawnProcess-1/6340,Reducing risk model from 4206 to 36 risk functions

2025-10-06T09:20:03.46,INFO,SpawnProcess-1/6340,minimum_asset_loss={'occupants': 0, 'structural': 0}

2025-10-06T09:20:03.47,INFO,SpawnProcess-1/6340,Storing risk model

2025-10-06T09:20:03.62,INFO,SpawnProcess-1/6340,Building 1_000 realizations

2025-10-06T09:20:03.64,INFO,SpawnProcess-1/6340,Building risk inputs from 1000 realization(s)

2025-10-06T09:20:03.72,INFO,SpawnProcess-1/6340,Building 1_000 realizations

2025-10-06T09:23:56.53,INFO,SpawnProcess-1/6340,Built 1129 risk inputs

2025-10-06T09:23:57.69,WARNING,SpawnProcess-1/6340,Sent 13 classical_risk tasks, 537.07 MB

2025-10-06T09:24:05.92,INFO,SpawnProcess-1/6340,Received 1 * 86 B in 8 seconds [unpik=0.00s] from classical_risk {'tot': '86 B'}

2025-10-06T09:24:06.03,ERROR,SpawnProcess-1/6340,Traceback (most recent call last): File "C:\Users\lisaj\anaconda3\envs\oqe\lib\site-packages\openquake\engine\engine.py", line 208, in run_calc calc.run(shutdown=True) File "C:\Users\lisaj\anaconda3\envs\oqe\lib\site-packages\openquake\calculators\base.py", line 339, in run raise exc from None File "C:\Users\lisaj\anaconda3\envs\oqe\lib\site-packages\openquake\calculators\base.py", line 326, in run self.result = self.execute() File "C:\Users\lisaj\anaconda3\envs\oqe\lib\site-packages\openquake\calculators\base.py", line 1272, in execute return smap.reduce(self.combine, self.acc) File "C:\Users\lisaj\anaconda3\envs\oqe\lib\site-packages\openquake\baselib\parallel.py", line 924, in reduce return self.submit_all().reduce(agg, acc) File "C:\Users\lisaj\anaconda3\envs\oqe\lib\site-packages\openquake\baselib\parallel.py", line 608, in reduce for result in self: File "C:\Users\lisaj\anaconda3\envs\oqe\lib\site-packages\openquake\baselib\parallel.py", line 594, in __iter__ yield from self._iter() File "C:\Users\lisaj\anaconda3\envs\oqe\lib\site-packages\openquake\baselib\parallel.py", line 584, in _iter out = result.get() File "C:\Users\lisaj\anaconda3\envs\oqe\lib\site-packages\openquake\baselib\parallel.py", line 414, in get raise etype(msg) ZeroDivisionError: File "C:\Users\lisaj\anaconda3\envs\oqe\lib\site-packages\openquake\baselib\parallel.py", line 437, in new val = func(*args) File "C:\Users\lisaj\anaconda3\envs\oqe\lib\site-packages\openquake\calculators\classical_risk.py", line 38, in classical_risk crmodel = monitor.read('crmodel') File "C:\Users\lisaj\anaconda3\envs\oqe\lib\site-packages\openquake\baselib\performance.py", line 398, in read return pickle.loads(dset[()]) File "C:\Users\lisaj\anaconda3\envs\oqe\lib\site-packages\openquake\risklib\scientific.py", line 405, in __setstate__ self.init() File "C:\Users\lisaj\anaconda3\envs\oqe\lib\site-packages\openquake\risklib\scientific.py", line 256, in init self._mlr_i1d = interpolate.interp1d(self.imls, self.mean_loss_ratios) File "C:\Users\lisaj\anaconda3\envs\oqe\lib\site-packages\scipy\interpolate\_interpolate.py", line 267, in __init__ _Interpolator1D.__init__(self, x, y, axis=axis) File "C:\Users\lisaj\anaconda3\envs\oqe\lib\site-packages\scipy\interpolate\_polyint.py", line 58, in __init__ self._set_yi(yi, xi=xi, axis=axis) File "C:\Users\lisaj\anaconda3\envs\oqe\lib\site-packages\scipy\interpolate\_polyint.py", line 131, in _set_yi self._y_axis = (axis % yi.ndim) ZeroDivisionError: integer division or modulo by zero

[Anirudh Rao]

I wonder if the job settings in the hazard job file are leading to some kind of saturation of the hazard curve values. Would it be possible to share your job_hazard.ini file too?

[Lisa Jusufi]

Dear Anirudh,

I have a problem attaching the hazard file directly here because its too large.
This is the link to it: https://we.tl/t-Fwggi2ekWH 

Note that I managed to run probabilistic damage with the same hazard file and with the fragility curves of GEM without a problem.

Kind regards,

Lisa

[Anirudh Rao]

Thanks Lisa, we'll try to see if we can reproduce the issue.

[Lisa Jusufi]

Thank you. Looking forward to your assistance.