diff --git a/tools/boundary/boundary.py b/tools/boundary/boundary.py
index 0dfe8f617..36adee86c 100644
--- a/tools/boundary/boundary.py
+++ b/tools/boundary/boundary.py
@@ -496,7 +496,8 @@ def add_coords(self, ds):
     def regrid_velocity(
                 self, usource, vsource, 
                 method='nearest_s2d', periodic=False, write=True, 
-                flood=False, fill='b', xdim='lon', ydim='lat', zdim='z', rotate=True, **kwargs):
+                flood=False, fill='b', xdim='lon', ydim='lat', zdim='z', rotate=True, 
+                time_attrs=None, time_encoding=None, **kwargs):
         """Interpolate velocity onto segment and (optionally) write to file.
 
         Args:
@@ -607,6 +608,12 @@ def regrid_velocity(
 
         ds_uv = self.rename_dims(ds_uv)
 
+        # Restore time attributes and encoding
+        if time_attrs:
+            ds_uv['time'].attrs = time_attrs
+        if time_encoding:
+            ds_uv['time'].encoding = time_encoding
+
         if write:
             self.to_netcdf(ds_uv, 'uv', **kwargs)
         
@@ -616,7 +623,8 @@ def regrid_tracer(
             self, tsource, 
             method='nearest_s2d', periodic=False, write=True, 
             flood=False, fill='b', xdim='lon', ydim='lat', zdim='z',
-            regrid_suffix='t', source_var=None, **kwargs):
+            regrid_suffix='t', source_var=None, 
+            time_attrs=None, time_encoding=None, **kwargs):
         """Regrid a tracer onto segment and (optionally) write to file.
 
         Args:
@@ -685,6 +693,12 @@ def regrid_tracer(
         
         tdest = self.rename_dims(tdest)
         tdest = tdest.rename({name: f'{name}_{self.segstr}'})
+
+        # Restore time attributes and encoding
+        if time_attrs:
+            tdest['time'].attrs = time_attrs
+        if time_encoding:
+            tdest['time'].encoding = time_encoding
         
         if write:
             self.to_netcdf(tdest, name, **kwargs)
diff --git a/tools/boundary/write_glorys_boundary_daily.py b/tools/boundary/write_glorys_boundary_daily.py
index 078cc9e83..8ca112740 100755
--- a/tools/boundary/write_glorys_boundary_daily.py
+++ b/tools/boundary/write_glorys_boundary_daily.py
@@ -31,6 +31,7 @@
 from os import path
 
 import xarray
+import numpy as np
 import yaml
 from boundary import Segment
 
@@ -53,18 +54,24 @@ def write_day(date, glorys_dir, segments, variables, output_prefix):
         return
 
     glorys = (
-        xarray.open_dataset(file_path)
+        xarray.open_dataset(file_path, decode_times=False)
         .rename({'latitude': 'lat', 'longitude': 'lon', 'depth': 'z'})
     )
 
+    # Capture time attributes and encoding
+    time_attrs = glorys['time'].attrs if 'time' in glorys.coords else None
+    time_encoding = glorys['time'].encoding if 'time' in glorys.coords else None
+
     for segment in segments:
         for variable in variables:
             if variable == 'uv':
                 print(f"Processing {segment.border} {variable}")
-                segment.regrid_velocity(glorys['uo'], glorys['vo'], suffix=f"{date:%Y%m%d}", flood=False)
+                segment.regrid_velocity(glorys['uo'], glorys['vo'], suffix=f"{date:%Y%m%d}", flood=False,
+                                        time_attrs=time_attrs, time_encoding=time_encoding )
             elif variable in ['thetao', 'so', 'zos']:
                 print(f"Processing {segment.border} {variable}")
-                segment.regrid_tracer(glorys[variable], suffix=f"{date:%Y%m%d}", flood=False)
+                segment.regrid_tracer(glorys[variable], suffix=f"{date:%Y%m%d}", flood=False,
+                                      time_attrs=time_attrs, time_encoding=time_encoding)
 
 def concatenate_files(nsegments, output_dir, variables, ncrcat_names, first_date, last_date, adjust_timestamps=False):
     """Concatenate annual files using ncrcat."""
@@ -97,7 +104,7 @@ def adjust_file_timestamps(file_path):
     Adjust timestamps for the first and last records in a file while preserving attributes and raw numerical format.
     """
     with xarray.open_dataset(file_path, decode_times=False) as ds:
-        # Explicitly load the dataset into memory to ovid lazy-loaded
+        # Explicitly load the dataset into memory if it's lazy-loaded
         ds.load()
 
         if 'time' in ds:
@@ -124,7 +131,7 @@ def adjust_file_timestamps(file_path):
             # Assign the new time variable back to the dataset
             ds = ds.assign_coords(time=new_time)
 
-            # Reapply the original encoding
+            # Reapply the original encoding to ensure consistency 
             ds['time'].encoding = time_encoding
 
             # Save the updated dataset