RB and TV ICs with small changes

e.g. TV T=3.0 and not 10.0

Author: ray-chew

RKLM_Python/inputs/rising_bubble.py

@@ -8,6 +8,7 @@
 from management.variable import States
 
 from physics.low_mach.second_projection import euler_backward_non_advective_impl_part
+from management import data
 
 class UserData(object):
     NSPEC = 1
@@ -85,6 +86,7 @@ def __init__(self):
         self.zmax =   1.0
 
         self.u_wind_speed = 0.0
+        self.v_wind_speed = 0.0
         self.w_wind_speed = 0.0
 
         self.bdry_type = np.empty((3), dtype=object)
@@ -101,8 +103,8 @@ def __init__(self):
         # self.dtfixed = 1.9 / self.t_ref
 
         # self.CFL = 1.0
-        # self.dtfixed0 = 20.0 / self.t_ref
-        # self.dtfixed = 20.0 / self.t_ref
+        # self.dtfixed0 = 21.69 / self.t_ref
+        # self.dtfixed = 21.69 / self.t_ref
         # self.dtfixed0 = 0.05
         # self.dtfixed = 0.05
         self.dtfixed0 = 100.0
@@ -146,7 +148,7 @@ def __init__(self):
         if self.continuous_blending == True:
             self.output_suffix = "_%i_%i_%.1f" %(self.inx-1,self.iny-1,self.tout[-1])
 
-        aux = 'debug_CFLfixed'
+        aux = 'debug_imbal_CFLfixed'
         self.aux = aux
         self.output_suffix = "_%i_%i_%.1f_%s" %(self.inx-1,self.iny-1,self.tout[-1],aux)
 

RKLM_Python/inputs/travelling_vortex_2D.py

@@ -123,7 +123,8 @@ def __init__(self):
         self.initial_projection = True
         self.initial_impl_Euler = False
 
-        self.tout = np.arange(0.0,10.05,0.05)[1:]
+        # self.tout = np.arange(0.0,10.05,0.05)[1:]
+        self.tout = np.arange(0.0,3.05,0.05)[1:]
         self.stepmax = 20000
 
         self.output_base_name = "_travelling_vortex"