- 积分
- 241
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2018-5-13
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
本帖最后由 biguo 于 2020-2-14 15:32 编辑
记录:筛选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
|
|