fortran多层次文件夹嵌套数据 【批处理】【文件】【多层文件夹】【系统命令】grads

3415519943

1. 
   

复制代码

在处理数据时经常会碰到文件夹嵌套文件夹的情况，如一个大文件夹包含a个高度文件夹，每个高度文件夹包含b个变量文件夹，每个变量文件夹包含c个日期文件，而日期文件中才包含真正的数据，共计a*b*c个文件，要对这些文件读取写出完整的命令太麻烦，所以调用系统命令读取文件夹名称和文件名。故分享实例代码如下，本代码读取嵌套文件数据并转换成grd格式，无法涵盖所有情况，仅供参考：

1. program read_files
   
2.       USE DFLIB
   
3.   implicit none
   
4.   character*100 folder_path,output,ii, file_name
   
5.   character*20 varstr,hstr(4),tstr
   
6.   integer varnum,hnum(4),tnum,maxhnum,locmax(1)
   
7.   integer vari,hi,ti
   
8.   integer i,numfiles,filenum,ios
   
9.   logical(4) res
   
10.   
    
11.   open(100,file='D:\grads\dataofzhenduan.grd',form='binary')
    
12.   output='none.txt'
    
13.   call outnone(output)!!
    
14.   
    
15.   folder_path = 'D:\grads\zhenduanshuju\data'
    
16.   i=1
    
17.   write(ii,'(i2)'),i
    
18.   output='dir'//trim(adjustl(ii))//'.txt'
    
19.   call ListFile(folder_path,outPut)
    
20.   varnum=numfiles(output)
    
