當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖

beefly

禁蟲

地溝油冶煉專家

專家經(jīng)驗(yàn): +458
應(yīng)助: 88 (初中生)
金幣: 13026.2
散金: 27307
紅花: 300
帖子: 4673
在線: 927.1小時
蟲號: 408372
注冊: 2007-06-21
性別: GG
專業(yè): 理論和計算化學(xué)
管轄: 計算模擬

[交流] 在Fortran程序中，自定義打印數(shù)據(jù)列數(shù)的三種方法已有7人參與

有時候?yàn)榱吮阌跀?shù)據(jù)顯示和后期的數(shù)據(jù)處理，我們希望能控制Fortran程序打印的數(shù)據(jù)列數(shù)，但是write語句不允許在格式控制符中加入變量。其實(shí)通過一些技巧，這是可以實(shí)現(xiàn)的。

示例代碼如下：

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

首先，由用戶定義輸出的列數(shù)：

CODE:

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

     3  用戶輸入

例1，打印數(shù)組a。這里利用了Fortran90引入的一個wreite選項(xiàng)advance='no'，表示打印之后不換行。注意此選項(xiàng)對Fortran 77編譯器（例如pgf77，g77）無效。每次只打印一個數(shù)組的元素，當(dāng)數(shù)據(jù)編號能被列數(shù)整除時，執(zhí)行一個空的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，打印數(shù)組b。這里利用字符串frmt存儲輸出格式，并利用write語句把自定義的整型列數(shù)作為字符打印到frmt，把frmt預(yù)留的第2到第5個字符位置改為用戶定義的列數(shù)。在上面的例子中，frmt改后成為( 3f11.3)。我們一般只輸出幾列到幾十列，因此預(yù)留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，打印數(shù)組c。這個例子可看作是例2的擴(kuò)展，輸出更靈活一些。不僅是列數(shù)，每個元素的輸出格式也都由用戶按照Fortran的格式定義。例如：

CODE:

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

       4d16.6    用戶輸入

然后根據(jù)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加上一個出錯判斷：若給定格式的輸出無法執(zhí)行，則按照默認(rèn)格式輸出。

回復(fù)此樓

» 本帖已獲得的紅花（最新10朵）

小范范1989

xk6891

» 猜你喜歡

271求調(diào)劑已經(jīng)有12人回復(fù)
322求調(diào)劑已經(jīng)有3人回復(fù)
環(huán)境工程 085701，267求調(diào)劑已經(jīng)有15人回復(fù)
材料0856 英一數(shù)二 323 求調(diào)劑已經(jīng)有3人回復(fù)
301求調(diào)劑已經(jīng)有3人回復(fù)
348求調(diào)劑已經(jīng)有3人回復(fù)
350求調(diào)劑已經(jīng)有3人回復(fù)
309分085801求調(diào)劑已經(jīng)有7人回復(fù)
生物學(xué)308分求調(diào)劑（一志愿華東師大）已經(jīng)有3人回復(fù)
求調(diào)劑：一志愿：南京大學(xué) 專業(yè)：0705 總分320 ，本科985，四六級已過已經(jīng)有3人回復(fù)

» 本主題相關(guān)價值貼推薦，對您同樣有幫助:

