introduce a minimum value for the scattering RMS theta0

1 files changed, 1 insertions, 1 deletions

diff --git a/src/path.c b/src/path.c
index 80bba50..5a87cf6 100644
--- a/src/path.c
+++ b/src/path.c
@@ -175,7 +175,7 @@ void path_eval(path *p, double s, double *pos, double *dir, double *T, double *t
     dir[2] = interp1d(s,p->s,p->dz,N);
 
     /* FIXME: This should be the *residual* scattering RMS. */
-    *theta0 = p->theta0*sqrt(s);
+    *theta0 = fmax(MIN_THETA0, p->theta0*sqrt(s));
 }
 
 void path_free(path *p)