三维模式输出的数据很大，传输起来很麻烦。
有时候我们只需要几个关键变量，因此把他们抽取出来写成nc文件比较方便。
例如，要单独保存megan输出的isoprene变量
isop为三维变量[nt,nx,ny], lat lon为cmaq的grid文件，二维变量
将其写入一个文件
import netCDF4 as nc
da=nc.Dataset('./BVOC_emission.nc','w',format='NETCDF4')
da.createDimension("lat", lat.shape[0])
da.createDimension("lon", lat.shape[1])
da.createDimension("time", cmaqt.shape[0])
da.createVariable("isoprene","f8",("time","lat","lon"))
da.createVariable("lat","f8",("lat","lon"))
da.createVariable("lon","f8",("lat","lon"))
da.variables["isoprene"][:,:,:] = isop
da.variables["lat"][:,:] = lat
da.variables["lon"][:,:] = lon
da.close()
利用ncview查看
image.png
