使用ERA5气象数据
使用ERA5气象数据
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站