Author Topic: openquake poe values  (Read 4939 times)

Anne Hulsey

  • Jr. Member
  • **
  • Posts: 70
    • View Profile
openquake poe values
« on: September 21, 2021, 12:24:58 AM »
Hi Kuanshi,

I just discovered that the poe terms in the oq .ini files is for the investigation time, not annual. Could you update line 238 in FetchOpenQuake.py to be (something like):

1 - np.exp(-n_years * (1/RP))

Thanks,
Anne

kuanshi

  • Jr. Member
  • **
  • Posts: 56
    • View Profile
Re: openquake poe values
« Reply #1 on: September 21, 2021, 12:28:46 AM »
Hi Anne,

Sure - thank you so much for catching this!  Yes, I'll modify the FetchOpenQuake with your correction.

Thanks,
Kuanshi

Anne Hulsey

  • Jr. Member
  • **
  • Posts: 70
    • View Profile
Re: openquake poe values
« Reply #2 on: September 21, 2021, 05:53:23 PM »
Great, thanks!

After a little more trial and error, I'm still not sure what the best option is. Here is the line I used for what I described above:

'poes': 1 - np.exp(-scen_info['EqRupture']['TimeSpan']*(1.0/scen_info['EqRupture'].get('ReturnPeriod', 100)))

But it gives a lot of significant digits for 500yr RP in 50 years instead of 0.1, as I would prefer. I can't think of a good way to generically round it to the appropriate number of decimal places (perhaps to one significant digit, e.g., 0.1 and 0.002)? For now I am hardcoding 0.1 into the backend script.