|
5 | 5 |
|
6 | 6 | import numpy as np |
7 | 7 | import xarray as xr |
8 | | -from parcels import FieldSet |
9 | 8 |
|
| 9 | +from parcels import FieldSet |
10 | 10 | from virtualship.instruments.drifter import Drifter, DrifterInstrument |
11 | 11 | from virtualship.models import Location, Spacetime |
12 | 12 | from virtualship.models.expedition import Waypoint |
@@ -111,15 +111,16 @@ def test_simulate_drifters(tmpdir) -> None: |
111 | 111 |
|
112 | 112 | def test_drifter_depths(tmpdir) -> None: |
113 | 113 | CONST_TEMPERATURE = 1.0 # constant temperature in fieldset |
| 114 | + DEPTH_FACTOR = 3.0 # factor to multiply surface values by at depth for test |
114 | 115 |
|
115 | 116 | v = np.full((2, 2, 2, 2), 1.0) |
116 | 117 | u = np.full((2, 2, 2, 2), 1.0) |
117 | 118 | t = np.full((2, 2, 2, 2), CONST_TEMPERATURE) |
118 | 119 |
|
119 | 120 | # different values at depth (random) |
120 | | - v[:, -1, :, :] = 1.0 * np.random.randint(0, 10) |
121 | | - u[:, -1, :, :] = 1.0 * np.random.randint(0, 10) |
122 | | - t[:, -1, :, :] = CONST_TEMPERATURE * np.random.randint(0, 10) |
| 121 | + v[:, -1, :, :] = 1.0 * DEPTH_FACTOR |
| 122 | + u[:, -1, :, :] = 1.0 * DEPTH_FACTOR |
| 123 | + t[:, -1, :, :] = CONST_TEMPERATURE * DEPTH_FACTOR |
123 | 124 |
|
124 | 125 | fieldset = FieldSet.from_data( |
125 | 126 | {"V": v, "U": u, "T": t}, |
@@ -184,6 +185,9 @@ def test_drifter_depths(tmpdir) -> None: |
184 | 185 | "Depth drifter depth should be constant" |
185 | 186 | ) |
186 | 187 |
|
| 188 | + if not drifter_surface.temperature[0] != drifter_depth.temperature[0]: |
| 189 | + breakpoint() |
| 190 | + |
187 | 191 | assert drifter_surface.temperature[0] != drifter_depth.temperature[0], ( |
188 | 192 | "Surface and deeper drifter should have different temperature measurements" |
189 | 193 | ) |
0 commit comments