想请教一下怎么画最小值点

chanhui

这是我自己编写的gs脚本，不知道为什么运行时没有报错但是也没有画出位势的最小值点，参考了一些画台风路径的脚本，想知道到底是哪里出了问题
'reinit'
'Sdfopen D:\ruanjian\download-nc\2022-6-18-27-ece6b8baff59.nc'
'set grads off'
'set grid off'
'set lat 30 60'
'set lon 90 145'
'set z 1'
'set xlint 5'
'set ylint 5'
'set ylopts 1 5 0.12'
'set xlopts 1 5 0.12'
'set mpdset hires cnworld'
'd z/100'
'define locx=minloc(min(z,lat=30,lat=60),lon=90, lon=145)'
'define locy=minloc(min(z,lon=90, lon=145),lat=30,lat=60)'
'd locx'
locxx=subwrd(result,4)
'd locy'
locyy=subwrd(result,4)
'q gr2xy 'locxx' 'locyy''
x1=subwrd(result,3)
y1=subwrd(result,6)
'set strsiz 0.2'
'set line 2'
'draw mark 3 'x1' 'y1' 0.5'

tbag

我记得好像是需要sublin然后再subwrd

你先把'define locx=minloc(min(z,lat=30,lat=60),lon=90, lon=145)'之后所有语句去掉，然后'd locx'看看输出的是什么

river

如果你的GrADS是2.0.2以后的版本，建议使用最新的命令，更简单
amaxlocx()
amaxlocy()
aminlocx()
aminlocy()
具体用法自己去看官网介绍 http://cola.gmu.edu/grads/gadoc/gadocindex.html

chanhui

tbag 发表于 2022-12-2 04:41
我记得好像是需要sublin然后再subwrd

你先把'define locx=minloc(min(z,lat=30,lat=60),lon=90, lon=145 ...

试过，出来了一个数，可能是格点坐标，但是转换成经纬度却不成功

chanhui

river 发表于 2022-12-2 12:52
如果你的GrADS是2.0.2以后的版本，建议使用最新的命令，更简单
amaxlocx()
amaxlocy()

好的，谢谢您

栤唥de嶶笑

minloc()函数返回的是格点数，即第几个格点，实际定位出来的经纬度值，还要结合你的ctl里的分辨率和起点来计算。公式为