- 积分
 - 100
 
	- 贡献
 -  
 
	- 精华
 
	- 在线时间
 -  小时
 
	- 注册时间
 - 2016-9-29
 
	- 最后登录
 - 1970-1-1
 
 
 
 
 
 
 | 
	
 
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册 
 
 
 
x
 
program main 
implicit none 
integer,parameter::dy=365,nfile=16 
real,parameter::a=28.5,b=23.5,c=33.5 
real,dimension(dy,nfile)::var 
real,dimension(dy)::u,v 
integer x,y,l 
character*200 filename 
open(11,file='D:\grads\Year\1987\input.txt') 
do x=1,nfile,1 
read(11,*) filename 
 open(12,file='D:\grads\Year\1987\1987P01.dat',form='binary') 
 read(12,*) ((var(y,x),y = 1,dy),x=1,nfile)  
! print *,var  
enddo 
close(11) 
close(12) 
 
open(3,file='D:\grads\TXT\1989\u.txt',action='write',form='formatted',status='new') 
open(4,file='D:\grads\TXT\1989\v.txt',action='write',form='formatted',status='new')   
do y=1,nfile 
    do x=1,dy 
       if (var(x,y) == -9.96921e+36)then 
          l=x 
         var(l,y)=0 
     endif 
    enddo        
 enddo 
 
do x=1,dy 
       u(x) = 0.5*(var(x,12) + var(x,13) - var(x,4) - var(x,5)) 
       v(x) = 0.25/cos(a) * (var(x,5) + 2*var(x,9) + var(x,13) - var(x,4) - 2*var(x,8) - var(x,12)) 
  write(3,*) u(x) 
  write(4,*) v(x) 
enddo 
close(3) 
close(4) 
end program main 
 
各位叔叔阿姨哥哥姐姐和小伙伴们,上述是本人在批量处理文件时的程序,input.txt内为所有data数据的路径+文件名,一共16行,即每一年有16个data数据。读取过程中出现了一个问题,如下图 
。 
注:data文件转化成txt文件后为一列real型数据,无特殊字符和描述文件说明出现,写了ctl配合用nc也可以直接画出图。 
想了好几天,问了同学也没想明白问题出在哪里?拜托各位大神帮帮忙,小生感激不尽。 
 |   
- 
提示出现的问题 
 
 
 
 
 
 
 
 |