diff --git a/src/PhysicalModels/MagneticModels.jl b/src/PhysicalModels/MagneticModels.jl
index e60c328..78e7f2f 100644
--- a/src/PhysicalModels/MagneticModels.jl
+++ b/src/PhysicalModels/MagneticModels.jl
@@ -9,8 +9,8 @@ struct Magnetic <: Magneto
   αr::Ref{Float64}
   χe::Float64
   
-function Magnetic(; μ0::Float64, αr::Ref{Float64}, χe::Float64=0.0)
-  new(μ0, αr, χe)
+function Magnetic(; μ0::Float64, αr::Float64, χe::Float64=0.0)
+  new(μ0, Ref(αr), χe)
 end
 function (obj::Magnetic)(Λ::Float64=1.0)
   μ, αr, χe = obj.μ, obj.αr, obj.χe
diff --git a/test/TestConstitutiveModels/PhysicalModelTests.jl b/test/TestConstitutiveModels/PhysicalModelTests.jl
index eab5106..dea916f 100644
--- a/test/TestConstitutiveModels/PhysicalModelTests.jl
+++ b/test/TestConstitutiveModels/PhysicalModelTests.jl
@@ -666,8 +666,7 @@ end
 @testset "Magnetic3D" begin
   ∇φ = VectorValue(1.0, 2.0, 3.0)
   a=40e-3 
-  ra=Ref(a)
-  modelID = Magnetic(μ0=1.2566e-6, αr=ra ,χe=0.0)
+  modelID = Magnetic(μ0=1.2566e-6, αr=a ,χe=0.0)
   Ψ, ∂Ψφ, ∂Ψφφ = modelID()
 
   Km=Kinematics(Magneto,Solid)
@@ -690,8 +689,7 @@ end
 @testset "Magnetic2D" begin
   ∇φ = VectorValue(1.0, 2.0)
   a=40e-3 
-  ra=Ref(a)
-  modelID = Magnetic(μ0=1.2566e-6, αr=ra ,χe=0.0)
+  modelID = Magnetic(μ0=1.2566e-6, αr=a ,χe=0.0)
   Ψ, ∂Ψφ, ∂Ψφφ = modelID()
 
   Km=Kinematics(Magneto,Solid)