提取753站资料出现错误。。跪求大神相助。。

張顧.煒Zane

program swctotal
implicit none
character*8 stid(5)
real lat(5),lon(5),alt,tim
integer,parameter::yrn=49,n=5
integer startyr,startmo,startdy,dyn(12,yrn)
integer i,j,yr,nlev,flag,mo,dy
integer ap,t,t1,t2,h,r,w,s
integer stn
real rain2(17898,5),rain
!137站
character*8::name(n)=(/'59058','59046','59037','59023','59021'/)
tim=0.0
nlev=1
flag=1
do j=1,5
i=1
open(1,file='h:/753rain/753/'//trim(name(j))//'.txt',status='old')
open(2,file='h:/753rain/total.txt')
10 read(1,*,end=100)stid,lat,lon,alt,startyr,startmo,startdy,ap,t,t1,t2,h,r,w,s
do while(startyr>1959)
rain=r
goto 20
enddo
20 if(abs(rain-32700)<0.0001)then
rain=-9990
elseif(abs(rain-32766)<0.0001)then
rain=-9990
elseif(abs(rain-32744)<0.0001)then
rain=-9990
elseif(int(rain/1000.0)==30)then
rain=rain-30000
elseif(int(rain/1000.0)==31)then
rain=rain-31000
elseif(int(rain/1000.0)==32)then
rain=rain-32000
endif
rain2(i,j)=rain
i=i+1
goto 10
100 continue
close(1)
do 40 i=1,17898
rain2(i,j)=rain2(i,j)/10
write(2,*)stid(j),lat(j),lon(j),tim,nlev,flag,rain2(i,j)
print*,stid(j),lat(j),lon(j),tim,nlev,flag,rain2(i,j)
40 continue
nlev=0
write(2,*)stid,lat,lon,tim,nlev,flag
enddo
close(2)
end


程序是这样的。fortran没有报错，能运行。但是提取出的数据有问题。。贴上图。
10 read(1,*,end=100)stid,lat,lon,alt,startyr,startmo,startdy,ap,t,t1,t2,h,r,w,s
do while(startyr>1959)
rain=r
goto 20
enddo
应该是这边的错了，但是不知道该怎么改。我想提取出这几个站的1959年以后的数据。如果没有这部分，试提取一个站的所有年份的是可以的。

wuwei2163

有可能是fortan版本的问题 建议下载一个新版的试试

張顧.煒Zane

wuwei2163 发表于 2013-6-23 01:30

                               
登录/注册后可看大图

有可能是fortan版本的问题 建议下载一个新版的试试

会这样么。。我用的是论坛里下载的版本。。

張顧.煒Zane

fortran是6.6的版本。。

平流层的萝卜

this is resolved ,isn't it?

張顧.煒Zane

平流层的萝卜 发表于 2013-6-23 12:59

                               
登录/注册后可看大图

this is resolved ,isn't it?

你那只能搞一个站啊  搞多个站 就数组溢出了

平流层的萝卜

那再循环下不就可以了吗？

張顧.煒Zane

平流层的萝卜 发表于 2013-6-23 13:03

                               
登录/注册后可看大图

那再循环下不就可以了吗？

goodjob man！

麦田_smile

張顧.煒Zane 发表于 2013-6-23 13:42

                               
登录/注册后可看大图

goodjob man！

这资料很变态啊 每个站点独立一个txt，里头从51到05年 逐日的 我想要挑某年某一天的资料岂不是挖太行山的工作量

Snowyea

受教了，多谢大神！