- 积分
 - 1086
 
	- 贡献
 -  
 
	- 精华
 
	- 在线时间
 -  小时
 
	- 注册时间
 - 2011-6-22
 
	- 最后登录
 - 1970-1-1
 
 
 
 
 
 
 | 
	
 
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册 
 
 
 
x
 
下面是我的程序 
请问读取批量文件数据出了什么问题,加入这个批量读取后,程序就没有执行,不知道问题出在哪里了,也不理解批量处理的每一行是什么意思,请各位高手帮忙解决一下,这个批量的处理程序也是我从网上搜到的,稀里糊涂的贴进去发现行不通。 
 
渴求答案,希望各位能抽一分钟帮我解答。 
 
program l_bound2015 
implicit none 
real::date4=650.0 
integer,parameter::N=5282,nfile=100 
real*8 upar_temp(N,12) 
integer x,y,l 
real Hr,hra,hrn 
 
!!!!!!!!!!!!!批量读取文件!!!!!!!!!!!!!!! 
integer ::stat1 
character(11) ::filename 
 
open(13,file="input.list")  !The name of input file 
open(14,file="output.list",status="replace") !The name of output file 
 
do  
read(13,iostat=stat1) filename 
 
if(stat1/=0) exit 
open(15,file=filename) 
 
print*,2015 
read(15,*) 
read(15,*) 
read(15,*) 
read(15,*) 
read(15,*) 
read(15,*) 
read (15,*) ((upar_temp(x,y),y=1,12),x=7,N) 
 
do x=7,N 
 
!1 
if(upar_temp(x,12).lt.1000)then 
 
if(upar_temp(x,4).GE.86.and.upar_temp(x,4).le.94)then 
write(14,"(1x,F8.4,1x,F9.3,1x,F7.3)") upar_temp(x,4),upar_temp(x,7)+upar_temp(x,12) 
endif 
 
 
print*,"1" 
 
!2 
elseif (upar_temp(x,12).ge.1000.and.upar_temp(x,12).lt.2000)then 
Hr=-6.416*(upar_temp(x,12)/1000.0)+97.0 
print*,Hr 
hra=Hr+3.0 
hrn=Hr-5.0 
print*,'000' 
 
if(upar_temp(x,4).GE.hrn.and.upar_temp(x,4).le.hra)then 
print*,'111' 
write(14,"(1x,F8.4,1x,F9.3,1x,F7.3)") upar_temp(x,4),upar_temp(x,7)+upar_temp(x,12) 
endif 
print*,"2" 
 
 
!3 
elseif(upar_temp(x,12).ge.2000.and.upar_temp(x,12).lt.7562)then 
Hr=-1.223*(upar_temp(x,12)/1000)+87 
hra=Hr+3.0 
hrn=Hr-5.0 
 
if(upar_temp(x,4).GE.hrn.and.upar_temp(x,4).le.hra)then 
write(14,"(1x,F8.4,1x,F9.3,1x,F7.3)") upar_temp(x,4),upar_temp(x,7)+upar_temp(x,12) 
endif 
print*,"3" 
 
 
!4 
elseif(upar_temp(x,12).ge.7562.and.upar_temp(x,12).le.10000)then 
Hr=-4.0*(upar_temp(x,12)/1000)+108 
hra=Hr+3.0 
hrn=Hr-5.0 
 
if(upar_temp(x,4).GE.hrn.and.upar_temp(x,4).le.hra)then 
write(14,"(1x,F8.4,1x,F9.3,1x,F7.3)") upar_temp(x,4),upar_temp(x,7)+upar_temp(x,12) 
endif 
print*,"4" 
 
!5 
else 
if(upar_temp(x,4).GE.63.and.upar_temp(x,4).le.71)then 
write(14,"(1x,F8.4,1x,F9.3,1x,F7.3)") upar_temp(x,4),upar_temp(x,7)+upar_temp(x,12) 
endif 
print*,"5" 
 
endif 
enddo 
 
 
 
close(15) 
enddo 
end 
 
 |   
 
 
 
 |