为什么用minloc函数无法定位到低压中心的最小值所在位置

爱喝柠檬汁

'reinit'
'sdfopen E:\grads\era.nc'
'set grads off'
'set mpdset hires'
'set map 15 1 1'
'set lev 500'
'set lat 55 65'
'set lon 100 125'

'set t 12'                     
'd z'
*******计算低压中心的经纬度
'define loa=minloc(min(z,lat=59,lat=61),lon=114,lon=116)'
'define lab=minloc(min(z,lon=114,lon=116),lat=59,lat=61)-90'

*******提取低压中心的经纬度
'd loa'
a=sublin(result,1)
lo=subwrd(a,4)
'd lab'
c=sublin(result,1)
la=subwrd(c,4)

*****************
'q w2xy 'lo' 'la''
xx=subwrd(result,3)
yy=subwrd(result,6)
***转换到图中的xy坐标

'draw mark 3 'xx' 'yy' 0.1'
'draw string 'xx' 'yy' D'
***在低压中心标点

'printim E:\grads\123.png white'
;

爱喝柠檬汁

楼上是相应的gs和图，
我想计算的低压中心是在112-114E，59-61N范围内，不知道为什么定位到外面去了
很费解，求各位路过的大神帮我看看呀

guzhong

你再看看你的
'define loa=minloc(min(z,lat=59,lat=61),lon=114,lon=116)'
'define lab=minloc(min(z,lon=114,lon=116),lat=59,lat=61)-90'
再说一遍你要的位置是哪。。

爱喝柠檬汁

'reinit'
'sdfopen E:\grads\era.nc'
'set grads off'
'set mpdset hires'
'set map 15 1 1'
'set lev 500'
'set lat 55 65'
'set lon 105 120'

'set t 12'                     
'd z'
*******计算低压中心的经纬度
'define loa=minloc(min(z,lat=59,lat=61),lon=112,lon=114)'
'define lab=minloc(min(z,lon=112,lon=114),lat=59,lat=61)-90'
*******提取低压中心的经纬度
'd loa'
a=sublin(result,1)
lo=subwrd(a,4)
'd lab'
c=sublin(result,1)
la=subwrd(c,4)
******
'q w2xy 'lo' 'la''
xx=subwrd(result,3)
yy=subwrd(result,6)
***转换到图中的xy坐标

'draw mark 3 'xx' 'yy' 0.1'
'draw string 'xx' 'yy' D'
***在低压中心标点

'printim E:\grads\123.png white'
;

爱喝柠檬汁

guzhong 发表于 2017-11-25 10:15
你再看看你的
'define loa=minloc(min(z,lat=59,lat=61),lon=114,lon=116)'
'define lab=minloc(min(z,lo ...

不好意思，之前那个是我笔误，新改的图在楼上，
确实定位不到中心呀，定到的是边界上的角（114E,61N）
不知道为什么？我觉得是这个函数使用的不对，或者有其他定位的方法吗？

爱喝柠檬汁

求各位路过的大神帮帮我呀，

爱喝柠檬汁