Skip to content

Commit

Permalink
Fix typo in documentation
Browse files Browse the repository at this point in the history
  • Loading branch information
jialinl6 committed Mar 1, 2024
1 parent 4cfe9ad commit 39de957
Showing 1 changed file with 4 additions and 4 deletions.
8 changes: 4 additions & 4 deletions docs/src/dry_deposition.md
Original file line number Diff line number Diff line change
Expand Up @@ -6,7 +6,7 @@ This is an implementation of a box model used to calculate changes in gas specie
Here's an example of how concentration of different species, such as SO₂, O₃, NO₂, NO, H₂O₂ and CH₂O change due to dry deposition.

We can create an instance of the model in the following manner:
```julia @example 1
```julia
using AtmosphericDeposition
using ModelingToolkit
using DifferentialEquations
Expand All @@ -17,11 +17,12 @@ using Unitful
model = DrydepositionG(t)
```
Before running any simulations with the model we need to convert it into a system of differential equations.
```julia @example 1
```julia
sys = structural_simplify(get_mtk(model))
tspan = (0.0, 3600*24)
u0 = [2.0,10.0,5,5,2.34,0.15] # initial concentrations of SO₂, O₃, NO₂, NO, H₂O₂, CH₂O
sol = solve(ODEProblem(sys, u0, tspan, []),AutoTsit5(Rosenbrock23()), saveat=10.0) # default parameters
```

```@setup 1
using AtmosphericDeposition
Expand All @@ -39,7 +40,6 @@ u0 = [2.0,10.0,5,5,2.34,0.15] # initial concentrations of SO₂, O₃, NO₂, NO
sol = solve(ODEProblem(sys, u0, tspan, []),AutoTsit5(Rosenbrock23()), saveat=10.0) # default parameters
```

```
which we can plot as
```@example 1
using Plots
Expand All @@ -48,7 +48,7 @@ plot(sol, xlabel="Time (second)", ylabel="concentration (ppb)", legend=:outerrig

## Parameters
The parameters in the model are:
```julia @example 1
```julia
parameters(sys) # [z, z₀, u_star, L, ρA, G, T, θ]
```
where ```z``` is the top of the surface layer [m], ```z₀``` is the roughness length [m], ```u_star``` is friction velocity [m/s], and ```L``` is Monin-Obukhov length [m], ```ρA``` is air density [kg/m3], ```T``` is surface air temperature [K], ```G``` is solar irradiation [W m-2], ```Θ``` is the slope of the local terrain [radians].
Expand Down

0 comments on commit 39de957

Please sign in to comment.