如何用ncL循环读取.txt文件里的经纬度，并提取出该点的模拟结果的降水值

lll1234567

begin
a = addfile("/usr/test/QLwrfout/xiajiwrfout_d02_2017-08-23_08:00:00.nc","r")
diri="/home/syrus/桌面/"  
  fname = "shanxijiangshui.txt"
  lines = asciiread(diri+fname,-1,"string")
rain_exp1 = wrf_user_getvar(a,"RAINNC",-1)
rain_con1 = wrf_user_getvar(a,"RAINC",-1)
rain_tot1 = rain_exp1 + rain_con1
times = wrf_user_getvar(a,"times",-1)
xlat=wrf_user_getvar(a,"XLAT",0)
xlon=wrf_user_getvar(a,"XLONG",0)
lat= tofloat(str_get_field(lines(1:),1," "))
lon= tofloat(str_get_field(lines(1:),2," "))
locationij=wrf_user_ll_to_ij(a,lon,lat,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)
pre=rain_tot1(:,locationij(1)-1,locationij(0)-1)
print(pre)
asciiwrite("test.txt",pre)
end
为什么出来后全是
9.96921e+36
9.96921e+36
9.96921e+36
9.96921e+36
9.96921e+36
9.96921e+36
9.96921e+36
9.96921e+36
9.96921e+36
9.96921e+36
9.96921e+36
9.96921e+36
9.96921e+36
9.96921e+36
9.96921e+36
9.96921e+36
9.96921e+36
9.96921e+36
9.96921e+36
9.96921e+36
9.96921e+36
9.96921e+36
9.96921e+36
9.96921e+36
9.96921e+36
9.96921e+36
9.96921e+36
9.96921e+36
9.96921e+36
9.96921e+36
9.96921e+36
9.96921e+36
9.96921e+36
9.96921e+36
9.96921e+36
9.96921e+36
9.96921e+36
9.96921e+36
9.96921e+36
9.96921e+36
9.96921e+36
9.96921e+36
9.96921e+36
9.96921e+36
9.96921e+36
9.96921e+36

我的尾巴呢

您好，我想问一下您的陕西降水资料应该是实测资料吧，您是在哪儿下的呢，我在中国气象服务网下载的全球地面气象站定时观测资料中没有经纬度。

lll1234567

我的尾巴呢 发表于 2018-3-23 22:55
您好，我想问一下您的陕西降水资料应该是实测资料吧，您是在哪儿下的呢，我在中国气象服务网下载的全球地面 ...

有经纬度的，你仔细看一下，可能用写字板打开才能看到

放逐流年

请问解决了吗？