From 11be82fc82818ba9627de1b52737b6271fa8ce5d Mon Sep 17 00:00:00 2001 From: Jussi Aittoniemi Date: Sat, 18 Nov 2023 15:10:57 +0100 Subject: [PATCH] replace hard-coded heat capacity and seawater density with values from parameters --- warmth/mesh_model.py | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/warmth/mesh_model.py b/warmth/mesh_model.py index c9c87ca..80745e3 100644 --- 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 )