Added constants to restrict column density and emissivity ranges.

slcs-jsc · Jan 9, 2024 · 174e2ff · 174e2ff
1 parent f9eb43a
commit 174e2ff
Show file tree

Hide file tree

Showing 2 changed files with 16 additions and 1 deletion.
diff --git a/src/jurassic.c b/src/jurassic.c
@@ -5007,7 +5007,7 @@ void read_tbl(
 	    continue;
 
 	  /* Check ranges... */
-	  if (u < 0 || u > UMAX || eps < 0 || eps > 1) {
+	  if (u < UMIN || u > UMAX || eps < EPSMIN || eps > EPSMAX) {
 	    nrange++;
 	    continue;
 	  }

diff --git a/src/jurassic.h b/src/jurassic.h
@@ -168,6 +168,16 @@
 #define C2 1.43877506
 #endif
 
+/*! Minimum emissivity. */
+#ifndef EPSMIN
+#define EPSMIN 0
+#endif
+
+/*! Maximum emissivity. */
+#ifndef EPSMAX
+#define EPSMAX 1
+#endif
+
 /*! Standard gravity [m/s^2]. */
 #ifndef G0
 #define G0 9.80665
@@ -228,6 +238,11 @@
 #define TSUN 5780.
 #endif
 
+/*! Minimum column density [molecules/cm^2]. */
+#ifndef UMIN
+#define UMIN 0
+#endif
+
 /*! Maximum column density [molecules/cm^2]. */
 #ifndef UMAX
 #define UMAX 1e30