CORDEX气象数据下载与处理

1. 数据集介绍

动力降尺度介绍What is regional downscaling? – Cordex

CORDEX网站主页Regional Climate Change simulations for CORDEX domains – Cordex

我所用到的domain-ESA主页 Region 7: East Asia – Cordex

PS: 关于CMIP6的CORDEX

引自周天军 , 邹立维 , 陈晓龙 . 第六次国际耦合模式比较计划（CMIP6）评述 [J]. 气候变化研究进展 , 2019, 15 (5): 445-456 目前还没有发布由CMIP6 GCM驱动的CORDEX，只有CMIP5的，关于CORDEX-CMIP6要求介绍（模拟年限和内容）Requests to CMIP6 – Cordex，最新状态见CORDEX simulations status | simulation-status

1. 数据下载（现在只能下CMIP5）

参考努力努力再努力搬砖博客【IDM下载CIMAP6-CORDEX data】_cmip6数据下载-CSDN博客

与CMIP6数据直接下载不同，注意一定要注册登录之后才能下载，建议使用谷歌批量下载插件，每次批量之前检查是否已经登录网站

1. 数据预处理

3.1 查看数据

下载数据之后用matlab查看信息（或者panoply打开）

ncdisp(ncname,'/', 'full');#ncname是含路径文件名

输出：

 rotated_latitude_longitude
           Size:       1x1
           Dimensions: 
           Datatype:   int32
           Attributes:
                       grid_mapping_name         = 'rotated_latitude_longitude'
                       grid_north_pole_latitude  = 61
                       grid_north_pole_longitude = -63.7
                       north_pole_grid_longitude = 0

发现CORDEX数据不是按照常见地理坐标系存的，采用的是旋转极投影，关于各种投影的介绍可以参考（一）Cartopy中的地图投影 - 码农教程，旋转极投影在第17个

那么现在问题来了，因为我们要对模式数据进行偏差校正（bias correction)，而常见的观测数据（observed data）都是地理坐标系（-180~180，-90~90），那么怎么把旋转极坐标下的模式nc文件转为地理坐标系呢？

3.2 转换坐标系

这是官网回答，见最后一条：

FAQ Data and Access – Cordex

这部分我找了很多资料，网页列在下面，总结可以有三种方法

如何将旋转的NetCDF转换回正常纬度/经度？ - 优文库

3.2.1 数学关系转换（附matlab/python代码）

① 参考coordinate system - Manually transforming rotated lat/lon to regular lat/lon? - Geographic Information Systems Stack Exchange（回答内含matlab/python代码）

经过用TLC点测试，行不通，可能因为代码中让输入南极点在旋转网格中的坐标，但是nc文件信息只给出了北极点的，简单的纬度减180并不ok

最后，关于如何用北极点坐标计算，找到了MathWorks上大佬的回答Rotated grid transform - File Exchange - MATLAB Central (mathworks.cn) ：SP_lon = NP_lon - 180, SP_lat = -NP_lat.

然后我这时候发现，nc文件里面给的lat,lon二维矩阵其实就是读取出的变量的经纬度（与变量矩阵对应行列）

这时我人已经傻了，浪费一整天时间就这，就这？然后我用代码测试了几个点，TLC点和其他5个点都测试成功，误差控制在0.01°内

② 参考大佬卷心没有菜博客nc文件曲面坐标转等距离坐标并对NC文件求四季平均值（内涵极地投影坐标系（极射赤道投影））_卷心没有菜的博客-CSDN博客

由于python和linux我都不是很熟，因此还没有进行过尝试

3.2.2 软件

①CDO，气象软件一霸，虽然我没用过，软件安装与使用参考卷心没有菜博客https://blog.csdn.net/qq_44907989/category_11742932.html，或者cdo随笔一：安装、插值与提取 - 墨天轮

答疑网站cdo remap /remapgrid help - CDO - Project Management Service

②CDE，Rotation of Coordinates Based On CORDEX Domains (agrimetsoft.com)

③ARCGIS

想通过gis转换坐标系，失败

3.3 插值

4. 总结

1.在下载数据之前没有将网站的介绍看完，尤其是Q&A页面，其实我关心的CMIP6和坐标系转换问题都在里面

2.下载完一个数据就先跑一次，而不是先下完所有的数据才开始处理（浪费大量时间）