Differences

This shows you the differences between two versions of the page.

--- tornadobox [2025/03/19 12:54] – admin
+++ tornadobox [2025/04/03 17:13] (current) – admin
@@ Line 1: / Line 1: @@
-==== Tornado Box====
+==== 3-D Fortran CFD code for tornado-like vortices ====
-This page provides the Fortran program used for
+Grab the code: [[https://github.com/bfiedler/tornadobox | tornadobox at gitbub]].
+This page describes the Fortran program used for
 {{ :tornadobox:suction_vortices_and_spiral_breakdown_in_numerical.pdf | my 2009 QJ article}}.
 I also still offer [[https://12characters.net/suction/|a more ancient page]]
@@ Line 250: / Line 252: @@
 .1,.01
 </code>
-That should finish in a few minutes. To simulate out to $t=30$, I  needed about 11 hours.
+That should finish in a few minutes. To simulate out to $t=30$, I  needed about 11 hours and 100G of disk space.
 Your previous ''answers'' file should still be appropriate.
 <code>
@@ Line 277: / Line 279: @@
 </code>
 You should see you have a lot of ''.sav'' files. Now lets extract what we need into python pickle files.
-Such files are easily readable by our plotting programs.
+Such files are easily readable by our plotting programs. This may need about 20 minutes:
 <code>
 $ sav2pkl.py b*.sav
 </code>
+If your python does not have the ''pyvista'' module, you will need to install it.
+To make the times tidy and complete, let's rename a file:
+<code>
+$ mv bF10E015000.pkl bF10D015000.pkl
+$ pkl2png.py -s -hg -vg   bF10D01500
+</code>
+Note that ''pkl2png.py'' will do the glob of ''bF10D01500*.pkl'', in this case finding one ''.pkl'' file.
+The ''-s'' option should have shown the image do the screen, like the first one at the top of this page.
+Next we zoom in to 25% of the camera distance, and plot at $t=30$:
+<code>
+$ pkl2png.py -s -hg -c .25 bF10D030000
+</code>
+And add a $w$ isosurface:
+<code>
+$ pkl2png.py -s -hg -w -c .25 bF10D030000
+</code>
+$ Now we go for it, and plot the frames for movie.  No ''-s'' option:
+<code>
+$ mkdir pngs
+$ pkl2png.py -hg -w -c .25 bF10D0 pngs
+</code>
+That required 10 minutes on my computer to produce 301 ''.png'' files.
+<code>
+$ cd pngs
+$ eog bF10D015000.png
+</code>
+The ''eog'' image viewer allows for the right arrow key to rapidly step through the images.
+<code>
+$ cd ~/boxout
+$ maketitle.py pngs/bF10D015000.png
+$ pip install moviepy
+$ png2mp4.py -t title.png pngs/bF10D
+$ mpv video.mp4
+</code>
+**FINISHED!**