Differences

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

--- tornadobox [2025/03/19 13:11] – 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 281: / Line 283: @@
 $ 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
+$ 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!**