fix all)unitless deprication

Author: TorkelE

src/networkapi.jl

@@ -1633,9 +1633,10 @@ function validate(rs::ReactionSystem, info::String = "")
     timeunits = get_unit(get_iv(rs))
 
     # no units for species, time or parameters then assume validated
-    (specunits in (MT.unitless, nothing)) && (timeunits in (MT.unitless, nothing)) &&
-        MT.all_dimensionless(get_ps(rs)) && return true
-
+    if (specunits in (MT.unitless, nothing)) && (timeunits in (MT.unitless, nothing))
+        all(u == 1.0 for u in ModelingToolkit.get_unit(get_ps(rs))) || return true
+    end
+    
     rateunits = specunits / timeunits
     for rx in get_rxs(rs)
         rxunits = get_unit(rx.rate)

src/reactionsystem.jl

@@ -532,8 +532,10 @@ struct ReactionSystem{V <: NetworkProperties} <:
         end
 
         if isempty(sivs) && (checks == true || (checks & MT.CheckUnits) > 0)
-            nonrx_eqs = Equation[eq for eq in eqs if eq isa Equation]
-            MT.all_dimensionless([unknowns; ps; iv]) || check_units(nonrx_eqs)
+            if !all(u == 1.0 for u in ModelingToolkit.get_unit([unknowns; ps; iv]))
+                nonrx_eqs = Equation[eq for eq in eqs if eq isa Equation]
+                check_units(nonrx_eqs)
+            end
         end
 
         rs = new{typeof(nps)}(eqs, rxs, iv, sivs, unknowns, spcs, ps, var_to_name, observed,