[CIG-SHORT] a question
Brad Aagaard
baagaard at usgs.gov
Mon Sep 18 10:56:33 PDT 2017
Yasemin,
Your description indicates that you want a spatial variation in your
friction parameters and you have created a SimpleDB spatial database
file, friction.spatialdb. However, your pylithapp.cfg file is setting
the friction spatial database to a UniformDB. You need to use a SimpleDB.
See Sections 4.5.1 and 4.5.2 on the differences in settings between the
SimpleDB and the UniformDB.
You want:
friction.db_properties = spatialdata.spatialdb.SimpleDB
friction.db_properties.label = Slip-weakening friction parameters
friction.db_properties.query_type = linear
friction.db_properties.iohandler.filename = friction.spatialdb
Also for unitless values, you can use None for units, although m/m
should work.
Regards,
Brad
On 09/18/2017 09:03 AM, Yasemin Korkusuz wrote:
> Hello Prof Aagaard,
>
> As I have written you before, I am using PyLith in order to obtain a 3D
> dynamic earthquake rupture simulation for a piecewise planar fault
> geometry which has a heterogeneous stress field. My case is almost
> similar to SCEC TPV205.
>
> On the other hand, our geometry has not only locking but also creeping
> segments.
> Therefore, /I would like to modify initial friction coefficients for
> some segments. /To do this I created a data file, but I am doing
> something wrong or I am not available to do this. In order to be sure,
> I wanted to ask you. I have added my data file and pylithapp.cfg.
>
> Since friction coefficients are unitless, I have written "m/m" for them.
> Also, the x,y,z components are for my initiation patch (in my new test
> case), as I have tested. before.
>
> The error is like following;
> ......
> File
> "/home/yasemin/softwares/pylith/pylith-2.1.4-linux-x86_64/lib/python2.7/site-packages/spatialdata/spatialdb/UniformDB.py",
> line 122, in _validateParameters
> % (self.label, len(params.values), len(params.data))
> ValueError: Incompatible settings for uniform spatial database '<bound
> method UniformDB.label of <spatialdata.spatialdb.UniformDB.UniformDB;
> proxy of <Swig Object of type 'spatialdata::spatialdb::UniformDB *' at
> 0x2c2e4b0> >>'.
> 'values' and 'data' must be lists of the same size.
> 'values' has size of 4 but 'data' has size of 1.
>
> I will be looking forward to seeing your suggestions.
>
> Best Regards
> Yasemin
>
>
> *Yasemin Korkusuz Öztürk*
> Boğaziçi Üniversitesi
> Kandilli Ras. & Dep. Arş. Ens.
> Jeofizik Anabilim Dalı
> 34684, Çengelköy, İstanbul
> Ofis: +90 216 516 33 78
> Cep: +90 536 942 44 02
More information about the CIG-SHORT
mailing list