diff --git a/docs/src/dry_deposition.md b/docs/src/dry_deposition.md
index ece25144..8b944633 100644
--- a/docs/src/dry_deposition.md
+++ b/docs/src/dry_deposition.md
@@ -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
@@ -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
@@ -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
@@ -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].