如何选取wrfout里一定范围的数据？

qingcheng

我用官网里的台风路径脚本作图，其中有一部分是提取最低气压作为台风的坐标位置，但是由于台风起点的中心气压可能小于陆地的某个低压，台风起点跑到了陆地上（其他点的位置都没有问题），我想问问如何直接选取某一范围的气压，谢谢


用大括号截取的方法我有尝试过但是出现错误

sslp = wrf_user_getvar(a,"slp",0)
   slp=sslp({18:23},:)
  dims = dimsizes(slp)

fatal:Dimension (south_north) of (sslp) does not have an associated coordinate variable

Wetter

为什么要用中括号？直接slp=sslp(18:23,:)不行么

qingcheng

Wetter 发表于 2019-3-29 16:10
为什么要用中括号？直接slp=sslp(18:23,:)不行么

因为我看到加括号表示截取，我把括号去掉以后也出现了错误
fatal:Subscript out of range, error in subscript #0
fatal:An error occurred reading lon2d

Wetter

qingcheng 发表于 2019-3-29 16:31
因为我看到加括号表示截取，我把括号去掉以后也出现了错误
fatal:Subscript out of range, error in sub ...

这报的错误又不是这儿的吧，不是lon2d的问题么？

qingcheng

Wetter 发表于 2019-3-29 17:30
这报的错误又不是这儿的吧，不是lon2d的问题么？

是这的错误，应该是因为我没加上经度的范围，后来加上以后还是有错

sslp = wrf_user_getvar(a,"slp",0)
   slp=sslp(18:23,109:126)
  dims = dimsizes(slp)

fatal:Subscript out of range, error in subscript #1
fatal:An error occurred reading sslp
fatal:["Execute.c":8578]:Execute: Error occurred at or near line 27 in file track.ncl

Wetter

你能不能printVarSummary(sslp)一下啊

qingcheng

Wetter 发表于 2019-3-30 13:36
你能不能printVarSummary(sslp)一下啊

是这个嘛？

Variable: slp
Type: float
Total Size: 42000 bytes
            10500 values
Number of Dimensions: 2
Dimensions and sizes:   [south_north | 84] x [west_east | 125]
Coordinates:
Number Of Attributes: 6
  description : Sea Level Pressure
  units :       hPa
  FieldType :   104
  MemoryOrder : XYZ
  stagger :
  coordinates : XLONG XLAT

Wetter

qingcheng 发表于 2019-3-30 16:17
是这个嘛？

Variable: slp

这不就很明显了么，Subscript out of range, error in subscript #1
第二个纬度只有125行，最多到124（ncl是从0开始的），而你。

qingcheng

Wetter 发表于 2019-3-30 17:27
这不就很明显了么，Subscript out of range, error in subscript #1
第二个纬度只有125行，最多到124（n ...

这样啊，但是84x125这个是格点，我设的范围是按照经纬度的，我试着在格点范围内截取范围，数据就会乱套

Wetter

qingcheng 发表于 2019-3-30 18:15
这样啊，但是84x125这个是格点，我设的范围是按照经纬度的，我试着在格点范围内截取范围，数据就会乱套{: ...

你得找到你的经纬度最相邻的格点坐标啊