fortran文件读取技巧
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
最近专业课海洋要素计算作业需要用fortran编程读取验潮站资料(如下,是txt,其中数据是有62行,限于篇幅,我只截取了前5行数据外加验潮站的4行头信息)
!要注意的是最后一列数据后面没有空格
program ex1
!**********************************************
!* written by 潘海东, 2015,3,20
!*一次性读取744个数据,放入一维数组中
!*(如果一次性读取太多的话可能会对系统造成某些影响)
!**********************************************
INTEGER :: b(62*12),i
character(len=20) time,location,u1,u2
open(10, file = 'C:\Users\Administrator\Desktop\sjs_t_97aug_m(harmonic).txt', status = 'old')
!********读取头信息
read(10,"(A3)") u1
read(10,"(A18)") time
read(10,"(A13)") location
read(10,"(A2)") u2
!********读取数据
read(10,100) b(:)
100 format(61(11(I4,1x),I4,/),11(I4,1x),I4) !在fortran的输入输出中'/'代表换行
close(10)
end
program ex2
!**********************************************
!* written by 潘海东, 2015,3,20
!*读取第3列至第11列
!**********************************************
INTEGER :: a(62,9),i
character(len=20) time,location,u1,u2
open(10, file = 'C:\Users\Administrator\Desktop\sjs_t_97aug_m(harmonic).txt', status = 'old')
!***读取头信息
read(10,"(A3)") u1
read(10,"(A18)") time
read(10,"(A13)") location
read(10,"(A2)") u2
!***读取验潮站数据
do i=1,62
read(10,100) a(i,:)
write(*,99) a(i,:)
end do
99 format(9(I4,1x))
100 format(10x,9(I4,1x),4x)
close(10)
end
!!!!这些技巧总结起来就是利用nX跳过不想要的东西。