don't recompute max speed for timestep

src/SRHD/rhs.jl

@@ -31,15 +31,15 @@ dtfactor(mesh::Mesh1d{SpectralElement}) = mesh.element.N
 function timestep(env, P::Project, hrsc::HRSC.AbstractArtificialViscosity, mesh::Mesh1d)
   @unpack cache, mesh = env
   @unpack equation    = P
+  @unpack max_v = get_static_variables(cache)
 
   # dg1d.new_broadcast_volume!(conservative_fixing, equation, mesh)
-  dg1d.new_broadcast_volume!(cons2prim_kastaun, equation, mesh)
-  dg1d.new_broadcast_volume!(maxspeed, equation, mesh)
-  max_v = get_variable(cache, :max_v)
+  # dg1d.new_broadcast_volume!(cons2prim_kastaun, equation, mesh)
+  # dg1d.new_broadcast_volume!(maxspeed, equation, mesh)
   vmax = dg1d.absolute_maximum(view(max_v,:))
-  vmax_limit = 0.99
+  vmax_limit = 0.9999999
   if vmax > vmax_limit
-    @warn "Limiting timestep due to maximum speed exceeding $vmax_limit"
+    @warn "Limiting timestep due to maximum speed $vmax exceeding $vmax_limit"
     vmax = vmax_limit
   end
   smoothed_mu = get_variable(cache, :smoothed_mu)