1

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

我编写了代码，但输出是一列值 f。

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

从： 1 2 4 5 ...

到：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
4

1 回答 1

0

这个说法

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

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

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

这应该可以很好地将数组的整行写入i输出文件中的单行。当然，这会使您的循环变得j多余，因此您可以将其删除。

于 2013-05-20T16:53:31.817 回答