use maximum of max_v instead of v for time step estimation

src/SRHD/rhs.jl

@@ -6,9 +6,9 @@ timestep(env, P, hrsc) = timestep(env, P, hrsc, env.mesh)
 
 
 function compute_maxspeed(mesh, equation, cache)
-  @unpack v = get_static_variables(cache)
+  @unpack max_v = get_static_variables(cache)
   dg1d.new_broadcast_volume!(maxspeed, equation, mesh)
-  vmax = dg1d.absolute_maximum(v)
+  vmax = dg1d.absolute_maximum(max_v)
   vmax_limit = 1e4
   if vmax > vmax_limit
     @warn "Limiting timestep due to maximum speed exceeding $vmax_limit"