diff --git a/code/scorpia_cpp/egd_progs/corello/demoIni/constBallDrag/README.md b/code/scorpia_cpp/egd_progs/corello/demoIni/constBallDrag/README.md
new file mode 100644
index 0000000000000000000000000000000000000000..7f50d19b48f88a77cf2bbc84c9b9792ebb206c90
--- /dev/null
+++ b/code/scorpia_cpp/egd_progs/corello/demoIni/constBallDrag/README.md
@@ -0,0 +1,20 @@
+# Info
+This is a similar scenario than the linear momentum test was.
+The main difference however is, that this time the ball is
+powered. Meaning that a special rehology is used to pull the
+ball with constant velocity. Such a setting could be used for
+measuring the drag.
+
+The viscosity of the surrounding media is not zero, but very low.
+This is important for measuring the drag properly, since we can
+not measure the viscus drag.
+
+## Integrators
+We will only use the jm integrator.
+
+
+
+
+
+
+
diff --git a/code/scorpia_cpp/egd_progs/corello/demoIni/constBallDrag/test_dragConstBall_jm.ini b/code/scorpia_cpp/egd_progs/corello/demoIni/constBallDrag/test_dragConstBall_jm.ini
new file mode 100644
index 0000000000000000000000000000000000000000..d728e13c60a122f9029bec32dbdc7b6dec1452ae
--- /dev/null
+++ b/code/scorpia_cpp/egd_progs/corello/demoIni/constBallDrag/test_dragConstBall_jm.ini
@@ -0,0 +1,197 @@
+; Test for conservation of the constant ball.
+;
+; A ball is pulled with constant velocity through a
+; low viscosity media. We can then measure
+; the drag. The szenario builds on the linearMomentum
+; test.
+;
+
+
+[Misc]
+hdf5_file=corello_constBallDrag_jm.hdf5
+
+typeOnlyFirst=yes
+NoStrainRateDump=yes
+NoStressDump=yes
+DumpOnlyRep=yes
+
+; Probability is larger than one, so dump all.
+MarkerSelectionProb=1.2
+
+; No density stabilization.
+RhoStab=false
+
+; Dump in intervals of one mili seconds
+StepBasedDumping=no
+DumpingStepInterval=5
+TimeBasedDumping=yes
+MinDumpTimeInterval=1ms
+
+
+
+[TimeStepping]
+
+FinalTime=600ms
+
+minDeltaTime=40us
+maxDeltaTime=0.5ms
+
+; Setup restriction
+nSuperSmallTimeSteps=10
+nRampSteps=40
+
+maxDisplacementX=0.1
+maxDisplacementY=0.1
+
+
+[Integrator]
+; Use the min mod lin integrator.
+type=jmMod
+
+
+[Domain]
+; This will lead of an equal grid spacing in each direction.
+;  A cell is 1cm in wide.
+Nx=301
+Ny=101
+
+nMarkerX=5
+nMarkerY=5
+MarkerRandDisplace=yes
+
+; Only use the full heatingterm
+WriteOnylFullHeatTerm=yes
+
+
+[SetUp]
+; Use the normal ball scenario
+type=constball
+
+; No gravity is acting
+gravY=0.0
+
+; The velocitoies that is read in
+;  negative x means moving to the left.
+vx=-4.0
+vy=0.0
+
+
+[StokesSolver]
+; Use the extended Navier-Stokes solver
+type=nsext
+
+
+[TemperatureSolver]
+; No not use one.
+type=nullext
+dumpSol=no
+
+
+[StokesBoundary]
+; In this section the condition for the stokes boundary conditions
+; this are velocities (x & y) and pressure are configured.
+; Note that not all options are aviable.
+
+; Here the conditions on the top can be configured.
+; This sets the x velocity condition at the top to
+; free slip.
+vel_x_top_type=freeSlip
+;vel_x_top_value=
+
+; Here we can controll the bottom condition for the x velocity..
+; We also set it to free slip.
+vel_x_bot_type=freeSlip
+;vel_x_bot_value=
+
+; Here we can controll the x velocity at the right boundary,
+; these are lower value of x. It is set to dirichlet.
+; We force it to zero.
+vel_x_rig_type=dirichlet
+vel_x_rig_value=0
+
+; Here we controll the left side of the domain, these are larger x values.
+vel_x_lef_type=dirichlet
+vel_x_lef_value=0
+
+; Here the conditions on the top can be configured.
+; This sets the y velocity condition at the top to
+; dirichlet and force it to zero.
+vel_y_top_type=dirichlet
+vel_y_top_value=0
+
+; Here we can controll the bottom condition for the y velocity..
+; We also set it to dirichlet with zero velocity.
+vel_y_bot_type=dirichlet
+vel_y_bot_value=0
+
+; Here we can controll the y velocity at the right boundary,
+; these are lower value of x. It is set to free slip.
+; We force it to zero.
+vel_y_rig_type=freeSlip
+;vel_y_rig_value=
+
+; Here we controll the y velocity at the left side of the domain,
+; these are larger x values.
+vel_y_lef_type=freeSlip
+;vel_y_lef_value=
+
+; Here we controll the pressure value. This is the linear offset.
+; This value is needed since we can only solve for the derivative
+; and we need an integration constant.
+pressure_value=1.0
+
+
+[TemperatureBoundary]
+; These section controles the temperature boundary of
+; the problem.
+
+; Here we control the temperature condition that is applied
+; at the top. We use a dirichlet condition andforce the value to 273K.
+; The value here is interpreted in Kelvin.
+temp_top_type=dirichlet
+temp_top_value=273
+
+; Here we controll the temperature that is applied to the
+; bottom of the domain. We use a dirichlet condition, the
+; value that is used is interpreted as Kelvin.
+temp_bot_type=dirichlet
+temp_bot_value=1573
+
+; Here we controll the temperature condition that is applied
+; to the left boundary, these is smaller x values.
+; It is set to insulation, this means that the derivation in
+; x direction is set to zero, and no temperature flux happens.
+temp_lef_type=insulation
+;temp_lef_value=
+
+; Here we controll the temperature condition that is applied
+; to the roight boundary. It is set to insulation meaning that
+; no temperature flux is allowed.
+temp_rig_type=insulation
+;temp_rig_value=
+
+
+[Rheology]
+; No not use a rheology, also constant marker properties.
+type=constball
+
+
+[MarkerToGrid]
+; Do not use navive interpolation, interpolate markers directly
+; to all grids that are needed.
+type=full
+
+
+[GridToMarker]
+; Only map velocity back to the marker.
+type=velo
+
+; Use change interpolation.
+JobType=c
+
+
+[ApplyBC]
+; The standard implementation uses the original conditions.
+type=pipe
+
+