stardis-test

commit e81ed18df517cf8b6a43399fb2f14dc6c9b2bd6e
parent 720d31e1242cceb174f04a212cff4935b642e942
Author: Vincent Forest <vincent.forest@meso-star.com>
Date:   Sun, 14 Apr 2024 14:19:06 +0200

Update the programmable 1D radiative environment

Add the stardis_t_range function, which returns the temperature range
that the radiative environment can have. In fact, this should be the
range of reference temperatures used to linearize radiative transfer,
not the temperatures themselves.

Diffstat:
M
src/sadist_conducto_radiative.c
 | 
2
+-
M
src/sadist_lib_radenv_1d.c
 | 
13
+++++++++++++

2 files changed, 14 insertions(+), 1 deletion(-)
diff --git a/src/sadist_conducto_radiative.c b/src/sadist_conducto_radiative.c
@@ -50,7 +50,7 @@
 #define FILENAME_SCENE "scene.txt"
 
 #define EMISSIVITY 1.0
-#define TREF 300.0 /* [K] */
+#define TREF 315.0 /* [K] */
 #define T0 280.0 /* [K] */
 #define T1 350.0 /* [K] */
 #define LAMBDA 0.1
diff --git a/src/sadist_lib_radenv_1d.c b/src/sadist_lib_radenv_1d.c
@@ -181,6 +181,19 @@ stardis_radiative_env_reference_temperature(const double dir[3], void* data)
   return radenv->reference_temperatures[dir[0] > 0]; /* [K] */
 }
 
+double*
+stardis_t_range(void* data, double t_range[2])
+{
+  struct radenv* radenv = data;
+  const double* trefs = NULL;
+  ASSERT(radenv);
+
+  trefs = radenv->reference_temperatures;
+  t_range[0] = MMIN(trefs[0], trefs[1]);
+  t_range[1] = MMAX(trefs[0], trefs[1]);
+  return t_range;
+}
+
 /*******************************************************************************
  * Legal notices
  ******************************************************************************/