- 积分
- 17
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2012-10-17
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
处理一个月能见度数据,一分钟一个数据,要求输出三列,一列时间,一列一小时的平均值,一列每小时的最小值,数据有缺测,所以要把数据线赋值为字符型,然后再转为实型,我都乱了,求帮助啊。
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
|
-
-
V5803807.011
444.96 KB, 下载次数: 17, 下载积分: 金钱 -5
2011年7月资料
-
-
read.f90
1.06 KB, 下载次数: 3, 下载积分: 金钱 -5
程序,好多错
|