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
        fc = dataset.createVariable('florida_current_transport',np.float64,(tdim.name,))
        fc.units = 'Sv'
        fc[:] = self.fc

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