Question

我在Fortran中是全新的，我需要编写一些相对简单的代码。

我有一些文件（各种文件，例如200个文件）;特定节点的每个文件，通过一些简化，每个文件包含：

T(i), X(i)

这些是我的输入，我希望输出文件包含：

T(i) X(i)1 X(i)2 ... X(i)n

问题在于我无法将输出文件的不同列中的数据分开，它们在1列中相继出现。

我的代码是：

PROGRAM Output

implicit none
integer ::nn,n,i,j,l
real,dimension(:),allocatable::t,x,y,z
character(len=10)::TD 

open(11,file='outputX.txt')
allocate (t(1000),x(1000),y(1000),z(1000)) 

n=5                  ! Number of Files
nn=50                ! Number of Rows in each File
 Do i=1,n            ! loop for opening different files 
   write(TD,10)i 
   write(*,*)TD 
   open(1,file=TD) 

      Do l=1,nn      ! loop for reading rows in each file
        read(1,*)t(j),x(j)
        write(11,*)x(j) !!!! This is my PROBLEM, all the data shows in
                               ! one column, I want each file in separately
      Enddo

 Enddo

    10  format('100',i3.3,'') 

    deallocate(x,y,z,t) 
    END PROGRAM Output

我得到的输出是这样的：

但实际上我想要：

11    21   31
12    22   32
13    23   33

Answer 1

您的代码有几个问题

 Do i=1,n            ! loop for opening different files 
   write(TD,10)i 
   write(*,*)TD 
   open(1,file=TD) 
   Do l=1,nn      ! loop for reading rows in each file
     read(1,*)t(j),x(j)
     write(11,*)x(j) !!!! This is my PROBLEM, all the data shows in
                               ! one column, I want each file in separately
  End do
End do

索引j完全未定义。您应该将j=1或0和j=j+1放在循环内的某处。

另一个问题是你的输出。您正在按顺序读取文件。将每个文件打印到单独的列中非常困难。每个文件的单独行很容易：

write(11,*) x(1:nn)

在内循环之后

。

或者更精细的控制并避免换行

write(11,'999(g0,1x)') x(1:nn)

（g0是一般编辑描述符，仅使用必要的宽度）。这只有在您修复我上面提到的j问题时才有效！

要将其放入单独的列中，您必须

同时打开所有文件，然后从每个文件中读取并在一个写入命令中打印读取数据。

或

将所有文件中的所有数据存储到2D数组中的单独列中，然后打印2D数组。

Answer 2

所以我得到的是x(1) belongs to 11，x(2) belongs to 12，x(3) belongs to 13等等，不是吗？你可以尝试这个：

write(11,100) x(j), x(j+3), x(j+6)
100 format(1X,11F20.4,11F20.4,11F20.4)

Answer 3

2D阵列......对于这个问题，它是向后的但它应该给你一些帮助......（希望）

PROGRAM Output
IMPLICIT NONE
INTEGER, PARAMETER              :: nfiles = 5
INTEGER, PARAMETER              :: nrows  = 50
integer                         :: iFile, iRow   !File and row counter/#
real,dimension(:,:),allocatable :: t,x,y,z       !The are now 2D
INTEGER,dimension(nFiles)       :: lFile         !logical file indexed by file#
LOGICAL,dimension(nFiles)       :: Open4Biz      !logical for closing
character(len=10)::TD                            !Unsure about this

open(11,file='outputX.txt')
allocate (t(iRow,iFile),x(iRow,iFile),y(iRow,iFile),z(iRow,iFile)) 

Open4Biz(:) = .FALSE.                            !Initialize

!n=5                  ! Number of Files              ^^Moved up/renamed^^
!nn=50                ! Number of Rows in each File  ^^Moved up/renamed^^
 Files_Loop1: Do iFile = 1, nFiles     !I think that the loop identifier is std f95 (std ifort anyhow)
   write(TD,10) iFile          !Unsure about this
   write(*,*)TD 
   lFile(iFile) = 10+iFile
   open(lFile(iFile),file=TD)  !Probably put in some logic if the file is not found
   Open4Biz(iFile) = .TRUE.
 ENDDO Files_Loop1

 Rows_Loop: Do iRow = 1, nn                !The loops are backwards from normal
   Files_Loop2: Do iFile = 1, nFiles
     read(lFile(iFile),*) t(iRow, iFile), x(iRow, iFile)
   Enddo Files_Loop2
   write(11,*) x(iRow,:)
 Enddo Rows_Loop

10  format('100',i3.3,'')

667  CONTINUE                         !This is label to 'jump to' from a bad open
Files_Loop3: Do iFile = 1, nFiles
  IF(Open4Biz(iFile) CLOSE(lFile(iFile))
ENDDO Files_Open_Loop

IF(ALLOCATED(X)) deallocate(x) 
IF(ALLOCATED(Y)) deallocate(y) 
IF(ALLOCATED(Z)) deallocate(z) 
IF(ALLOCATED(T)) deallocate(t)
END PROGRAM Output

如何从多个文件中读取数据并写入单个文件的列

3 个答案: