From c9eff6ec07f3c06a70e1e02f66440b5f59030cd7 Mon Sep 17 00:00:00 2001
From: =?UTF-8?q?Andr=C3=A9s=20Rueda-Ram=C3=ADrez?= <aruedara@uni-koeln.de>
Date: Tue, 6 Jun 2023 12:13:37 +0200
Subject: [PATCH] minor modifications to OT

---
 .../elixir_mhdmultiion_orszag_tang.jl         | 29 ++++++++++---------
 1 file changed, 15 insertions(+), 14 deletions(-)

diff --git a/examples/tree_2d_dgsem/elixir_mhdmultiion_orszag_tang.jl b/examples/tree_2d_dgsem/elixir_mhdmultiion_orszag_tang.jl
index 2dcb3a4ae3..a226d15822 100644
--- a/examples/tree_2d_dgsem/elixir_mhdmultiion_orszag_tang.jl
+++ b/examples/tree_2d_dgsem/elixir_mhdmultiion_orszag_tang.jl
@@ -32,24 +32,24 @@ function initial_condition_orszag_tang(x, t, equations::IdealMhdMultiIonEquation
   B3 = 0.0
   
   rho1 = zero(x[1])
-  if x[2] > 0.75
-    rho1 = 0.49 * (tanh(50 * (x[2] - 1.0)) + 1) + 0.02
-  elseif x[2] > 0.25
-    rho1 = 0.49 * (-tanh(50 * (x[2] - 0.5)) + 1) + 0.02
+  if x[1] > 0.75
+    rho1 = 0.49 * (tanh(50 * (x[1] - 1.0)) + 1) + 0.02
+  elseif x[1] > 0.25
+    rho1 = 0.49 * (-tanh(50 * (x[1] - 0.5)) + 1) + 0.02
   else
-    rho1 = 0.49 * (tanh(50 * (x[2])) + 1) + 0.02
+    rho1 = 0.49 * (tanh(50 * (x[1])) + 1) + 0.02
   end
 
-  if x[2] < 0.25
-    rho2 = 0.49 * (-tanh(50 * (x[2])) + 1) + 0.02
-  elseif x[2] < 0.75
-    rho2 = 0.49 * (tanh(50 * (x[2] - 0.5)) + 1) + 0.02
+  if x[1] < 0.25
+    rho2 = 0.49 * (-tanh(50 * (x[1])) + 1) + 0.02
+  elseif x[1] < 0.75
+    rho2 = 0.49 * (tanh(50 * (x[1] - 0.5)) + 1) + 0.02
   else
-    rho2 = 0.49 * (-tanh(50 * (x[2] - 1.0)) + 1) + 0.02
+    rho2 = 0.49 * (-tanh(50 * (x[1] - 1.0)) + 1) + 0.02
   end
 
-  p1 = T # * rho1
-  p2 = T # * rho2
+  p1 = T * rho1
+  p2 = T * rho2
 
   return prim2cons(SVector(B1, B2, B3, rho1, v1, v2, v3, p1, rho2, v1, v2, v3, p2), equations)
 
@@ -81,12 +81,13 @@ mesh = TreeMesh(coordinates_min, coordinates_max,
 
 semi = SemidiscretizationHyperbolic(mesh, equations, initial_condition, solver,
                                     source_terms=source_terms_collision_ion_ion)
+                                    #source_terms=source_terms_standard)
 
 
 ###############################################################################
 # ODE solvers, callbacks etc.
 
-tspan = (0.0, 1.0)
+tspan = (0.0, 0.5)
 ode = semidiscretize(semi, tspan)
 
 summary_callback = SummaryCallback()
@@ -95,7 +96,7 @@ analysis_interval = 100
 analysis_callback = AnalysisCallback(semi, interval=analysis_interval)
 alive_callback = AliveCallback(analysis_interval=analysis_interval)
 
-save_solution = SaveSolutionCallback(interval=10,
+save_solution = SaveSolutionCallback(interval=100,
                                      save_initial_solution=true,
                                      save_final_solution=true,
                                      solution_variables=cons2prim)