raise exception if wavelength distribution of multiple dust data file…

…s do not match

miex_app.py

@@ -5,9 +5,9 @@
 import streamlit as st
 
 
-# Converter for lambda/n/k database files
 def conv(line):
-    return line.replace("D", "e")
+    # Converter for lambda/n/k database files
+    return line.replace(b"D", b"e")
 
 
 st.set_page_config(page_title="MIEX", page_icon=None)
@@ -79,7 +79,7 @@ def conv(line):
     abun = np.array([1.0])
 else:
     st.info(
-        "all data files have to contain the refractive index and the same wavelength distribution"
+        "All data files have to contain the refractive index and the same wavelength distribution. See https://github.com/mlietzow/MIEX-Python/tree/main/ri-data for exemplary files."
     )
     with col1:
         ncomp = st.number_input(
@@ -252,6 +252,10 @@ def conv(line):
                     st.stop()
             else:
                 if nlam <= len(w):
+                    if icomp > 0 and not np.array_equal(wavelength, w[:nlam]):
+                        st.error("Wavelength distribution in dust data files do not match.")
+                        st.stop()
+
                     wavelength = w[:nlam]
                     ri_real[icomp] = n[:nlam]
                     ri_imag[icomp] = k[:nlam]

run_miex.py

@@ -47,7 +47,7 @@
 
 def conv(line):
     # Converter for lambda/n/k database files
-    return line.replace("D", "e")
+    return line.replace(b"D", b"e")
 
 
 def main(input_filename):
@@ -166,12 +166,17 @@ def main(input_filename):
     # read lambda/n/k database
     for icomp in range(ncomp):
         w, n, k = np.loadtxt(
-            "ri-data/" + fnames[icomp], comments="#", unpack=True, converters=conv
+            "ri-data/" + fnames[icomp],
+            comments="#",
+            unpack=True,
+            converters=conv,
         )
+        if icomp > 0 and not np.array_equal(wavelength, w[:nlam]):
+            raise Exception("Wavelength distribution in dust data files do not match.")
+
         wavelength = w[:nlam]
         n_real[icomp] = n[:nlam]
         k_imag[icomp] = k[:nlam]
-        # first two lines of file give information about the content of the file; no need to save in a variable
 
     # define radial step width
     radminlog = np.log10(radmin)