使用ERA5气象数据

使用ERA5气象数据

---

采用的是best track形式的飓风数据，鉴于NetCDF格式的数据产品应用更为广泛，此处以ERA5格网数据(0.25°×0.25°)为例进行演示。

一、ERA5数据下载

1. ERA5数据的时间跨度为1940年至今，可从ECMWF上下载

2. ERA5数据提供网页下载方式和API方式，这里因为数据量较小，采用网页操作

图 2. ERA5数据下载选项（针对Ike）

3. 从下载好的数据中可以看到ERA5在存储经纬度时用的是一维，ADCIRC读取NetCDF文件时需要二维，不然程序会broken

图 3. ERA5数据维度

二、ERA5数据格式转换

1. 目前在官网当中没有对NWS=14选项的描述，在wiki上找到了蛛丝马迹，但是描述很少

https://wiki.adcirc.org/NWS#Value_Seen_in_fort.15_File

图 4. NWS=14选项（V55以上）

直接读ERA5数据需要三个文件，其中fort.22有明确的描述。另外两个文件测试之后发现需对下载的ERA5经纬度坐标进行重新编辑

图 5. NWS=14时fort.22文件格式
针对ERA5数据，fort.22内容为：
time
time
hours
longitude
longitude
latitude
latitude
sp
u10
v10
2. 采用如下代码读取原始ERA5数据文件，将经纬度转换为二维

三、运行ADCIRC

1. 在fort.15里面将NWS对应的值改为14，WTIMINC修改为3600.0（因为这里用的ERA5数据为小时间隔，3600秒），并在末尾加上NCDATE的具体数值

2. 将fort.221.nc、fort.222.nc和fort.22三个文件同fort.14和fort.15放到一起运行即可

图 6. NWS=14运行成功
四、结果展示

此处仅展示同上一篇中提到的验潮站的结果对比，之前因为best track数据第6列标记错了一个地方，导致成图存在时延，这次改了过来。整体变化趋势上都比较接近，模式结果同实测的差距主要还是来自参考面以及控制文件的参数。

---

参考来源：B站