ncl读取某个经纬度附近的温度值

hahahahaha

但是一直报错，所以代码和报错贴出来，希望有大神能给点意见
代码：begin
a = addfile("/usr/local/wrf/WRFV3/test/em_real/wrfout_d01_2015-11-01_00:00:00.nc","r")
T2=wrf_user_getvar(a,"T2",0)
times = wrf_user_getvar(a,"times",-1)
xlat=wrf_user_getvar(a,"XLAT",0)
xlon=wrf_user_getvar(a,"XLONG",0)
locationij=wrf_user_ll_to_ij(a,121.7854194,31.19645278,True);
;print(locationij)
lonlat=wrf_user_ij_to_ll(a,locationij(0),locationij(1),True);
print(lonlat)
;xx=xlat(locationij(1)-1,locationij(0)-1)
;print(xx)
t=T2(:,locationij(1)-1,locationij(0)-1)
print(t)
asciiwrite("test.txt",t)
end

hahahahaha

报错：
Copyright (C) 1995-2013 - All Rights Reserved
University Corporation for Atmospheric Research
NCAR Command Language Version 6.1.2
The use of this software is governed by a License Agreement.
See http://www.ncl.ucar.edu/ for more details.


Variable: lonlat
Type: float
Total Size: 8 bytes
            2 values
Number of Dimensions: 1
Dimensions and sizes:   [lon_lat_location | 2]
Coordinates:
(0)     121.6453
(1)     31.39188
fatal:["Execute.c":7432]:performASSIGN_VAR_VAR_OP: Number of subscripts on rhs do not match
                        number of dimensions of variable, (3) Subscripts used, (2)

fatal:["Execute.c":8128]:Execute: Error occurred at or near line 17 in file test.ncl


提示是t=T2(:,locationij(1)-1,locationij(0)-1)这一句有问题

freekiller

改为 t=T2(locationij(1)-1,locationij(0)-1)

你只读了一个时次

hahahahaha

freekiller 发表于 2016-6-15 17:27
改为 t=T2(locationij(1)-1,locationij(0)-1)

你只读了一个时次

出现了新的问题：读取T2是发生错误
Copyright (C) 1995-2013 - All Rights Reserved
University Corporation for Atmospheric Research
NCAR Command Language Version 6.1.2
The use of this software is governed by a License Agreement.
See http://www.ncl.ucar.edu/ for more details.


Variable: lonlat
Type: float
Total Size: 8 bytes
            2 values
Number of Dimensions: 1
Dimensions and sizes:   [lon_lat_location | 2]
Coordinates:
(0)     121.6453
(1)     31.39188
fatal:Subscript out of range, error in subscript #0
fatal:An error occurred reading T2
fatal:["Execute.c":8128]:Execute: Error occurred at or near line 17 in file test.ncl

Juliustan

没看懂，不过还是mark一下，日后好相见

hahahahaha

Juliustan 发表于 2016-6-16 08:56
没看懂，不过还是mark一下，日后好相见

我也不懂。。。

freekiller

hahahahaha 发表于 2016-6-16 07:33
出现了新的问题：读取T2是发生错误
Copyright (C) 1995-2013 - All Rights Reserved
University Corp ...

经纬度写反了， 不在模式区域了               

hahahahaha

freekiller 发表于 2016-6-16 11:08
经纬度写反了， 不在模式区域了

但是我用了WPS中namelist.wps里设置的经纬度，还是同样的报错，好奇怪

freekiller

hahahahaha 发表于 2016-6-16 12:29
但是我用了WPS中namelist.wps里设置的经纬度，还是同样的报错，好奇怪

你再试试，否则这么简单的脚本，不应该出错

hahahahaha

freekiller 发表于 2016-6-16 14:50
你再试试，否则这么简单的脚本，不应该出错

试了，确实还是那样。不管加不加：，