求助：我这有个提取最大风速的脚本，如何得到其对应的坐标点

15918702650

求助：我这有个提取最大风速的脚本，只能输出最大风速值，请问大神们如何添加命令可以得到最大风速对应的坐标呢（经纬度）。
a  = addfile("wrfout_d03_2017-08-22_06:00:00.nc","r")
times = wrf_user_getvar(a,"times",-1) ;    get all times in the file
ntimes = dimsizes(times)              ;number of times in the file
do i=0, ntimes-1
u = wrf_user_getvar(a,"U10",i)
v = wrf_user_getvar(a,"V10",i)
uv2d = sqrt(u*u+v*v)
print( i +" : " + max(uv2d))
end do
end

御风者

用maxind可以找到极大值点。

15918702650

御风者 发表于 2017-12-3 14:22
用maxind可以找到极大值点。

感谢解答，请问例子中X是指一个两维阵列，这个阵列是怎么表示呢
dims = dimsizes(X)
  x1d = ndtooned(X)      ; convert 2D array to 1D for use in maxind
  inds = ind_resolve(maxind (x1d), dims)    ; convert 1D array back to 2D
  ilat = inds(0,0)        ; select the latitude index where the X array is at its' maximum  
  ilon = inds(0,1)        ; select the longitude index where the X array is at its' maximum
  lat_max = X&lat(ilat)   ; insert the latitude index into the lat coordinate variable
  lon_max = X&lon(ilon)   ; insert the longitude index into the lon coordinate variable
  print("Maximum value located at "+lat_max+", "+lon_max)