Hi,

I am trying to initial stress in a 2D subduction zone with a triangle mesh. I put the gravity (using the example) in my subduction (same example than in the manual) and now I am trying to do the step8. I already re-run the step3 with commenting out the CellFilterAvg settings and run pylith gravity_initstress.cfg gravity_vardensity.cfg. But when I did ./generate_statedb.py, I have this error:

“Traceback (most recent call last):

File “./generate_statedb.py”, line 63, in

quadCoords = calcQuadCoords(vertices, cells, qpts)

File “./generate_statedb.py”, line 42, in calcQuadCoords

quadCoords[:,:,iDim] = numpy.dot(cellCoords[:,:,iDim], qpts.transpose())

ValueError: shapes (523,3) and (4,4) not aligned: 3 (dim 1) != 4 (dim 0)”

I understood that in the generate_statedb.py file the qpts array are the basis function evaluated at the quadrature points but I still don’t understand how to calculate this value? Furthermore, the example of 2D gravity is for lines and my model is for triangle so I need to say in generate_statedb.py the good size. But I don’t understand how?

Can you help me?

Thank you in advance,

Lauriane Baylé

Phd Student Institut des Sciences de la Terre de Paris

Sorbonne University