ncl提取一条等高线上的坐标

要提取一条等高线上的坐标，可以使用NCL中的`gsn_csm_contour_line`函数来绘制等高线，并且该函数会返回等高线上的坐标信息。具体操作步骤如下：

1. 读取数据文件，例如：

   f = addfile("filename.nc", "r")
   z = f->z

2. 定义绘图区域和绘图选项，例如：

   res = True
   res@cnFillOn = False
   res@cnLinesOn = True
   res@cnLineLabelsOn = False
   res@cnLevelSelectionMode = "ManualLevels"
   res@cnMinLevelValF = 10.0
   res@cnMaxLevelValF = 50.0

3. 调用`gsn_csm_contour_line`函数绘制等高线，并且返回等高线上的坐标信息，例如：

   contours = gsn_csm_contour_line(wks, z, res)
   x = contours(0)(0)(:,0)
   y = contours(0)(0)(:,1)

上述代码中，`contours`变量是一个三维数组，第一维表示等高线的数量，第二维表示每条等高线的段数，第三维表示每条等高线上的点数。因此，`contours(0)(0)`表示第一条等高线的第一段，`(:,0)`表示提取该段上所有点的x坐标，`(:,1)`表示提取该段上所有点的y坐标。

如果需要提取所有等高线上的坐标信息，可以使用类似以下代码：

   x = new((/0/), float)
   y = new((/0/), float)
   do i=0,dim_contours-1
      do j=0,dim_contours(i)-1
         x = fspan(x, contours(i)(j)(:,0))
         y = fspan(y, contours(i)(j)(:,1))
      end do
   end do

上述代码中，`dim_contours`是一个一维数组，表示每条等高线的段数。`fspan`函数用于将多个一维数组合并成一个一维数组。因此，`x`和`y`变量最终分别表示所有等高线上的点的x坐标和y坐标。