vvt 发表于 2019-10-24 16:48
[mw_shl_code=fortran,true]    open(fid,file=trim(fpath)//'/p1.dat')
    write(fid,120) x(1:n)
    wr ...

刚刚也找到原因了，十分感谢！

[Fortran] 纯文本查看 复制代码

    open(fid,file=trim(fpath)//'/p1.dat')
    write(fid,120) x(1:n)
    write(fid,120) h0
    do i = 1,n
        write(fid,120) p(i,1:n)
    end do
    do i = 1,n
        write(fid,120) h(i,1:n)
    end do
    close(fid)
120 format(512(es16.8e3,1x))