此前采用的是best track形式的飓风数据,鉴于NetCDF格式的数据产品应用更为广泛,此处以ERA5格网数据(0.25°×0.25°)为例进行演示。
一、ERA5数据下载
(相关资料图)
1. ERA5数据的时间跨度为1940年至今,可从ECMWF上下载
https://cds.climate.copernicus.eu/cdsapp#!/dataset/reanalysis-era5-single-levels?tab=form
2. ERA5数据提供网页下载方式和API方式,这里因为数据量较小,采用网页操作
3. 从下载好的数据中可以看到ERA5在存储经纬度时用的是一维,ADCIRC读取NetCDF文件时需要二维,不然程序会broken
二、ERA5数据格式转换
1. 目前在官网当中没有对NWS=14选项的描述,在wiki上找到了蛛丝马迹,但是描述很少
https://wiki.adcirc.org/NWS#Value_Seen_in_fort.15_File
直接读ERA5数据需要三个文件,其中fort.22有明确的描述。另外两个文件测试之后发现需对下载的ERA5经纬度坐标进行重新编辑
针对ERA5数据,fort.22内容为:
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放到一起运行即可
四、结果展示
此处仅展示同上一篇中提到的验潮站的结果对比,之前因为best track数据第6列标记错了一个地方,导致成图存在时延,这次改了过来。整体变化趋势上都比较接近,模式结果同实测的差距主要还是来自参考面以及控制文件的参数。
关键词: