use t = default_t()

Author: TorkelE

docs/src/api/catalyst_api.md

@@ -36,7 +36,7 @@ models, and stochastic chemical kinetics jump process models.
 
 ```@example ex1
 using Catalyst, DifferentialEquations, Plots
-import Catalyst: t_nounits as t
+t = default_t()
 @parameters β γ
 @species S(t) I(t) R(t)
 

docs/src/catalyst_applications/simulation_structure_interfacing.md

@@ -104,7 +104,7 @@ Finally, we note that we cannot change the values of solution unknowns or parame
 Catalyst is built on an *intermediary representation* implemented by (ModelingToolkit.jl)[https://github.com/SciML/ModelingToolkit.jl]. ModelingToolkit is a modelling framework where one first declares a set of symbolic variables and parameters using e.g.
 ```@example ex2
 using ModelingToolkit
-import ModelingToolkit: t_nounits as t
+t = default_t()
 @parameters σ ρ β
 @variables x(t) y(t) z(t)
 nothing # hide

docs/src/catalyst_functionality/chemistry_related_functionality.md

@@ -10,7 +10,7 @@ While Catalyst has primarily been designed around the modelling of biological sy
 We will first show how to create compound species through [programmatic model construction](@ref programmatic_CRN_construction), and then demonstrate using the DSL. To create a compound species, use the `@compound` macro, first designating the compound, followed by its components (and their stoichiometries). In this example, we will create a CO₂ molecule, consisting of one C atom and two O atoms. First, we create species corresponding to the components:
 ```@example chem1
 using Catalyst
-import Catalyst: t_nounits as t
+t = default_t()
 @species C(t) O(t) 
 ```
 Next, we create the `CO2` compound species:
@@ -97,7 +97,7 @@ In all of these cases, the side to the left of the `~` must be enclosed within `
 ### Compounds with multiple independent variables
 While we generally do not need to specify independent variables for compound, if the components (together) have more than one independent variable, this have to be done:
 ```@example chem1
-import Catalyst: t_nounits as t
+t = default_t()
 @variables s
 @species N(s) O(t) 
 @compound NO2(t,s) ~ N + 2O
@@ -118,7 +118,7 @@ which correctly finds the (rather trivial) solution `C + 2O --> CO2`. Here we no
 Let us consider a more elaborate example, the reaction between ammonia (NH₃) and oxygen (O₂) to form nitrogen monoxide (NO) and water (H₂O). Let us first create the components and the unbalanced reaction:
 ```@example chem2
 using Catalyst # hide
-import Catalyst: t_nounits as t
+t = default_t()
 @species N(t) H(t) O(t) 
 @compounds begin
     NH3 ~ N + 3H

docs/src/catalyst_functionality/compositional_modeling.md

@@ -51,7 +51,7 @@ plot(TreePlot(rn), method=:tree, fontsize=12, nodeshape=:ellipse)
 We could also have directly constructed `rn` using the same reaction as in
 `basern` as
 ```@example ex1
-import Catalyst: t_nounits as t
+t = default_t()
 @parameters k
 @species A(t), B(t), C(t)
 rxs = [Reaction(k, [A,B], [C])]
@@ -114,7 +114,7 @@ ability to substitute the value of these variables into the DSL (see
 [Interpolation of Julia Variables](@ref dsl_description_interpolation_of_variables)). To make the repressilator we now make
 three genes, and then compose them together
 ```@example ex1
-import Catalyst: t_nounits as t
+t = default_t()
 @species G3₊P(t)
 @named G1 = repressed_gene(; R=ParentScope(G3₊P))
 @named G2 = repressed_gene(; R=ParentScope(G1.P))

docs/src/catalyst_functionality/constraint_equations.md

@@ -29,7 +29,7 @@ Here, to create differentials with respect to time (for our differential equatio
 
 ```@example ceq1
 using Catalyst, DifferentialEquations, Plots
-import Catalyst: t_nounits as t, D_nounits as D
+t = default_t(), D_nounits as D
 
 # set the growth rate to 1.0
 @parameters λ = 1.0
@@ -72,7 +72,7 @@ As an alternative to the previous approach, we could have constructed our
 `ReactionSystem` all at once by directly using the symbolic interface:
 ```@example ceq2
 using Catalyst, DifferentialEquations, Plots
-import Catalyst: t_nounits as t, D_nounits as D
+t = default_t(), D_nounits as D
 
 @parameters λ = 1.0
 @variables V(t) = 1.0
@@ -104,7 +104,7 @@ advanced_simulations) tutorial.
 Let's first create our equations and unknowns/species again
 ```@example ceq3
 using Catalyst, DifferentialEquations, Plots
-import Catalyst: t_nounits as t, D_nounits as D
+t = default_t(), D_nounits as D
 
 @parameters λ = 1.0
 @variables V(t) = 1.0
@@ -128,7 +128,7 @@ plot(sol; plotdensity = 1000)
 ```
 We can also model discrete events. Similar to our example with continuous events, we start by creating reaction equations, parameters, variables, and unknowns. 
 ```@example ceq3
