Add input.prm and reference it.

git-svn-id: https://svn.dealii.org/trunk@16063 0785d39b-7218-0410-832d-ea1e28bc413d

diff --git a/deal.II/examples/step-33/doc/results.dox b/deal.II/examples/step-33/doc/results.dox
index 586a91170e09913fa5280e2ba80c2870a7e337b0..6bdcb6c4c46c002075a7d89fe844affbcd3c4d41 100644 (file)
--- a/deal.II/examples/step-33/doc/results.dox
+++ b/deal.II/examples/step-33/doc/results.dox
@@ -3,7 +3,7 @@
 
 We run the problem with the mesh <code>slide.inp</code> (this file is in the
 same directory as the source code for this program) and the following input
-deck: 
+deck (available as <code>input.prm</code> in the same directory): 
 @verbatim
 # Listing of Parameters
 # ---------------------

diff --git a/deal.II/examples/step-33/input.prm b/deal.II/examples/step-33/input.prm
new file mode 100644 (file)
index 0000000..f2f5e6e
--- /dev/null
+++ b/deal.II/examples/step-33/input.prm
@@ -0,0 +1,86 @@
+# Listing of Parameters
+# ---------------------
+
+# The input grid 
+set mesh = slide.inp
+
+# Stabilization parameter
+set diffusion power = 2.0
+
+# Scaled value for gravity.  Positive means gravity points down.
+set gravity = 1.0
+
+# --------------------------------------------------
+# Boundary conditions
+# We may specify boundary conditions for up to MAX_BD boundaries.
+# Your .inp file should have these boundaries designated.
+subsection boundary_1
+  set no penetration = true # reflective boundary condition
+end
+
+subsection boundary_2
+  # outflow boundary
+  # set w_2 = pressure
+  # set w_2 value = 1.5 - y
+end
+
+subsection boundary_3
+  set no penetration = true # reflective
+  # set w_3 = pressure
+  # set w_3 value = 1.0
+end
+
+subsection boundary_4
+  set no penetration = true #reflective
+end
+
+# --------------------------------------------------
+# Initial Conditions
+# We set the initial conditions of the conservative variables.  These lines
+# are passed to the expression parsing function.  You should use x,y,z for
+# the coordinate variables.
+
+subsection initial condition
+  set w_0 value = 0
+  set w_1 value = 0
+  set w_2 value = 10*(x<-0.7)*(y> 0.3)*(y< 0.45) + (1-(x<-0.7)*(y> 0.3)*(y< 0.45))*1.0
+  set w_3 value = (1.5-(1.0*1.0*y))/0.4
+end
+
+# --------------------------------------------------
+# Time stepping control
+subsection time stepping
+  set final time = 10.0 # simulation end time
+  set time step  = 0.02 # simulation time step
+end
+
+subsection linear solver
+  set output         = quiet
+  set method         = gmres
+  set ilut fill      = 1.5
+  set ilut drop tolerance = 1e-6
+  set ilut absolute tolerance = 1e-6
+  set ilut relative tolerance = 1.0
+end
+
+# --------------------------------------------------
+# Output frequency.
+# You may wish to set this > time step if you dont want output at every step
+subsection output
+  set step = 0.01
+end
+
+# --------------------------------------------------
+# Refinement control
+subsection refinement
+  set refinement = shock # none only other option
+  set shock value = 1.5
+  set shock levels = 1 # how many levels of refinement to allow
+end
+
+# --------------------------------------------------
+# Flux parameters
+subsection flux
+ set stab = mesh
+ #set stab value = 1.0
+end