Fortran 循环进入文件夹，然后循环处理文件夹内的文件（批处理文件）

 1 program main
 2 implicit none
 3 
 4   character(len = 3) :: cTemp
 5   character(len = 100) :: fileplace
 6   integer, parameter :: time_sequence=27    ! number of time sequences = np
 7   integer, parameter :: np=27    ! corresponding pressure = time_sequence
 8   integer, parameter :: ntap=15    ! number of taps
 9   integer :: tap(ntap)    ! I index of taps (places)
10   integer :: k, i, j, m,n
11   real :: x,y,u,v,w,rho
12   real :: p(np,ntap)
13   
14   data tap /603,643,666,685,700,715,727,740,751,762,773,799,826,857,895/ ! I index of the tap
15 
16 
17   do j = 1, ntap
18 
19     write(*,"('Begin to figure out No.' i2 ' tap!')") j
20 
21     do k = 584, 584+time_sequence-1
22 
23       write(cTemp, '(i3)') k
24       fileplace = '/home/postproc'//trim(adjustl(cTemp))//"/" ！ Name of the folder, not include the name of the target file
25       open(k,file=trim(adjustl(fileplace)) // "results_grid_1_tecplot.dat", form='formatted', status='old')
26 
27       do i = 1, tap(j)+1
28         read(k,*)
29       end do
30 
31       read(k,'(1x,7(1pe15.6))') x,y,u,v,w,rho,p(k-583,j)
32       close(k)
33 
34     end do
35 
36   end do
37 
38   open(10,file='time_pressure_1-7.dat', status='new')
39 
40   do m = 1, time_sequence
41     write(10,'(1x,I3,7(1pe15.6))') m, p(m,1), p(m,2), p(m,3), p(m,4), p(m,5), p(m,6), p(m,7)
42   end do
43 
44   close(10)
45 
46   open(20,file='time_pressure_8-15.dat', status='new')
47 
48   do m = 1, time_sequence
49     write(20,'(1x,I3,8(1pe15.6))') m, p(m,8), p(m,9), p(m,10), p(m,11), p(m,12), p(m,13), p(m,14), p(m,15)
50   end do
51 
52   close(20)
53 
54 end