... | ... | @@ -89,10 +89,10 @@ eqTracer = slim.ShallowWaterTracer2d(domain, "implicit",eq, name="tracer", offli |
|
|
- ```initial_time``` initial time for the temporal solver
|
|
|
- ```final_time``` final time for the final solver
|
|
|
|
|
|
Import the tracer initial condition from the file created in pre-processing
|
|
|
Import the tracer initial condition from the file created in preprocessing
|
|
|
|
|
|
```python
|
|
|
eqTracer.set_initial_condition(pre_data_dir+"/tracer/tracer.idx")
|
|
|
eqTracer.set_initial_condition(data_dir+"/tracer/tracer.idx")
|
|
|
```
|
|
|
|
|
|
Set the diffusivity model used
|
... | ... | @@ -100,7 +100,7 @@ Set the diffusivity model used |
|
|
```python
|
|
|
eqTracer.set_diffusivity("okubo", okubo_coefficient=0.04)
|
|
|
```
|
|
|
- ```"okubo"``` diffusivity model used: either a constant diffusivity or the okubo scheme
|
|
|
- ```"okubo"``` diffusivity model used: either a constant diffusivity or the Okubo scheme
|
|
|
- ```okubo_coefficient``` coefficient for Okubo scheme [in m^0.85/s] (default: 0.03 m^0.85/s)
|
|
|
|
|
|
Add the boundary condition
|
... | ... | @@ -108,8 +108,8 @@ Add the boundary condition |
|
|
```python
|
|
|
eqTracer.set_boundary_coast("coast")
|
|
|
eqTracer.set_boundary_coast("island")
|
|
|
eqTracer.set_boundary_open("mexico", (pre_data_dir + "concentration_atlantic.nc",'tracer'))
|
|
|
eqTracer.set_boundary_open("atlantic", (pre_data_dir + "concentration_atlantic.nc",'tracer'))
|
|
|
eqTracer.set_boundary_open("mexico", (data_dir + "concentration_atlantic.nc",'tracer'))
|
|
|
eqTracer.set_boundary_open("atlantic", (data_dir + "concentration_atlantic.nc",'tracer'))
|
|
|
```
|
|
|
- ```boundary_coast``` condition for coast (no flux condition)
|
|
|
- ```boundary_open``` impose the external condition
|
... | ... | @@ -144,9 +144,9 @@ Create a ```Domain``` object, the computational domain defined by the mesh and t |
|
|
```python
|
|
|
domain = slim.Domain(data_dir+'mesh.msh', data_dir_base+'bathymetry_smooth/bathymetry_smooth.idx')
|
|
|
```
|
|
|
- ```(pre_data_dir+"gravity.nc","gravity")``` is a pre-processing file with the gravity over the domain
|
|
|
- ```(data_dir+"gravity.nc","gravity")``` is a preprocessing file with the gravity over the domain
|
|
|
- ```rho``` is the water density
|
|
|
- ```solve_on_sphere``` ```True``` if you solve the equation on a sphere, ```False``` if it is in a projected doamin
|
|
|
- ```solve_on_sphere``` ```True``` if you solve the equation on a sphere, ```False``` if it is in a projected domain
|
|
|
|
|
|
Create the Tracer equation
|
|
|
|
... | ... | @@ -161,10 +161,10 @@ eqTracer = slim.ShallowWaterTracer2d(domain, "implicit","output_offline/", name= |
|
|
- ```final_time``` final time for the final solver
|
|
|
|
|
|
|
|
|
Import the tracer initial condition from the file created in pre-processing
|
|
|
Import the tracer initial condition from the file created in preprocessing
|
|
|
|
|
|
```python
|
|
|
eqTracer.set_initial_condition(pre_data_dir+"/tracer/tracer.idx")
|
|
|
eqTracer.set_initial_condition(data_dir+"/tracer/tracer.idx")
|
|
|
'''
|
|
|
|
|
|
Set the diffusivity model used
|
... | ... | @@ -179,8 +179,8 @@ Add the boundary condition |
|
|
```python
|
|
|
eqTracer.set_boundary_coast("coast")
|
|
|
eqTracer.set_boundary_coast("island")
|
|
|
eqTracer.set_boundary_open("mexico", (pre_data_dir + "concentration_atlantic.nc",'tracer'))
|
|
|
eqTracer.set_boundary_open("atlantic", (pre_data_dir + "concentration_atlantic.nc",'tracer'))
|
|
|
eqTracer.set_boundary_open("mexico", (data_dir + "concentration_atlantic.nc",'tracer'))
|
|
|
eqTracer.set_boundary_open("atlantic", (data_dir + "concentration_atlantic.nc",'tracer'))
|
|
|
```
|
|
|
- ```boundary_coast``` condition for coast (no flux condition)
|
|
|
- ```boundary_open``` impose the external condition
|
... | ... | |