From ae90fe0e0c47e39a803c453e565e0eebea2ecd65 Mon Sep 17 00:00:00 2001
From: tlatorre <tlatorre@uchicago.edu>
Date: Thu, 6 Sep 2018 10:13:52 -0500
Subject: introduce a minimum value for the scattering RMS theta0

---
 src/path.c | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

(limited to 'src/path.c')

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)
-- 
cgit