diff --git a/Documentation/Examples/cascade_example.cc b/Documentation/Examples/cascade_example.cc
index 231f60308db573aa6e9d88dc54b2191e6ac2cb60..aa70f3d62dbb880951985bc54754c2017a8ec0e9 100644
--- a/Documentation/Examples/cascade_example.cc
+++ b/Documentation/Examples/cascade_example.cc
@@ -236,6 +236,9 @@ public:
 // The example main program for a particle cascade
 //
 int main() {
+
+  const LengthType height_atmosphere = 112.8_km;
+  
   feenableexcept(FE_INVALID);
   // initialize random number sequence(s)
   random::RNGManager::GetInstance().RegisterRandomStream("cascade");
@@ -283,7 +286,7 @@ int main() {
   process::sibyll::NuclearInteraction sibyllNuc(env, sibyll);
   process::sibyll::Decay decay(trackedHadrons);
   ProcessCut cut(20_GeV);
-  ObservationLevel obsLevel(1400_m);
+  ObservationLevel obsLevel(height_atmosphere - 2000_m);//1400_m);
 
   process::TrackWriter::TrackWriter trackWriter("tracks.dat");
   process::EnergyLoss::EnergyLoss eLoss;
@@ -318,7 +321,7 @@ int main() {
     cout << "input angles: theta=" << theta << " phi=" << phi << endl;
     cout << "input momentum: " << plab.GetComponents() / 1_GeV << endl;
     Point pos(rootCS, 0_m, 0_m,
-              112.8_km); // this is the CORSIKA 7 start of atmosphere/universe
+              height_atmosphere); // this is the CORSIKA 7 start of atmosphere/universe
     stack.AddParticle(std::tuple<particles::Code, units::si::HEPEnergyType,
                                  corsika::stack::MomentumVector, geometry::Point,
                                  units::si::TimeType, unsigned short, unsigned short>{