Differences

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

--- dustdevils [2025/04/01 09:31] – [Some interactive images at $t=29.8$] admin
+++ dustdevils [2025/04/12 19:29] (current) – admin
@@ Line 1: / Line 1: @@
-===== Dustdevils redux =====
+===== Dust devils redux: using numerical python =====
-This page shows a revisit to {{ :dustdevils:fiedlerkanak2001.pdf | Rayleigh-BeÂnard convection as a tool for studying dust devils}}, but with the computations done with numerical python.
+This page shows a revisit to {{ :dustdevils:fiedlerkanak2001.pdf | Rayleigh-Benard convection as a tool for studying dust devils}}, but with the computations done with numerical python.
+It is a substantial upgrade to [[https://12characters.net/dustdevil/ |my old dust devil page]].
 To reproduce these simulations, you may use my ipython notebook:
-[[https://github.com/bfiedler/metr4323/blob/main/N140_convection3d.ipynb|3D convection, periodic]]
+[[https://github.com/bfiedler/metr4323/blob/main/N140_convection3d.ipynb|3D convection, periodic]].
-In cell number 28, set ''nexp'' to 91:
+To reproduce the visualizations, you will need {{https://github.com/bfiedler/dustdevils/blob/main/dustdevils.tar.gz | my companion python programs}}.
+In the notebook, in cell number 28, set ''nexp'' to 91:
 <code>
 # choose an integer, see choices for nexp below
@@ Line 48: / Line 52: @@
 Try
 <code>
-./trace.py R4096g0w4 1 400
+./hist.py R4096g0w4 1 400
 </code>
 if all went okay,
 <code>
-./trace.py R4096g0w4 1 400 > R4096g0w4.hist
+./hist.py R4096g0w4 1 40 > R4096g0w4.hist
-</code>
+</code>0
 And then
 <code>
-./plottrace.py -s R4096g0w4.hist
+./plothist.py -s R4096g0w4.hist
 </code>
 That should produce the following image, which shows a dust devil event just prior to t=30.
@@ Line 71: / Line 74: @@
 <code>
 mkdir pngs
-./hcontour.py -o pngs R4096g0w4 000001 000400
+./npy2contour.py -o pngs R4096g0w4 000001 000400
 ./maketitle.py -o titleh.png
 ./png2mp4.py pngs/0 -t titleh.png
@@ Line 86: / Line 89: @@
 <code>
 mkdir vista
-./looknpy.py  R4096g0w4 000270 000320 -o vista
+./npy2vista.py  R4096g0w4 000270 000320 -o vista
 ./maketitlevista.py
 ./png2mp4.py vista/0 -t titlevista.png
@@ Line 96: / Line 99: @@
 <code>
-./looknpy.py -s R4096g0w4 000298
+./npy2vista.py -s R4096g0w4 000298
 open scene.html
 </code>
@@ Line 105: / Line 108: @@
 $\zeta=-40$ (purple)   , $p=-2$ (green), and $s$ at the surface
-I did not finish implementing the command line options for ''looknpy.py'',
+I did not finish implementing the command line options for ''npy2vista.py'',
 but the source code can be easily edited to select the desired isosurfaces.
+(You may want to copy ''npy2vista.py'' to a file of some other name, and edit that...)
 Change these lines:
@@ Line 142: / Line 146: @@
 <code>
-./looknpy.py -s R4096g0w4 000298
+./npy2vista.py -s R4096g0w4 000298
 open scene.html
 </code>
@@ Line 149: / Line 153: @@
 $w=.5$ (green), $w=-.5$ (grey) , $c=6$ (red)
+The suitably renamed ''scene.html'' can then be uploaded to a website, as done here.