diff --git a/modules/util/split_nc4.py b/modules/util/split_nc4.py
index 9b1bf38d8100ee15d22faad63249848564cba9f0..f2636a24679ddd4e8c585e3dc4dca4f6cf1c8c79 100644
--- a/modules/util/split_nc4.py
+++ b/modules/util/split_nc4.py
@@ -1,5 +1,6 @@
 import netCDF4 as nc
 import numpy as np
+import xarray as xr
 
 
 def split_dataset(input_file, output_pattern, dim_name, chunk_size):
@@ -55,5 +56,12 @@ def split_dataset(input_file, output_pattern, dim_name, chunk_size):
 
     ds.close()
 
+
+def concatenate_nc4_files(nc_files, output_file, concat_dim_name='time'):
+    datasets = [xr.open_dataset(nc_file) for nc_file in nc_files]
+    combined = xr.concat(datasets, dim=concat_dim_name)
+    combined.to_netcdf(output_file)
+    print(f"All files combined and saved to {output_file}")
+
 # Call the function
 # split_dataset('input.nc', 'output_{}.nc', 'time', 10)
\ No newline at end of file