Я вычисляю несжимаемый 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 файл и .exe файл Я сделал на основе слайдов лекции