Test the radiative transfer

AWehrhahn · Mar 9, 2021 · 6c9aff4 · 6c9aff4
1 parent 0c79c40
commit 6c9aff4
Show file tree

Hide file tree

Showing 3 changed files with 600 additions and 15 deletions.
diff --git a/.gitignore b/.gitignore
@@ -85,3 +85,4 @@ src/eos/eos_math_special.c
 src/eos/eos.c
 src/sme/hlinop.c
 src/sme/hlinprof.c
+debug_radiative_transfer.npz
diff --git a/test/sme_synth.py b/test/sme_synth.py
@@ -176,12 +176,12 @@ def InputLineList(self, linelist):
             names of the elements (with Ionization level)
         """
         try:
-            atomic = linelist.atomic.T
-            species = linelist.species
+            atomic = linelist["atomic"].T
+            species = linelist["species"]
         except AttributeError:
             raise TypeError("linelist has to be a LineList type")
 
-        nlines = len(linelist)
+        nlines = 1
         species = np.asarray(species, "U8")
 
         assert (
@@ -265,16 +265,16 @@ def InputModel(self, teff, grav, vturb, atmo):
             raise ValueError("Turbulence velocity must be positive or zero")
 
         try:
-            motype = atmo.depth
+            motype = atmo["depth"]
             depth = atmo[motype]
             ndepth = len(depth)
-            t = atmo.temp
-            xne = atmo.xne
-            xna = atmo.xna
-            rho = atmo.rho
+            t = atmo["temp"]
+            xne = atmo["xne"]
+            xna = atmo["xna"]
+            rho = atmo["rho"]
             vt = np.full(ndepth, vturb) if np.size(vturb) == 1 else vturb
-            wlstd = atmo.wlstd
-            opflag = atmo.opflag
+            wlstd = atmo["wlstd"]
+            opflag = atmo["opflag"]
             args = [
                 ndepth,
                 teff,
@@ -291,9 +291,9 @@ def InputModel(self, teff, grav, vturb, atmo):
             ]
             type = "sdddusdddddd"  # s : short, d: double, u: unicode (string)
 
-            if atmo.geom == "SPH":
-                radius = atmo.radius
-                height = atmo.height
+            if atmo["geom"] == "SPH":
+                radius = atmo["radius"]
+                height = atmo["height"]
                 motype = "SPH"
                 args = args[:5] + [radius] + args[5:] + [height]
                 type = type[:5] + "d" + type[5:] + "d"