From e3e8004e09b6a1a2cc0941d979ea65b213581b4a Mon Sep 17 00:00:00 2001 From: Ian Hinder Date: Sat, 4 Jun 2011 12:02:57 +0200 Subject: EulerSR.m: Change initial data to 2D shock tube --- Examples/EulerSR.m | 12 +++++++----- 1 file changed, 7 insertions(+), 5 deletions(-) (limited to 'Examples') diff --git a/Examples/EulerSR.m b/Examples/EulerSR.m index e868742..e477de3 100644 --- a/Examples/EulerSR.m +++ b/Examples/EulerSR.m @@ -34,13 +34,15 @@ initialShockCalc = Name -> "eulersr_initial_shock", Schedule -> {"at CCTK_INITIAL as eulersr_initial"}, ConditionalOnKeyword -> {"initial_data", "shock"}, + Shorthands -> {X}, Equations -> { - rho -> rhoR0 StepFunction[x-0.5] + rhoL0 (1-StepFunction[x-0.5]), - v[1] -> vR0 StepFunction[x-0.5] + vL0 (1-StepFunction[x-0.5]), - v[2] -> 0, - v[3] -> 0, - epsi -> epsiR0 StepFunction[x-0.5] + epsiL0 (1-StepFunction[x-0.5]) + X -> x+y (*+y+z*), + rho -> rhoR0 StepFunction[X] + rhoL0 (1-StepFunction[X]), + v[1] -> vR0 StepFunction[X] + vL0 (1-StepFunction[X]), + v[2] -> vR0 StepFunction[X] + vL0 (1-StepFunction[X]), + v[3] -> vR0 StepFunction[X] + vL0 (1-StepFunction[X]), + epsi -> epsiR0 StepFunction[X] + epsiL0 (1-StepFunction[X]) } }; -- cgit v1.2.3