如何在R中使用ncdf4包正确读取气候研究数据（CRU）？

Question

CRU netCDF数据具有月平均值。我的文件包含1901-1910期间的数据。当我读取变量（tmp）时，它显示10年内120个月中只有12个等级。

library(ncdf4)
nc<- nc_open(file.choose())
lat<- ncvar_get(nc, "latitude")
lon<- ncvar_get(nc, "longitude")
tm<- ncvar_get(nc, "tmp")

dim(tm)
[1] 720 360  12

数据文件的链接在这里： - https://drive.google.com/open?id=0B9ACBRSUXCjnWWxZclR3alNIdEE

Answer 1

如果查看文件的标题

ncdump -h cru_tmp_clim_1901-1910.nc 

你会看到该文件只有12次。

netcdf cru_tmp_clim_1901-1910 {
dimensions:
time = 12 ;
latitude = 360 ;
longitude = 720 ;
bounds = 2 ;

如果你使用cdo：

cdo showdate cru_tmp_clim_1901-1910.nc 
1905-01-16  1905-02-15  1905-03-16  1905-04-16  1905-05-16  1905-06-16  1905-07-16  1905-08-16  1905-09-16  1905-10-16  1905-11-16  1905-12-16

你看到的数据仅适用于1905年。

我认为有人通过了一个已被删减的数据文件。我建议从源代码重新下载原始CRU数据。请记住，如果您有疑问，最好使用ncdump来示例文件头和元数据，以查看其中的内容。