# zhufei 2019-3-24
rm(list = ls())
library(ncdf4)
setwd('C:/Users/zhufe/Desktop/R_spatial') # 设置nc文件所在的路径
nc <- nc_open("E4E81D21036EAC6D84FEA2521AA84E3D_ferret_listing.nc")
pcp <- ncvar_get( nc = nc, varid = 'pcp')
lon <- ncvar_get( nc = nc, varid = 'lon')
lat <- ncvar_get( nc = nc, varid = 'lat')
#  建立索引
id <- rep(1:10,each = 8) # 根据pcp第三纬的长度，然后除以8，确定1:10这个变量
#id <- id[1:60623]      #  如果不能整除，变成整除然后取子集
e <- data.frame(id)     # 建立数据框，第一列是id
for (i in 1:80) {
# 数据的重新整理  变成一列
    pre <- as.matrix(pcp[,,i])  # 第三纬每一纬变成矩阵
    dim(pre) <- c(100*70,1)     #  按照经纬度重塑数据
    e[i,2:7001] <- pre          #  保存数据
}
# 每八列求和
a <- data.frame()
for (j in 1:10) {  # pcp第三维的长度确定j/8的长度
    a[j,1:7001] <- colSums(subset(e,id == j),
                           na.rm = T,dims = 1) # 每八行计算一天的降水量
}