diff --git a/tools/jobsystem/run_dg_benchmarks.py b/tools/jobsystem/run_dg_benchmarks.py
index 9f5a149395f330b918b903bcd3c52b03aeb8702e..197f8600b95fbb8c2dfd180721cca20585e52d57 100644
--- a/tools/jobsystem/run_dg_benchmarks.py
+++ b/tools/jobsystem/run_dg_benchmarks.py
@@ -84,7 +84,7 @@ def modifyJobskript(name, degree, logfolder,penalty, args):
     if mat=='Holzapfel' and degree > 1:
         args.update({'ReferenceLevel': '3'})
     else:
-        args.update({'ReferenceLevel': '4'})
+        args.update({'ReferenceLevel': '3'})
         args.update({'DGPenalty': penalty})
     filename = name + '_' + mat.lower() + '_p' + str(degree) + 'DG' + str(penalty)
     jobname = args['MechProblem'] + mat + 'P' + str(degree) + 'DG' + str(penalty)
@@ -131,6 +131,18 @@ def elliposid_jobs(mesh, run_jobs=False):
             'NewtonDamping': '0.75'
         }
     ]
+    arglists3 = [
+        {
+            'ActiveMaterial': 'Holzapfel',
+            'QuasiCompressiblePenalty': 'Ciarlet',
+            'DGSign': '-1',
+            'DGPenalty': '700',
+            'Overlap': 'dG1',
+            'Overlap_Distribution': '1',
+            'InterpolateStartVector': 'false',
+            'NewtonDamping': '0.75'
+        }
+    ]
     penalty1 = {'20'}
     for d in range(1,3):
         for arglist1 in arglists1:
@@ -150,6 +162,15 @@ def elliposid_jobs(mesh, run_jobs=False):
                 if run_jobs:
                     print('running job', jobname)
                     os.system('sbatch ' + jobname)
+    penalty3 = {'700'}
+    for d in range(1,3):
+        for arglist3 in arglists3:
+            arglist3.update({"MechProblem": mesh})
+            for p3 in penalty3:
+                jobname = modifyJobskript('run_'+mesh.lower() + 'DG', d, '../data/benchmarks/', p3, arglist3)
+                if run_jobs:
+                    print('running job', jobname)
+                    os.system('sbatch ' + jobname)