FORTRAN读MICAPS数据报错 求大神

海soso

小白一枚 借鉴了平流层的萝卜大大读MICAPS数据的帖子http://bbs.06climate.com/forum.p ... 13422&extra=&page=1我想把年份和月份也输出 就报错了 求大神们帮看一下
程序：
parameter(nt=1460)
         integer::n,num
       real,allocatable::lon(:),lat(:),vi(:)
                                                !动态数组分别用来储存站点的经纬度和能见度值
       real:: temp(14)
       character*8,allocatable::sta(:),yr(:,:),mon(:,:)  
                                                 !动态数组，储存站号，站号必须是8个字符，5个的话会stnmap出不出                     
       character*12 filename(nt)        !用于储存1460个时次的文件名

       open(11,file='d:/wu/data1/1994/1994.txt')
        do i=1,nt
         read(11,*) filename(i)
         print*,filename(i)
        enddo
       close(11)               
       open(12,file='d:/wu/vis1994.grd',form='unformatted')
         open(13,file='d:/wu/vis1994.txt')         !储存1460个时次的能见度场
        do k=1,nt
         print*,'文件数',k
           open(11,file='d:/wu/data1/1994/'//filename(k))
         read(11,*)
         read(11,*) n,n,n,n,n        !将该时次的站点数赋值于n
         allocate(lat(n))
         allocate(lon(n))
         allocate(vi(n))
         allocate(sta(n))
         allocate(yr(nt,n))
         allocate(mon(nt,n))
         !print*,n

         close(11)
         open(11,file='d:/wu/data1/1994/'//filename(k)) !获知n后，重新读取。
         read(11,*)
         read(11,*)
          do i=1,n            
                                    !只读能见度，风速暂略     此循环为该时次的站点循环
           read(11,*) sta(i),lon(i),lat(i),(temp(j),j=1,14),vi(i)
           yr(k,i)=filename(k)(1:2)
           mon(k,i)=filename(k)(3:4)

           !print*,vi(i)

            write(12)   sta(i),lat(i),lon(i),yr(k,i),mon(k,i),vi(i)
              write(13,*) sta(i),lat(i),lon(i),yr(k,i),mon(k,i),vi(i)  

          enddo

           close(12)
         close(13)
        !一个时次的能见度场输入完毕
        deallocate(lat)
        deallocate(lon)
        deallocate(vi)
        deallocate(sta)
        close(11)
        enddo
        end

海soso

附上FORTRAN报错提示 应该是动态数据的问题。。。不知道怎么改==

四叶草

你的数组分配后没有释放。   
   allocate(lat(n))
   allocate(lon(n))
   allocate(vi(n))
   allocate(sta(n))
   allocate(yr(nt,n))
   allocate(mon(nt,n))
与下面不对应，下面应再加两个语句
  deallocate(yr)
   deallocate(mon)

海soso

四叶草 发表于 2016-1-10 18:43
你的数组分配后没有释放。   
   allocate(lat(n))
   allocate(lon(n))

谢谢啦 我加上这二句后 发现不报错了 但是只生成了包含一个时次的文件就结束了。。。应该运行1460个时次的呀。。。这是怎么回事呢 是加的位置不对吗
附上加了以后的程序：
  parameter(nt=1460)
         integer::n,num
       real,allocatable::lon(:),lat(:),vi(:)
                                                !动态数组分别用来储存站点的经纬度和能见度值
       real:: temp(14)
       character*8,allocatable::sta(:),yr(:,:),mon(:,:)  
                                                 !动态数组，储存站号，站号必须是8个字符，5个的话会stnmap出不出                     
       character*12 filename(nt)        !用于储存1460个时次的文件名
      
       open(11,file='d:/wu/data1/1994/1994.txt')
        do i=1,nt
         read(11,*) filename(i)
         !print*,filename(i)
        enddo
       close(11)               
       open(12,file='d:/wu/vis1994ym.grd',form='unformatted')
         open(13,file='d:/wu/vis1994ym.txt')         !储存1460个时次的能见度场
        do k=1,nt
         print*,'文件数',k
           open(11,file='d:/wu/data1/1994/'//filename(k))
         read(11,*)
         read(11,*) n,n,n,n,n        !将该时次的站点数赋值于n
         allocate(lat(n))
         allocate(lon(n))
         allocate(vi(n))
         allocate(sta(n))
         allocate(yr(nt,n))
         allocate(mon(nt,n))
         !print*,n
         
         close(11)
         open(11,file='d:/wu/data1/1994/'//filename(k)) !获知n后，重新读取。
         read(11,*)
         read(11,*)
          do i=1,n            
                                    !只读能见度，风速暂略     此循环为该时次的站点循环
           read(11,*) sta(i),lon(i),lat(i),(temp(j),j=1,14),vi(i)
           yr(k,i)=filename(k)(1:2)
           mon(k,i)=filename(k)(3:4)
          if(vi(i)==9999) then !若为缺测站点，直接pass
           continue
            else        !不是缺测站点
             num=1
             endif
      
           !print*,vi(i)
               
           if(num==1) then   
            tim=0.0;nlev=1;nflag=1
            !write(12)   sta(i),lat(i),lon(i),tim,nlev,nflag,vi(i)
              !write(13,*) sta(i),lat(i),lon(i),tim,nlev,nflag,vi(i)
            write(12)   sta(i),lat(i),lon(i),yr(k,i),mon(k,i),vi(i)
              write(13,*) sta(i),lat(i),lon(i),yr(k,i),mon(k,i),vi(i)  
           endif
          enddo
       nlev=0
       write(12) sta(n),lat(n),lon(n),tim,nlev,nflag
           close(12)
         close(13)        !一个时次的能见度场输入完毕
        deallocate(lat)
        deallocate(lon)
        deallocate(vi)
        deallocate(sta)
        deallocate(yr)
        deallocate(mon)
        close(11)
        enddo
        
        end

海soso

四叶草 发表于 2016-1-10 18:43
你的数组分配后没有释放。   
   allocate(lat(n))
   allocate(lon(n))

恩 问题已经都解决了！谢谢~

自导自演

最后怎么解决的呢

double_fight

飘过~~大神快快出现

海soso

自导自演 发表于 2016-1-10 22:41
最后怎么解决的呢

删掉close(12)和close(13)就好了 外层的循环就可以了