1. 环境准备与数据获取在开始处理ERA5风场数据之前我们需要准备好Python环境和相关数据。首先确保你安装了最新版的mikeio库1.0版本这个版本相比旧版有很多API变化这也是很多工程师遇到问题的根源。安装依赖库很简单使用pip命令即可pip install mikeio netCDF4 numpyERA5数据可以从欧洲中期天气预报中心(ECMWF)获取。这里有个小技巧如果你只需要特定区域和时间段的数据可以在下载时直接指定这些参数这样能大大减小文件体积。我通常会选择下载包含u10东西向风速、v10南北向风速和sp表面气压变量的数据。下载后的数据是NetCDF格式我们可以用以下代码快速检查数据内容import netCDF4 as nc file nc.Dataset(era5_wind_and_pressure.nc, r) print(file.variables.keys()) # 查看包含哪些变量 print(file.variables[time]) # 查看时间维度信息2. 数据预处理关键步骤ERA5数据的预处理有几个容易踩坑的地方我花了很长时间才搞明白。首先是纬度顺序问题 - ERA5的纬度是从大到小排列的90°到-90°而MIKE21默认是从小到大。如果不处理这个问题后续计算会出现各种奇怪的结果。解决方法是用numpy的flipud函数翻转纬度轴lat np.flipud(np.array(file.variables[latitude][:])) u np.flip(np.array(file.variables[u10][:]), axis1) v np.flip(np.array(file.variables[v10][:]), axis1)另一个常见问题是单位转换。ERA5的风速单位已经是m/s但有时需要根据实际情况进行调整。比如在某些项目中我发现将风速放大1.1倍能获得更好的模拟结果u u * 1.1 v v * 1.1时间处理也很关键。ERA5使用hours since 1900-01-01这样的相对时间表示法需要转换为Python的datetime对象time_dt [] tstart dt.datetime(1900, 1, 1, 0) for i in time: time_dt.append(tstart dt.timedelta(hoursint(i)))3. 创建Dfs2文件结构mikeio 1.0版本在创建Dfs2文件时有了很大变化。首先需要定义网格信息这里最容易出错的是投影参数设置。对于经纬度坐标应该使用LONG/LAT投影geometry mikeio.Grid2D( x0lon[0], nxlen(lon), dx0.25, y0lat[0], nylen(lat), dy0.25, projectionLONG/LAT )接下来定义变量信息。这里要注意EUMType和EUMUnit的准确使用错误的类型会导致MIKE21无法正确识别数据items [ ItemInfo(u, EUMType.Wind_Velocity, EUMUnit.meter_per_sec), ItemInfo(v, EUMType.Wind_Velocity, EUMUnit.meter_per_sec), ItemInfo(p, EUMType.Pressure, EUMUnit.pascal) ]4. 数据写入与验证将所有处理好的数据打包成Dataset对象时新版本的API变得更加简洁d [np.array(u), np.array(v), np.array(p)] dataset mikeio.Dataset(datad, timetime_dt, itemsitems, geometrygeometry)最后写入Dfs2文件。这里有个重要参数是时间间隔dtERA5数据通常是每小时一次所以设置为3600秒dfs Dfs2() dfs.write( datadataset, filenameoutput_wind.dfs2, dt3600, titleERA5 Wind Data )写入完成后建议用mikeio快速检查一下生成的文件ds mikeio.read(output_wind.dfs2) print(ds)5. 常见问题排查在实际项目中我遇到过几个典型问题。首先是内存不足错误处理大范围高分辨率ERA5数据时特别常见。解决方法是可以分时间段处理数据或者降低空间分辨率。另一个常见问题是时间对齐错误。有时候ERA5数据的时间戳不是整点这会导致MIKE21计算时出现问题。我的经验是强制对齐到整点时间time_dt [t.replace(minute0, second0) for t in time_dt]投影参数错误也很常见。如果你处理的是局部区域数据确保使用正确的投影字符串。我曾经花了两天时间才找出一个因为投影参数错误导致的计算偏差问题。最后记得检查生成的文件是否能被MIKE21正常读取。有时候Python端一切正常但MIKE21会报错这通常是因为变量定义不匹配或者单位设置有问题。