def getLayeredModel(self):
''' Get a layered model, see :class:`pyrocko.cake.LayeredModel`. '''
def iterLines():
for il, m in enumerate(self.iterLayers()):
yield self.d[il], m, ''
return LayeredModel.from_scanlines(iterLines())
def getLayeredModel(self):
''' Get a layered model, see :class:`pyrocko.cake.LayeredModel`. '''
def iterLines():
for il, m in enumerate(self.iterLayers()):
yield self.d[il], m, ''
return LayeredModel.from_scanlines(iterLines())
def PsGrnArray2LayeredModel(psgrn_input_path):
"""
Read PsGrn Input file and return velocity model.
Parameters
----------
psgrn_input_path : str
Absolute path to the psgrn input file.
Returns
-------
:class:`LayeredModel`
"""
a = num.loadtxt(psgrn_input_path, skiprows=136)
b = a[:, 1: -1]
b[:, 3] /= 1000.
return LayeredModel.from_scanlines(
read_nd_model_str(
re.sub('[\[\]]', '', num.array2string(
b, precision=4,
formatter={'float_kind': lambda x: "%.3f" % x}))))