关于ERA-Interim的问题

[Nianqitongs]

关于ERA-Interim的问题

[xuguangning1218]

具体步骤如下：
1、收集好CMA数据中台风中心的经纬度，由于CMA的台风中心经纬度是0.1度的，因此需要通过除以10转化为正常的经纬度
2、由于ERA-Interim目前已经暂停更新，因此建议使用CDS中的ERA5 Pressure Level数据集，里面包含了ERA-Interim的所有数据及相应更新数据
3、一般CDS数据下载有两种方式：网站下载和代码下载。其中网站下载无法选择相应的grid精度，建议使用代码下载
4、注册CDS，https://cds.climate.copernicus.eu/user/register?destination=%2F%23!%2Fhome
5、您可以登录帐号后，查看这个链接https://cds.climate.copernicus.eu/api-how-to配置您的机器，以获得代码下载权限
6、如果您不熟悉怎么配置下载参数，可以访问链接，通过网页选择您要下载的数据然后点击最下面的“Show API Request”按钮可以查看对应的生成的下载代码（通过代码可以增加grid参数选择对应的经纬度分辨率，如grid:[0.5,0.5]即表示经纬度为0.5度0.5度的空间分辨率，默认为0.25）或通过这个链接学习使用。下面为一个示例代码：
import cdsapi
c = cdsapi.Client()
c.retrieve(
'reanalysis-era5-pressure-levels',
{
'product_type': 'reanalysis',
'format': 'grib',
'variable': 'geopotential',
'pressure_level': '1000',
'year': '1979',
'month': '02',
'day': '02',
'time': [
'00:00', '01:00', '02:00',
],
'area': [
20, 80, -20,
120,
],
},
'download.grib')
7、了解如何使用代码下载CDS数据后，您可以通过第1步寻找到的台风中心经纬度，生成您想要下载的再分析数据的范围、对应时间等，并填入area及相应参数中，以实现个性化下载。示例代码中的area表示为一个20S20N,80E100E的4040的范围，其中经纬度分辨率为默认的0.25*0.25
8、关于下载的小技巧，建议您不要一次性下载所有需要的数据，分批逐个下载最佳，因为下载需要构建，生成下载。而且生成下载有效期，构建太多数据可能有效期内无法下载完成

[Nianqitongs]

作者您好，关于这个代码示例我还有个疑问，在选取下载数据的经纬度时无法跨越太平洋上180度经度线，是否如果要下载180度经度线两边的数据需要分开下载。比如说下面这张图上代表的代码范围就是东经95度到东经180度，不知道我理解的对不对。 念 ***@***.***  

…

------------------ 原始邮件 ------------------ 发件人: "xuguangning1218/Typhoon_Path" ***@***.***>; 发送时间: 2023年5月8日(星期一) 中午1:32 ***@***.***>; ***@***.******@***.***>; 主题: Re: [xuguangning1218/Typhoon_Path] 关于ERA-Interim的问题 (Issue #4) 具体步骤如下： 1、收集好CMA数据中台风中心的经纬度，由于CMA的台风中心经纬度是0.1度的，因此需要通过除以10转化为正常的经纬度 2、由于ERA-Interim目前已经暂停更新，因此建议使用CDS中的ERA5 Pressure Level数据集，里面包含了ERA-Interim的所有数据及相应更新数据 3、一般CDS数据下载有两种方式：网站下载和代码下载。其中网站下载无法选择相应的grid精度，建议使用代码下载 4、注册CDS，https://cds.climate.copernicus.eu/user/register?destination=%2F%23!%2Fhome 5、您可以登录帐号后，查看这个链接https://cds.climate.copernicus.eu/api-how-to配置您的机器，以获得代码下载权限 6、如果您不熟悉怎么配置下载参数，可以访问链接，通过网页选择您要下载的数据然后点击最下面的“Show API Request”按钮可以查看对应的生成的下载代码（通过代码可以增加grid参数选择对应的经纬度分辨率，如grid:[0.5,0.5]即表示经纬度为0.5度0.5度的空间分辨率，默认为0.25）或通过这个链接学习使用。下面为一个示例代码： import cdsapi c = cdsapi.Client() c.retrieve( 'reanalysis-era5-pressure-levels', { 'product_type': 'reanalysis', 'format': 'grib', 'variable': 'geopotential', 'pressure_level': '1000', 'year': '1979', 'month': '02', 'day': '02', 'time': [ '00:00', '01:00', '02:00', ], 'area': [ 20, 80, -20, 120, ], }, 'download.grib') 7、了解如何使用代码下载CDS数据后，您可以通过第1步寻找到的台风中心经纬度，生成您想要下载的再分析数据的范围、对应时间等，并填入area及相应参数中，以实现个性化下载。示例代码中的area表示为一个20S20N,80E100E的4040的范围，其中经纬度分辨率为默认的0.25*0.25 8、关于下载的小技巧，建议您不要一次性下载所有需要的数据，分批逐个下载最佳，因为下载需要构建，生成下载。而且生成下载有效期，构建太多数据可能有效期内无法下载完成 — Reply to this email directly, view it on GitHub, or unsubscribe. You are receiving this because you authored the thread.Message ID: ***@***.***>

[xuguangning1218]

作者您好，关于这个代码示例我还有个疑问，在选取下载数据的经纬度时无法跨越太平洋上180度经度线，是否如果要下载180度经度线两边的数据需要分开下载。比如说下面这张图上代表的代码范围就是东经95度到东经180度，不知道我理解的对不对。 念 @.***  
…

[是的，如果您因为选择的格点很大，出现这个情况，则需要下载两份再拼接再一起，如下图所示：

[Nianqitongs]

好的，谢谢您！真是帮了大忙了。 念 ***@***.***  

…

------------------ 原始邮件 ------------------ 发件人: "xuguangning1218/Typhoon_Path" ***@***.***>; 发送时间: 2023年5月8日(星期一) 下午3:20 ***@***.***>; ***@***.******@***.***>; 主题: Re: [xuguangning1218/Typhoon_Path] 关于ERA-Interim的问题 (Issue #4) 作者您好，关于这个代码示例我还有个疑问，在选取下载数据的经纬度时无法跨越太平洋上180度经度线，是否如果要下载180度经度线两边的数据需要分开下载。比如说下面这张图上代表的代码范围就是东经95度到东经180度，不知道我理解的对不对。 念 @.***   … [是的，如果您因为选择的格点很大，出现这个情况，则需要下载两份再拼接再一起，如下图所示： — Reply to this email directly, view it on GitHub, or unsubscribe. You are receiving this because you authored the thread.Message ID: ***@***.***>