- 积分
- 40
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2014-10-16
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
- program tiqu
- integer,parameter::m=1,day=22647,yr=63,dd=92*63,yrd=92
- character*5 id(m)
- character*4 yrc
- integer year(m,day),mon(m,day),dy(m,day),temp,tmax,tmin,rain(m,day)
- integer year1(m,dd),mon1(m,dd),dy1(m,dd),rain1(m,dd)
- real rainyr(1951:2012)
- integer a(1951:2012)
- integer i,yrn
- !!!!!!读入站点名!!!!!!
- open(10,file='f:\fortran\shuju\52sta.txt')
- read(10,*) (id(i),i=1,m)
- do i=1,m
- j=1
- !!!!!!读入51-13年数据!!!!!!
- open(33,file='f:\fortran\shuju\51-13\'//id(i)//'.txt')
- 100 read(33,*) id(i),year(i,j),mon(i,j),dy(i,j),rain(i,j)
- print *, id(i),year(i,j),mon(i,j),dy(i,j),rain(i,j)
- j=j+1
- if(.not.eof(33)) goto 100
- write(33,*)
- close(33)
- enddo
- !!!!!!提取3-5月降水!!!!!!
- do i=1,m
- open(44,file='f:\fortran\shuju\51-13\'//id(i)//'.txt')
- do j=1,day
- if(mon(i,j).ge.3.and.mon(i,j).le.5) then
- if(rain(i,j)>=30000) then
- rain(i,j)=0.0
- end if
- if(rain(i,j)==-32744.or.rain(i,j)==32744.or.rain(i,j)==32766.or.rain(i,j)==-32766) then
- rain(i,j)=0.0
- end if
- write(44,*) id(i),year(i,j),mon(i,j),dy(i,j),rain(i,j) !按台站存放!
- endif
- enddo
- close(44)
- enddo
- end
复制代码 运行出来提取的是气温,要怎么改才能提取降水
文件格式如下
51053 1974 7 25 220 301 145 0
51053 1974 7 26 234 296 179 0
51053 1974 7 27 246 325 139 0
51053 1974 7 28 266 343 166 0
51053 1974 7 29 284 359 168 0
51053 1974 7 30 248 352 194 32700
51053 1974 7 31 228 305 141 13
51053 1974 8 1 249 327 166 0
51053 1974 8 2 231 298 171 32700
51053 1974 8 3 243 329 158 32700
51053 1974 8 4 243 324 155 0
51053 1974 8 5 246 332 153 32700
51053 1974 8 6 261 332 185 0
51053 1974 8 7 254 321 199 32700
51053 1974 8 8 215 278 176 2
51053 1974 8 9 207 282 140 8
51053 1974 8 10 224 297 143 32700
51053 1974 8 11 213 281 144 0
51053 1974 8 12 195 268 114 0
求大神解答!!!
|
|