catch failing get_horizons_coord()

jgieseler · Feb 29, 2024 · f9947af · f9947af
1 parent 982ad2e
commit f9947af
Showing 1 changed file with 5 additions and 1 deletion.
diff --git a/solarmach/__init__.py b/solarmach/__init__.py
@@ -218,7 +218,11 @@ def __init__(self, date, body_list, vsw_list=[], reference_long=None, reference_
         self.reference_lat = reference_lat
         self.coord_sys = coord_sys
 
-        pos_E = get_horizons_coord(399, self.date, None)  # (lon, lat, radius) in (deg, deg, AU)
+        try:
+            pos_E = get_horizons_coord(399, self.date, None)  # (lon, lat, radius) in (deg, deg, AU)
+        except ValueError:
+            print('')
+            print('!!! No ephemeris found for Earth for date {self.date} - there probably is a problem with JPL HORIZONS.')
         if coord_sys=='Carrington':
             self.pos_E = pos_E.transform_to(frames.HeliographicCarrington(observer='Sun'))
         elif coord_sys=='Stonyhurst':