-import Catalyst: t_nounits as t
+t = default_t()
 @parameters k_on switch_time k_off
 @species A(t) B(t)
 

docs/src/catalyst_functionality/dsl_description.md

@@ -54,7 +54,7 @@ the model. We do this by creating a mapping from each symbolic variable
 representing a chemical species to its initial value
 ```@example tut2
 # define the symbolic variables
-import Catalyst: t_nounits as t
+t = default_t()
 @species X(t) Y(t) Z(t) XY(t) Z1(t) Z2(t)
 
 # create the mapping
@@ -610,7 +610,7 @@ parameters outside of the macro, which can then be used within expressions in
 the DSL (see the [Programmatic Construction of Symbolic Reaction Systems](@ref programmatic_CRN_construction)
 tutorial for details on the lower-level symbolic interface). For example,
 ```@example tut2
-import Catalyst: t_nounits as t
+t = default_t()
 @parameters k α
 @species A(t)
 spec = A

docs/src/catalyst_functionality/example_networks/hodgkin_huxley_equation.md

@@ -15,7 +15,7 @@ We begin by importing some necessary packages.
 using ModelingToolkit, Catalyst, NonlinearSolve
 using DifferentialEquations, Symbolics
 using Plots
-import Catalyst: t_nounits as t, D_nounits as D
+t = default_t(), D_nounits as D
 ```
 
 We'll build a simple Hodgkin-Huxley model for a single neuron, with the voltage,

docs/src/catalyst_functionality/example_networks/smoluchowski_coagulation_equation.md

@@ -52,7 +52,7 @@ end
 We'll store the reaction rates in `pars` as `Pair`s, and set the initial condition that only monomers are present at ``t=0`` in `u₀map`.
 ```julia
 # unknown variables are X, pars stores rate parameters for each rx
-import Catalyst: t_nounits as t
+t = default_t()
 @species k[1:nr] (X(t))[1:N]
 pars = Pair.(collect(k), kv)
 

docs/src/catalyst_functionality/parametric_stoichiometry.md

@@ -35,7 +35,7 @@ rx.substrates[1],rx.substoich[1]
 We could have equivalently specified our systems directly via the Catalyst
 API. For example, for `revsys` we would could use
 ```@example s1
-import Catalyst: t_nounits as t
+t = default_t()
 @parameters k₊, k₋, m, n
 @species A(t), B(t)
 rxs = [Reaction(k₊, [A], [B], [m], [m*n]),
@@ -177,7 +177,7 @@ calculate and plot the average amount of protein (which is also plotted in the
 MomentClosure.jl
 [tutorial](https://augustinas1.github.io/MomentClosure.jl/dev/tutorials/geometric_reactions+conditional_closures/)).
 ```@example s1
-import Catalyst: t_nounits as t
+t = default_t()
 function getmean(jprob, Nsims, tv)
     Pmean = zeros(length(tv))
     @variables P(t)

docs/src/catalyst_functionality/programmatic_CRN_construction.md

@@ -15,12 +15,12 @@ and then define symbolic variables for each parameter and species in the system
 (the latter corresponding to a `variable` or `unknown` in ModelingToolkit
 terminology)
 ```@example ex
-import Catalyst: t_nounits as t
+t = default_t()
 @parameters α K n δ γ β μ
 @species m₁(t) m₂(t) m₃(t) P₁(t) P₂(t) P₃(t)
 nothing    # hide
 ```
-Note: Each species is declared as a function of time. Here, we first import the *time independent variable using `import Catalyst: t_nounits as t`, and then use it to declare out species.
+Note: Each species is declared as a function of time. Here, we first import the *time independent variable using `t = default_t()`, and then use it to declare out species.
 
 !!! note
        For users familiar with ModelingToolkit, chemical species must be declared
@@ -117,7 +117,7 @@ Reaction(rate, nothing, [P₁,...,Pₙ], nothing, [β₁,...,βₙ])
 Finally, we note that the rate constant, `rate` above, does not need to be a
 constant or fixed function, but can be a general symbolic expression:
 ```julia
-import Catalyst: t_nounits as t
+t = default_t()
 @parameters α, β
 @species A(t), B(t)
 rx = Reaction(α + β*t*A, [A], [B])
@@ -133,7 +133,7 @@ reactions using the [`@reaction`](@ref) macro.
 
 For example, the repressilator reactions could also have been constructed like
 ```julia
-import Catalyst: t_nounits as t
+t = default_t()
 @species P₁(t) P₂(t) P₃(t)
 rxs = [(@reaction hillr($P₃,α,K,n), ∅ --> m₁),
        (@reaction hillr($P₁,α,K,n), ∅ --> m₂),
@@ -162,7 +162,7 @@ rx = @reaction hillr(P,α,K,n), A --> B
 ```
 is equivalent to
 ```julia
-import Catalyst: t_nounits as t
+t = default_t()
 @parameters P α K n
 @variables A(t) B(t)
 rx = Reaction(hillr(P,α,K,n), [A], [B])