From 1bc7fd14efcb5167066c2fa0d00a17b5eccb65e0 Mon Sep 17 00:00:00 2001
From: hrocho <hrocho@vodacionline.sk>
Date: Thu, 8 Feb 2018 00:59:31 +0100
Subject: [PATCH] deff zcheck

---
 calculations.cpp | 6 ++++--
 1 file changed, 4 insertions(+), 2 deletions(-)

diff --git a/calculations.cpp b/calculations.cpp
index 5cca39c..c60ed2f 100644
--- a/calculations.cpp
+++ b/calculations.cpp
@@ -134,12 +134,14 @@ double bethek_barkas(double zp_eff,double eta, double zt){
 double bethek_density_effect(double beta, int zt){
     double gamma = 1/sqrt(1-(beta*beta));
     double x = log(beta * gamma) / 2.3025851;
-    int i = zt-1;
+    int i;
     double del = 0;
     
     if(zt>97){  // check if data are available, if not take highest z data
-        i=96;
+        zt=97;
         }
+    i = zt-1;
+    
     if (x < density_effect::x0[i] ){
         if(density_effect::del_0[i] > 0.)del = density_effect::del_0[i] * pow(10.0,(2.*(x-density_effect::x0[i])));
         }