【已解决】如何给这个程序加个循环

倔强的葱头

topmad 发表于 2013-4-25 09:30

                               
登录/注册后可看大图

你这样  每次组合完了文件名以后  print一下 然后pause  看看对不对 在进行下面的 一步步print 就知道错在 ...

嗯 好主意谢谢

倔强的葱头

topmad 发表于 2013-4-24 16:32

                               
登录/注册后可看大图

葱头，你的失败代码是啥
个人感觉，一般这里错误，都是在文件名的地方

real dt(360,640)
character(len=20)nouse,form,name(12),outname(12)
do iyear=0,1
do i=1,12
write(form,"(i1)") lm
write(name(lm),*) lm= iyear*100+i +19600,".txt"
write(outname(lm),*) lm= iyear*100+i +19600,".xls"
enddo
enddo

open(1,file=name(lm))
do i=1,6
read(1,*)nouse
enddo
do i=1,360
read(1,*)(dt(i,j),j=1,640)
enddo
close(1)

open(2,file=outname(lm))
do i=1,360
do j=1,640
if(dt(i,j).ne.-9999.0) write(2,"(3f8.1)") 72+(j-1)*0.1,18+360*0.1-i*0.1,dt(i,j)
enddo
enddo
close(2)

enddo
end
这是我失败的。。。

纳川

文件名有问题。write(name(lm),*) lm= iyear*100+i +196000,".txt"和write(outname(lm),*) lm= iyear*100+i +196000,".xls"中 lm= iyear*100+i +196000、 lm= iyear*100+i +196000为整型变量，不能直接赋值给字符型变量name(lm)，应先转换成字符型变量。

纳川

改成如下试试：
do iyear=0,100
do i=1,12
lm= iyear*100+i +196000
write(form,"(i1)") lm
write(name(lm),*)trim(form),".txt"
write(outname(lm),*) trim(form),".xls"

倔强的葱头

纳川 发表于 2013-4-26 11:26

                               
登录/注册后可看大图

改成如下试试：
do iyear=0,100
do i=1,12

forrtl: severe (161): Program Exception - array bounds exceeded
Image              PC        Routine            Line        Source
fpiliang.exe       004010E2  Unknown               Unknown  Unknown

Incrementally linked image--PC correlation disabled.


这个错误是什么意思

倔强的葱头

纳川 发表于 2013-4-26 11:26

                               
登录/注册后可看大图

改成如下试试：
do iyear=0,100
do i=1,12

问题已解决。。。是文件的问题

倔强的葱头

纳川 发表于 2013-4-26 11:26

                               
登录/注册后可看大图

改成如下试试：
do iyear=0,100
do i=1,12

real dt(72,128)
character(len=20)nouse,form,name(600),outname(600)

icount=0
do iyear=1,50 !年
do lm=1,12    !月
   
    ifile= 196000+ iyear*100+ lm
    icount = icount +1

   write(form,"(i10)") ifile
   write(name(icount),*) trim(form),".txt"
   write(outname(icount),*) trim(form),".xls"
print*,"正在处理",trim(form)

open(1,file=name(icount))
do i=1,6
read(1,*)nouse
enddo
do i=1,72
read(1,*)(dt(i,j),j=1,128)
enddo
close(1)


open(2,file=outname(icount))
do i=1,72
do j=1,128
if(dt(i,j).ne.-9999.0) write(2,"(3f8.1)") 72+(j-1)*0.5,18+72*0.5-i*0.5,dt(i,j)
enddo
enddo
close(2)

enddo
enddo  !年
end