wip

README.md

@@ -1,4 +1,4 @@
-# Athena-nbody: an N-body code based on Athena++
+# Canoe: Comprehensive Atmosphere N' Ocean Engine
 
 [![build](https://github.com/chengcli/canoe/actions/workflows/main.yml/badge.svg)](https://github.com/chengcli/canoe/actions/workflows/main.yml)
 [![build](https://github.com/chengcli/canoe/actions/workflows/mac.yml/badge.svg)](https://github.com/chengcli/canoe/actions/workflows/mac.yml)

cmake/ktable_test.cmake

@@ -23,3 +23,4 @@ set(HYDROSTATIC ON)
 set(RFM ON)
 set(NETCDF ON)
 set(DISORT ON)
+set(PYTHON_BINDINGS ON)

doc/doxygen/configure

@@ -5,7 +5,7 @@
 PROJECT_NAME           = "Canoe"
 #PROJECT_NUMBER         = "v19"
 PROJECT_LOGO           = logo1.png
-PROJECT_BRIEF          = Planetary Atmosphere Simulator
+PROJECT_BRIEF          = Comprehensive Atmosphere N' Ocean Engine
 OUTPUT_DIRECTORY       = .
 
 # The doxygen documentation says this about the following flag:

doc/examples/robert/intro.dox

@@ -0,0 +1,19 @@
+<a name="Intro"></a>
+<h1> Introduction </h1>
+
+The second example shows the the rising of a warm bubble, which to some extent,
+resembles the mushroom cloud of an atomic bomb. This test case was proposed by
+@cite robert1993. The original setting was is 2D and we extended it to 3D.
+The size of the domain is 1 km $\times$ 1 km $\times$ 2.5km. The model resolution is
+uniformlly 5 m in each dimension. The bubble is initially centered at
+$x_0 = y_0 = 500 \text{m}, z_0 = 260 \text{m}$ with a Gaussian-shaped potential
+temperature anomaly as:
+
+@f{eqnarray*}
+    \Delta \theta &=& 0.5\;\text{K}, \qquad\quad\quad r \leq a \\
+                  &=& 0.5e^{-(r-a)^2/s^2}\;\text{K}, \quad r > a
+@f}
+
+where $r$ is the distance to the center, $a=50$ m and $s=100$ m. The background
+atmosphere is isentropic with surface pressure being 1 bar and surface temperature
+303.15 K. The bubble is then released and evolve to 1200 s.

doc/examples/robert/results.dox

@@ -0,0 +1,14 @@
+<h1> Results </h1>
+
+<p>Here we show the animation of the rising bubble</p>
+
+<div>
+<img src="mov/robert.gif" alt="Animation of a rising air bubble">
+</div>
+
+<p>The temperature anomaly of the bubble is colored in blue/red with blue being colder
+and red being warmer. Initially, the warm part of the bubble is hiden at its center
+and you only see the colder outer part. Multiple convective overturnings develop seen during the
+rising motion to deliver the warm gas upward and outward. Unlike the 2D solution shown in @cite robert1993,
+the 3D evolution of the bubble features a vortex ring surounding the core of the rising
+center and the head of the rising bubble is minimum compared to the ring.

doc/examples/robert/robert.cpp

@@ -0,0 +1 @@
+../../../examples/2019-Li-snap/robert.cpp