M0 M0const*nojac(M0factor) Mobility turning parameter
M0factor phiA^2*phiB^2+terpf.phiC^2*phiA^2+terpf.phiC^2*phiB^2 Help variable
M0const 2e-5[m^3/s] Help variable