Source code for openquake.sub.tests.slab.rupture_test_pai

import os
import unittest
import shutil

from openquake.sub.slab.rupture import calculate_ruptures
from openquake.sub.create_inslab_nrml import create
from openquake.sub.build_complex_surface import build_complex_surface

BASE_DATA_PATH = os.path.dirname(__file__)




@unittest.skipUnless('OQ_RUN_SLOW_TESTS' in os.environ, 'slow')
class RuptureCreationPAITest(unittest.TestCase):



    def setUp(self):
        """
        """
        relpath = '../data/ini/test_kt_z1.ini'
        self.ini_fname = os.path.join(BASE_DATA_PATH, relpath)
        #
        # prepare the input folder and the output folder
        in_path = os.path.join(BASE_DATA_PATH, '../data/profiles/pai_kt_z1/')
        out_path = os.path.join(BASE_DATA_PATH, '../data/tmp/')
        #
        # cleaning the tmp directory
        if os.path.exists(out_path):
            shutil.rmtree(out_path)
        #
        # first we create the complex surface. We use the profiles used for
        # the subduction in the model for the Pacific Islands
        max_sampl_dist = 10.
        build_complex_surface(in_path, max_sampl_dist, out_path,
                              upper_depth=50, lower_depth=750)




    def test_create(self):
        """
        Test rupture creation
        """
        reff = os.path.join(BASE_DATA_PATH, '../data/ini/')
        calculate_ruptures(self.ini_fname, False, reff)

        label = 'test'
        tmps = '../data/tmp/ruptures.hdf5'
        rupture_hdf5_fname = os.path.abspath(os.path.join(BASE_DATA_PATH,
                                                          tmps))
        output_folder = os.path.join(BASE_DATA_PATH, '../data/tmp/')
        investigation_t = '1.'
        create(label, rupture_hdf5_fname, output_folder, investigation_t)