Question

我正在围绕rectangle计算不可压缩的Navier-Stokes 2D流。

通过使用FORTRAN 77中的类中提供的示例代码，我设法解决并获取了.data文件的输出，但是我在绘制它时遇到了麻烦。我的主程序如下：

  open(unit=60,file='flow.data',form='unformatted')
  open(unit=70,file='hist.data',form='formatted') 
  call setflw
  call setgrd
  call slvflw
  close(unit=60)
  close(unit=70)
  stop
  end

并且我有一个子例程，可将流数据输出到.dat文件中，例如：

  write(60) re,cfl,dt,nlast,time
  write(60) mx,i1,i2,my,j1,j2
  write(60)((x(i,j),y(i,j),i=1,mx),j=1,my)
  write(60)((u(i,j),v(i,j),p(i,j),i=1,mx),j=1,my)

所以我确定我已经在flow.data文件中找到了它的数据。但是（我认为）是二进制文件，并且由于我还是个初学者，所以我不知道如何绘制它们。我想知道在FORTRAN 77或Python或gnuplot中（如果可能的话）分别绘制流线/压力场的方法

这是我根据讲座幻灯片制作的.for file和.exe file

在Python / FORTRAN中绘制Karman涡流的流线

0 个答案: