fortran处理能见度数据的问题

不能用真名924

处理一个月能见度数据，一分钟一个数据，要求输出三列，一列时间，一列一小时的平均值，一列每小时的最小值，数据有缺测，所以要把数据线赋值为字符型，然后再转为实型，我都乱了，求帮助啊。



program visibility_assic
implicit none
integer,parameter::x=60,y=2000,z
CHARACTER v2(2000,2000)*5,v1(2000)*4
integer,dimension(y)::v_min
real,dimension(y)::v_ave
character(len=10)::header
character(4),dimension(y)::t
integer::i,j
                                    
open(1,file='E:\1\2011\V5803807.011',status='old')
open(2,file='E:\1\out\V5803807.011.txt',status='unknown')
   read(1,*) header
   write(2,*) 'time  v_min  v_ave'
do j=1,y
   read(1,'(a4,60i5)') t(j),v1(1:x,j)
end do
!!!!
do j=1,y
   do i=1,x
  v2=v1
   z=x
   if v1(i:x,j)="/////" then
z=x-1
v1(i,j)=0
v2(i,j)=99999
endif
enddo
!!
   real::v11(1,z),v21(1,x)
   read(v1,*) v11
   read(v2,*) v21
   end
!v11(1:x)="v1(1:x)"
!v21(1:x)="v2(1:x)"
end
!!!!
   v_min(j)=minval(v21(1:x,j))
   v_ave(j)=sum(v11(1:z,j))/z   
enddo
   write(2,'(a4,a2,i5,a2,f8.1)') t(j),'  ',v_min(j),'  ',v_ave(j)
end do

!print*,header
!print*,t(1)
!print*, v1(1:10,1)
!print*,v_min(1)
!print*,v_ave(1)
!stop
close (1)
close (2)

end program visibility_assic

寒江雪（王训

其实除了fortran，你的初始数据有规律的话，用linux的sed，awk，处理就行。

又是那隻貓

不会用Fortran。。

kongfeng0824

楼主的能见度数据是那一套呢？

topmad

数据很乱，导致程序不好懂

大概说一下你的问题呢  有缺测  是要舍去还是插补订正？ 自己要形成一个完整的系统的算法  编程有具体疑问大家来讨论 放上一个数据一个代码，留作业的提问方式不太好

新生活运动

通过read和write将进行字符型变量和其他类型的转换

不能用真名924

topmad 发表于 2013-7-16 20:28

                               
登录/注册后可看大图

数据很乱，导致程序不好懂

大概说一下你的问题呢  有缺测  是要舍去还是插补订正？ 自己要形成一个完整的 ...

嗯，多谢批评指正

ibeing

想问楼主问题解决了吗？正好我也要处理能见度的数据

Candy_Girl

我也在弄这个最近，头疼

yuyuxuan

Candy_Girl 发表于 2017-3-27 20:16
我也在弄这个最近，头疼

楼主能见度数据那里下载的可以告诉一下吗