ncl 160站数据读取记录

biguo

记录：筛选160个站点中某个区域内的站点数据。zd.txt是指160站数据的帮助文件，即包含站点序号，站点经纬度的文件。

;筛选出经纬度范围内的站点序号
begin
f=asciiread("/cygdrive/c/Users/admin/Desktop/tt/zd.txt",(/160,3/),"float")
j=0
seri_number=new((/1,43/),integer)
do i=0,159
    if(f(i,2).ge.105..and.f(i,2).le.130..and.f(i,1).ge.16..and.f(i,1).le.30.) then
    print(i+0)
    seri_number(:,j)=i
    j=j+1
end if
end do
  print(j)
  print(seri_number)
;---------------------保存筛选出的站点序号，后面分析使用-------------------
opt=True
filename="105-130Eseri_number.txt"
opt@fout=filename
SYS_RM="/bin/rm -f"+filename
system(SYS_RM)
write_matrix(seri_number, "10i7", opt)   ;10表示10列的意思   

;-----------------------读取筛选站点的降水数据-----------------------------
f2=asciiread("/cygdrive/f/yjs/sj/160pre/r1602.txt",(/69,160/),"float")
f3=asciiread("/cygdrive/f/yjs/sj/160pre/r1603.txt",(/69,160/),"float")
seri_num=asciiread("105-130Eseri_number.txt",(/43/),"integer")
print(seri_num)
n=dimsizes(seri_num)                       ;站点序号的个数
year=ispan(1951,2018,1)                    ;选择研究时间段
ny=dimsizes(year)                          ;研究的年份的个数
pre=new((/ny,n/),"float")                  ;定义新数组pre
do i=0,n-1
pre(:,i)=f2(:ny-1,seri_num(i))             ;将筛选后的站点降水赋给pre
end do
printVarSummary(pre)
end

sxc123

楼主可以发一下zd.txt文件吗？ 3314274740@qq.com

biguo

来太晚了..当时就是在气候中心的网站上下载的..