diff --git a/simpful/simpful.py b/simpful/simpful.py
index df11616..23b9011 100644
--- a/simpful/simpful.py
+++ b/simpful/simpful.py
@@ -1023,7 +1023,7 @@ def plot_surface(self, variables, output, detail=40, color_map="plasma", outputf
             C.append(temp)
         C = array(C)
 
-        A,B = meshgrid(A,B)
+        A,B = meshgrid(A,B, indexing="ij")
 
         fig = figure(figsize=(8,6))
         ax = axes(projection='3d')