У меня есть файл netcdf здесь , который я прочитал с ncdf4_1.16 и версиями R 3.5.0 и 3.5.1 (пробовал оба).
Я читаю переменную tsl с var=ncvar_get(nc,"tsl"). tsl - это двумерная переменная с размерами [3,87648]. Всякий раз, когда я читаю данные, последние 10 временных шагов в измерении 2 для каждого уровня измерения 1 каждый раз различаются, но всегда представляют собой очень маленькие числа (например, 1e-37). Когда я открываю тот же файл с помощью ncdump, он читает те же числа, что и NA, и я ожидаю, что эти последние временные шаги будут.
Я думал, что эти небольшие числа могут быть NaN, поэтому, на всякий случай, я попытался var[is.nan(var)] = NA, но это не сработало, так что я думаю, что нет. Я ничего не понимаю и волнуюсь, что это может случиться снова, когда я думаю, что с NA обращаются, но получается, что числа очень близки к 0 ...
Итак, короткий вопрос: что происходит?