- 积分
- 1502
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2012-4-12
- 最后登录
- 1970-1-1
|
Fortran
系统平台: |
fortran |
问题概况: |
编了个小程序,想把特定时间段的数据提出来,调试成功后,运行也没有问题。但是没有生成我想要的文件 |
问题截图: |
- |
我看过提问的智慧: |
看过 |
自己思考时长(天): |
2 |
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
本帖最后由 黄小仙儿 于 2014-7-5 17:51 编辑
编了个小程序,想把特定时间段的数据提出来,调试成功后,运行也没有问题。但是没有生成我想要的文件,即文件time.txt
这是我的程序:
我是想把20:00,20:20,20:40,21:00,21:20,21:40......这几个时间点,前后5分钟发生闪电的经纬度数提出来(比如20点的,我需要19点57分30秒——20点02分30秒之间的闪电经纬度),写入新的的文件中
program incorpfile
integer,parameter::n=11644
integer :: hour(n),mini(n),sec(n)
integer i,j
logical alive
real(8) :: lat(n),lon(n),m
open(10,file="D:/Fortran/Console11/25yuan.txt",status='old')
do i=1,n
read(10,*)hour(i),mini(i),sec(i),lat(i),lon(i)
!print*,i
enddo
do i=1,n
m=hour(i)*3600+mini(i)*60+sec(i)
if((hour(i)*3600+mini(i)*60+sec(i))>=(70650.+i*1200.).and.(hour(i)*3600+mini(i)*60+sec(i))<=(70950.+i*1200.)) then
open(12,file="D:/Fortran/Console11/time.txt",form="formatted")
write(12,100)m,lat(i),lon(i)
endif
100 format(f10.4,f10.4,f10.4)
enddo
close(12)
end
原始文件25yuan.txt的数据我截了一段:第一列是闪电发生的小时,第二列是分钟,第三列是秒数,第四列是纬度,第五列是经度
18 7 33.3271784 31.1279 106.6388
18 7 33.603658 31.1241 106.5709
18 5 9.6678056 31.1581 106.4682
18 5 10.3223859 31.1619 106.4464
18 0 6.8301053 31.7269 106.1537
18 35 30.8234903 31.3522 105.8865
18 43 49.5188434 31.7886 105.758
18 56 55.5062223 31.2948 105.4081
18 0 36.4096751 31.1789 105.3401
18 58 58.8170708 31.6236 105.2973
18 46 42.6784351 31.305 105.2613
18 58 58.9550511 31.447 105.2417
求大神帮忙~~
|
|