write time as days

aprofiles/io/write_profiles.py

@@ -19,11 +19,11 @@ def write(profiles, base_dir, verbose):
         - verbose (bool): Verbose mode. Defaults to False.
     """
 
-    def _convert_time_after_epoch(ds, resolution='ms'):
-        time_attrs = ds["time"].attrs
-        ds = ds.assign_coords(time=ds.time.data.astype(f"datetime64[{resolution}]").astype(int))
-        ds["time"] = ds["time"].assign_attrs(time_attrs)
-        ds["time"].attrs['units'] = 'milliseconds since 1970-01-01T00:00:00'
+    def _convert_time_after_epoch(ds):
+        ds = ds.assign_coords(time=ds.time.data.astype(f"datetime64[ms]").astype(f'float32'))
+        # milliseconds to days
+        ds['time'] = ds['time'] / (1000 * 60 * 60 * 24) 
+        ds["time"].attrs['units'] = 'days since 1970-01-01T00:00:00'
         return ds
 
     def _classify_scene(ds):
@@ -131,7 +131,7 @@ def _classify_retrieval_scene(ds):
         ds_towrite = ds_towrite.drop(nodim_var)
 
     # converts time
-    ds_towrite = _convert_time_after_epoch(ds_towrite, resolution='ms')
+    ds_towrite = _convert_time_after_epoch(ds_towrite)
 
     # add altitude direction
     ds_towrite["altitude"] = ds_towrite["altitude"].assign_attrs({

docs/changelog.rst

@@ -4,6 +4,12 @@ Changelog
 .. image:: _static/images/history-solid.svg
    :class: awesome-svg
 
+0.6.5
+^^^^^^^
+Sep 20, 2024
+
+- Work on CF compliance (write time as days).
+
 0.6.4
 ^^^^^^^
 Sep 16, 2024

pyproject.toml

@@ -1,6 +1,6 @@
 [tool.poetry]
 name = "aprofiles"
-version = "0.6.4"
+version = "0.6.5"
 description = "Analysis of atmospheric profilers measurements"
 authors = ["augustinm <[email protected]>"]
 license = "GPL-3.0"