python创建nc

### 如何使用Python创建NetCDF（.nc）文件

为了创建NetCDF文件，需先安装`netCDF4`库。可以通过pip安装此库：

pip install netcdf4

下面是一个完整的例子展示如何利用`netCDF4`和`numpy`模块创建一个新的NetCDF文件，并向其添加一些基本信息。

#### 导入库与初始化设置

import netCDF4 as nc  # 导入netCDF4用于操作NetCDF文件
import numpy as np    # 导入numpy用于数值计算

fn = 'example_file.nc'  # 定义目标路径及文件名为example_file.nc
ds = nc.Dataset(fn, 'w', format='NETCDF4')  # 创建新的NetCDF文件对象，指定写入模式('w')

#### 设置全局属性

ds.description = "这是一个通过netCDF4-Python创建的例子"
ds.history = "Created for demonstration purposes only."
ds.source = "netCDF4 Python module tutorial"

#### 增加维度

time_dim = ds.createDimension('time', None)       # 创建时间维度，长度设为None意味着无限大
lat_dim = ds.createDimension('latitude', 73)      # 创建纬度维度，固定大小73
lon_dim = ds.createDimension('longitude', 144)    # 创建经度维度，固定大小144

#### 添加变量

times = ds.createVariable('time', 'f8', ('time',))        # 时间轴作为浮点数存储
lats = ds.createVariable('latitude', 'f4', ('latitude',)) # 纬度数组作为单精度浮点数保存
lons = ds.createVariable('longitude', 'f4', ('longitude',))
temperature = ds.createVariable('temp', 'f4', ('time', 'latitude', 'longitude'))

#### 赋予单位和其他元数据给变量

lats.units = "degrees_north"                          # 给定纬度的单位
lons.units = "degrees_east"                           # 经度的单位
times.units = 'hours since 0001-01-01 00:00:00.0'
times.calendar = 'gregorian'

temperature.standard_name = "air_temperature"         # 温度的标准名称
temperature.long_name = "Air Temperature (instantaneous)"
temperature.units = "K"                               # 单位开尔文(K)

#### 关闭文件连接

完成所有的编辑之后，记得关闭文件以确保所有更改都被正确保存到磁盘上[^2]。

ds.close()