bugfix 2

JuliaGeodynamics · Oct 28, 2021 · b8abea8 · b8abea8
1 parent f0b4e16
commit b8abea8
Show file tree

Hide file tree

Showing 3 changed files with 7 additions and 4 deletions.
diff --git a/src/Density/Density.jl b/src/Density/Density.jl
@@ -35,8 +35,11 @@ end
 # Calculation routine
 function ComputeDensity(P,T, s::ConstantDensity)
     @unpack ρ   = s
-
-    return ρ*1.0
+    if length(T)>1
+        return Value(ρ).*ones(size(T))
+    else
+        return ρ*1.0
+    end
 end
 
 # Print info 

diff --git a/src/Energy/Conductivity.jl b/src/Energy/Conductivity.jl
@@ -38,7 +38,7 @@ function ComputeConductivity(P, T, s::ConstantConductivity)
     @unpack k   = s
 
     if length(T)>1
-        return k.*ones(size(T))
+        return Value(k).*ones(size(T))
     else
         return Value(k)
     end

diff --git a/src/Energy/HeatCapacity.jl b/src/Energy/HeatCapacity.jl
@@ -35,7 +35,7 @@ end
 function ComputeHeatCapacity(P, T, s::ConstantHeatCapacity)
     @unpack cp   = s
     if length(T)>1
-        return cp*ones(size(T))
+        return Value(cp)*ones(size(T))
     else
         return Value(cp)
     end