diff --git a/src/test/java/org/orekit/rugged/atmosphericrefraction/MultiLayerModelTest.java b/src/test/java/org/orekit/rugged/atmosphericrefraction/MultiLayerModelTest.java
index 1ac7b056249fa7d5ea069eb59b5d589e12161802..5b3d7cd4fd7bf5c2c3be22d5eea3ff19b5589867 100644
--- a/src/test/java/org/orekit/rugged/atmosphericrefraction/MultiLayerModelTest.java
+++ b/src/test/java/org/orekit/rugged/atmosphericrefraction/MultiLayerModelTest.java
@@ -17,6 +17,7 @@
 package org.orekit.rugged.atmosphericrefraction;
 
 import org.apache.commons.math3.geometry.euclidean.threed.Vector3D;
+import org.apache.commons.math3.util.FastMath;
 import org.junit.Test;
 import org.orekit.bodies.GeodeticPoint;
 import org.orekit.bodies.OneAxisEllipsoid;
@@ -25,7 +26,9 @@ import org.orekit.frames.FramesFactory;
 import org.orekit.rugged.errors.RuggedException;
 import org.orekit.rugged.intersection.AbstractAlgorithmTest;
 import org.orekit.rugged.intersection.IntersectionAlgorithm;
+import org.orekit.rugged.intersection.duvenhage.DuvenhageAlgorithm;
 import org.orekit.rugged.intersection.duvenhage.MinMaxTreeTile;
+import org.orekit.rugged.intersection.duvenhage.MinMaxTreeTileFactory;
 import org.orekit.rugged.raster.TileUpdater;
 import org.orekit.rugged.utils.ExtendedEllipsoid;
 import org.orekit.rugged.utils.NormalizedGeodeticPoint;
@@ -54,8 +57,8 @@ public class MultiLayerModelTest extends AbstractAlgorithmTest {
         // check mayon volcano tiles
 
 
-
         setUpMayonVolcanoContext();
+
         final IntersectionAlgorithm algorithm = createAlgorithm(updater, 8);
         Vector3D position = new Vector3D(-3787079.6453602533, 5856784.405679551, 1655869.0582939098);
         Vector3D los = new Vector3D( 0.5127552821932051, -0.8254313129088879, -0.2361041470463311);
@@ -64,10 +67,21 @@ public class MultiLayerModelTest extends AbstractAlgorithmTest {
 
 
 
+
         MultiLayerModel model = new MultiLayerModel();
 
 
 
+        final double latitude  = FastMath.toRadians(13.27);
+        final double longitude = FastMath.toRadians(123.68);
+        MinMaxTreeTile tile = new MinMaxTreeTileFactory().createTile();
+        updater.updateTile(latitude, longitude, tile);
+        double altitude = tile.interpolateElevation(latitude, longitude);
+        final GeodeticPoint groundGP = new GeodeticPoint(latitude, longitude, altitude);
+        Vector3D groundP = earth.transform(groundGP);
+
+
+
         // Vector3D initialPos, Vector3D initialLos, Vector3D initialZenith, double altitude, Tile tile
         // model.getPointOnGround();
 
@@ -77,6 +91,6 @@ public class MultiLayerModelTest extends AbstractAlgorithmTest {
 
     @Override
     protected IntersectionAlgorithm createAlgorithm(TileUpdater updater, int maxCachedTiles) {
-        return null;
+        return new DuvenhageAlgorithm(updater, maxCachedTiles, false);
     }
 }