fix coriolis precomputed for order > 1

dgpy/scripts/slim.py

@@ -532,7 +532,7 @@ class ShallowWater2d:
             netcdf or .msh file containing the coriolis term for the whole domain. 
         """
         self._coriolis = slim_private._load_function(coriolis, self._domain._groups)
-        self._coriolis_PC = dgpy.functionPrecomputed(self._domain._groups, 3)
+        self._coriolis_PC = dgpy.functionPrecomputed(self._domain._groups, self._domain._order*2+1)
         self._coriolis_PC.compute(self._coriolis)
         self._equation.setCoriolisFactor(self._coriolis_PC)