版块导航: 正在加载中...

登录注册

应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

beefly

专家顾问 (职业作家)

地沟油冶炼专家

专家经验: +457
应助: 86 (初中生)
金币: 12778.4
散金: 27307
红花: 300
帖子: 4664
在线: 925.5小时
虫号: 408372
注册: 2007-06-21
性别: GG
专业: 理论和计算化学
管辖: 计算模拟

[交流] 在Fortran程序中，自定义打印数据列数的三种方法已有7人参与

有时候为了便于数据显示和后期的数据处理，我们希望能控制Fortran程序打印的数据列数，但是write语句不允许在格式控制符中加入变量。其实通过一些技巧，这是可以实现的。

示例代码如下：

CODE:

    IMPLICIT DOUBLE PRECISION (A-H,O-Z)

    dimension a(20),b(20),c(20)

    character*100 frmt

    a=1.d0

    b=2.d0

    c=3.d0

    write(*,*)"Please specify the number of columns (default: 5):"

    read(*,*) ncol

    if(ncol .lt. 1) ncol=5

    write(*,"(/,'a=')")

    do i=1,20

      write(*,"(f20.12)",advance='no')a(i)

      if(mod(i,ncol) .eq. 0) write(*,*)

    end do

    write(*,*)

    frmt='(    f11.3)'

    write(frmt(2:5),"(i4)")ncol

    write(*,"(/,'b=')")

    write(*,trim(frmt))b

    write(*,*)

    write(*,*)"Please specify the format (for example, 5f12.6):"

    read(*,*)frmt

    frmt="(" // trim(frmt) // ")"

    write(*,"(/,'c=')")

    write(*,trim(frmt))c

    end

首先，由用户定义输出的列数：

CODE:

    Please specify the number of columns (default: 5):

     3  用户输入

例1，打印数组a。这里利用了Fortran90引入的一个wreite选项advance='no'，表示打印之后不换行。注意此选项对Fortran 77编译器（例如pgf77，g77）无效。每次只打印一个数组的元素，当数据编号能被列数整除时，执行一个空的write命令，起到换行作用。输出为：

CODE:

    a=

          1.000000000000      1.000000000000      1.000000000000

          1.000000000000      1.000000000000      1.000000000000

          1.000000000000      1.000000000000      1.000000000000

          1.000000000000      1.000000000000      1.000000000000

          1.000000000000      1.000000000000      1.000000000000

          1.000000000000      1.000000000000      1.000000000000

          1.000000000000      1.000000000000

例2，打印数组b。这里利用字符串frmt存储输出格式，并利用write语句把自定义的整型列数作为字符打印到frmt，把frmt预留的第2到第5个字符位置改为用户定义的列数。在上面的例子中，frmt改后成为( 3f11.3)。我们一般只输出几列到几十列，因此预留4个字符位置一般足够了。输出为：

CODE:

    b=

          2.000      2.000      2.000

          2.000      2.000      2.000

          2.000      2.000      2.000

          2.000      2.000      2.000

          2.000      2.000      2.000

          2.000      2.000      2.000

          2.000      2.000

例3，打印数组c。这个例子可看作是例2的扩展，输出更灵活一些。不仅是列数，每个元素的输出格式也都由用户按照Fortran的格式定义。例如：

CODE:

    Please specify the format (for example, 5f12.6):

       4d16.6    用户输入

然后根据Fortran格式，还要把前后的括号加上。输出如下：

CODE:

    c=

        0.300000D+01    0.300000D+01    0.300000D+01    0.300000D+01

        0.300000D+01    0.300000D+01    0.300000D+01    0.300000D+01

        0.300000D+01    0.300000D+01    0.300000D+01    0.300000D+01

        0.300000D+01    0.300000D+01    0.300000D+01    0.300000D+01

        0.300000D+01    0.300000D+01    0.300000D+01    0.300000D+01

用户自定义的格式未必总是符合语法的。如果作为正式的程序，最好给write加上一个出错判断：若给定格式的输出无法执行，则按照默认格式输出。

回复此楼

» 本帖已获得的红花（最新10朵）

小范范1989

xk6891

» 猜你喜欢

职称评审没过，求安慰已经有51人回复
毕业后当辅导员了，天天各种学生超烦已经有5人回复
26申博自荐已经有3人回复
A期刊撤稿已经有4人回复
垃圾破二本职称评审标准已经有17人回复
投稿Elsevier的Neoplasia杂志，到最后选publishing options时页面空白，不能完成投稿已经有22人回复
EST投稿状态问题已经有7人回复

» 本主题相关价值贴推荐，对您同样有帮助:

急：求用Fortran或者其他语言编写一个程序，处理下面的数据！！！！！！！！！！已经有37人回复
Fortran怎么从一个dat文件中读取不同类型数据赋值给不同的变量？已经有8人回复
Fortran如何读取不确定行数的数据已经有3人回复
高分求将我的MATLAB程序改写为FORTRAN程序（1000BB）已经有13人回复
Fortran输出二进制数据，用matlab读出，结果出错，求助问题何在？已经有4人回复
fortran编译中$符号的使用已经有4人回复
fortran程序在SSH计算中心上运行的问题！已经有13人回复
请教一下CFX如何导入FORTRAN用户自定义程序已经有5人回复
关于复数在Fortran中如何定义的问题已经有3人回复
用Fortran编写程序输出的数据一般用什么程序作图比较好？已经有28人回复
fortran，带字符的数据读取与存储！已经有11人回复
一个主程序，一个subroutine程序，要如何将两个程序组合在fortran里运行已经有4人回复
fortran中如何处理字符和数据混合的文件已经有6人回复
Fortran数据读入格式的修改已经有5人回复
请教 fortran 运行错误的原因已经有13人回复
关于fortran中的rewind 语言已经有5人回复
Fortran中子程序与函数的区别已经有7人回复
如何改变fortran的浮点精度已经有28人回复
请教一个fortran小程序编译出错的问题，谢谢已经有9人回复
写了一个fortran90的小程序，编译通不过，请大侠帮忙已经有59人回复
【求助】fortran中是如何实现将两个变量的值在同一行进行打印出来的已经有7人回复
【求助】fortran中怎么存储不确定长度的数据【已完结】已经有6人回复

beefly《西太平洋大学现代英汉词典》[bi:fli]牛肉一般地

1楼 2015-02-01 04:11:34

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

beefly

专家顾问 (职业作家)

地沟油冶炼专家

专家经验: +457
应助: 86 (初中生)
金币: 12778.4
散金: 27307
红花: 300
帖子: 4664
在线: 925.5小时
虫号: 408372
注册: 2007-06-21
性别: GG
专业: 理论和计算化学
管辖: 计算模拟

引用回帖:

7楼: Originally posted by pippi6 at 2015-02-02 22:40:40
“但是write语句不允许在格式控制符中加入变量” lz 的这个论断是不正确的。比如 n 是一个正整数，我希望打印 n 用合适的整数位，就可以写
write(6,"(a,I<int(log10(n+0d0))+1>