今天试了一下用PyNIO写nc文件，难用的一腿

算了，还是用netcdf模块写了。

代码如下：

# -*- coding: utf-8 -*-

import numpy as np
import sys
import os
from netCDF4 import Dataset

def creatspinc(value, filename):

gridspi = Dataset(filename, 'w', format='NETCDF4')

# dimensions

gridspi.createDimension('time', None)

gridspi.createDimension('lat', 1000)   #len(lat)

gridspi.createDimension('lon', 2000)

# Create coordinate variables for dimensions

times = gridspi.createVariable('time', np.float64, ('time',))

latitudes = gridspi.createVariable('lat', np.float32, ('lat',))

longitudes = gridspi.createVariable('lon', np.float32, ('lon',))

# Create the actual variable

var = gridspi.createVariable('var', np.float32, ('time', 'lat', 'lon',))

# Global Attributes

import time

gridspi.description = 'var'

gridspi.history = 'Created ' + time.ctime(time.time())

gridspi.source = 'netCDF4 python module tutorial'

# Variable Attributes

latitudes.units = 'degree_north'

longitudes.units = 'degree_east'

times.units = 'days since 2019-01-01 00:00:00'

times.calendar = 'gregorian'

# data

lats = np.linspace(25.,35.,1000) #notice: the last numb is not included

lons = np.linspace(95.,115.,2000) #notice: the last numb is not included

latitudes[:] = lats

longitudes[:] = lons

#--Fill in values

var[0:value.shape[0],:,:] = value

# Fill in times

from datetime import datetime, timedelta

from netCDF4 import num2date, date2num

dates = []

for n in range(var.shape[0]):

dates.append(datetime(2019, 1, 1) + n * timedelta(days=1))

times[:] = date2num(dates, units = times.units,calendar = times.calendar)

print 'time values (in units %s): ' % times.units +'\n', times[:]

dates = num2date(times[:], units=times.units, calendar=times.calendar)

gridspi.close()

return

out_dir='~/'
out_file_nc = out_dir+'out.nc'
data_nc = np.zeros(shape=(61,1000,2000))

creatspinc(data_nc, out_file_nc)

最后这一段是我随便写了个变量调用写.nc文件的函数来测试的。

最近有小伙伴问我为什么公众号文章不能留言，我百度了一下，说这几个月开通的公众号账户文章都没有留言功能，也不知道以后会不会给留言。不过可以赞赏的哟！