diff --git a/demo/gptools_mds_ne_JRW.py b/demo/gptools_mds_ne_JRW.py
index <HASH>..<HASH> 100644
--- a/demo/gptools_mds_ne_JRW.py
+++ b/demo/gptools_mds_ne_JRW.py
@@ -181,7 +181,7 @@ k = gptools.GibbsKernel1dTanh(
 #                                        enforce_bounds=True)
 
 # Set noise kernel:
-nk = gptools.DiagonalNoiseKernel(1, n=0, initial_noise=0.0, fixed_noise=True, noise_bound=(0.0, 10.0))
+nk = gptools.DiagonalNoiseKernel(1, n=0, initial_noise=0.001, fixed_noise=False, noise_bound=(0.0, 10.0))
 
 # Create and populate GP:
 gp = gptools.GaussianProcess(k, noise_k=nk)
Try adding some noise for JRW demo.