'''

Created on 02/02/2017

@author: ramato
'''
import iris
from iris.util import unify_time_units

DIR = '/project/ccrs/ramato/monitoring/anglo'

def main(): 
    for rid in ['ao554', 'ap163', 'ap448', 'ap430']:
        dir = '{}/mi-{}/apd'.format(DIR, rid)
        cubes = iris.load('{}/{}.????s.pr.nc'.format(dir, rid))
        unify_time_units(cubes)
        newcube = cubes.concatenate_cube()
        try:
            newcube.shape
            iris.save(newcube, '{}/{}.all.pr.nc'.format(dir, rid))
            print 'saved concatenated cube for', rid
        except AttributeError:
            print 'concatenation has failed for ', rid



if __name__ == "__main__":
    main()