As illustrated in #359, bottomDepth >= 0 for all cells in order to ensure model stability. This is also true for the drying slope case and a general implementation of MPAS-O will require the capability to simulate on a mesh with bottomDepth < 0, indicating land above sea level.

cc @mark-petersen and @sbrus89