IMPORTANT: Code error in calculating advective flux in bug fix release 0.7.2

[csyhuang]

The equation in #83 was correct but the fortran code implementation was unfortunately wrong. This shall be fixed in release 1.0.0. Thanks Christopher Polster for spotting that.

As mentioned in #83 , the correct expression for $u_e$ in the integral should be:

$$u_e(\lambda, \phi, \phi^\prime, z, t) \equiv \omega_e a \cos\phi = \frac{u(\lambda, \phi + \phi^\prime, z, t)\cos\phi}{\cos(\phi+\phi^\prime)} - u_{\text{REF}}(\phi, z, t)$$

where $\omega = \frac{u}{a\cos\phi}$ is the angular velocity (and $\phi^\prime \in [0, \Delta\phi]$).

The integral, which should have been computed as

in release v0.7.2 was computed by mistake as

Release v1.0.0 has this mistake corrected.