Hi @justinquinn, that's great! I've updated to the new version (revision 25847) and managed to recompile ISSM. But there was still an error when I try to use frictionschoof
:
[0] ??? Error using ==> ./classes/IoModel.cpp:2866
[0] SetFilePointerToData error message: could not find data with name "md.friction.coefficient" in binary file
The frictionschoof()
class doesn't actually use md.friction.coefficient
(compared to the other sliding laws e.g. Weertman), the friction coefficient is stored in md.friction.C
instead. Either there's something hardcoded in the code to use md.friction.coefficient
, or I missed compiling something. Could you assist in debugging this?
Another thing I noticed was that using the extrude function (e.g. model.extrude(md, 3, 3)
) caused the Schoof m
exponent md.friction.m
numpy array to revert from a shape like (1234, 1)
to (1234,)
. I had to use md.friction.m = md.friction.m[..., np.newaxis]
to change it back before the model would run (but then there's the SetFilePointerToData
error above).