snapshot...

modules/util/split_nc4.py

 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