... | ... | @@ -58,7 +58,7 @@ mesh = data_dir+"mesh.msh" |
|
|
path_bathy = data_dir_base+"bathymetry_smooth/bathymetry_smooth.idx"
|
|
|
```
|
|
|
|
|
|
Create the reef map with the following arguments
|
|
|
Create reef map with the following arguments
|
|
|
- .msh file with the mesh
|
|
|
- .msh file whose field value is 1 at reefs and 0 elsewhere
|
|
|
- Mesh projection
|
... | ... | @@ -96,7 +96,7 @@ Create a ```Particle_tracker``` object with the following arguments |
|
|
tracker = slimPost.Particle_tracker(mesh, "output", path_bathy, 50, 118 ,1, 1800, 3600, diffusivity = 'okubo_map.msh', path_output = "output/particles/")
|
|
|
```
|
|
|
|
|
|
Set the open boundaries
|
|
|
Set open boundaries
|
|
|
|
|
|
```python
|
|
|
tracker.set_open_boundary_condition(["mexico","atlantic_east","atlantic_north"])
|
... | ... | @@ -129,9 +129,9 @@ for i in range(100): |
|
|
tracker.add_particle_at_point(623320, 3.01958e6, 50,i, i+2, 0, i+2)
|
|
|
```
|
|
|
|
|
|
Set the settling properties
|
|
|
Set settling properties
|
|
|
- settle over shallow (0), deep (1) or both (2) reefs (here ```2```)
|
|
|
- number of shallow reefs (taken following ID order, here ```10```)
|
|
|
- number of shallow reefs (X first reefs following ID order, here ```10```)
|
|
|
|
|
|
```python
|
|
|
tracker.set_settling_over_which_reefs(2,10)
|
... | ... | |