diff --git a/CHT/flow-over-plate/buoyantPimpleFoam-fenics/Solid/heat.py b/CHT/flow-over-plate/buoyantPimpleFoam-fenics/Solid/heat.py index 57c44b15e..6eed16305 100644 --- a/CHT/flow-over-plate/buoyantPimpleFoam-fenics/Solid/heat.py +++ b/CHT/flow-over-plate/buoyantPimpleFoam-fenics/Solid/heat.py @@ -178,17 +178,17 @@ def determine_heat_flux(V_g, u, k, flux): n = n_cp else: # update solution u_n.assign(u_np1) - t += dt + t += float(dt) n += 1 if precice.is_time_window_complete(): tol = 10e-5 # we need some tolerance, since otherwise output might be skipped. - if abs((float(t)+tol) % dt_out) < 2*tol: # output if t is a multiple of dt_out + if abs((t + tol) % dt_out) < 2*tol: # output if t is a multiple of dt_out print("output vtk for time = {}".format(float(t))) file_out << u_n # Update dirichlet BC - u_D.t = t + dt(0) + u_D.t = t + float(dt) # Hold plot precice.finalize() diff --git a/FSI/cylinderFlap/OpenFOAM-FEniCS/Solid/cyl-flap.py b/FSI/cylinderFlap/OpenFOAM-FEniCS/Solid/cyl-flap.py index 7675931dc..6157e9bd6 100644 --- a/FSI/cylinderFlap/OpenFOAM-FEniCS/Solid/cyl-flap.py +++ b/FSI/cylinderFlap/OpenFOAM-FEniCS/Solid/cyl-flap.py @@ -229,7 +229,7 @@ def avg(x_old, x_new, alpha): n = n_cp else: u_n.assign(u_np1) - t += dt + t += float(dt) n += 1 if precice.is_time_window_complete(): diff --git a/FSI/flap_perp/OpenFOAM-FEniCS/Solid/perp-flap.py b/FSI/flap_perp/OpenFOAM-FEniCS/Solid/perp-flap.py index 7977ea0e3..008ffb082 100644 --- a/FSI/flap_perp/OpenFOAM-FEniCS/Solid/perp-flap.py +++ b/FSI/flap_perp/OpenFOAM-FEniCS/Solid/perp-flap.py @@ -222,7 +222,7 @@ def avg(x_old, x_new, alpha): n = n_cp else: u_n.assign(u_np1) - t += dt + t += float(dt) n += 1 if precice.is_time_window_complete(): diff --git a/HT/partitioned-heat/fenics-fenics/heat.py b/HT/partitioned-heat/fenics-fenics/heat.py index 7a4c408e2..21c3742b0 100644 --- a/HT/partitioned-heat/fenics-fenics/heat.py +++ b/HT/partitioned-heat/fenics-fenics/heat.py @@ -228,7 +228,7 @@ def determine_gradient(V_g, u, flux): n = n_cp else: # update solution u_n.assign(u_np1) - t += dt + t += float(dt) n += 1 if precice.is_time_window_complete(): @@ -243,8 +243,8 @@ def determine_gradient(V_g, u, flux): error_out << error_pointwise # Update Dirichlet BC - u_D.t = t + dt(0) - f.t = t + dt(0) + u_D.t = t + float(dt) + f.t = t + float(dt) # Hold plot precice.finalize()