... | ... | @@ -43,3 +43,15 @@ $` v_r (r,\theta) = v_{\theta}(r) \psi(r,\theta)\exp-\frac {(r-r_c)^2}{10\sigma |
|
|
|
|
|
where $`\psi(r,\theta) \in [-10^{-2},10^{-2}]`$ is a uniformly distributed random variable drawn for each grid cell.
|
|
|
|
|
|
2. The disk aspect ratio is quite unusual: $`h_0=0.09`$.
|
|
|
This choice is made to trigger a vortex trough Rossby Wave Instability (Lovelace 1999 criterium).
|
|
|
|
|
|
3. Output we only output the disk at 5 and 10 steps of 6.28 (physical time in code units).
|
|
|
|
|
|
The plot of the density field in disk10.h5 is obtained running:
|
|
|
|
|
|
```
|
|
|
<build>/utils/view_horizontal_slice.py disk10.h5 density -g polar
|
|
|
```
|
|
|
|
|
|
![cavity](uploads/71cf84621098b3af4866c03ca2c58249/cavity.png) |