#!/bin/usr/python

import netCDF4 as net
import sys

# Open base file to write merged dataset into
dg = net.Dataset('base.nc', 'a')
prec = dg.variables['PREC'][:]


for i in range(1,69):
#for i in range(67,69):
    si = str(i).zfill(2)
    pfile = 'PREC/PREC_p{}.nc'.format(si)
    print pfile
    pc = net.Dataset(pfile)
    tprec = pc.variables['PREC'][:,0,:]
    prec[:,i-1,:] = tprec
    
print sum(sum(prec))

# Write and close file
dg.variables['PREC'][:] = prec
dg.close()