[流程控制] 关于通过循环来优化程序的。

  if(t==60000.)then
    if(thita>0..and.thita<=0.005)then
    write(11,101)t,Exc,y(4),y(1),y(6),y(3)
    else if(thita>0.005.and.thita<=0.01)then
    write(12,102)t,Exc,y(4),y(1),y(6),y(3)
    else if(thita>0.01.and.thita<=0.015)then
    write(13,103)t,Exc,y(4),y(1),y(6),y(3)
    else if(thita>0.015.and.thita<=0.02)then
    write(14,104)t,Exc,y(4),y(1),y(6),y(3)
    else if(thita>0.02.and.thita<=0.025)then
    write(15,105)t,Exc,y(4),y(1),y(6),y(3)
    else if(thita>0.025.and.thita<=0.03)then
    write(16,106)t,Exc,y(4),y(1),y(6),y(3)

    end if
    end if

    101 format(1X,F15.6,6F15.5)
    102 format(1X,F15.6,6F15.5)
    103 format(1X,F15.6,6F15.5)
    104 format(1X,F15.6,6F15.5)
    105 format(1X,F15.6,6F15.5)
    106 format(1X,F15.6,6F15.5)


    open(unit=11,file='0.005.txt')
    open(unit=12,file='0.01.txt')
    open(unit=13,file='0.015.txt')
    open(unit=14,file='0.02.txt')
    open(unit=15,file='0.025.txt')
    open(unit=16,file='0.03.txt')

	integer  i
	real v 
	character(10) :: filename='0.005.txt'
	v = 0.0
	do i=11, 16
		v = v + 0.005
		write(filename(1:5),'(f5.3)') v
		open(unit=i,file=filename)
	end do
	if(t==60000.0 .and. thita>0 .and. thita<=0.03)then
		i = 11 + int( (thita-tiny(thita)) / 0.005)
		write(i,101) t,Exc,y(4),y(1),y(6),y(3)
	end if
	101 format(1X,F15.6,6F15.5)