问题标签 [cdo-climate]

我有两个 netcdf 文件：rsds.nc名为rsns.nc. rsds.nc包含一个名为 rsds 的变量和rsns.nc一个名为 rsns 的变量。现在我想rsus.nc通过分别减去rsds.nc和中的变量来获得上升辐射rsns.nc。

我尝试了以下方法：

他们都产生了一个rsus.nc，但这个文件中的变量 rsus 丢失了。知道为什么会这样吗？

正如标题所示，我想计算一个.nc名为sfcWind_1999.nc. 在每个网格点（纬度/经度组合）上，我想计算time维度中的中位数。是否有任何 nco 属性可以做到这一点？

我有一个我创建的多变量 ncdf，并想为每个变量添加额外的数据。纬度和经度维度将保持不变，但我想通过将新数据附加到每个变量来扩展时间维度。新数据集具有相同数量的变量，dim1 和 dim2，但其 dim3 开始于第一个数据集的 dim3 时间结束的位置。现有的 ncdf（有 42 个变量）：这是启动代码：

这是我之前的问题的后续，因为示例和数据更加具体。我有一个 304 x 448 的冰数据阵列（在极地立体投影中），其中的 netcdf 可以在这里获得。每个网格框对应的纬度/经度点可以在这个文件中找到。我想将这些数据转换为一个常规的纬度/经度网格（比如 180x360），其中每个网格框都是一个度数。有没有什么python方法可以做到这一点？等使用底图，或者使用 CDO（气候数据运营商）？

提前谢谢了。詹姆士

我有一个包含空气质量数据的 netcdf 文件。有4个维度：

我想制作一个只保留级别 1 的新 netcdf 文件，删除其他级别，因为我只需要表面。

我怎样才能做到这一点 ？

我尝试使用nccopy命令行执行此操作，但它复制了数据集。我不知道如何删除其他级别。

我正在尝试从由纬度/经度边界定义的用户指定（非矩形）域中的每个网格单元中提取数据。我的输入文件位于曲线网格上。我尝试了各种方法python，cdo，ncks，但我仍然无法弄清楚。我只想要输入 ncfile 的多边形域子集中的每个网格单元的时间序列信息。我的输入 ncfile 信息在这里给出：

我尝试过的一些事情是

我认为问题是我的输入文件的 x,y 维度以米为单位（没有负号，可能是“位置 1 处的无效字符”），我要求 cdo 提取 lat/lon 维度度。我的 ncfile 中有变量“映射”，它可能有助于将米转换为纬度/经度，但我不知道该怎么做。

并且 ncks 在这里对我不起作用，可能是由于我在 cdo 中看到的相同的仪表 <-> lon/lat 问题。

虽然我想要提取一个多边形，但我尝试过的这些示例只是矩形子集，我想得到多边形我可以做多个矩形子集来实现我的最终多边形形状，但是如果有更好的方法来做到这一点，任何建议都是赞赏。

谢谢

我有一个 netcdf 文件，其中的数据是 lon、lat 和 time 的函数。我想计算在时间维度上求和的每个网格单元中缺失条目的总数，最好使用 CDO 或 NCO，因此我不需要调用 R、python 等。

我知道如何获得缺失值的总数

正如我对这个相关问题的回答： count number of missing values in netcdf file - R

CDO 可以告诉我在空间上的总和

但我不知道如何随着时间的推移求和。有没有办法例如在 ncap2 中指定要与 number_miss 相加的维度？

我想从 cru 全局数据中对一个区域进行子集化。“cmsaf”包、box_mergetime 函数可以对 CMIP5 和 CORDEX .nc 数据进行子集化，但在 CRU .nc 数据中会出现以下错误。

数据文件的链接。

_FillValue 或 missing_value 是否仍然占用存储空间？

如果有一个包含一些空值的 2-dimission 数组，我如何将它写入 netcdf 文件以节省存储空间？

我想从 MERRA2 netcdf 文件中提取网格单元，这些文件仅使用 R 标记为海：

然而，由于土地区域下边界的弯曲性质，部分土地区域在我的选择范围内。我想选择最接近陆地的海域区域，但我的选择中有一些陆地区域我想过滤掉。

我有另一个 nc 文件，其中包含我理解的陆地和海域的一部分，可以从中提取我需要的网格单元：

所以我必须定义一个掩码，我只能使用它从我的第一个文件中选择具有一定阈值的海洋点，因为陆地和海洋的覆盖范围很小。请问有没有办法用 R 或者在 ArcGIS 中做到这一点