[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