Abstract The intent of this paper is to encourage improved numerical implementation land models. Our contributions in are two-fold. First, we present a unified framework formulate and implement model equations. We separate the representation physical processes from their solution, enabling use established robust methods solve Second, introduce set synthetic test cases (the laugh tests) evaluate...