日气温偏热天数答案

【问题标题】：Number of hot days from daily temperature日气温偏热天数
【发布时间】：2022-01-11 01:00:06
【问题描述】：

我正在关注answer，根据每日 tmax 数据计算一年中炎热的天数（温度超过 35 摄氏度）。

我在 2000 年至 2016 年期间使用 CHELSA 的 tmax，并根据我的边界框要求对其进行了裁剪。

这是我已经完成的步骤（使用 2001 年数据的示例 - 1 个月 1 nc 文件）：

将月度数据合并到年度：cdo mergetime chelsa_daily_2001*.nc chelsa_annual_2001.nc
计算热天：cdo gec,308.15 chelsa_annual_2001.nc chelsa_hotdays_2001.nc Chelsa 的温度是开尔文，所以热天的阈值是 308.15
一年中的总天数：cdo yearsum chelsa_hotdays_2001.nc chelsa_hotdays_yearsum_2001.nc

下面是结果，不幸的是不是我所期望的。

为什么天数不是整数？我是否遗漏了脚本中的某些内容？

UPDATE1（根据 Adrian 的回复）

我已经通过自制软件安装了 ncview，但不幸的是无法打开。出现以下错误：

Note: could not open file /Users/xxx/.ncviewrc for reading
Error: Can't open display:

我尝试使用QGIS打开nc输出，结果还是float。

更新2

好的，我设法使用 ncdump 对其进行了检查，这是包含该值的第一行。有点困惑，因为我尝试使用1年的数据，总数超过365。这是怎么发生的？

【问题讨论】：

【解决方案1】：

我强烈怀疑 panoply 正在对导入的数据执行某种空间插值。

请像这样直接使用ncdump查看原始字段

ncdump chelsa_hotdays_yearsum_2001.nc | less

（我将管道调到 less 以便您可以浏览数据）。或者您也可以打开ncview 中的文件并将光标移到数据上，您将看到对话框中显示的字段值。

ncview chelsa_hotdays_yearsum_2001.nc

【讨论】：

谢谢阿德里安，我尝试安装 ncview 但失败了。并尝试使用 qgis 打开，结果仍处于浮动状态。请参阅更新的问题。你觉得我写的剧本对吗？
脚本看起来不错，如果你这样运行它，那么 chelsa_hotdays_2001.nc 必须只包含 1 和 0 以及最后一个文件的整数天数。 QGIS 图没有告诉您任何信息，它只是将范围划分为一组浮点颜色。您需要以某种方式“查看”原始值，或者使用 ncdump，或者您也可以将其读入 python 或 R 中。如果您更新您的问题以让我们知道您在哪个系统上，这将有所帮助，然后我们可以就安装提供建议。或者，如果你在 Windows 中，我们可以提供一个 3 行脚本来打印内容......但我相信它很好
我正在使用 macOS 10.5 并设法使用 ncdump chelsa_hotdays_yearsum_2001.nc | less 并滚动到包含值的第一行。你是对的，结果是整数，但令人惊讶的是总数超过 365，而我只使用 1 年的数据。查看更新的问题。