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)