update TP-one-patch-purely-postive

TP-one-patch/TP-one-patch-purely-postive-pc.py

@@ -14,16 +14,16 @@ import helpers as hlp
 sym.init_printing()
 
 
-solver_tol = 5E-6
+solver_tol = 1E-6
 
 ############ GRID #######################
-use_case="TP-one-patch-purely-positive-pc"
-mesh_resolution = 30
-timestep_size = 0.0005
-number_of_timesteps = 2500
+
+mesh_resolution = 40
+timestep_size = 0.001
+number_of_timesteps = 10
 # 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
+number_of_timesteps_to_analyse = 1
 starttime = 0
 
 Lw = 0.25  #/timestep_size
@@ -32,18 +32,19 @@ Lnw=Lw
 l_param_w = 40
 l_param_nw = 40
 
-include_gravity = True
-debugflag = False
+include_gravity = False
+debugflag = True
 analyse_condition = True
 
-output_string = "./output/nondirichlet_pure_positive_pc_number_of_timesteps{}_".format(number_of_timesteps)
+use_case="TP-one-patch-purely-positive-pc"
+output_string = "./output/pure_positive_pc_number_of_timesteps{}_".format(number_of_timesteps)
 
 ##### Domain and Interface ####
 # global simulation domain domain
-sub_domain0_vertices = [df.Point(-1.0,-1.0), #
-                        df.Point(1.0,-1.0),#
-                        df.Point(1.0,1.0),#
-                        df.Point(-1.0,1.0)]
+sub_domain0_vertices = [df.Point(-1.0, -1.0),  #
+                        df.Point(1.0, -1.0),  #
+                        df.Point(1.0, 1.0),  #
+                        df.Point(-1.0, 1.0)]
 
 subdomain0_outer_boundary_verts = {
     0: [sub_domain0_vertices[0],