search

GeoscienceAustralia / eqrm

Automatically exported from code.google.com/p/eqrm
Other
5 stars 4 forks source link

Problem with using the characteristic recurrence model distribution #105

Open GoogleCodeExporter opened 9 years ago

GoogleCodeExporter commented 9 years ago 
What steps will reproduce the problem?
1.run_type =hazard

2.fault source with:
recurrence_model distribution="characteristic" 

3.Have recurrence min and max mag exactly 0.5 apart.
eg.  recurrence_min_mag="4.5" recurrence_max_mag="5.00" 

What is the expected output? What do you see instead?

If you look at the hazrd curves they are not as expected. see attached image 
"aceh_exceedance"

calc_activities_Characteristic(magnitude,b,m0,mMax):

recurrence_functions.py:265: RuntimeWarning: divide by zero encountered in 
double_scalars
  C=((beta*exp(-beta*(mMax-m0-m1-m2)))*m2)/(1-exp(-beta*(mMax-m0-m2)))
recurrence_functions.py:271: RuntimeWarning: invalid value encountered in 
double_scalars
  ((1-1*exp(-beta*(mMax-m0-m2)))*(1+C)))
recurrence_functions.py:275: RuntimeWarning: invalid value encountered in 
double_scalars
  ((1-1*exp(-beta*(mMax-m0-m2)))*(1+C)))

Please use labels and text to provide additional information.

Original issue reported on code.google.com by Vanessa....@gmail.com on 23 Dec 2012 at 10:55

Attachments:

GoogleCodeExporter commented 9 years ago 
It is not appropriate to use the characteristic recurrence model distribution 
for a recurrence magnitude range of 0.5 or less.

Add a warning message to the log file if a user specifies the recurrence_model 
distribution="characteristic" and sets a recurrence magnitude range of 0.5 or 
less.

Original comment by Vanessa....@gmail.com on 23 Dec 2012 at 11:02

GoogleCodeExporter commented 9 years ago 
Hi Vanessa,
Thanks for submitting the ticket.  I have a question.
Why exactly is a range less than 0.5 a problem?

Original comment by duncan.g...@gmail.com on 6 Jan 2013 at 10:36