From 260d6c5e0b8b54ff64ba4222a91524a1ffe41633 Mon Sep 17 00:00:00 2001
From: Stephen Fegan <sfegan@gmail.com>
Date: Mon, 1 Jul 2024 11:47:59 +0200
Subject: [PATCH] Small fix

---
 include/simulation/vcl_iact_ground_map.hpp | 11 ++---------
 1 file changed, 2 insertions(+), 9 deletions(-)

diff --git a/include/simulation/vcl_iact_ground_map.hpp b/include/simulation/vcl_iact_ground_map.hpp
index 39ee2ed2..7c4c7e48 100644
--- a/include/simulation/vcl_iact_ground_map.hpp
+++ b/include/simulation/vcl_iact_ground_map.hpp
@@ -84,8 +84,6 @@ template<typename VCLArchitecture> class alignas(VCLArchitecture::vec_bytes) VCL
 
   unsigned iobs_ = 0;
   double zobs_ = 0.0;
-  double ref_index_ = 1.0;
-  double ref_index_correction_ = 0.0;
 #endif
 };
 
@@ -101,9 +99,6 @@ VCLIACTTrackVisitor<VCLArchitecture>(atm, config, rng, adopt_atm, adopt_rng)
 {
   iobs_ = 0;
   zobs_ = this->atm_->zobs(iobs_);
-  double n_minus_one = this->atm_->n_minus_one(zobs_);
-  ref_index_ = 1.0 + n_minus_one;
-  ref_index_correction_ = -n_minus_one;
 }
 
 template<typename VCLArchitecture> VCLIACTGroundMap<VCLArchitecture>::
@@ -140,13 +135,11 @@ propagate_rays(calin::math::ray::VCLRay<double_real> ray, double_bvt ray_mask,
 
   ray_mask &= (ray.z()>zobs_) & (ray.uz()<0);
 
-  double_vt dz = ray.z()-zobs_;
   ray.mutable_ct() -= vcl::select(ray_mask,
-    (dz*ref_index_correction_ + this->atm_->template vcl_propagation_ct_correction_to_iobs<VCLArchitecture>(ray.z(), 0))/ray.uz(),
+    (this->atm_->template vcl_propagation_ct_correction_to_iobs<VCLArchitecture>(ray.z(), 0))/ray.uz(),
     0);
 
-  ray_mask = ray.propagate_to_z_plane_with_mask(ray_mask,
-    VCLIACTTrackVisitor<VCLArchitecture>::atm_->zobs(0), false);
+  ray_mask = ray.propagate_to_z_plane_with_mask(ray_mask, zobs_, false);
 
   double_at t;
   double_at x;