#!/bin/usr/python3
import xarray as xr
from glob import glob
import sys

args = sys.argv
gcm = args[1]

data = sorted(glob('{}/*.nc'.format(gcm)))
if len(data) == 0:
    sys.exit()

print("RUnning: {}".format(gcm))
ds = xr.open_mfdataset(data, combine='nested', concat_dim='times')
out = "../merge/{}.nc".format(gcm)
print(out)
ds.to_netcdf(out)