Question

我想知道是否有人可以提供建议，

我正在尝试使用RNCEP包和函数NCEP.interp从RNCEP数据库中提取特定样本的降水率。

首先，我从.xlsx导入一个工作表，其中一列包含SampleYate，格式为％Y-％m-％d％H：％M：％S。

其次，我使用NCEP.interp，当它出现以下错误并停止时，它似乎运行到18％：

Error in if (format(dt.f, "%m-%d %H:%M:%S") > "12-31 17:59:59") { : 
missing value where TRUE/FALSE needed

还有一系列警告涉及下载长度，例如：

50: In download.file(paste("http://www.esrl.noaa.gov/psd/thredds/dodsC/Datasets/ncep.reanalysis",  ... : downloaded length 522 != reported length 200

但是，我不认为警告是一个问题（包括完整性）。我的完整代码是：

install.packages("RNCEP", dependencies =TRUE)
library(RNCEP)

SampPre6 <- read.xlsx("C:....Rainfall.xlsx",
                sheetName = "1")
SampPre6$Sample_Date <- strptime(as.character(SampPre6$Sample_Date), "%Y-%m-%d %H:%M:%S")
format(SampPre6$Sample_Date, "%Y-%m-%d %H:%M:%S")

PratePre6 <- NCEP.interp(variable='prate.sfc', level='gaussian', 
                   lat=SampPre6$Latitude,lon=SampPre6$Longitude,
                   dt=SampPre6$Sample_Date,reanalysis2=TRUE, 
                   keep.unpacking.info = TRUE)

数据的一个例子是：

   Sample_ID         Sample_Date Latitude Longitude
1        151 2013-03-15 23:29:59 31.39725  121.5068
2        177 2013-04-13 00:14:59 31.39695  121.5067
3        500 2013-04-13 06:59:59 23.26420  113.5147
4        501 2013-04-13 08:15:00 23.26168  113.5214
5        502 2013-04-13 09:00:00 23.27464  113.5529
6        247 2013-04-26 19:44:59 22.45018  114.1511
7        246 2013-04-26 21:14:59 22.45174  114.1346
8        245 2013-04-26 22:30:00 22.43341  114.1168
9        300 2013-05-18 00:59:59 31.39733  121.5067
10       584 2013-06-04 01:15:00 23.18805  113.4910

总之，任何人都可以了解导致进程停止的错误消息吗？

提前致谢，伊恩

Answer 1

dt参数应该是帮助文件中指定的“字符”类型？NCEP.interp

在运行NCEP.interp（）之前，请尝试包含以下内容......

SampPre6 $ Sample_Date＆lt; - format（SampPre6 $ Sample_Date，“％Y-％m-％d％H：％M：％S”）

NCEP.interp失败＆＃39;错误if（格式（dt.f，＆＃34;％m-％d％H：％M：％S＆＃34;）

1 个答案: