minor

src/main/java/org/orekit/rugged/atmosphericrefraction/AtmosphericRefraction.java

@@ -18,6 +18,7 @@ package org.orekit.rugged.atmosphericrefraction;
 
 
 import org.apache.commons.math3.geometry.euclidean.threed.Vector3D;
+import org.orekit.rugged.raster.Tile;
 
 /**
  * Interface for atmospheric refraction.
@@ -25,6 +26,6 @@ import org.apache.commons.math3.geometry.euclidean.threed.Vector3D;
  */
 public interface AtmosphericRefraction {
 
-    double getDeviation(Vector3D pos, Vector3D los, Vector3D zenith, double altitude);
+    double getDeviation(Vector3D pos, Vector3D los, Vector3D zenith, double altitude, Tile tile);
 
 }

src/main/java/org/orekit/rugged/atmosphericrefraction/MultiLayerModel.java

@@ -18,6 +18,9 @@ package org.orekit.rugged.atmosphericrefraction;
 
 import org.apache.commons.math3.geometry.euclidean.threed.Vector3D;
 import org.apache.commons.math3.util.FastMath;
+import org.orekit.rugged.raster.Tile;
+import org.orekit.rugged.utils.ExtendedEllipsoid;
+import org.orekit.rugged.utils.NormalizedGeodeticPoint;
 
 import java.util.Collections;
 import java.util.Map;
@@ -53,7 +56,10 @@ public class MultiLayerModel implements AtmosphericRefraction {
     }
 
     @Override
-    public double getDeviation(Vector3D pos, Vector3D los, Vector3D zenith, double altitude) {
+    public double getDeviation(Vector3D pos, Vector3D los, Vector3D zenith, double altitude, Tile tile) {
+
+        new ExtendedEllipsoid(ellipsoid.getEquatorialRadius(), ellipsoid.getFlattening(),
+                ellipsoid.getBodyFrame());
 
         double incidenceAngleSin = FastMath.sin(Vector3D.angle(los, zenith));
         double previousRefractionIndex = -1;
@@ -73,6 +79,8 @@ public class MultiLayerModel implements AtmosphericRefraction {
             previousRefractionIndex = entry.getValue();
         }
 
+        NormalizedGeodeticPoint geodeticPoint = tile.cellIntersection(pos, los, 0, 0);
+
         return pos.getX() + xDistance;
     }
 }