使用nco计算nc文件的中位数

Question

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

Answer 1

如manual中所述，NCO可以报告单个变量的中位数，但不能一步报告文件中的所有变量。因此，您必须自己创建循环变量，无论是在Bash中还是在ncap2中。

ncap2 -O -v -s "foo=gsl_stats_median_from_sorted_data(var_nm.sort());print(foo)" in.nc out.nc

要么

ncmdn var_nm in.nc

Answer 2

中位数是第50个百分点，因此您可以使用CDO中的百分位函数来执行此操作：

cdo timpctl,50 in.nc -timmin in.nc -timmax in.nc median.nc

该函数需要知道数据边界值以计算直方图，因此管道min和max命令。