dat2nc.py 文件源码

def create_netcdf(savef, dts, dat):
    """ Write Florida current data to netcdf file """
    dataset = Dataset(savef, 'w', format='NETCDF4_CLASSIC') 

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

    # Create data variable
    fc = dataset.createVariable('florida_current_transport',np.float64,(tdim.name),fill_value=1.e20)
    fc.units = 'Sv'
    fc[:] = dat

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