ie

TP-one-patch/mesh_study/R-one-patch-mesh-study.py

@@ -27,26 +27,26 @@ mesh_study = True
 # resolutions = {128: 1e-7} #[1,2,3,4,5,10,20,40,75,100]
 resolutions = {
                 # 1: 1e-7,
-                # 2: 1e-7,
-                # 4: 1e-7,
-                # 8: 1e-7,
-                # 16: 1e-7,
+                2: 1e-7,
+                4: 1e-7,
+                8: 1e-7,
+                16: 1e-7,
                 32: 1e-7,
-                # 64: 1e-7,
-                # 128: 1e-7,
-                # 256: 1e-7,
-                # 512: 1e-7,
+                64: 1e-7,
+                128: 1e-7,
+                256: 1e-7,
+                512: 1e-7,
                 }
 
 ############ GRID #######################
 # mesh_resolution = 20
-timestep_size = 0.00001
-number_of_timesteps = 10
+timestep_size = 0.01
+number_of_timesteps = 70
 plot_timestep_every = 1
 # decide how many timesteps you want analysed. Analysed means, that we write out
 # subsequent errors of the L-iteration within the timestep.
 number_of_timesteps_to_analyse = 5
-starttimes = [0.5]
+starttimes = [0.0,0.25,0.5]
 # starttimes = [0.0, 0.05]
 
 # starttimes = {
@@ -66,9 +66,9 @@ Lnw=Lw
 lambda_w = 0
 lambda_nw = 0
 
-include_gravity = False
+include_gravity = True
 debugflag = False
-analyse_condition = True
+analyse_condition = False
 
 if mesh_study:
     output_string = "./output/{}-{}_timesteps{}_P{}".format(datestr, use_case, number_of_timesteps, FEM_Lagrange_degree)