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 !