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

涛声依旧了

赞一个，以前的台风数据比较齐全，登陆地段，路径趋势都有的，现在变成这种最佳台风路径数据集后，就简单多了，唉。。。。

kongfeng0824

!========================批处理读取文件==========================  
  do k=1949,2009
    write(cTemp,'(I4)') k
filename='CH'//cTemp//'BST.txt'  !或filename='CH'//trim(adjustl(cTemp))//'BST.txt'
    write(*,*) "正在处理文件:",filename,"..."
    do i=1,12
      n(i)=0
    end do
    NUM=1

请教兰溪版主，这段中的filename那一行的意思。最近我也在批量处理降水的数据，我把目录下的文件名都读进了一个txt文件中，我应该如何打开这些数据处理呢？请指点一二。

兰溪之水

kongfeng0824 发表于 2013-2-27 23:33

                               
登录/注册后可看大图

!========================批处理读取文件==========================  
  do k=1949,2009
    write(cTe ...

看看这个http://bbs.06climate.com/forum.php?mod=viewthread&tid=3079
写得很清楚了~

liuchao0120

一定支持原创，兰溪加油，论坛加油啦

胆小鬼yy

传说中的谁 发表于 2011-9-29 20:42
抢沙发之后慢慢看继续编辑~~~~~~~看起来像是处理台风路径的？兰溪辛苦啦，果断收藏，以后肯定用得 ...

感谢分享

Figo

版主，我编译了你的程序，运行后提示数组重复定义。请看一下，谢谢！

兰溪之水

Figo 发表于 2015-7-8 00:15
版主，我编译了你的程序，运行后提示数组重复定义。请看一下，谢谢！

可以在end do前面把动态数组释放

学而无术

一定支持原创，兰溪加油，论坛加油啦

葱儿

真是帮了我大忙，谢谢

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



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