分享一个写得很烂的Fortran批处理CMA_TCBestTrack数据

yesyesnotion

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,"(A11)",iostat=stat1) filename
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



以上是我的程序，为什么不执行命令呢，什么结果都没有出

葱儿

太好了，一直都想要这样的程序！！！谢谢！！！

zdcgl

感谢楼主的无私分享

liuxunx

版主你好，在使用过程中出现了之前有人提到过的数组重复定义的问题，根据你的回复在enddo之前释放数组的话，会出现数组没有声明的错误