diff --git a/wrinkling_testsuite/cosserat-continuum.parset b/wrinkling_testsuite/cosserat-continuum.parset
new file mode 100644
index 0000000000000000000000000000000000000000..d02e5fe48431b59c6992ee00e30ed1bd7c5bbf91
--- /dev/null
+++ b/wrinkling_testsuite/cosserat-continuum.parset
@@ -0,0 +1,77 @@
+# Number of grid levels
+numLevels = 3
+
+# Number of homotopy steps for the Dirichlet boundary conditions
+numHomotopySteps = 1
+
+# Tolerance of the trust region solver
+tolerance = 1e-6
+
+# Max number of steps of the trust region solver
+maxTrustRegionSteps = 1000
+
+# Initial trust-region radius
+initialTrustRegionRadius = 1
+
+# Number of multigrid iterations per trust-region step
+numIt = 200
+
+# Number of presmoothing steps
+nu1 = 3
+
+# Number of postsmoothing steps
+nu2 = 3
+
+# Number of coarse grid corrections
+mu = 1
+
+# Number of base solver iterations
+baseIt = 100
+
+# Tolerance of the multigrid solver
+mgTolerance = 1e-8
+
+# Tolerance of the base grid solver
+baseTolerance = 1e-8
+
+# Measure convergence
+instrumented = 0
+
+############################
+#   Problem specifications
+############################
+
+[materialParameters]
+
+# shell thickness
+thickness = 0.05
+
+# Lame parameters
+# corresponds to E = 1e6, nu=0.3
+mu = 3.8462e+05
+lambda = 2.7149e+05
+#mu = 1
+#lambda = 1
+
+# Cosserat couple modulus
+mu_c = 0
+#mu_c = 3.8462e+05
+
+# Length scale parameter
+L_c = 0.02
+
+# Curvature exponent
+q = 2
+
+# Shear correction factor
+kappa = 1
+
+[]
+
+# 2d problem
+#path = /home/haile/sander/data/richards/twosquares/
+#gridFile = twosquares0.grid
+
+# 3d problem
+path = /home/haile/sander/data/contact/tetracubes/
+gridFile = tetracube0.grid