FORTRAN批处理

1368984969

大家好，新年快乐。在程序的调试过程中出现了一下问题，已经尝试了好长时间还是不行，所以请求大神可以帮帮我。
程序是用来处理全球分辨率为1°的AOD栅格数据，是想让数据最后按照经度，纬度，value三列来输出，之前的数据中没有经纬度。在处理单个文件时没有问题，是在后面加上批处理的程序才出现问题。
问题：编译、链接都没有出现错误，但在最后的运行中出现问题。想请教一下是怎么回事？麻烦了。

下面是我的代码

1. PROGRAM aod_process
   
2. IMPLICIT NONE
   
3. INTEGER,PARAMETER::row=180,col=360,n=11
   
4. INTEGER::i,j,k
   
5. REAL,DIMENSION(row)::lat
   
6. REAL,DIMENSION(col)::lon
   
7. INTEGER,DIMENSION(row,col)::aod
   
8. CHARACTER(len=605)a(n)
   
9. CHARACTER(len=55)::filename!每个文件名有55个字符，55*11=605
   
10. OPEN(20,file='D:\AOD\Terra_AOD_C6\2000\T.txt')!将filename.txt文档中的文件名读到数组a(n)中
    
11. DO i=1,n
    
12. READ(20,*)a(i)
    
13. PRINT*,a(i)
    
14. END DO
    
15. Close(20)
    
16. DO k=1,n
    
17. WRITE(filename,'(a55)')a(k)!利用循环将a(k)赋给变量filename
    
18. OPEN(15,file='D:\AOD\Terra_AOD_C6\2000\'//filename//'',status='old',form='formatted')
    
19. read(15,*) ((aod(i,j),j=1,col),i=1,row)
    
20. OPEN(10,position='Append',file='D:\AOD\Terra_AOD_C6\2000\2000.txt')
    
21. DO i=1,row
    
22.   DO j=1,col
    
23.   lat(i)=89.5-(i-1)
    
24.   lon(j)=0.5+(j-1)
    
25. WRITE(10,1000)lon(j),lat(i),aod(i,j)
    
26. 1000 FORMAT(F7.2,F7.2,I6)
    
27.   END DO
    
28. END DO
    
29. END DO
    
30. END PROGRAM aod_process
    

复制代码

topmad

显示你找不到15号文件
看看循环中打开了15号文件 没有关闭 又继续循环了
加上就好了

1368984969

topmad 发表于 2016-2-21 19:37
显示你找不到15号文件
看看循环中打开了15号文件 没有关闭 又继续循环了
加上就好了

谢谢你，在你的指导下问题解决了，