Assign time when creating coords

Author: havardhhaugen

bris/outputs/verif.py

@@ -222,8 +222,11 @@ def create_nan_array(shape, dtype=np.float32) -> np.ndarray:
     def finalize(self) -> None:
         """Write forecasts and observations to file"""
 
+        frts = self.intermediate.get_forecast_reference_times()
+        frts_unix = utils.datetime_to_unixtime(frts).astype(np.double)
+
         coords = {}
-        coords["time"] = (["time"], [], cf.get_attributes("time"))
+        coords["time"] = (["time"], frts_unix, cf.get_attributes("time"))
         coords["leadtime"] = (
             ["leadtime"],
             self.intermediate.pm.leadtimes.astype(np.float32) / 3600,
@@ -270,9 +273,6 @@ def finalize(self) -> None:
 
         self.ds = xr.Dataset(coords=coords)
 
-        frts = self.intermediate.get_forecast_reference_times()
-        self.ds["time"] = utils.datetime_to_unixtime(frts).astype(np.double)
-
         # Load forecasts
         fcst_shape = (
             len(frts),