急：求用Fortran或者其他語言編寫一個程序，處理下面的數(shù)據(jù)�。。。。。。。。�！已經(jīng)有37人回復(fù)
Fortran怎么從一個dat文件中讀取不同類型數(shù)據(jù)賦值給不同的變量？已經(jīng)有8人回復(fù)
Fortran如何讀取不確定行數(shù)的數(shù)據(jù) 已經(jīng)有3人回復(fù)
高分求將我的MATLAB程序改寫為FORTRAN程序（1000BB）已經(jīng)有13人回復(fù)
Fortran輸出二進(jìn)制數(shù)據(jù)，用matlab讀出，結(jié)果出錯，求助問題何在？已經(jīng)有4人回復(fù)
fortran編譯中$符號的使用已經(jīng)有4人回復(fù)
fortran程序在SSH計算中心上運(yùn)行的問題！已經(jīng)有13人回復(fù)
請教一下CFX如何導(dǎo)入FORTRAN用戶自定義程序已經(jīng)有5人回復(fù)
關(guān)于復(fù)數(shù)在Fortran中如何定義的問題已經(jīng)有3人回復(fù)
用Fortran編寫程序輸出的數(shù)據(jù)一般用什么程序作圖比較好？已經(jīng)有28人回復(fù)
fortran，帶字符的數(shù)據(jù)讀取與存儲！已經(jīng)有11人回復(fù)
一個主程序，一個subroutine程序，要如何將兩個程序組合在fortran里運(yùn)行已經(jīng)有4人回復(fù)
fortran中如何處理字符和數(shù)據(jù)混合的文件已經(jīng)有6人回復(fù)
Fortran數(shù)據(jù)讀入格式的修改已經(jīng)有5人回復(fù)
請教 fortran 運(yùn)行錯誤的原因已經(jīng)有13人回復(fù)
關(guān)于fortran中的rewind 語言已經(jīng)有5人回復(fù)
Fortran中子程序與函數(shù)的區(qū)別已經(jīng)有7人回復(fù)
如何改變fortran的浮點(diǎn)精度已經(jīng)有28人回復(fù)
請教一個fortran小程序編譯出錯的問題，謝謝已經(jīng)有9人回復(fù)
寫了一個fortran90的小程序，編譯通不過，請大俠幫忙已經(jīng)有59人回復(fù)
【求助】fortran中是如何實(shí)現(xiàn)將兩個變量的值在同一行進(jìn)行打印出來的已經(jīng)有7人回復(fù)
【求助】fortran中怎么存儲不確定長度的數(shù)據(jù)【已完結(jié)】已經(jīng)有6人回復(fù)

beefly《西太平洋大學(xué)現(xiàn)代英漢詞典》[bi:fli]牛肉一般地

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

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

beefly

禁蟲

地溝油冶煉專家

專家經(jīng)驗(yàn): +458
應(yīng)助: 88 (初中生)
金幣: 13026.2
散金: 27307
紅花: 300
帖子: 4673
在線: 927.1小時
蟲號: 408372
注冊: 2007-06-21
性別: GG
專業(yè): 理論和計算化學(xué)
管轄: 計算模擬

引用回帖:

