import xarray as xr
merged = xr.open_dataset("path.nc")
laguar = ['spatial_ref_laguar', 'ndvi_laguar','geo_laguar','soil_laguar','meteo_laguar','lon_0_laguar','lat_0_laguar','time_0_laguar', 'lon_1_laguar', 'lat_1_laguar','geo_property_laguar', 'lon_2_laguar', 'lat_2_laguar','lon_3_laguar', 'lat_3_laguar', 'time_3_laguar']
negrar = ['spatial_ref_negrar','ndvi_negrar', 'geo_negrar', 'soil_negrar', 'meteo_negrar', 'lon_0_negrar','lat_0_negrar','time_0_negrar', 'lon_1_negrar','lat_1_negrar','geo_property_negrar','lon_2_negrar','lat_2_negrar','lon_3_negrar', 'lat_3_negrar', 'time_3_negrar']
L = merged[laguar]
L.attrs={"datasets_names":['laguar']}
N = merged[negrar]
N.attrs={"datasets_names":['negrar']}
N.to_netcdf("negrar_filtered.nc")
L.to_netcdf("laguar_filtered.nc")