Define the SPDE model

The smoothness parameter was chosen as $\nu = 1$, which means $\alpha = 2$ in a 2D plane.

spde <- inla.spde2.matern(mesh = mesh, alpha = 2, constr = TRUE)