7樓: Originally posted by pippi6 at 2015-02-02 22:40:40
“但是write語句不允許在格式控制符中加入變量” lz 的這個論斷是不正確的。比如 n 是一個正整數(shù)，我希望打印 n 用合適的整數(shù)位，就可以寫
write(6,"(a,I<int(log10(n+0d0))+1>

' n=',n
...

謝謝！我以前試過括號，引號，等等，都不行。10年前在smth的計算版中問過，所有人都告訴我是不可能的，讓我斷了此念...

贊一下

回復(fù)此樓

beefly《西太平洋大學(xué)現(xiàn)代英漢詞典》[bi:fli]牛肉一般地

8樓2015-02-02 22:58:17

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

查看全部 9 個回答

qchem

超級版主

應(yīng)助: 280 (大學(xué)生)
金幣: 6871.6
散金: 1411
紅花: 43
帖子: 2586
在線: 669.2小時
蟲號: 183011
注冊: 2006-02-12
專業(yè): 理論和計算化學(xué)

★
小木蟲: 金幣+0.5, 給個紅包，謝謝回帖

還可以考慮用內(nèi)部文件來實(shí)現(xiàn)

[ 發(fā)自小木蟲客戶端 ]

贊一下

回復(fù)此樓

4樓2015-02-01 23:50:26

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

馬墨跡

專家顧問

應(yīng)助: 1 (幼兒園)
金幣: 806
散金: 60
紅花: 2
帖子: 228
在線: 204.8小時
蟲號: 1981113
注冊: 2012-09-06
性別: GG
專業(yè): 防災(zāi)工程

馬克

[ 發(fā)自小木蟲客戶端 ]

回復(fù)此樓

墨跡不磨嘰

5樓2015-02-02 00:23:04

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

yebin2006

兌換貴賓

應(yīng)助: 1 (幼兒園)
金幣: 5042
散金: 19
帖子: 78
在線: 128.7小時
蟲號: 1269098
注冊: 2011-04-18
性別: GG
專業(yè): 電力系統(tǒng)

好帖，頂！

[ 發(fā)自小木蟲客戶端 ]

回復(fù)此樓

comeon~!

6樓2015-02-02 00:28:53

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

最具人氣熱帖推薦 [查看全部]	作者	回/看	最后發(fā)表

[考研] 材料0856 英一數(shù)二 323 求調(diào)劑 +3	袁sy 2026-04-01	3/150	2026-04-01 07:00 by 無際的草原
[考研] 考研調(diào)劑 +6	Amber00 2026-03-31	6/300	2026-04-01 00:42 by fmesaito
[考研] 求調(diào)劑 +9	akdhjs 2026-03-31	10/500	2026-03-31 17:14 by zhyzzh
[考研] 315求調(diào)劑 +6	akie... 2026-03-28	7/350	2026-03-31 16:48 by asdfzly
[考研] 362求調(diào)劑 +9	西南交材料專碩3 2026-03-31	9/450	2026-03-31 16:25 by w蟲蟲123
[考研] 一志愿中海洋材料357 +4	麥恩莉. 2026-03-30	4/200	2026-03-31 14:35 by 記事本2026
[考研] 085404 22408 315分 +5	zhuangyan123 2026-03-31	6/300	2026-03-31 13:48 by limeifeng
[考研] 266分，求材料冶金能源化工等調(diào)劑 +8	哇呼哼呼哼 2026-03-27	10/500	2026-03-31 13:35 by Huaxue_Wang
[考研] 一志愿大連理工大學(xué)，機(jī)械工程學(xué)碩，341 +3	西瓜田的守望者 2026-03-30	3/150	2026-03-31 11:08 by asdfzly
[考研] 085600，專業(yè)課化工原理，320分求調(diào)劑 +6	大饞小子 2026-03-29	6/300	2026-03-31 10:03 by 氯化亞硝酰
[考研] 本科211總分289，08工學(xué)真心求調(diào)劑 +3	utopiaE 2026-03-30	3/150	2026-03-30 23:42 by ms629
[考研] 一志愿北京化工大學(xué)材料與化工（085600）296求調(diào)劑 +25	稻妻小編 2026-03-26	25/1250	2026-03-30 20:11 by 滴滴上岸呀
[考研] 085600 286分材料求調(diào)劑 +11	麻辣魷魚 2026-03-27	12/600	2026-03-30 19:33 by Wang200018
[考研] 抱歉 +3	田洪有 2026-03-30	3/150	2026-03-30 19:11 by 迷糊CCPs
[考研] 334分一志愿武理材料求調(diào)劑 +16	李李不服輸 2026-03-26	16/800	2026-03-30 17:47 by wangjy2002
[考研] 本科雙非材料，跨考一志愿華電085801電氣，283求調(diào)劑，任何專業(yè)都可以 +6	芝士雪baoo 2026-03-28	8/400	2026-03-29 08:16 by 松花缸1201
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +6	崔wj 2026-03-26	6/300	2026-03-29 01:11 by hanserlol
[考研] 352分化工與材料 +5	海納百川Ly 2026-03-27	5/250	2026-03-28 03:39 by fmesaito
[考研] 一志愿鄭大085600，310分求調(diào)劑 +5	李瀟可 2026-03-26	5/250	2026-03-27 11:14 by 不吃魚的貓
[考研] 一志愿天津大學(xué)339材料與化工求調(diào)劑 +3	江往賣魚 2026-03-26	3/150	2026-03-26 09:42 by 王小欠i