Question

我尝试使用 nxn 矩阵格式编写输出file.dat。

我编写代码，但输出是值为f的列。

现在的问题是：如何更改要写入的文件的输出格式？

自： 1 2 4 五 ...

到：1,2,3,4 // 5,6,8，.. //

program eccen
    implicit none
    integer, parameter:: grid=800
    integer::i,j,k,n,m
    real*8,allocatable::f(:,:)
    real*8::xx(grid),yy(grid),mval,Mxval
    real*8,allocatable::x(:),y(:)

    open(10,file='3d_disk.txt')
    n=0
        DO
                READ(10,*,END=100)
                n=n+1
        END DO

 100     continue
        rewind(10)

    allocate(x(n),y(n))

    do i=1, n
        read(10,*) x(i),y(i)
    end do


    mval=-20.
    Mxval=20.
    do i=1, grid
        xx(i) = mval + ((Mxval - mval)*(i-1))/(grid-1)
        yy(i) = mval + ((Mxval - mval)*(i-1))/(grid-1)
    end do

    open(20,file='3d_map.dat')

    allocate(f(n,n))
    f=0
    do i=1,grid
        do j=1,grid
            m=0.
            do k=1, n
                if (x(k) > xx(i) .and. x(k) < xx(i+1) .and. &
                 & y(k) > yy(j) .and. y(k) < yy(j+1)) then  
                    m=m+1 ! CONTA IL NUMERO DI PARTICELLE
                end if
            end do
            f(i,j)=float(m+1)

我认为修改必须在这里：

                write(20,*) f(i,j)
            end do
            write(20,*)
     print *,i
        end do    
end program eccen

为：

    do i=1,grid
      do j=1,grid
        write(20,*) f(i,j)
      end do
    end do

end do
            write(20,*)
     print *,i
        end do

 end program eccen

Answer 1

本声明

write(20,*) f(i,j)

将写入f(i,j)的值然后移动到下一行，因此您获得的文件正是您的代码指定的文件。如果您希望文件包含n行，每行包含n值，请尝试

write(20,*) f(i,:)

应该很好地将数组的整行i写入输出文件中的一行。当然，这会使您的循环超过j，因此您可以删除它。

在fortran中以数组格式写

1 个答案: