Manual data entry gone wrong

CliMA · Aug 28, 2020 · c4f67a8 · c4f67a8
1 parent 64c3d3d
commit c4f67a8
Showing 1 changed file with 6 additions and 6 deletions.
diff --git a/verification/lid_driven_cavity/plot_lid_driven_cavity.py b/verification/lid_driven_cavity/plot_lid_driven_cavity.py
@@ -44,11 +44,11 @@
 ζ_Ghia = {
       100: [nan, 40.0110, 22.5378, 16.2862, 12.7844, 10.4199,  8.69628, 7.43218, 6.57451, 6.13973, 6.18946, 6.82674,   8.22110, 10.7414,  15.6591,  30.7923, nan],
       400: [nan, 53.6863, 34.6351, 26.5825, 21.0985, 16.8900, 13.7040, 11.4537, 10.0545,  9.38889, 9.34599, 9.88979,  11.2018,  13.9068,  19.6859,  35.0773, nan],
-     1000: [nan, 42.1124, 23.8707, 18.3120, 16.0458, 14.8061, 14.1374, 14.0928, 14.8901, 16.8350, 20.2666, 25.4341,  32.2953,   40.5437,  51.0557,  75.5980, nan],
-     3200: [nan, 49.9664, 34.2327, 30.4779, 27.9514, 25.8572, 24.4639, 24.1457, 25.3889, 28.9413, 35.8795, 47.1443,  61.7864,   75.6401,  89.3391, 126.670,  nan],
-     5000: [nan, 56.7091, 41.3394, 38.0436, 35.3504, 33.0486, 31.5791, 31.3793, 33.0115, 37.3609, 45.8622, 60.0065,  77.9509,   91.5682, 103.436,  146.702,  nan],
-     7500: [nan, 61.4046, 50.0769, 46.8901, 43.5641, 40.6123, 38.6951, 38.3834, 40.3982, 45.9128, 56.9345, 75.6334,  98.2364,  111.115,  125.131,  180.927,  nan],
-    10000: [nan, 66.0352, 57.7756, 54.3725, 50.3792, 46.8672, 44.6303, 44.3287, 46.8271, 53.5905, 67.1400, 90.0231, 116.275,   127.928,  145.073,  209.452,  nan]
+     1000: [nan, 42.1124, 23.8707, 18.3120, 16.0458, 14.8061, 14.1374, 14.0928, 14.8901, 16.8350, 20.2666, 25.4341,  32.2953,   40.5437,  51.0557,  75.5980, nan][::-1],
+     3200: [nan, 49.9664, 34.2327, 30.4779, 27.9514, 25.8572, 24.4639, 24.1457, 25.3889, 28.9413, 35.8795, 47.1443,  61.7864,   75.6401,  89.3391, 126.670,  nan][::-1],
+     5000: [nan, 56.7091, 41.3394, 38.0436, 35.3504, 33.0486, 31.5791, 31.3793, 33.0115, 37.3609, 45.8622, 60.0065,  77.9509,   91.5682, 103.436,  146.702,  nan][::-1],
+     7500: [nan, 61.4046, 50.0769, 46.8901, 43.5641, 40.6123, 38.6951, 38.3834, 40.3982, 45.9128, 56.9345, 75.6334,  98.2364,  111.115,  125.131,  180.927,  nan][::-1],
+    10000: [nan, 66.0352, 57.7756, 54.3725, 50.3792, 46.8672, 44.6303, 44.3287, 46.8271, 53.5905, 67.1400, 90.0231, 116.275,   127.928,  145.073,  209.452,  nan][::-1]
 }
 
 def plot_lid_driven_cavity_frame(Re, n):
@@ -122,7 +122,7 @@ def plot_lid_driven_cavity_frame(Re, n):
     plt.savefig(f"lid_driven_cavity_Re{Re}_{n:05d}.png")
     plt.close("all")
 
-Re = 100
+Re = 10000
 ds = xr.open_dataset(f"lid_driven_cavity_Re{Re}.nc")
 
 joblib.Parallel(n_jobs=-1)(