set decode_times to false

AugustinMortier · Sep 19, 2024 · 898c203 · 898c203
1 parent ed2533a
commit 898c203
Show file tree

Hide file tree

Showing 2 changed files with 11 additions and 4 deletions.
diff --git a/aprofiles/cli/aprocess.py b/aprofiles/cli/aprocess.py
@@ -125,7 +125,6 @@ def main(
 
         # list all files in out directory
         datepath = Path(basedir_out, yyyy, mm, dd)
-        onlyfiles = [str(e) for e in datepath.iterdir() if e.is_file()]
 
         if update_calendar:
             # create calendar
@@ -134,6 +133,8 @@ def main(
             if not path.is_file():
                 utils.json_calendar.make_calendar(basedir_out, yyyy, mm, calname)
 
+            # list all files in out directory
+            onlyfiles = [str(e) for e in datepath.iterdir() if e.is_file()]
             # add to calendar
             for file in tqdm(onlyfiles, desc="calendar     ", disable=disable_progress_bar):
                 utils.json_calendar.add_to_calendar(file, basedir_out, yyyy, mm, dd, calname)
@@ -146,11 +147,15 @@ def main(
             if not path.is_file():
                 utils.json_map.make_map(basedir_out, yyyy, mm, mapname)
 
+            # list all files in out directory
+            onlyfiles = [str(e) for e in datepath.iterdir() if e.is_file()]
             # add to map
             for file in tqdm(onlyfiles, desc="map          ", disable=disable_progress_bar):
                 utils.json_map.add_to_map(file, basedir_out, yyyy, mm, dd, mapname)
 
     if update_climatology:
+        # list all files in out directory
+        onlyfiles = [str(e) for e in datepath.iterdir() if e.is_file()]
         # get station id from file name
         stations_id = ["-".join(onlyfile.split("/")[-1].split("AP_")[1].split("-", 5)[:5]) for onlyfile in onlyfiles]
         # exclude moving stations

diff --git a/aprofiles/cli/utils/json_climatology.py b/aprofiles/cli/utils/json_climatology.py
@@ -28,9 +28,11 @@ def compute_climatology(basedir, station_id, variables, aerosols_only):
             if station_id in file and file.endswith(".nc"):
                 station_files.append(os.path.join(root, file))
 
+    print(station_files)
+
     try:
         # open dataset with xarray
-        ds = xr.open_mfdataset(station_files, parallel=False)
+        ds = xr.open_mfdataset(station_files, parallel=False, decode_times=False)
 
         # convert time index
         ds = convert_time_int_to_datetime(ds)
@@ -51,9 +53,9 @@ def compute_climatology(basedir, station_id, variables, aerosols_only):
         attrs["today"] = datetime.today().strftime("%Y-%m-%d")
 
         # seasonal resampling
-        Qds = ds.resample(time="Q").mean()
+        Qds = ds.resample(time="QE").mean()
         # add number of days per season as a new variable
-        Qds["ndays"] = ds.scene.resample(time="D").count().resample(time="Q").count()
+        Qds["ndays"] = ds.scene.resample(time="D").count().resample(time="QE").count()
 
         # work with selected variable