21.   open(31,file=''//trim(output)//'')
    
22.   
    
23.   do i=1,varnum
    
24.       read(31,*),varstr
    
25.       folder_path = 'D:\grads\zhenduanshuju\data\'//trim(varstr)
    
26.       write(ii,'(i2)'),10+i
    
27.       output='dir'//trim(adjustl(ii))//'.txt'
    
28.       call ListFile_n(folder_path,outPut)
    
29.       hnum(i)=numfiles(output)
    
30.       open(i,file=''//trim(output)//'')
    
31.   enddo
    
32.   
    
33.   rewind(31)
    
34.   maxhnum=maxval(hnum)
    
35.   locmax=maxloc(hnum)
    
36.   i=locmax(1)!!
    
37.   close(i)!!
    
38.   !print*, locmax
    
39.   write(ii,'(i2)'),10+locmax
    
40.   output='dir'//trim(adjustl(ii))//'.txt'
    
41.   open(32,file=''//trim(output)//'')
    
42.   
    
43.       
    
44.   
    
45.   folder_path = 'D:\grads\zhenduanshuju\data\air\100'
    
46.   i=3
    
47.   write(ii,'(i2)'),i
    
48.   output='dir'//trim(adjustl(ii))//'.txt'
    
49.   call ListFile(folder_path,outPut)
    
50.   tnum=numfiles(output)
    
51.   open(33,file=''//trim(output)//'')
    
52.   !print*,varnum,hnum,tnum,maxhnum
    
53.   
    
54.   do ti=1,tnum
    
55.       read(33,*),tstr
    
56.       !print*,tstr
    
57.       do vari=1,varnum-1
    
58.           read(31,*),varstr
    
59.           do hi=1,maxhnum
    
60.               read(32,*),hstr(4)!!
    
61.               read(vari,*,iostat=ios),hstr(vari)
    
62.               if(hstr(4)/=hstr(vari).OR.ios /= 0)then
    
63.                   BACKSPACE(vari)
    
64.               endif
    
65.               if(hstr(4)/=hstr(vari))then
    
66.                   file_name='none.txt'
    
67.                   call writefile(file_name)
    
68.                   print*,file_name
    
69.                   cycle
    
70.               endif
    
71.               file_name='D:\grads\zhenduanshuju\data\'//trim(varstr)//'\'//trim(hstr(4))//'\'//trim(tstr)
    
72.               call writefile(file_name)
    
73.               print*,file_name
    
74.           enddo
    
75.           rewind(32)
    
76.       enddo
    
77.       file_name='uv'
    
78.       call writefileuv(file_name,tstr)
    
79.       rewind(31)
    
80.       rewind(1)
    
81.       rewind(2)
    
82.       rewind(3)
    
83.   enddo
    
84.   do i=1,varnum
    
85.       close(i)
    
86.   enddo
    
87.   close(31)
    
88.   close(32)
    
89.   close(33)
    
90.   close(100)
    
91.   pause
    
92.     end program read_files
    
93. 
    
94.     subroutine ListFile(fPath,outPut)
    
95.     character*(*),InTent(In):: fPath,outPut
    
96.     character*100 CMD
    
97.     LOGICAL(4) res
    
98.     CMD="dir /on/b "//trim(fPath)//" >"//trim(outPut)
    
99.     res=SYSTEMQQ(CMD)
    
100.     end subroutine
     
101.    
     
102.     subroutine ListFile_n(fPath,outPut)
     
103.     character*(*),InTent(In):: fPath,outPut
     
104.     character*100 CMD
     
105.     LOGICAL(4) res
     
106.     CMD="dir /o-n/b "//trim(fPath)//" >"//trim(outPut)
     
107.     res=SYSTEMQQ(CMD)
     
108.     end subroutine
     
109.    
     
110.     ! Get the number of files in the folder
     
111.     function numfiles(folder)
     
112.     integer numfiles,status
     
113.     character*100 folder
     
114.     numfiles=0
     
115.     open(unit=10, file=''//trim(folder)//'', status='old')
     
116.     do i = 1, 1000
     
117.     read(10, '(A)', iostat=status) file_name
     
118.     if (status /= 0) exit
     
119.     numfiles = numfiles + 1
     
120.     end do
     
121.     close(10)
     
122.     end function
     
123.    
     
124.     subroutine writefile(filename)
     
125.     real filedata(71,101)
     
126.     integer::i=1,aa
     
127.     character*100 filename
     
128.     open(90,file=''//trim(filename)//'')
     
129.     !if(aa==4)then
     
130.      !   read(90,*)
     
131.       !  read(90,*)
     
132.        ! read(90,*)
     
133.         !do i=1,71
     
134.          !   read(90,*),(filedata(i,j),j=1,101)
     
135.           !  write(100),(filedata(i,j),j=1,101)
     
136.         !enddo
     
137.         !i=1
     
138.         !do i=1,71
     
139.         !    read(90,*),(filedata(i,j),j=1,101)
     
140.         !    write(100),(filedata(i,j),j=1,101)
     
141.         !enddo
     
142.     !else
     
143.         read(90,*)
     
144.         read(90,*)
     
145.         read(90,*)
     
146.         read(90,*)
     
147.         do i=1,71
     
148.             read(90,*),(filedata(i,j),j=1,101)
     
149.             write(100),(filedata(i,j),j=1,101)
     
150.         enddo
     
151.     !endif
     
152.     close(90)
     
153.     end
     
154.    
     
155.     subroutine outnone(outPut)
     
156.     real::none(71,101)=9999.990
     
157.     integer i
     
158.     character*100 output
     
159.     open(unit=11, file=''//trim(output)//'')
     
160.     write(11,*),9999.990
     
161.     write(11,*),9999.990
     
162.     write(11,*),9999.990
     
163.     write(11,*),9999.990
     
164.     do i=1,71
     
165.     write(11,*),(none(i,j),j=1,101)
     
166.     enddo
     
167.     close(11)
     
168.     end subroutine
     
169.    
     
170.     subroutine writefileuv(varstr,tstr)
     
171.     real filedata(142,101,12)
     
172.     integer::i,j,k
     
173.     character*100 varstr,filefile
     
174.     character*20 tstr,hstr
     
175.     do k=1,12
     
176.         read(32,*),hstr
     
177.         filefile='D:\grads\zhenduanshuju\data\'//trim(varstr)//'\'//trim(hstr)//'\'//trim(tstr)
     
178.         open(90,file=''//trim(filefile)//'')
     
179.         read(90,*)
     
180.         read(90,*)
     
181.         read(90,*)
     
182.         do i=1,142
     
183.             read(90,*),(filedata(i,j,k),j=1,101)
     
184.         enddo
     
185.     enddo
     
186.     rewind(32)
     
187.     do k=1,12
     
188.         do i=1,71
     
189.             write(100),(filedata(i,j,k),j=1,101)
     
190.         enddo
     
191.         print*,k
     
192.     enddo
     
193.     do k=1,12
     
194.         do i=72,142
     
195.             write(100),(filedata(i,j,k),j=1,101)
     
196.         enddo
     
197.         print*,k
     
198.     enddo
     
199.     close(90)
     
200.     end subroutine

复制代码

以下有逐行注释的代码，本代码泛用性不强，其实可以通过宏定义提高泛用性，等以后有时间了再来修改吧，如果想要了解思路，可以下载代码文件
Source1.f90 (7.15 KB, 下载次数: 0, 售价: 1 贡献)

2024-4-3 16:50 上传

点击文件名下载附件
逐行注释


数据说明.txt (1.19 KB, 下载次数: 1)

2024-4-3 16:56 上传

点击文件名下载附件

data.zip (38.98 MB, 下载次数: 50)

2024-4-3 16:58 上传

点击文件名下载附件