determine initial viscosity using all dynamic variables as indicators

src/SRHD/initialdata.jl

@@ -430,7 +430,7 @@ function initialdata_hrsc!(env, P::Project, hrsc::Union{HRSC.EntropyViscosity,HR
 
   @unpack mu, cellmax_v = get_cell_variables(env.cache)
   @unpack max_v         = get_static_variables(env.cache)
-  @unpack D             = get_dynamic_variables(env.cache)
+  @unpack D, S, tau     = get_dynamic_variables(env.cache)
   dg1d.new_broadcast_volume!(maxspeed, P.equation, env.mesh)
   Npts, K = layout(env.mesh)
   mat_max_v = reshape(view(max_v, :), (Npts, K))
@@ -438,8 +438,8 @@ function initialdata_hrsc!(env, P::Project, hrsc::Union{HRSC.EntropyViscosity,HR
     cellmax_v[k] = maximum(view(mat_max_v, :, k))
   end
   HRSC.compute_viscosity!(
-      mu,
-      D, cellmax_v,
+      mu, deepcopy(mu),
+      (D,S,tau), cellmax_v,
       mda_hrsc.av)
 end