P0 P1 P2 P3 P4 P1dc P2dc P3dc P4dc RT0 RT1 RT2 RT3 RT4 macro plot int1d int2d solve movemesh adaptmesh trunc checkmovemesh on func buildmesh square Eigenvalue min max imag exec LinearCG NLCG Newton BFGS LinearGMRES catch try intalledges jump average mean load savemesh convect abs si...
+int3d(th1)(-1/dt*convect([u1,u2,u3],-dt,siold[mpirank])*qi1) ) // macro SDR(si2,qi2) (int3d(th1)(1/dt*(si2*qi2)) +int3d(th1)(Df*(Grad(si2)’*Grad(qi2))) +int3d(th1)(Nu[mpirank]*siold[mpirank]/(Ks[mpirank]+siold[mpirank])*IpH* ...