hstress -solid.pm "Hydrostatic stress"
Cveq Cv0*exp((1-f)*Va*hstress/(k_B_const*T)) "Equilibrium vacancy concentration"
G -((c-Cveq)/taus) "Rosenberg-Ohring recombination"
strate (G-cTIME)*Va*f+(1-f)*G*Va "Strain rate"
PFluxX -e_const*Znum*(Dv/(k_B_const*T))*ec.Ex*c "Flux due to electric field, x direction"
PFluxY -e_const*Znum*(Dv/(k_B_const*T))*ec.Ey*c "Flux due to electric field, y direction"
SFluxX -(Dv/(k_B_const*T))*d(hstress,x)*Va*f*c "Flux due to stress gradient, x direction"
SFluxY -(Dv/(k_B_const*T))*d(hstress,y)*Va*f*c "Flux due to stress gradient, y direction"
TFluxX -(Dv/(k_B_const*T^2))*Tx*Qstar*c "Flux due to temperature gradient, x direction"
TFluxY -(Dv/(k_B_const*T^2))*Ty*Qstar*c "Flux due to temperature gradient, y direction"