updates from Rene Gassmoeller

maxrudolph · Oct 18, 2017 · 6dae0a0 · 6dae0a0
1 parent 3289b33
commit 6dae0a0
Show file tree

Hide file tree

Showing 3 changed files with 570 additions and 416 deletions.
diff --git a/benchmarks/vankeken_subduction/plot_first_timestep.py b/benchmarks/vankeken_subduction/plot_first_timestep.py
@@ -38,7 +38,7 @@
 vk_T = np.flipud(np.loadtxt(vk_T_file))
 
 #output_dir = sys.argv[1]
-output_dir = 'output_case1'
+output_dir = 'output_case6'
 print('Reading output from ',output_dir);
 
 # In[2]:
@@ -71,14 +71,14 @@
 mesh_files = sorted(glob.glob(output_dir + '/solution/mesh*.h5'))
 
 # Load the statistics file
-stats_file = output_dir + '/statistics'
-stats = np.loadtxt(stats_file,usecols=(0,1,2,4,12,15))
-time = stats[:,1]
-tavg = stats[:,4]
-ncell = stats[:,3]
-vrms = stats[:,5]
+#stats_file = output_dir + '/statistics'
+#stats = np.loadtxt(stats_file,usecols=(0,1,2,4,12,15))
+#time = stats[:,1]
+#tavg = stats[:,4]
+#ncell = stats[:,3]
+#vrms = stats[:,5]
 # load the output filenames from the statistics file
-output_files=np.loadtxt(stats_file,usecols=(10),dtype=str)
+output_files=np.loadtxt(stats_file,usecols=(10),type=str)
 output_times=[]
 for i in range(len(output_files)):
     if output_files[i] != '""':