use nearest not linear interpolation to prevent nan in extrapolated p…

…roperties
desihub · Jan 12, 2022 · 36b3742 · 36b3742
1 parent 0aaba73
commit 36b3742
Show file tree

Hide file tree

Showing 2 changed files with 258 additions and 258 deletions.
diff --git a/py/desisim/templates.py b/py/desisim/templates.py
@@ -1506,10 +1506,10 @@ def make_star_templates(self, nmodel=100, vrad_meansig=(0.0, 200.0),
         else:
             from scipy.interpolate import griddata
             baseflux = griddata(base_properties, self.baseflux,
-                                input_properties, method='linear')
+                                input_properties, method='nearest')
             interptemplateid = griddata(base_properties, np.arange(nbase),
-                                        input_properties, method='linear')
-            interptemplateid = np.array([int(tempid) for tempid in interptemplateid])
+                                        input_properties, method='nearest')
+            #interptemplateid = np.array([int(tempid) for tempid in interptemplateid])
 
         # Build each spectrum in turn.
         if restframe: