Merge branch 'release/0.16.2'

climetlab/plotting/backends/magics/backend.py

@@ -109,7 +109,17 @@ def _push_layer(self, data):
         self._layers.append(Layer(data))
 
     def plot_grib(self, path: str, offset: int):
-        assert offset < 2**31, "grib_field_position must be less than 2**31"
+        if offset >= 2**31:
+            # Because Magics (version <4.14) does not support input
+            # grib files whihc are > 2**31 Bytes.
+            from climetlab.readers.grib.codes import CodesReader
+
+            tmp = self.temporary_file(".grib")
+            reader = CodesReader.from_cache(path)
+            handle = reader.at_offset(int(offset))
+            handle.save(tmp)
+            path = tmp
+            offset = 0
 
         self._push_layer(
             mgrib(

climetlab/sources/constants.py

@@ -136,9 +136,33 @@ def sin_local_time(self, date):
     def insolation(self, date):
         from ecmwf_mlkit.variables.insolation import insolation
 
+        date = to_datetime(date)
         result = insolation([date], self.latitude_(), self.longitude_())
         return result.flatten()
 
+    def toa_incident_solar_radiation(self, date):
+        from earthkit.meteo.solar import toa_incident_solar_radiation
+
+        date = to_datetime(date)
+        result = toa_incident_solar_radiation(
+            date - datetime.timedelta(hours=1),
+            date,
+            self.latitude_(),
+            self.longitude_(),
+        )
+        return result.flatten()
+
+    def cos_solar_zenith_angle(self, date):
+        from earthkit.meteo.solar import cos_solar_zenith_angle
+
+        date = to_datetime(date)
+        result = cos_solar_zenith_angle(
+            date,
+            self.latitude_(),
+            self.longitude_(),
+        )
+        return result.flatten()
+
 
 class ConstantField:
     def __init__(self, date, param, proc, shape):

climetlab/version

@@ -1 +1 @@
-0.16.1
+0.16.2