How to use pylith_powerlaw_gendb.py?

Hello,

I would like to use the pylith_powerlaw_gendb.py function to create a spatial database for the power law rheology. In the doc, it says that we can make all the information in a file called powerlaw_gendb.cfg and that we need input spatial databases with the laboratory parameters. But I haven’t found any examples of how to use this function and I have no idea what these files look like.

Thank you for your help!

We do not have an example parameter file for the utility pylith_powerlaw_gendb.py in the PyLith code base. This has been on our todo list and we haven’t gotten around to it. It looks like the utility needs a little updating to work with the power-law bulk rheology implementation in PyLith v3. I will try to create a working example in the next week.

I have updated Step 8 of examples/reverse-2d to demonstrate the use of the pylith_powerlaw_gendb Python script. The example parameter file is in the main branch of the PyLith repository https://github.com/geodynamics/pylith/blob/main/examples/reverse-2d/powerlaw_gendb.cfg. The updated documentation for the example is online at Step 8: Slip on Two Faults and Power-law Viscoelastic Materials — PyLith 3.0.3 documentation. The next PyLith release will include these updates.

Thank you very much to have created an example so quickly ! I have dowloaded the different files needed but when I try the utility function pylith_powerlaw_gendb.py as it is written in the documentation (without changing anything just to check), I get this error message :

– db_state_variables ← ‘spatialdata.spatialdb.UniformDB’
– unrecognized property ‘db_state_variables’

I have tried with the SimpleDB spatial parameter, but I still have the same error …
Do I need to activate a specific thing to use the ‘pylith_powerlaw_gendb.py’ function ?

Thank you in advance !

In order to generate.a spatial databsae file compatible with PyLith v3, I had to update the pylith_powerlaw_gendb script. Use the one in the main branch of the PyLith repository https://github.com/geodynamics/pylith/blob/main/applications/pylith_powerlaw_gendb.

Ok, thank you a lot for your help !