added ou correlation time k scaling

cpp/full_code/ornstein_uhlenbeck_process.hpp

@@ -52,8 +52,10 @@ class ornstein_uhlenbeck_process{
 
         inline double gamma(double kabs)
         {
-            return this->kolmogorov_constant*sqrt(this->kolmogorov_constant/
-                            this->ou_energy_amplitude);
+            //return this->kolmogorov_constant*sqrt(this->kolmogorov_constant/
+             //               this->ou_energy_amplitude);
+            return 10.*pow(kabs,2./3.)*this->kolmogorov_constant*sqrt(this->kolmogorov_constant/
+                           this->ou_energy_amplitude); //TODO take out 10
         }