Question

我有一个每日降水数据的3维数组[lon，lat，365]和一个3维的播种日期文件[lon，lat，30]。每年在网格上播种一次。我有一个代码（如下所示）提取每个网格单元上的季节总降水量（播种至收获）。

但是，现在我只想提取每个网格单元播种日期后30-60天的降水量值（即，播种= Kharif季节的x，y网格单元上的157天，而305天是狂犬病作物的播种日期） x，y网格单元）。结果将是播种日期后30天[lon，lat，30days]的3维数组。播种后，即30-60天，80-100天，我需要不同的天数窗口。每个网格单元的播种日期可能不同，因此不同网格单元的播种后30天可能不同，但天数的长度（30-60天）应该相同。

非常感谢您根据播种（日期）和收获（日期）文件选择季节内特定天数的任何帮助，请

仅作为示例：

#precipitation aray
prec<-array(0,dim=nyear)

for (year in c(2:30)) #first year not, because growing season might have started in year 0
{
  print(year+1980)
  inputprec.file<-paste("E:/Paper_2018/Prec/","prec_",1980+year,".nc",sep="")
  inputprec.file.previous<-paste("E:/Paper_2018/Prec/","prec_",1980+year-1,".nc",sep="")

  nc<-nc_open(inputprec.file)
  lons<-ncvar_get(nc,"longitude")
  lats<-ncvar_get(nc,"latitude")
  val<-ncvar_get(nc,"pr")
  nc_close(nc)

  x<-which(abs(lons-lon)<0.01)
  y<-which(abs(lats-lat)<0.01)
  #for kharif season harvest date (260) is > sowing (157)

  if(hdate[year]>sdate[year])
  {
    prec[year]<-sum(val[x,y,sdate[year]:hdate[year]]) # season total     
    }else

  {
#for Rabi sowing start in November and harvets occurs in March-April  next year
    prec[year]<-sum(val[x,y,1:hdate[year]]) # code to extract winter season total precipitation from sdate: hdate spans over two years (1982-1983)
    # nc<-nc_open(inputprec.file.previous)
    # val<-ncvar_get(nc,"pr")
    # nc_close(nc)
    prec[year]<-prec[year]+sum(val[x,y,sdate[year-1]:365])

  }
}

使用功能从网格数据中提取亚季节性温度数据

0 个答案: