diff --git a/examples/vertical_EAS.cpp b/examples/vertical_EAS.cpp
index ba69ffe2814e599a3f7bf6eb582a864beffc62c5..6e12bb9c404ea434ffb071d30cc134faa8b70b79 100644
--- a/examples/vertical_EAS.cpp
+++ b/examples/vertical_EAS.cpp
@@ -42,6 +42,7 @@
 #include <corsika/media/UniformMagneticField.hpp>
 #include <corsika/media/ShowerAxis.hpp>
 #include <corsika/media/CORSIKA7Atmospheres.hpp>
+#include <corsika/media/WMM.hpp>
 
 #include <corsika/modules/BetheBlochPDG.hpp>
 #include <corsika/modules/LongitudinalProfile.hpp>
@@ -129,7 +130,7 @@ int main(int argc, char** argv) {
   // build a Linsley US Standard atmosphere into `env`
   create_5layer_atmosphere<setup::EnvironmentInterface, MyExtraEnv>(
       env, AtmosphereId::LinsleyUSStd, center, Medium::AirDry1Atm,
-      MagneticFieldVector{rootCS, 0_T, 50_uT, 0_T});
+      get_wmm(rootCS, 2022.5, 10_km, 49, 8.4));
 
   // pre-setup particle stack
   unsigned short const A = std::stoi(std::string(argv[1]));