|
50 | 50 | }, |
51 | 51 | { |
52 | 52 | "cell_type": "code", |
53 | | - "execution_count": 3, |
| 53 | + "execution_count": null, |
54 | 54 | "metadata": { |
55 | 55 | "tags": [ |
56 | 56 | "hide-input" |
|
59 | 59 | "outputs": [], |
60 | 60 | "source": [ |
61 | 61 | "R0 = 50e3\n", |
62 | | - "ETA0 = .2\n", |
| 62 | + "ETA0 = -.2\n", |
63 | 63 | "LAT0 = 36\n", |
64 | 64 | "\n", |
65 | 65 | "dxy = 5e3" |
|
74 | 74 | }, |
75 | 75 | { |
76 | 76 | "cell_type": "code", |
77 | | - "execution_count": 4, |
| 77 | + "execution_count": null, |
78 | 78 | "metadata": { |
79 | 79 | "tags": [ |
80 | 80 | "hide-input" |
81 | 81 | ] |
82 | 82 | }, |
83 | 83 | "outputs": [], |
84 | 84 | "source": [ |
85 | | - "def simulate_gaussian_ssh(R0, eta0, R):\n", |
86 | | - " return -eta0 * jnp.exp(-(R / R0)**2)\n", |
| 85 | + "def simulate_gaussian_ssh(R0, eta0, r):\n", |
| 86 | + " return -eta0 * jnp.exp((r / R0)**2)\n", |
87 | 87 | "\n", |
88 | 88 | "\n", |
89 | 89 | "def simulate_gaussian_geos(R0, X, Y, ssh, coriolis_factor):\n", |
|
94 | 94 | " return u_geos, v_geos\n", |
95 | 95 | "\n", |
96 | 96 | "\n", |
97 | | - "def simulate_gaussian_cyclo(R0, X, Y, R, ssh, coriolis_factor):\n", |
98 | | - " azim_geos = -(2 * geometry.GRAVITY * R * ssh / (coriolis_factor * R0 ** 2))\n", |
99 | | - " azim_cyclo = 2 * azim_geos / (1 + jnp.sqrt(1 + 4 * azim_geos / (coriolis_factor * R)))\n", |
100 | | - " u_cyclo = -azim_cyclo * Y / R\n", |
101 | | - " v_cyclo = azim_cyclo * X / R\n", |
| 97 | + "def simulate_gaussian_cyclo(R0, X, Y, r, ssh, coriolis_factor):\n", |
| 98 | + " azim_geos = -(2 * geometry.GRAVITY * r * ssh / (coriolis_factor * R0 ** 2))\n", |
| 99 | + " azim_cyclo = 2 * azim_geos / (1 + jnp.sqrt(1 + 4 * azim_geos / (coriolis_factor * r)))\n", |
| 100 | + " u_cyclo = -azim_cyclo * Y / r\n", |
| 101 | + " v_cyclo = azim_cyclo * X / r\n", |
102 | 102 | " return u_cyclo, v_cyclo\n", |
103 | 103 | "\n", |
104 | 104 | "\n", |
|
0 commit comments