search

ydluo / qdyn

A Quasi-DYNamic earthquake simulator
56 stars 28 forks source link

Normal stress coupling needs to be added to partial diffs #73

Open martijnende opened 1 year ago

martijnende commented 1 year ago

In derivs_all.f90:151 the time derivative of the normal stress needs to be added, i.e.:

    dmain_var = ( dtau_per + dtau_dt - sigma*dmu_dtheta*dth_dt - dtau_dP*dP_dt ) &
                     / ( sigma*dmu_dv + pb%zimpedance )

should become

    dmain_var = ( dtau_per + dtau_dt - sigma*dmu_dtheta*dth_dt - dtau_dP*(dP_dt - dsigma_dt) &
                     / ( sigma*dmu_dv + pb%zimpedance )

with dtau_dP = -mu (rename variable to make more sense...).

This only affects 3D simulations with FEAT_STRESS_COUPL = 1.