Backstress

nika

I am trying to calculate the backstress, applying md=backstressfrominversion(md);

but I'm getting an error: Reference to non-existent field 'strainrate'.

Error in thomasparams (line 67)

e1=averaging(md,md.results.strainrate.principalvalue1,smoothing)/md.constants.yts; %

convert to s^-1

How can I solve for strain rate to get the field md.results.strainrate?

Thanks

chrisb

In order to calculate backstress, the model must contain strain rates either from observed or modeled velocities. Strain rates can be calculated, for example, by calling md=mechanicalproperties(md,vx,vy), where vx and vy are the velocity components you want to use for the strain rate calculation. This will save results for stress and strainrate in md.results.stress and md.results.strainrate, respectively. Then you can call the backstressfrominversion function.

alain

Hello，I use the function md=mechanicalproperties(md,vx,vy)，but it show mistakes about the input of vx or md.inversion.vx_obs. Could you offer a concrete example？Thanks

mathieumorlighem

Hi Alain,

you can type

help mechanicalproperties

you will see some examples. Basically, you need to provide horizontal velocities:


%   Example:
%      md=mechanicalproperties(md,md.initialization.vx,md.initialization.vy);
%      md=mechanicalproperties(md,md.inversion.vx_obs,md.inversion.vy_obs);

I realized on your other post that there is a bug in the current public branch that will be fixed soon.

Best,

Mathieu