Defining Deterministic Scenario with varying fault rupture parameters

[johanna pollet]

Hello all,

I am running deterministic scenarios and I was wondering if there is a way to define and run one scenario with varying fault rupture parameters. 

For example, given the following fault rupture scenario:

<?xml version="1.0" encoding="utf-8"?>

<nrml xmlns:gml="http://www.opengis.net/gml" xmlns="http://openquake.org/xmlns/nrml/0.4">

    <singlePlaneRupture>

        <magnitude>5.8</magnitude> (instead specify a range from 5.5 to 6.5 with a step of 0.1)

        <rake>90</rake> (instead specify a range from 80 to 100 with a step of 5)

        <hypocenter lat="45.5" lon="-73.6" depth="14"/>

        <planarSurface strike="130" dip="90">

            <topLeft lon="-73.64322" lat="45.52540" depth="11.06931"/>

            <topRight lon="-73.55682" lat="45.47458" depth="11.06931"/>

            <bottomLeft lon="-73.64322" lat="45.52540" depth="16.93069"/>

            <bottomRight lon="-73.55682" lat="45.47458" depth="16.93069"/>

        </planarSurface>

    </singlePlaneRupture>

 

</nrml>

Thanks in advance for your help.

Best,

Johanna

[anirudh.rao]

Hi Johanna,

Something like this has been possible since engine-3.10. The alternate rupture models can be provided in a single csv file to run a scenario hazard/damage/risk calculation, as shown in the example case here: https://github.com/gem/oq-engine/tree/engine-3.12/openquake/qa_tests_data/scenario/case_13. Unfortunately, documentation of this feature is somewhat lagging.

Prior discussion on Github is here: https://github.com/gem/oq-engine/issues/5603

Best regards,

Anirudh

[johanna pollet]

Hi Anirudh, 

Thank you for your answer, it helps a lot! 

I've been trying to run a scenario similar to case 13 on the GitHub but haven't been successful so far. It seems there is something wrong with the mesh definition and I keep getting the following error. 

I'm attaching the ruptures.csv file I'm using. Do you have any idea what I am doing wrong? Or is there another thread explaining how to input a fault rupture as a csv file?

Traceback (most recent call last):

  File "c:\program files\openquake engine\python3.6\lib\site-packages\openquake\engine\engine.py", line 346, in run_calc

    calc.run(exports=exports, **kw)

  File "c:\program files\openquake engine\python3.6\lib\site-packages\openquake\calculators\base.py", line 252, in run

    self.result = self.execute()

  File "c:\program files\openquake engine\python3.6\lib\site-packages\openquake\calculators\event_based.py", line 302, in execute

    self._read_scenario_ruptures()

  File "c:\program files\openquake engine\python3.6\lib\site-packages\openquake\calculators\event_based.py", line 257, in _read_scenario_ruptures

    aw = readinput.get_ruptures(oq.inputs['rupture_model'])

  File "c:\program files\openquake engine\python3.6\lib\site-packages\openquake\commonlib\readinput.py", line 600, in get_ruptures

    s1, s2 = mesh.shape[1:]

ValueError: too many values to unpack (expected 2)

Thanks again.

Best

Johanna