- 积分
- 1084
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 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
|
|