我有一个netcdf文件here,我使用ncdf4_1.16和R版本3.5.0和3.5.1(均尝试过)读取。
我正在用tsl
读取变量var=ncvar_get(nc,"tsl")
。 tsl
是一个二维变量,其维度为[3,87648]
。每当我读取数据时,维度2的每个级别在维度2中的最后10个时间步每次都会出现不同的值,但总是很小的数字(例如1e-37)。当我使用ncdump打开相同的文件时,它将读取与NA相同的数字,这正是我期望这些最后的步骤。
我认为这些较小的数字可能是NaN,因此,以防万一,我尝试了var[is.nan(var)] = NA
,但是它没有用,所以我想没有。我很无能为力,担心这种情况还会再次发生,因为我认为要解决NA问题,但结果是数字非常接近0 ...
那么,简短的问题是:发生了什么事?