Problem changing controlfile parameter csm_hysteretic_damping

The EQRM seems to bomb whenever
csm_hysteretic_damping = None 
OR
csm_hysteretic_damping = trapezoidal

I need to change these parameters to reproduce the Newcastle 1989 scenario work 
we have done in the past. See below for error when using the 
setdata_ScenRisk.py (attached) from the demo  

ERROR MESSAGE: 
Davids-iMac:test8 davidrobinson$ python setdata_ScenRisk.py
Logfile is './EQRM.log' with logging level of DEBUG, console logging level is 
INFO
WARNING: max_width term in EQRM control file is deprecated. Replaced with 
scenario_max_width=15.
Logfile is './output/scen_risk/log-0.txt' with logging level of DEBUG, console 
logging level is INFO
host name: Davids-iMac.local
SVN version: 1241
event_set_handler = generate
P0: Generating event set
P0: Saving event set to ./output/scen_risk/newc_event_set
P0: Event set created. Number of events=1002
P0: Sites set created. Number of sites=6305
P0: do site 1 of 6305
csm_hysteretic_damping None

Traceback (most recent call last):
  File "setdata_ScenRisk.py", line 97, in <module>
    main(locals())
  File "/Users/davidrobinson/Documents/datafiles/workstuff/svn_sandpit/eqrm/eqrm_code/analysis.py", line 441, in main
    damage) = sites.calc_total_loss(SA, eqrm_flags, overloaded_MW)
  File "/Users/davidrobinson/Documents/datafiles/workstuff/svn_sandpit/eqrm/eqrm_code/structures.py", line 269, in calc_total_loss
    loss_aus_contents=eqrm_flags.loss_aus_contents)
  File "/Users/davidrobinson/Documents/datafiles/workstuff/svn_sandpit/eqrm/eqrm_code/damage_model.py", line 184, in aggregated_building_loss
    self.building_loss(ci=ci, loss_aus_contents=loss_aus_contents)
  File "/Users/davidrobinson/Documents/datafiles/workstuff/svn_sandpit/eqrm/eqrm_code/damage_model.py", line 150, in building_loss
    damage_states = self.get_building_states()
  File "/Users/davidrobinson/Documents/datafiles/workstuff/svn_sandpit/eqrm/eqrm_code/damage_model.py", line 111, in get_building_states
    (SA, SD) = self.get_building_displacement()
  File "/Users/davidrobinson/Documents/datafiles/workstuff/svn_sandpit/eqrm/eqrm_code/damage_model.py", line 145, in get_building_displacement
    point = self.capacity_spectrum_model.building_response(self.SA)
  File "/Users/davidrobinson/Documents/datafiles/workstuff/svn_sandpit/eqrm/eqrm_code/capacity_spectrum_model.py", line 137, in building_response
    rtol=rtol,maxits=maxits)
  File "/Users/davidrobinson/Documents/datafiles/workstuff/svn_sandpit/eqrm/eqrm_code/equivalent_linear_solver.py", line 48, in solve
    SA,SD,SAcap,exit_flag=update_function(intersection_x) # update curves
  File "/Users/davidrobinson/Documents/datafiles/workstuff/svn_sandpit/eqrm/eqrm_code/capacity_spectrum_model.py", line 155, in updated_response
    non_linear_damping=self._non_linear_damping(displacement)
  File "/Users/davidrobinson/Documents/datafiles/workstuff/svn_sandpit/eqrm/eqrm_code/capacity_spectrum_model.py", line 176, in _non_linear_damping
    self.kappa,SA,SD,self.csm_hysteretic_damping)
  File "/Users/davidrobinson/Documents/datafiles/workstuff/svn_sandpit/eqrm/eqrm_code/capacity_spectrum_functions.py", line 53, in nonlin_damp
    csm_hysteretic_damping)
  File "/Users/davidrobinson/Documents/datafiles/workstuff/svn_sandpit/eqrm/eqrm_code/capacity_spectrum_functions.py", line 161, in hyst_area_rand
    return Harea
UnboundLocalError: local variable 'Harea' referenced before assignment
Original issue reported on code.google.com by RobinsonDavidJ0@gmail.com on 4 Sep 2012 at 12:23
Attachments:
setdata_ScenRisk.py
GeoscienceAustralia / eqrm

Problem changing controlfile parameter csm_hysteretic_damping #84
