replace hard-coded heat capacity and seawater density with values fro…

…m parameters
equinor · Nov 18, 2023 · 11be82f · 11be82f
1 parent 587b059
commit 11be82f
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/warmth/mesh_model.py b/warmth/mesh_model.py
@@ -650,7 +650,7 @@ def boundary(x):
 
             self.thermalCond.x.array[i] = conductivity_effective
             self.mean_porosity.x.array[i] = mean_porosity
-            self.c_rho.x.array[i] = 1000*((self.c_rho0.x.array[i]/1000) * (1-mean_porosity) + mean_porosity*1000)
+            self.c_rho.x.array[i] = self._parameters.cp*((self.c_rho0.x.array[i]/self._parameters.cp) * (1-mean_porosity) + mean_porosity*self._parameters.rhowater)
 
         # self.rhpFcn.x.array[:] = np.multiply( self.rhp0.x.array[:], (1.0-self.mean_porosity.x.array[:]) )
         self.rhpFcn.x.array[:] = np.multiply( self.rhp0.x.array[:], 1.0 )