From 63955c258f3e4114ef2c894aa2a2423508611b93 Mon Sep 17 00:00:00 2001 From: Ansgar Wehrhahn <31626864+AWehrhahn@users.noreply.github.com> Date: Fri, 3 Jun 2022 16:35:52 +0200 Subject: [PATCH] fix the lower abundance limit --- src/pysme/sme_synth.py | 4 +++- 1 file changed, 3 insertions(+), 1 deletion(-) diff --git a/src/pysme/sme_synth.py b/src/pysme/sme_synth.py index cb9edb29..617fef2f 100644 --- a/src/pysme/sme_synth.py +++ b/src/pysme/sme_synth.py @@ -245,7 +245,9 @@ def InputAbund(self, abund): """ abund = abund("sme", raw=True) # TODO: does this make sense or should we fix the C code? - abund[np.isnan(abund)] = -12 + abund[np.isnan(abund)] = -13 + abund[abund < -13] = -13 + _smelib.InputAbund(abund) def Opacity(self):