observations.py 文件源码

def write_to_netcdf(self, ncfile):
        """ Write observation data to netcdf file """

        # Open ncfile and create coords
        dataset = Dataset(ncfile, 'w', format='NETCDF4_CLASSIC')
        tdim = dataset.createDimension('time', None)

        # Create time coordinate
        time = dataset.createVariable('time',np.float64,(tdim.name,))
        time.units = 'hours since 0001-01-01 00:00:00.0'
        time.calendar = 'gregorian'
        time[:] = date2num(self.dates, time.units, calendar=time.calendar)

        # Create variables
        ek = dataset.createVariable('t_ek10',np.float64,(tdim.name,))
        ek.units = 'Sv'
        ek[:] = self.ekman

        umo = dataset.createVariable('t_umo10',np.float64,(tdim.name,))
        umo.units = 'Sv'
        umo[:] = self.umo

        fc = dataset.createVariable('t_gs10',np.float64,(tdim.name,))
        fc.units = 'Sv'
        fc[:] = self.fc

        moc = dataset.createVariable('t_moc_mar_hc10',np.float64,(tdim.name,))
        moc.units = 'Sv'
        moc[:] = self.moc

        # Close file
        print 'SAVING: %s' % ncfile
        dataset.close()