2012-09-17 109 views
1

我是新來的Fortran,我想寫一個程序來讀取.txt文件,其中我有24480行和〜6000列。 在每行(作爲個體)中,我有基因型如1和2所示,例如,如果在第一行中,我有204種基因型,這種基因型的前半部分(= 102)屬於個體的父系,基因型屬於個體的大壩。另外每行的值不相等。那麼,我該如何定義Fortran以逐行讀取此文件,並將每行分爲兩部分,並將每個元素(i)放在均值+ i旁邊。行閱讀和編輯在Fortran 90

例如我展示我的文件的兩行簡短: ROW1:112122121112122111112121111211122121111121和 2行:21112111112112222121112121211121221212121111121112 1212

等。 任何幫助將提前感激。

+0

嘗試用正確的標籤。你將有更大的回答機會。 –

+1

對不起,你能清楚地指出你想要的結果。最好展示一個例子結果。 –

+0

首先我想讀取並輸入我的數據從文本文件到fortran並選擇感興趣的基因型,其次,選擇感興趣的個體,最後爲選定的基因型和個體生成一個基因型文件,例如作爲2維矩陣。 – user1677032

回答

2

你需要閱讀每行一個字符串,獲取其長度與LEN_TRIM(),然後讀取每個半到整型數組:

 character*10000 line 

    read(unit,'(a)')line 
    len=len_trim(line) 
    read(line,'(10000i1)')a(:len/2),b(:len/2)