Comment lire correctement les données de recherche sur le climat (CRU) en utilisant le paquetage ncdf4 dans R?

Question

Les données CRU netCDF ont des valeurs moyennes mensuelles. Mon fichier contient des données pour la période 1901-1910. Quand je lis la variable (tmp), elle montre seulement 12 niveaux sur 120 mois en 10 ans.

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 

lien de fichier de données est ici: - https://drive.google.com/open?id=0B9ACBRSUXCjnWWxZclR3alNIdEE

Answer 1

Si vous regardez l'en-tête du fichier

ncdump -h cru_tmp_clim_1901-1910.nc 

vous verrez que le fichier a seulement 12 tranches de temps.

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

Si vous utilisez 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 

vous voyez que les données sont seulement pour 1905.

Je pense que quelqu'un vous a passé un fichier de données qui a été coupé vers le bas. Je vous conseille de re-télécharger les données CRU d'origine de la source. Rappelez-vous, si vous avez des doutes, toujours bon d'utiliser ncdump pour exemple l'en-tête du fichier et les métadonnées pour voir ce qu'il y a dedans.