问题标签 [ncl]

有没有一种快速简单的方法可以从 bash 的命令行将 HDF5 文件转换为 netcdf(4)？或者一个简单的脚本在 R、NCL 或 python 中自动处理这种转换？

我想用TTF字体在NCL中绘制中文单词。我已将 TTF 字体转换为 svg 路径。如果我能找到一种在 NCL 中绘制 svg 路径的方法，我想我们可以为所有语言绘制每个单词，没有任何限制。

svg 路径数据如下：

d="M1204 1679h162v-677q273 140 510 307l100 -144q-216 -136 -610 -333v-735q0 -136 150 -136h140q144 0 162 136q18 126 32 328q88 -36 170 -64q-20 -190 -44 -324q-38 -228 -282 -228h-216q-274 0 -274 270v1600zM26 99q261 90 618 236v630h-576v148h576v562h162v-1922 h-162v415q-192 -78 -554 -227z"

我相信标签 M、h、v、z 很简单。但是如何处理标签 q,c？我们如何绘制贝塞尔曲线？

我正在使用打开其他 NCL 脚本的 NCL 脚本。似乎在同一目录中的 NCL 脚本可以正常打开，但是如果它们应该从另一个目录（具有相同的父目录）中获取，则无法打开它们。似乎 NCL 没有设置正确的“主”目录。我该怎么做呢？我通过 conda 安装程序 ( https://www.ncl.ucar.edu/Download/conda.shtml ) 获得了 NCL，并且那里提到的测试运行良好。

我有一个全年累积降水的 3D nc 变量 PP = (time, lat, lon)，我想计算每日值。

我想在示例中做一些链接，但找不到任何示例。

我试图做一个循环，其中减法应该在 [index+1] 减去值 [index] 的值之间（如下所示）

但这种方式会产生错误。

我需要每天的值，并且只将累积的降水量作为输入信息，并认为这样的循环应该是最好的方法。

如何从变量中进行计算？？？

我正在合作，ncl但如果你有cdo或有nco提议，欢迎他们。

我正在尝试将代码从 NCL 转换为 PyNGL/PyNio。正如我在过渡指南中所读到的，为绘图例程处理资源文件的方式应该与在 NCL 中的方式相同。但是，在我的脚本中，python 不包含资源文件中的任何资源，即使我遵循了转换指南中的基本示例。这可能是因为我在 conda 环境中工作吗？还是 PyNGL 中存在错误？

我在 Conda 环境中使用 Python 3.7。我已经用我自己的代码以及下面包含的一个非常简单的示例代码进行了尝试。

文件 title_app.png 中出现的文本既不显示为红色也不显示为粗体。

我有以下shapefile和netcdf 文件。

我想从包含在 shapefile 边界内的 netcdf 文件中提取数据。

您对我如何实现这一目标有什么建议吗？

shapefile 对应于 SREX 区域 11 North Europe (NEU)，netcdf 文件是 CMIP6 气候模型数据输出的示例（ua 变量）。我想要的输出必须是 netcdf 格式。

更新

到目前为止，我尝试使用 NCL 和 CDO 创建一个 netcdf 掩码，并将此掩码应用于原始 netcdf 数据集。下面是步骤（和NCL 脚本）：

输出几乎是正确的。见下图。但是 shapefile 的南部边界（SREX 11，NEU）没有被正确捕获。所以我想生成 netcdf 掩码的 NCL 脚本有问题。

我有很多 netCDF 文件，它们是使用 python gdal 库从 tif 文件创建的。我想将这些文件中的每一个合并到一个大型 nc 文件中，但它们没有时间变量。

我尝试了此处描述的方法，但是当我要查看生成的大型 nc 文件的每个波段的时间戳时，所有文件都为零。

此外，结果ncdump -h name-of-the-file如下。

我想要做的是创建一个时间变量，它与全局变量相同：GDAL_TEMPORAL_STOP 并将每个文件合并到一个大的 nc 文件中。理想情况下，我更喜欢在 python 中执行此操作。我认为可能有某种方法可以复制并将其global variable:GDAL_TEMPORAL_STOP = "2017-01-10" ;复制为每个文件的时间变量，但我找不到方法，所以这就是我需要你帮助的地方。

我在 python 方面没有经验，我正在将 NCL 脚本转换为 python，希望 python 运行得更快。环顾四周，我没有找到我认为是 NCL 脚本中最简单计算的答案。看看如何完成更艰难的计算，我也没有找到如何在 python 中完成这些的答案。

大部分计算是在将 3 维变量转换为 1 维变量并查询它们在数组空间中的值和位置之后完成的。知道了 t 个变量在数组空间中的位置，我们可以得到对应于 t 个变量整数值的 p 个变量值。

计算如下：

• 将 p 变量中的值设置为其默认值 _FillValue，
• 计算网格点的数量（体积） t 变量中每个可能的整数值都会出现一个值（时间和空间的总和），
• 计算 t 变量中每个可能的整数值的开始和结束时间索引，
• 将持续时间计算为 t 变量中结束时间和开始时间之间的差值（+1，因为数字），
• 计算 t 变量中每个可能的整数值的平均（时空）纬度和经度，
• 计算 t 变量中每个可能的整数值的面积（体积/持续时间），
• 从 p 变量计算平均 p，其中它在空间时间上对应于 t 变量中每个可能的整数值，并且
• 从 p 变量计算 p 百分位数，其中它在时空中对应于 t 变量中的每个可能的整数值。

所有这些计算都将值保存在一维数组中，维度大小等于 t 变量中的最大整数值。例如，at 变量可能有 0 到 100 之间的整数。0 整数值被忽略，因此示例中的每个一维数组应该有 100 个值；（100 卷、100 个开始时间、100 个结束时间等）。

最后，将所有一维数组写入（制表符分隔的）文本文件，每一列都是一维数组。

我需要用 NetCDF 文件做一些计算。所以我有两个具有以下尺寸和大小的变量：

我需要的是生成一个新数组 ，C它的形状是(1,12,84,228)内容B传播到A.

通常，这很容易在 NCL 中使用 conform 函数完成。我不确定在 Python 中这个等价物是什么。

谢谢你。

我有两个 .nc 文件：zmksrf_lai.nc和zmksrf_pft.nc。我想将 0.5 度的 LAI 文件插入 0.05 度的 PFT 文件，这样我就有了高分辨率的 LAI 数据。我尝试了 CDO，也尝试了 ncremap，但没有奏效。非常感谢您的建议。

和

但是当我这样做时

使用 NCO，我尝试了 ncremap 并出现以下错误。

所有必要的模块都已加载。两个文件的元数据/维度是：

和